A059889 a(n) = |{m : multiplicative order of 7 mod m=n}|.
4, 6, 8, 26, 4, 42, 12, 48, 52, 66, 12, 778, 4, 138, 80, 300, 12, 528, 12, 1430, 72, 138, 28, 15216, 24, 66, 1216, 966, 28, 3630, 28, 1344, 360, 58, 108, 16988, 28, 138, 176, 12752, 28, 7398, 12, 4422, 1900, 122, 12, 131028, 240, 536, 744, 1046, 28, 23744, 44
Offset: 1
Keywords
Links
- Max Alekseyev, Table of n, a(n) for n = 1..388
Crossrefs
Programs
-
Maple
with(numtheory): a:= n-> add(mobius(n/d)*tau(7^d-1), d=divisors(n)): seq(a(n), n=1..40); # Alois P. Heinz, Oct 12 2012
-
Mathematica
a[n_] := DivisorSum[n, MoebiusMu[n/#] * DivisorSigma[0, 7^#-1] &]; Array[a, 60] (* Amiram Eldar, Jan 25 2025 *)
-
PARI
a(n) = sumdiv(n, d, moebius(n/d) * numdiv(7^d-1)); \\ Amiram Eldar, Jan 25 2025
Comments