（PARI）为A292544（n）=Mod（2，n）^ eulerphi（n）==eulerphi（n）；
isok（n）=isA292544（n）&&isA29254（欧拉法（n））；
（PARI）{ZK（m）=my（z，k）；z=znorder（Mod（2，m））；k=znlog（eulerphi（m），Mod（2,m）））；if（type（k）！=“t_INT”，return（））；[z，k]；}
{获得功率2（m）=我的（m2，t，zk，zk2，r）；
m2=欧拉比（m）；
t=估价（m2，2）；
平方米=2^t；
如果（m2==1，返回（0））；
zk=zk（m）；
zk2=ZK（m2）；
如果（！zk||！zk2，return（））；
r=[zk[1]，zk2]，zk[2]-t-zk2[2]+1]；\\求解r[1]*i=r[2]*j+r[3]
r/=含量（r）；
如果（gcd（r[1]，r[2]）>1，return（））；
（（r[2]*升力（Mod（-r[3]/r[2]，r[1]））+r[3]）/r[1]+r[2]*x）*zk[1]-zk[2]+1；}\\getpower2（m）返回带x参数的z*i-k+1（参见公式部分），即getpower2*（13）返回12*x+7，即13*2^（12*x%7）是所有x>=0的项。
