|
|
|
|
1, 2, 2, 4, 2, 4, 2, 6, 3, 4, 2, 8, 2, 4, 4, 8, 2, 6, 2, 8, 4, 4, 2, 12, 3, 4, 4, 8, 2, 8, 2, 10, 4, 4, 4, 12, 2, 4, 4, 12, 2, 8, 2, 8, 6, 4, 2, 16, 3, 6, 4, 8, 2, 8, 4, 12, 4, 4, 2, 16, 2, 4, 6, 12, 4, 8, 2, 8, 4, 8, 2, 18, 2, 4, 6, 8, 4, 8, 2, 16, 5, 4, 2, 16, 4, 4, 4, 12, 2, 12, 4, 8, 4, 4, 4, 20, 2, 6, 6, 12, 2, 8, 2, 12, 8
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Related to Menon's identity. See Conclusions and further work section of the arXiv file linked.
|
|
LINKS
|
|
|
MATHEMATICA
|
a239611[n_] := Sum[If[GCD[x^2 + y^2, n] == 1, GCD[x^2 + y^2 - 1, n], 0], {x, 1, n}, {y, 1, n}];
a079458[n_] := Product[{p, e} = pe; Which[p==2, 2^(2e-1), Mod[p, 4]==3, (p^2-1)p^(2e-2), Mod[p, 4]==1, (p-1)^2 p^(2e-2)], {pe, FactorInteger[n]}];
a[1] = 1; a[n_] := a239611[n]/a079458[n];
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|