cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A177976 Square array T(n,k) read by antidiagonals up. Cumulative column sums of A177975.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 6, 8, 4, 1, 1, 10, 15, 13, 5, 1, 1, 12, 29, 29, 19, 6, 1, 1, 18, 42, 63, 49, 26, 7, 1, 1, 22, 69, 106, 118, 76, 34, 8, 1, 1, 28, 95, 189, 225, 201, 111, 43, 9, 1, 1, 32, 134, 289, 434, 427, 320, 155, 53, 10, 1, 1, 42, 172, 444, 729, 888, 748, 484, 209, 64, 11, 1
Offset: 1

Views

Author

Mats Granvik, May 16 2010

Keywords

Comments

Each row is described by both a binomial expression and a closed form polynomial. The closed form polynomials given in A177977 extends this table to the left. For example the 0th column is A002321 and the -1st column is A092149.
Also number of ordered k-tuples of integers from [ 1..n ] with no global factor. - Seiichi Manyama, Jun 12 2021

Examples

			Table begins:
  1..1...1....1.....1.....1......1......1.......1.......1.......1
  1..2...3....4.....5.....6......7......8.......9......10......11
  1..4...8...13....19....26.....34.....43......53......64......76
  1..6..15...29....49....76....111....155.....209.....274.....351
  1.10..29...63...118...201....320....484.....703.....988....1351
  1.12..42..106...225...427....748...1233....1937....2926....4278
  1.18..69..189...434...888...1671...2948....4939....7930...12285
  1.22..95..289...729..1624...3303...6260...11209...19150...31447
  1.28.134..444..1209..2890...6278..12659...24034...43405...75139
  1.32.172..626..1850..4761..11067..23762...47841...91301..166506
  1.42.237..911..2850..7763..19074..43209...91598..183678..351261
  1.46.287.1203..4059.11829..30911..74129..165737..349426..700699
  1.58.377.1657..5878.18016..49474.124516..291706..643355.1347344
  1.64.452.2130..8044.26117..75676.200313..492185.1135761.2483392
  1.72.552.2766.11020.37599.114199.316228..811416.1952182.4443582
  1.80.652.3462.14566.52311.166747.483340.1295295.3248246.7692894

Links

Crossrefs

Column k=1..5 gives A000012, A002088, A015631, A015634, A015650.
Cf. A177975, A177977, A344527, A345229.

Programs

  • PARI
    T(n, k) = sum(j=1, n, sumdiv(j, d, moebius(j/d)*binomial(d+k-2, d-1))); \\ Seiichi Manyama, Jun 12 2021
  • PARI
    T(n, k) = binomial(n+k-1, k)-sum(j=2, n, T(n\j, k)); \\ Seiichi Manyama, Jun 12 2021

Formula

From Seiichi Manyama, Jun 12 2021: (Start)
G.f. of column k: (1/(1 - x)) * Sum_{j>=1} mu(j) * x^j/(1 - x^j)^k.
T(n,k) = Sum_{j=1..n} Sum_{d|j} mu(j/d) * binomial(d+k-2,d-1).
T(n,k) = binomial(n+k-1,k) - Sum_{j=2..n} T(floor(n/j),k). (End)

A344521 a(n) = Sum_{1 <= i <= j <= k <= n} gcd(i,j,k).

Original entry on oeis.org

1, 5, 13, 28, 47, 82, 116, 172, 235, 321, 397, 538, 641, 798, 980, 1192, 1361, 1655, 1863, 2218, 2553, 2912, 3210, 3766, 4171, 4661, 5183, 5840, 6303, 7168, 7694, 8510, 9283, 10095, 10951, 12190, 12929, 13932, 14990, 16414, 17315, 18925, 19913, 21438, 23055, 24500, 25674, 27862
Offset: 1

Views

Author

Seiichi Manyama, May 22 2021

Keywords

Links

Crossrefs

Column k=3 of A345229.
Partial sums of A309322.
Cf. A272718, A344132, A344522, A344992.

Programs

  • Mathematica
    a[n_] := Sum[Sum[Sum[GCD[i, j, k], {i, 1, j}], {j, 1, k}], {k, 1, n}]; Array[a, 50] (* Amiram Eldar, May 25 2021 *)
