A154300 Primes of the form (1+2+...+m)/57 = A000217(m)/57.
3, 13, 29, 113
Offset: 1
Programs
-
Mathematica
lst={};s=0;Do[s+=n/57;If[Floor[s]==s,If[PrimeQ[s],AppendTo[lst,s]]],{n,0,2*9!}];lst Select[Accumulate[Range[1000]]/57,PrimeQ] (* Harvey P. Dale, Jun 24 2015 *)
-
PARI
d=57*2;for(m=1,999,(m^2+m)%d==0&isprime((m^2+m)/d)&print1(m",")) \\ print the m-values(!) - use A154304(57) to get A154300 as a vector. \\ - M. F. Hasler, Jan 06 2013
Extensions
Keywords fini,full added by R. J. Mathar, Aug 15 2010
Edited by M. F. Hasler, Jan 06 2013
Comments