|
%I #23 Oct 28 2018 03:30:25
%S 12,88,207
%N Values of m associated with A082813.
%C No more terms below 822. - _W. Edwin Clark_, May 24 2003
%C No more terms below 4000. - _Hugo Pfoertner_, May 28 2003
%C No more terms < 25500. - _Martin Renner_, Jan 31 2013
%H Hugo Pfoertner, <a href="http://www.randomwalk.de/sequences/a082813.zip">Search for further terms</a>
%p A:=[]:
%p for n from 4 to 210 do
%p if not(isprime(n)) then
%p F:=ifactors(n)[2];
%p if nops(F)>1 then
%p s:=sum(F['i'][2]*F['i'][1]^n,'i'=1..nops(F));
%p if isprime(s) then A:=[op(A),n]; print(n): fi;
%p fi:
%p fi:
%p od: # _Martin Renner_, Jan 31 2013
%o (PARI) for(m=1,+oo,f=factor(m);ispseudoprime(sum(i=1,matsize(f)[1],f[i,1]^m*f[i,2])) && print1(m,", ")) \\ _Jeppe Stig Nielsen_, Oct 27 2018
%Y Cf. A082813.
%K nonn,bref,more,hard
%O 1,1
%A _N. J. A. Sloane_, May 24 2003
%E a(2) and a(3) from _Joshua Zucker_, May 24 2003
| 