A054611 -id:A054611 - cp's OEIS frontend

A185651 A(n,k) = Sum_{d|n} phi(d)*k^(n/d); square array A(n,k), n>=0, k>=0, read by antidiagonals.

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 6, 3, 0, 0, 4, 12, 12, 4, 0, 0, 5, 20, 33, 24, 5, 0, 0, 6, 30, 72, 96, 40, 6, 0, 0, 7, 42, 135, 280, 255, 84, 7, 0, 0, 8, 56, 228, 660, 1040, 780, 140, 8, 0, 0, 9, 72, 357, 1344, 3145, 4200, 2205, 288, 9, 0

Offset: 0

Alois P. Heinz, Aug 29 2013

Dirichlet convolution of phi(n) and k^n. - Richard L. Ollerton, May 07 2021

			Square array A(n,k) begins:
  0, 0,  0,   0,    0,     0,     0, ...
  0, 1,  2,   3,    4,     5,     6, ...
  0, 2,  6,  12,   20,    30,    42, ...
  0, 3, 12,  33,   72,   135,   228, ...
  0, 4, 24,  96,  280,   660,  1344, ...
  0, 5, 40, 255, 1040,  3145,  7800, ...
  0, 6, 84, 780, 4200, 15810, 46956, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0..10 give A000004, A001477, A053635, A054610, A054611, A054612, A054613, A054614, A054615, A054616, A054617.
Rows n=0..10 give A000004, A001477, A002378, A054602, A054603, A054604, A054605, A054606, A054607, A054608, A054609.
Main diagonal gives A228640.
Cf. A000010, A075195, A258170.

with(numtheory):
A:= (n, k)-> add(phi(d)*k^(n/d), d=divisors(n)):
seq(seq(A(n, d-n), n=0..d), d=0..12);

a[, 0] = a[0, ] = 0; a[n_, k_] := Sum[EulerPhi[d]*k^(n/d), {d, Divisors[n]}]; Table[a[n - k, k], {n, 0, 12}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 06 2013 *)

A(n,k) = Sum_{d|n} phi(d)*k^(n/d).
A(n,k) = Sum_{i=0..min(n,k)} C(k,i) * i! * A258170(n,i). - Alois P. Heinz, May 22 2015
G.f. for column k: Sum_{n>=1} phi(n)*k*x^n/(1-k*x^n) for k >= 0. - Petros Hadjicostas, Nov 06 2017
From Richard L. Ollerton, May 07 2021: (Start)
A(n,k) = Sum_{i=1..n} k^gcd(n,i).
A(n,k) = Sum_{i=1..n} k^(n/gcd(n,i))*phi(gcd(n,i))/phi(n/gcd(n,i)).
A(n,k) = A075195(n,k)*n for n >= 1, k >= 1.  (End)

A001868 Number of n-bead necklaces with 4 colors.

Original entry on oeis.org

1, 4, 10, 24, 70, 208, 700, 2344, 8230, 29144, 104968, 381304, 1398500, 5162224, 19175140, 71582944, 268439590, 1010580544, 3817763740, 14467258264, 54975633976, 209430787824, 799645010860, 3059510616424, 11728124734500, 45035996273872, 173215372864600, 667199944815064

Offset: 0

N. J. A. Sloane

From Richard L. Ollerton, May 07 2021: (Start)
Here, as in A000031, turning over is not allowed.
(1/n) * Dirichlet convolution of phi(n) and 4^n, n>0. (End)

J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 162.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 7.112(a).

T. D. Noe, Table of n, a(n) for n=0..200
Joscha Diehl, Rosa Preiß, and Jeremy Reizenstein, Conjugation, loop and closure invariants of the iterated-integrals signature, arXiv:2412.19670 [math.RA], 2024. See p. 21.
E. N. Gilbert and J. Riordan, Symmetry types of periodic sequences, Illinois J. Math., 5 (1961), 657-665.
Yi Hu, Numerical Transfer Matrix Method of Next-nearest-neighbor Ising Models, Master's Thesis, Duke Univ. (2021).
Yi Hu and Patrick Charbonneau, Numerical transfer matrix study of frustrated next-nearest-neighbor Ising models on square lattices, arXiv:2106.08442 [cond-mat.stat-mech], 2021.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 4
Juhani Karhumäki, S. Puzynina, M. Rao, and M. A. Whiteland, On cardinalities of k-abelian equivalence classes, arXiv preprint arXiv:1605.03319 [math.CO], 2016.
J. Riordan, Letter to N. J. A. Sloane, Jul. 1978
Index entries for sequences related to necklaces

