|
|
A061599
|
|
Primes p such that the greatest prime divisor of p-1 is 5.
|
|
3
|
|
|
11, 31, 41, 61, 101, 151, 181, 241, 251, 271, 401, 541, 601, 641, 751, 811, 1201, 1601, 1621, 1801, 2161, 2251, 3001, 4001, 4051, 4801, 4861, 6481, 7681, 8101, 8641, 9001, 9601, 9721, 11251, 14401, 15361, 16001, 19441, 21601, 21871, 22501, 23041, 24001
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Prime numbers n for which cos(2Pi/n) is an algebraic number of 5th degree. - Artur Jasinski, Dec 13 2006
|
|
LINKS
|
|
|
FORMULA
|
Primes of the form 2^a*3^b*5^c + 1 with a and c > 0.
|
|
MATHEMATICA
|
Do[If[Take[FactorInteger[EulerPhi[2n + 1]][[ -1]], 1] == {5} && PrimeQ[2n + 1], Print[2n + 1]], {n, 1, 10000}] (* Artur Jasinski, Dec 13 2006 *)
|
|
PROG
|
(PARI) { default(primelimit, 167772161); n=0; forprime (p=3, 167772161, f=factor(p - 1)~; if (f[1, length(f)]==5, write("b061599.txt", n++, " ", p)) ) } \\ Harry J. Smith, Jul 25 2009
(PARI) list(lim)=my(v=List(), s, t); lim\=1; lim--; for(i=1, logint(lim\2, 5), t=2*5^i; for(j=0, logint(lim\t, 3), s=t*3^j; while(s<=lim, if(isprime(s+1), listput(v, s+1)); s<<=1))); Set(v) \\ Charles R Greathouse IV, Oct 29 2018
|
|
CROSSREFS
|
The 3rd in a family of sequences after A019434(=Fermat-primes) and A058383.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|