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 コードで生成された数値が 43 の倍数である場合、それらが 0 または 344 mod (559) に一致するかどうかを確認するにはどうすればよいでしょうか?