OFFSET
1,1
COMMENTS
For a(n) = m the prime p = A123571(n) is the least prime dividing m^(m+1)+1, where m = (p-2n-1)/(2n).
FORMULA
a(n) = (A123571(n) - 2n - 1)/(2n).
PROG
(PARI) A123571(n)={ local(m) ; forprime(p=1, 500000, if( (p-2*n-1) % (2*n) ==0, m=(p-2*n-1)/(2*n) ; if( (m^(m+1)+1) % p == 0, return(p) ; ) ; ) ; ) ; return(-1) ; }
{ for(n=1, 100, print1(A177714(n)", ") ; ) ; } \\ R. J. Mathar, May 28 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, May 11 2010
EXTENSIONS
More terms from R. J. Mathar, May 28 2010
STATUS
approved