A295574 -id:A295574 - cp's OEIS frontend

A295575 a(n) = Sum_{1 <= j <= n/2, gcd(j,n)=1} j^3.

0, 1, 1, 1, 9, 1, 36, 28, 73, 28, 225, 126, 441, 153, 416, 496, 1296, 469, 2025, 1100, 1710, 1225, 4356, 1800, 4959, 2556, 5581, 4410, 11025, 3872, 14400, 8128, 11090, 8128, 15822, 8910, 29241, 13041, 21996, 16400, 44100, 15426, 53361, 27830, 33716, 29161, 76176, 27936, 77652, 37828

Offset: 1

N. J. A. Sloane, Dec 08 2017

nonn

If p is an odd prime, a(p) = (n^2-1)^2/64. - Robert Israel, Dec 10 2017

Seiichi Manyama, Table of n, a(n) for n = 1..10000
John D. Baum, A Number-Theoretic Sum, Mathematics Magazine 55.2 (1982): 111-113.

In the Baum (1982) paper, S_1, S_2, S_3, S_4 are A023896, A053818, A053819, A053820, and S'_1, S'_2, S'_3, S'_4 are A066840, A295574, A295575, A295576.

f:= n -> add(t^3, t = select(t->igcd(t,n)=1, [$1..n/2])) :
map(f, [$1..100]); # Robert Israel, Dec 10 2017

f[n_] := Plus @@ (Select[Range[n/2], GCD[#, n] == 1 &]^3); Array[f, 50] (* Robert G. Wilson v, Dec 10 2017 *)

a(n) = sum(j=1, n\2, (gcd(j, n)==1)*j^3); \\ Michel Marcus, Dec 10 2017

A295576 a(n) = Sum_{1 <= j <= n/2, gcd(j,n)=1} j^4.

Original entry on oeis.org

0, 1, 1, 1, 17, 1, 98, 82, 273, 82, 979, 626, 2275, 707, 2674, 3108, 8772, 3027, 15333, 9044, 14994, 9669, 39974, 17668, 50085, 24310, 60597, 50470, 127687, 45604, 178312, 103496, 149908, 103496, 225302, 129750, 432345, 187017, 349830, 266088, 722666

Offset: 1

N. J. A. Sloane, Dec 08 2017

nonn

If p is an odd prime, a(p) = p*(p^2-1)*(3*p^2-7)/480. - Robert Israel, Dec 10 2017

Seiichi Manyama, Table of n, a(n) for n = 1..10000
John D. Baum, A Number-Theoretic Sum, Mathematics Magazine 55.2 (1982): 111-113.

In the Baum (1982) paper, S_1, S_2, S_3, S_4 are A023896, A053818, A053819, A053820, and S'_1, S'_2, S'_3, S'_4 are A066840, A295574, A295575, A295576.

f:= n -> add(t^4, t = select(t->igcd(t,n)=1, [$1..n/2])):
map(f, [$1..100]); # Robert Israel, Dec 10 2017

f[n_] := Plus @@ (Select[Range[n/2], GCD[#, n] == 1 &]^4); Array[f, 41] (* Robert G. Wilson v, Dec 10 2017 *)

a(n) = sum(j=1, n\2, (gcd(j, n)==1)*j^4); \\ Michel Marcus, Dec 10 2017

cp's OEIS Frontend

A295575 a(n) = Sum_{1 <= j <= n/2, gcd(j,n)=1} j^3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI