选择(n->numtheory:-issqrfree(二项式(2*n-1,n)),[1..2000])#罗伯特·伊斯雷尔2015年12月9日
N: =300:#查找所有术语<=2^N
进位:=进程(n,m,p)
#在基数p中添加n+m时的进位数。
局部A、B、C、j、nc、t;
A: =换算(m,基数,p);
B: =换算(n,基数,p);
C: =0;nc:=0;
如果nops(A)<nops(B),则A=[op(A),0$(nops(B)-nops(A))]
elif nops(A)>nops(B)然后B:=[操作(B),0$(nops(A)-nops(B))]
fi;
对于j从1到nops(A)do
t: =C+A[j]+B[j];
如果t>=p,则
nc:=nc+1;
C: =1;
其他的
C: =0
fi(菲涅耳)
日期:
数控;
结束进程:
Cands:={seq(2^j,j=0..N),seq(seq(2 ^j+2^k,k=0..j-1),j=1..N-1)}:
我从2点到10点
Cands:=选择(n->进位(n-1,n,ithprime(i))<=1,Cands)
日期:
选择(n->numtheory:-issqrfree(二项式(2*n-1,n)),Cands)#罗伯特·伊斯雷尔2015年12月10日