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