OFFSET
1,1
COMMENTS
A133870(n) : primes of form 32n+1.
EXAMPLE
a(5) = 10 because 10^16+1 = 10000000000000001 = 353 * 449 * 641 * 1409 * 69857 with A133870(5) = 449.
MATHEMATICA
aa = {}; Do[p = Prime[n]; If[Mod[p, 32] == 1, k = 1; While[ ! Mod[k^16 + 1, p] == 0, k++ ]; AppendTo[aa, k]], {n, 300}]; aa
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Oct 22 2012
STATUS
approved