OFFSET
1,2
COMMENTS
a(n) <= (p_n-1)/2, where p_n is the n-th odd prime (A065091).
MATHEMATICA
terms = 100;
a137576[n_] := Module[{t}, (t = MultiplicativeOrder[2, 2 n + 1])* DivisorSum[2 n + 1, EulerPhi[#]/MultiplicativeOrder[2, #] &] - t + 1];
A137576 = Array[a137576, 3 terms];
FirstPosition[A137576, #][[1]]& /@ Prime[Range[2, terms+1]] (* Jean-François Alcover, Jan 11 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, May 10 2008
EXTENSIONS
Corrected and extended by Ray Chandler, May 19 2008
STATUS
approved