OFFSET
1,2
COMMENTS
a(n) = m where A080076(n) = k*2^m + 1, k odd.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
N:= 10^6: # for all Proth primes <= N
Proth:= sort(convert(select(isprime, {seq(seq(k*2^m+1, k = 1 .. min(2^m, (N-1)/2^m), 2), m=1..ilog2(N-1))}), list)):
map(t -> padic:-ordp(t-1, 2), Proth);
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Feb 02 2016
STATUS
approved