Как изменить код Pari для проверки гипотезы

Question

liste(lim,mul,base=34,step=5590,offset=588)=my(v=List(),X=mul*base); lim\=1; while(X<lim,forstep(n=X+offset+step,lim,step, listput(v,n)); X*=base); Set(v)
list(lim)=setunion(liste(lim,23),liste(lim,223))
v=list(1e6)

Это код Пари.

Теперь, как я могу изменить эту программу, чтобы увидеть, правда ли, что если числа, сгенерированные с помощью этого кода Пари, кратны 43, то они совпадают либо с 0, либо с 344 модулями (559)?

Answer 1

Вы можете использовать функцию select:

subset1 = select((x) -> (x % 43 == 0), v);
subset2 = select((x) -> (x % 43 == 0) && (x % 559 == 0 || x % 559 == 344), v);
vecsort(subset1) == vecsort(subset2)
> 1

Таким образом, гипотеза работает для чисел, сгенерированных вашим кодом PARI / GP.