(平价)
A003415(n) =如果(n<=1,0,my(f=因子(n));n*和(i=1,#f~,f[i,2]/f[i,1]);
A327860(n) ={my(m=1,i=0,s=0,pr=1,nextpr);而((n>0),i=i+1;nextpr=prime(i)*pr;if((n%nextpr),my(e=((n%nextpr)/pr));m*=(素数(i)^e);s+=(e/prime(i));n-=(n%nextr));pr=nextpr);(s*m);};
isA328240(n)=isprime(A003415(A327860(n) );
|