|
|
PROG
| (PARI) diffpowers(n, m) =
{
local(a, c=0, c2=0, j, k, y);
a=vector(floor(n^2/log(n^2)));
for(j=1, n,
for(k=1, n,
y=j^m-k^(m-1);
if(ispseudoprime(y),
c++;
\\ print(j", "k", "y);
a[c]=y;
);
);
);
a=vecsort(a);
for(j=2, length(a),
if(a[j]!=a[j-1]&&a[j]!=0,
c2++;
print1(a[j]", ");
if(c2>100, break);
);
);
}
|
