如何修改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)

这是一个Pari代码。

现在，我如何修改这个程序，看看它是否真的，如果用Pari代码生成的数字是43的倍数，那么它们是否与0或344 mod (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代码生成的数字。