OFFSET
0,2
LINKS
Muniru A Asiru, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = A053447(2*n) for n >= 0. - Jianing Song, Oct 03 2022
MAPLE
a := n -> `if`(n=0, 1, numtheory:-order(4, 4*n+1)): seq(a(n), n = 0..68);
MATHEMATICA
Table[MultiplicativeOrder[4, 4*n + 1], {n, 0, 70}] (* Arkadiusz Wesolowski, Nov 27 2012 *)
PROG
(Sage)
def A003573(n):
s, m, N = 0, 1, 4*n + 1
while True:
k = N + m
v = valuation(k, 4)
s += v
m = k // 4^v
if m == 1: break
return s
print([A003573(n) for n in (0..70)]) # Peter Luschny, Oct 07 2017
(GAP) List([0..70], n->OrderMod(4, 4*n+1)); # Muniru A Asiru, Feb 16 2019
(PARI) a(n) = znorder(Mod(4, 4*n+1)); \\ Michel Marcus, Feb 16 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
a(0) = 1 added by Peter Luschny, Oct 07 2017
STATUS
approved