A154301 Primes of the form (1+2+...+m)/75 = A000217(m)/75.
17, 37, 149, 151
Offset: 1
Programs
-
Mathematica
lst={};s=0;Do[s+=n/75;If[Floor[s]==s,If[PrimeQ[s],AppendTo[lst,s]]],{n,0,3*9!}];lst
-
PARI
d=75*2;for(m=1,999,(m^2+m)%d==0&isprime((m^2+m)/d)&print1(m",")) \\ print the m-values(!) - use A154304(75) to get A154301 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