Column 4 of A075195.

Cf. A054611.

A001868 := proc(n) local d,s; if n = 0 then RETURN(1); else s := 0; for d in divisors(n) do s := s+phi(d)*4^(n/d); od; RETURN(s/n); fi; end;

a[n_] := (1/n)*Total[ EulerPhi[#]*4^(n/#) &  /@ Divisors[n]]; a[0] = 1; Table[a[n], {n, 0, 24}] (* Jean-François Alcover, Oct 21 2011 *)
mx=40;CoefficientList[Series[1-Sum[EulerPhi[i] Log[1-4*x^i]/i,{i,1,mx}],{x,0,mx}],x] (* Herbert Kociemba, Nov 01 2016 *)
k=4; Prepend[Table[DivisorSum[n, EulerPhi[#] k^(n/#) &]/n, {n, 1, 30}], 1] (* Robert A. Russell, Sep 21 2018 *)

a(n) = if (n, sumdiv(n, d, eulerphi(d)*4^(n/d))/n, 1); \\ Michel Marcus, Nov 01 2016

a(n) = (1/n)*Sum_{d|n} phi(d)*4^(n/d) = A054611(n)/n, n>0.
G.f.: 1 - Sum_{n>=1} phi(n)*log(1 - 4*x^n)/n. - Herbert Kociemba, Nov 01 2016
a(0) = 1; a(n) = (1/n) * Sum_{k=1..n} 4^gcd(n,k). - Ilya Gutkovskiy, Apr 17 2021
a(0) = 1; a(n) = (1/n)*Sum_{k=1..n} 4^(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 07 2021

A343498 a(n) = Sum_{k=1..n} gcd(k, n)^4.

Original entry on oeis.org

1, 17, 83, 274, 629, 1411, 2407, 4388, 6729, 10693, 14651, 22742, 28573, 40919, 52207, 70216, 83537, 114393, 130339, 172346, 199781, 249067, 279863, 364204, 393145, 485741, 545067, 659518, 707309, 887519, 923551, 1123472, 1216033, 1420129, 1514003, 1843746, 1874197

Offset: 1

Seiichi Manyama, Apr 17 2021

Seiichi Manyama, Table of n, a(n) for n = 1..10000

Column 4 of A343510.

Cf. A000010, A001158 (sigma_3(n)), A018804, A054611, A069097, A332517, A342423, A342433, A343497, A343499, A343514.

A343498:= func< n | (&+[d^4*EulerPhi(Floor(n/d)): d in Divisors(n)]) >;
[A343498(n): n in [1..50]]; // G. C. Greubel, Jun 24 2024

a[n_] := Sum[GCD[k, n]^4, {k, 1, n}]; Array[a, 50] (* Amiram Eldar, Apr 18 2021 *)
f[p_, e_] := p^(e-1)*(p^(3*e+4) - p^(3*e) - p + 1)/(p^3-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 50] (* Amiram Eldar, Nov 22 2022 *)

PARI
```
a(n) = sum(k=1, n, gcd(k, n)^4);
```

a(n) = sumdiv(n, d, eulerphi(n/d)*d^4);

a(n) = sumdiv(n, d, moebius(n/d)*d*sigma(d, 3));

my(N=40, x='x+O('x^N)); Vec(sum(k=1, N, eulerphi(k)*x^k*(1+11*x^k+11*x^(2*k)+x^(3*k))/(1-x^k)^5))

def A343498(n): return sum(k^4*euler_phi(n/k) for k in (1..n) if (k).divides(n))
[A343498(n) for n in range(1,51)] # G. C. Greubel, Jun 24 2024

a(n) = Sum_{d|n} phi(n/d) * d^4.
a(n) = Sum_{d|n} mu(n/d) * d * sigma_3(d).
G.f.: Sum_{k >= 1} phi(k) * x^k * (1 + 11*x^k + 11*x^(2*k) + x^(3*k))/(1 - x^k)^5.
Dirichlet g.f.: zeta(s-1) * zeta(s-4) / zeta(s). - Ilya Gutkovskiy, Apr 18 2021
Sum_{k=1..n} a(k) ~ Pi^4 * n^5 / (450*zeta(5)). - Vaclav Kotesovec, May 20 2021
Multiplicative with a(p^e) = p^(e-1)*(p^(3*e+4) - p^(3*e) - p + 1)/(p^3-1). - Amiram Eldar, Nov 22 2022
a(n) = Sum_{1 <= i, j, k, l <= n} gcd(i, j, k, l, n) = Sum_{d divides n} d * J_4(n/d), where the Jordan totient function J_4(n) = A059377(n). - Peter Bala, Jan 18 2024

A123045 Number of frieze patterns of length n under a certain group (see Pisanski et al. for precise definition).

Original entry on oeis.org

0, 2, 6, 12, 39, 104, 366, 1172, 4179, 14572, 52740, 190652, 700274, 2581112, 9591666, 35791472, 134236179, 505290272, 1908947406, 7233629132, 27488079132, 104715393912, 399823554006, 1529755308212, 5864066561554, 22517998136936, 86607703209516

Offset: 0

N. J. A. Sloane, Nov 11 2006

nonn

Alois P. Heinz, Table of n, a(n) for n = 0..1000
Shinsaku Fujita, alpha-beta Itemized Enumeration of Inositol Derivatives and m-Gonal Homologs by Extending Fujita's Proligand Method, Bull. Chem. Soc. Jpn. 2017, 90, 343-366; doi:10.1246/bcsj.20160369. See Table 8.
T. Pisanski, D. Schattschneider and B. Servatius, Applying Burnside's lemma to a one-dimensional Escher problem, Math. Mag., 79 (2006), 167-180. See F(n).

Cf. A054611, A018215.

The 8 sequences in Table 8 of Fujita (2017) are A053656, A000011, A256216, A256217, A123045, A283846, A283847, A283848.

with(numtheory):
V:=proc(n) local k, t1; t1:=0; for k in divisors(n) do t1 := t1+phi(k)*4^(n/k); od: t1; end; # A054611
H:=n-> if n mod 2 = 0 then (n/2)*4^(n/2); else 0; fi; # this is A018215 interleaved with 0's
A123045:=n-> `if`(n=0,0, (V(n)+H(n))/(2*n));

V[n_] := Module[{t1 = 0}, Do[t1 = t1 + EulerPhi[k] 4^(n/k), {k, Divisors[n]}]; t1];
H[n_] := If[Mod[n, 2] == 0, (n/2) 4^(n/2), 0];
a[n_] := If[n == 0, 0, (V[n] + H[n])/(2n)];
a /@ Range[0, 26] (* Jean-François Alcover, Mar 20 2020, from Maple *)

See Maple program.

A343490 a(n) = Sum_{k=1..n} 4^(gcd(k, n) - 1).

Original entry on oeis.org

1, 5, 18, 70, 260, 1050, 4102, 16460, 65574, 262420, 1048586, 4195500, 16777228, 67112990, 268436040, 1073758360, 4294967312, 17179936830, 68719476754, 274878169880, 1099511636076, 4398047559730, 17592186044438, 70368748407000, 281474976711700, 1125899923619900

Offset: 1

Seiichi Manyama, Apr 17 2021

nonn

Robert Israel, Table of n, a(n) for n = 1..1660

Column 4 of A343489.

Cf. A000010, A054611.

N:= 30: # for a(1)..a(N)
G:= add(numtheory:-phi(k)*x^k/(1-4*x^k),k=1..N):
S:= series(G,x,N+1):
seq(coeff(S,x,j),j=1..N); # Robert Israel, Sep 11 2023

a[n_] := Sum[4^(GCD[k, n] - 1), {k, 1, n}]; Array[a, 26] (* Amiram Eldar, Apr 17 2021 *)

PARI
```
a(n) = sum(k=1, n, 4^(gcd(k, n)-1));
```

a(n) = sumdiv(n, d, eulerphi(n/d)*4^(d-1));

my(N=40, x='x+O('x^N)); Vec(sum(k=1, N, eulerphi(k)*x^k/(1-4*x^k)))

a(n) = Sum_{d|n} phi(n/d)*4^(d - 1) = A054611(n)/4.

G.f.: Sum_{k>=1} phi(k) * x^k / (1 - 4*x^k).

cp's OEIS Frontend

A185651 A(n,k) = Sum_{d|n} phi(d)*k^(n/d); square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A001868 Number of n-bead necklaces with 4 colors.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A343498 a(n) = Sum_{k=1..n} gcd(k, n)^4.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

PARI

PARI

SageMath

Formula

A123045 Number of frieze patterns of length n under a certain group (see Pisanski et al. for precise definition).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A343490 a(n) = Sum_{k=1..n} 4^(gcd(k, n) - 1).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

Formula