nmax = 100; Rest[CoefficientList[Series[1/(1 - x)*Sum[EulerPhi[k]*x^k/(1 - x^k)^3, {k, 1, nmax}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, Jun 05 2021 *)
Accumulate[Table[Sum[EulerPhi[n/d] * d*(d+1)/2, {d, Divisors[n]}], {n, 1, 100}]] (* Vaclav Kotesovec, Jun 05 2021 *)
  • PARI
    a(n) = sum(i=1, n, sum(j=i, n, sum(k=j, n, gcd([i, j, k]))));

Formula

From Vaclav Kotesovec, Jun 05 2021: (Start)
a(n) ~ Pi^2 * n^3 / (36*zeta(3)).
G.f.: 1/(1-x) * Sum_{k>=1} phi(k) * x^k/(1 - x^k)^3, where phi is the Euler totient function (A000010).
a(n) = Sum_{k=1..n} Sum_{d|k} phi(k/d) * d*(d+1)/2. (End)
a(n) = Sum_{k=1..n} phi(k) * binomial(floor(n/k)+2,3). - Seiichi Manyama, Sep 13 2024

A344992 a(n) = Sum_{1 <= i <= j <= k <= m <= n} gcd(i,j,k,m).

Original entry on oeis.org

1, 6, 18, 44, 83, 159, 249, 401, 592, 867, 1163, 1655, 2122, 2796, 3594, 4594, 5579, 7046, 8394, 10328, 12339, 14699, 17021, 20441, 23526, 27317, 31379, 36323, 40846, 47300, 52786, 59954, 67191, 75380, 83720, 94662, 103837, 115137, 126851, 141059, 153440
Offset: 1

Views

Author

Vaclav Kotesovec, Jun 05 2021

Keywords

Comments

In general, if g.f.: 1/(1-x) * Sum_{j>=1} phi(j) * x^j/(1 - x^j)^k, where k > 2 and phi is the Euler totient function (A000010), then a(n) ~ zeta(k-1) * n^k / (k! * zeta(k)).

Links

Crossrefs

Column k=4 of A345229.
Partial sums of A309323.
Cf. A000010, A059358, A272718, A309322, A344521.

Programs

  • Mathematica
    Table[Sum[Sum[Sum[Sum[GCD[i, j, k, m], {i, 1, j}], {j, 1, k}], {k, 1, m}], {m, 1, n}], {n, 1, 100}]
nmax = 100; Rest[CoefficientList[Series[1/(1-x) * Sum[EulerPhi[k]*x^k/(1 - x^k)^4, {k, 1, nmax}], {x, 0, nmax}], x]]
Accumulate[Table[Sum[EulerPhi[n/d] * d*(d+1)*(d+2)/6, {d, Divisors[n]}], {n, 1, 100}]] (* faster *)
  • PARI
    a(n) = sum(i=1, n, sum(j=i, n, sum(k=j, n, sum(m=k, n, gcd([i, j, k, m]))))); \\ Michel Marcus, Jun 06 2021

Formula

G.f.: 1/(1-x) * Sum_{k>=1} phi(k) * x^k/(1 - x^k)^4, where phi is the Euler totient function (A000010).
a(n) = Sum_{k=1..n} Sum_{d|k} phi(k/d) * d*(d+1)*(d+2)/6.
a(n) ~ 15 * zeta(3) * n^4 / (4*Pi^4).
a(n) = Sum_{k=1..n} phi(k) * binomial(floor(n/k)+3,4). - Seiichi Manyama, Sep 13 2024

A345230 a(n) = Sum_{1 <= x_1 <= x_2 <= ... <= x_n <= n} gcd(x_1, x_2, ..., x_n).

Original entry on oeis.org

0, 1, 4, 13, 44, 140, 512, 1782, 6652, 24682, 93599, 354341, 1359470, 5210328, 20098886, 77621774, 300797854, 1167164438, 4539201401, 17674941735, 68933414989, 269143872226, 1052114789548, 4116808923486, 16124224585644, 63205911146740, 247961982954952
Offset: 0

Views

Author

Seiichi Manyama, Jun 11 2021

Keywords

Links

Crossrefs

Main diagonal of A345229.
Cf. A343517, A343553, A344525.

Programs

  • Maple
    a:= n-> coeff(series((1/(1-x))* add(numtheory[phi](k)
         *x^k/(1-x^k)^n, k=1..n), x, n+1), x, n):
seq(a(n), n=0..26);  # Alois P. Heinz, Jun 11 2021
  • Mathematica
    a[n_] := Sum[DivisorSum[k, EulerPhi[k/#] * Binomial[n + # - 2, n - 1] &], {k, 1, n}]; Array[a, 30, 0] (* Amiram Eldar, Jun 11 2021 *)
  • PARI
    a(n) = sum(k=1, n, sumdiv(k, d, eulerphi(k/d)*binomial(d+n-2, n-1)));
  • PARI
    a(n) = sum(k=1, n, eulerphi(k)*binomial(n\k+n-1, n)); \\ Seiichi Manyama, Sep 13 2024

Formula

a(n) = Sum_{k=1..n} Sum_{d|k} phi(k/d) * binomial(d+n-2, n-1).
a(n) = [x^n] (1/(1 - x)) * Sum_{k >= 1} phi(k) * x^k/(1 - x^k)^n.
a(n) ~ 2^(2*n-1) / sqrt(Pi*n). - Vaclav Kotesovec, Jun 11 2021
a(n) = Sum_{k=1..n} phi(k) * binomial(floor(n/k)+n-1,n). - Seiichi Manyama, Sep 13 2024
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.