OFFSET
0,3
COMMENTS
a(n) = 1 <=> n is in A000225 <=> n = 2^k - 1 with k = 0, 1, 2, ... - M. F. Hasler, Sep 20 2017
LINKS
Robert Israel, Table of n, a(n) for n = 0..1672
MAPLE
A002326 := proc(n) if n = 0 then 1 ; else numtheory[order](2, 2*n+1) ; end if ; end proc:
seq(A165781(n), n=0..60) ; # R. J. Mathar, Nov 16 2009
MATHEMATICA
a[n_] := (2^MultiplicativeOrder[2, 2n+1]-1)/(2n+1);
a /@ Range[0, 40] (* Jean-François Alcover, Jun 04 2020 *)
PROG
(PARI) a(n)=(2^znorder(Mod(2, n=2*n+1))-1)/n \\ M. F. Hasler, Sep 20 2017
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Ctibor O. Zizka, Sep 26 2009
EXTENSIONS
Sign in definition and offset corrected by R. J. Mathar, Nov 16 2009
STATUS
approved