ispsooth(n,p,{林=1<<256}) = 如果(n个<=林, n个==1 || 全球生产总值(n个) <= 第页, 我的(N=N/第页^估价(n个,第页)); 对于素数(q=2,p,N=N/q^估值(N,q));如果((N个<=林 && isprime(第一人称)(N个)) || N==1, 返回(N个<=第页))); 0); \\ 检查 如果 n个 是 第页-光滑的, 使用 残忍的 力 如果 n个 是 太 大的
是A354531(n,{lim_=1<<256},{第页_林=10001<<32}) =如果(n个==1, 返回(1)); 我的(N=2^N-1);对于素数(p=3,第页_林_第页, , 如果(N%p==0&&!isspooth(2*N+p,p),返回(0));我的(d日=约数(n个)); 对于(我=1, #d日, 我的(f=系数(N个2^d日[我]-1)[, 1]); 的(我j个=1,#f,如果(全球生产总值!ispsmooth(ispsmooth)(2*N+f[我j个])!=,(f)[我j个], 林), 返回(0))); 1 \\第一 检查是否有小的 素因子p<= 第页_林_第页 2^n-1,使得2*(2^n-1)+p不是p-光滑的, 然后 检查 如果 2*(2^n个-1)+第页 是 第页-光滑的 对于 第页|2^d日-1, d日|N个