OFFSET
1,1
COMMENTS
The primes in this sequence have the property that with the exception of 7 they are congruent to 1 mod 18 and with the exception of 7, 19, 73 are congruent to 1 mod 54.
MATHEMATICA
a = {}; Do[k = ((2^EulerPhi[3^(w + 1)] - 1)/3^(w + 1))/((2^EulerPhi[3^w] - 1)/3^w); Do[If[Mod[k, Prime[n]] == 0, AppendTo[a, Prime[n]]; Print[Prime[n]]], {n, PrimePi[2], PrimePi[379081]}], {w, 1, 20}]; Union[a] (*Artur Jasinski*)
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
Artur Jasinski, Nov 19 2008
EXTENSIONS
Edited by N. J. A. Sloane, Nov 26 2008
STATUS
approved