A372881 a(n) = Sum_{j=1..n} Sum_{k=1..n} gcd(j*k,n).
1, 7, 19, 44, 61, 133, 127, 256, 297, 427, 331, 836, 469, 889, 1159, 1408, 817, 2079, 1027, 2684, 2413, 2317, 1519, 4864, 2825, 3283, 4131, 5588, 2437, 8113, 2791, 7424, 6289, 5719, 7747, 13068, 3997, 7189, 8911, 15616, 4921, 16891, 5419, 14564, 18117, 10633, 6487, 26752
Offset: 1
Links
- Seiichi Manyama, Table of n, a(n) for n = 1..1000
Programs
-
PARI
a(n) = sum(j=1, n, sum(k=1, n, gcd(j*k, n)));