A154302 Primes of the form (1+2+...+m)/87 = A000217(m)/87.
5, 19, 43, 173
Offset: 1
Programs
-
Mathematica
lst={};s=0;Do[s+=n/87;If[Floor[s]==s,If[PrimeQ[s],AppendTo[lst,s]]],{n,0,4*9!}];lst
-
PARI
d=87*2;for(m=1,999,(m^2+m)%d==0&isprime((m^2+m)/d)&print1(m",")) \\ print the m-values(!) - use A154304(87) to get A154302 as a vector. \\ - M. F. Hasler, Jan 06 2013
Extensions
Edited by M. F. Hasler, Jan 06 2013
Comments