A154298 Primes of the form (1+...+m)/33 = A000217(m)/33, for some m.
2, 7, 17, 67
Offset: 1
Programs
-
Mathematica
lst={};s=0;Do[s+=n/33;If[Floor[s]==s,If[PrimeQ[s],AppendTo[lst,s]]],{n,0,9!}];lst Select[Accumulate[Range[200]]/33,PrimeQ] (* Harvey P. Dale, Aug 11 2025 *)
-
PARI
select(x->denominator(x)==1 & isprime(x), vector(66, m, m^2+m)/66) \\ - M. F. Hasler, Dec 31 2012
Extensions
Edited by M. F. Hasler, Dec 31 2012
Comments