(平价)
A_simple(n)=局部(m=2);n*=2;直到(numdiv(m)==n&numdiv)(m+n)==n,m++);米
A_try_pair(p,q,n,极限)=
{
/*A_prime()的助手*/
/*寻找0 mod p^(n-1)和-n*2 mod q^(n-1)的解*/
本地(m=中文(Mod(0,p^(n-1)),Mod(-n*2,q^(n-1)));
forstep(x=升力(m),极限,分量(m,1),
如果(isprime(x\p^(n-1))和isprime;
限制
}
A _ try _ above _ below(m,n)=
{
/*A_prime()的助手*/
/*函数假定numdiv(m)==n*2*/
如果(numdiv(m-n*2)==n*2,极限=m-n*3,
如果(numdiv(m+n*2)==n*2,极限=m,
0))
}
A_素数(n,极限,pairmax=30)=
{
如果(n%2==0||!isprime(n),错误(“仅适用于奇数素数”);
if(default(primelimit)<limit\nextprime(pairmax+1)^(n-1),
默认值(primelimit,limit\nextprime(pairmax+1)^(n-1));
);
/*numdiv==n*2的Even是{2^(n*2-1)}u{2*p^(n-1)}u}2^*/
/*潜在的解决方案必须来自不同的集合*/
/*尝试前两组的上下*/
A _ try _ above _ below(2^(n*2-1),n);
对于素数(p=3,(极限\2)^(1/(n-1)),
如果(A_try_above_below(2*p^(n-1),n),break));
/*numdiv==n*2的奇数是{p^(n*2-1)}u{p^(n-1)*q}*/
/*尝试a(n)和a(n”)+n*2是(小素数)^(n-1)*(大素数)*/
对于素数(p=3,pairmax,对于素数,
如果(p!=q,limit=A_try_pair(p,q,n,limit));
/*尝试使用numdiv==n*2的所有其他奇数*/
对于素数(p=pairmax+1,(极限\3)^(1/(n-1)),
对于素数(q=3,极限\p^(n-1),
如果(p!=q&A_try_above_blower(p^(n-1)*q,n),break));
对于素数(p=3,极限^(1/21),
如果(A_try_abover_below(p^21,n),break));
限制