OFFSET
1,1
COMMENTS
a(n) = 2*A065813(n) + 1, n > 1.
LINKS
MATHEMATICA
Do[p = Prime[n]; k = 1; While[ !PrimeQ[ (p^Prime[k] - 1)/(p - 1)], k++ ]; Print[ Prime[k]], {n, 1, 56} ]
PROG
(PARI) { allocatemem(932245000); for (n=1, 100, p=prime(n); q=2; while (!isprime((p^q - 1)/(p - 1)), q=nextprime(q + 1)); write("b065854.txt", n, " ", q) ) } \\ Harry J. Smith, Nov 01 2009
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
Vladeta Jovovic, Nov 26 2001
STATUS
approved