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.

Previous Showing 11-20 of 23 results. Next

A261050 Expansion of Product_{k>=1} (1+x^k)^(Fibonacci(k)).

Original entry on oeis.org

1, 1, 1, 3, 5, 10, 19, 36, 67, 127, 236, 438, 811, 1496, 2750, 5046, 9224, 16827, 30630, 55623, 100803, 182342, 329205, 593326, 1067591, 1917885, 3440207, 6162004, 11021921, 19688757, 35126020, 62590629, 111398910, 198044551, 351700332, 623918086, 1105715149
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 08 2015

Keywords

Crossrefs

Programs

  • Maple
    f:= proc(n) option remember; (<<1|1>, <1|0>>^n)[1, 2] end:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
           add(binomial(f(i), j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..50);  # Alois P. Heinz, Aug 08 2015
  • Mathematica
    nmax=40; CoefficientList[Series[Product[(1+x^k)^Fibonacci[k],{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ phi^n / (2 * sqrt(Pi) * 5^(1/8) * n^(3/4)) * exp(-1/10 + 2*5^(-1/4)*sqrt(n) + s), where s = Sum_{k>=2} (-1)^(k+1) * phi^k / ((phi^(2*k) - phi^k - 1)*k) = -0.3237251774053525012502809827680337358578568068831886835557918847... and phi = A001622 = (1+sqrt(5))/2 is the golden ratio.
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k/(k*(1 - x^k - x^(2*k)))). - Ilya Gutkovskiy, May 29 2018

A258342 Expansion of Product_{k>=1} (1+x^k)^(k*(k+1)*(k+2)).

Original entry on oeis.org

1, 6, 39, 224, 1131, 5412, 24411, 105078, 435048, 1740312, 6755877, 25533330, 94205738, 340064322, 1203313782, 4180514846, 14279610417, 48013553310, 159086287869, 519912616614, 1677331973910, 5345927500226, 16843574682291, 52494817082952, 161923200857711
Offset: 0

Views

Author

Vaclav Kotesovec, May 27 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=30; CoefficientList[Series[Product[(1+x^k)^(k*(k+1)*(k+2)),{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ 3^(1/5) * Zeta(5)^(1/10) / (2^(91/120) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401*Pi^16 / (1749600000000 * Zeta(5)^3) + 49*Pi^8 * Zeta(3) / (2700000 * Zeta(5)^2) - Zeta(3)^2 / (25*Zeta(5)) + (343*Pi^12/(405000000 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(11/5)) - 7*Pi^4 * Zeta(3) / (750 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (180000 * 2^(3/5) * 3^(4/5) * 5^(2/5) * Zeta(5)^(7/5)) + 3^(1/5) * Zeta(3) / (2^(3/5) * (5*Zeta(5))^(2/5))) * n^(2/5) + 7*Pi^4 / (180 * 2^(2/5) * 3^(1/5) * (5*Zeta(5))^(3/5)) * n^(3/5) + 5*3^(2/5) * (5*Zeta(5)/2)^(1/5)/4 * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.

A258345 Expansion of Product_{k>=1} (1+x^k)^(k*(k-1)*(k-2)).

Original entry on oeis.org

1, 0, 0, 6, 24, 60, 135, 354, 972, 2684, 6990, 17802, 44627, 111582, 277329, 684164, 1671984, 4050096, 9735209, 23238480, 55120950, 129940442, 304502583, 709464798, 1643920584, 3789158988, 8690016942, 19833550266, 45056952957, 101900481462, 229462378987
Offset: 0

Views

Author

Vaclav Kotesovec, May 27 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=40; CoefficientList[Series[Product[(1+x^k)^(k*(k-1)*(k-2)),{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ 3^(1/5) * Zeta(5)^(1/10) / (2^(91/120) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401 * Pi^16 / (1749600000000*Zeta(5)^3) + 49 * Pi^8 * Zeta(3) / (2700000 * Zeta(5)^2) - Zeta(3)^2 / (25*Zeta(5)) + (-343 * Pi^12 / (405000000 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(11/5)) + 7*Pi^4 * Zeta(3) / (750 * 2^(4/5) * 3^(2/5) * 5^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (180000 * 2^(3/5) * 3^(4/5) * 5^(2/5) * Zeta(5)^(7/5)) + 3^(1/5) * Zeta(3) / (2^(3/5) * (5*Zeta(5))^(2/5))) * n^(2/5) - 7*Pi^4 / (180 * 2^(2/5) * 3^(1/5) * (5*Zeta(5))^(3/5)) * n^(3/5) + 5*3^(2/5) * ((5*Zeta(5))/2)^(1/5)/4 * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.

A258346 Expansion of Product_{k>=1} (1+x^k)^(k*(k-1)*(k-2)/6).

Original entry on oeis.org

1, 0, 0, 1, 4, 10, 20, 39, 72, 144, 280, 567, 1112, 2187, 4204, 8073, 15309, 28986, 54548, 102286, 190881, 354717, 656194, 1208712, 2217624, 4052633, 7379630, 13390098, 24215587, 43649482, 78435884, 140513905, 250988186, 447037367, 794031641, 1406585604
Offset: 0

Views

Author

Vaclav Kotesovec, May 27 2015

Keywords

Crossrefs

Programs

  • Mathematica
    nmax=50; CoefficientList[Series[Product[(1+x^k)^(k*(k-1)*(k-2)/6),{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ (3*Zeta(5))^(1/10) / (2^(523/720) * 5^(2/5) * sqrt(Pi) * n^(3/5)) * exp(-2401 * Pi^16 / (10497600000000 * Zeta(5)^3) + 49*Pi^8 * Zeta(3) / (16200000 * Zeta(5)^2) - Zeta(3)^2 / (150*Zeta(5)) + (-343*Pi^12 / (2430000000 * 2^(3/5) * 15^(1/5) * Zeta(5)^(11/5)) + 7*Pi^4 * Zeta(3) / (4500 * 2^(3/5) * 15^(1/5) * Zeta(5)^(6/5))) * n^(1/5) + (-49*Pi^8 / (1080000 * 2^(1/5) * 15^(2/5) * Zeta(5)^(7/5)) + Zeta(3) / (2^(6/5) * (15*Zeta(5))^(2/5))) * n^(2/5) - 7*Pi^4 / (180 * 2^(4/5) * (15*Zeta(5))^(3/5)) * n^(3/5) + 5*(15*Zeta(5))^(1/5) / 2^(12/5) * n^(4/5)), where Zeta(3) = A002117, Zeta(5) = A013663.

A284992 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of Product_{j>=1} (1+x^j)^(j^k) in powers of x.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 4, 5, 2, 1, 1, 8, 13, 8, 3, 1, 1, 16, 35, 31, 16, 4, 1, 1, 32, 97, 119, 83, 28, 5, 1, 1, 64, 275, 457, 433, 201, 49, 6, 1, 1, 128, 793, 1763, 2297, 1476, 487, 83, 8, 1, 1, 256, 2315, 6841, 12421, 11113, 4962, 1141, 142, 10, 1, 1
Offset: 0

Views

Author

Seiichi Manyama, Apr 07 2017

Keywords

Examples

			Square array begins:
  1,  1,   1,    1,     1,      1,       1,        1, ...
  1,  1,   1,    1,     1,      1,       1,        1, ...
  1,  2,   4,    8,    16,     32,      64,      128, ...
  2,  5,  13,   35,    97,    275,     793,     2315, ...
  2,  8,  31,  119,   457,   1763,    6841,    26699, ...
  3, 16,  83,  433,  2297,  12421,   68393,   382573, ...
  4, 28, 201, 1476, 11113,  85808,  678101,  5466916, ...
  5, 49, 487, 4962, 52049, 561074, 6189117, 69540142, ...
		

Crossrefs

Columns k=0-5 give A000009, A026007, A027998, A248882, A248883, A248884.
Rows (0+1),2-3 give: A000012, A000079, A007689.
Main diagonal gives A270917.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(b(n-i*j, i-1, k)*binomial(i^k, j), j=0..n/i)))
        end:
    A:= (n, k)-> b(n$2, k):
    seq(seq(A(n, d-n), n=0..d), d=0..14);  # Alois P. Heinz, Oct 16 2017
  • Mathematica
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0,
         Sum[b[n - i*j, i - 1, k]*Binomial[i^k, j], {j, 0, n/i}]]];
    A[n_, k_] := b[n, n, k];
    Table[Table[A[n, d - n], {n, 0, d}], {d, 0, 14}] // Flatten (* Jean-François Alcover, Feb 10 2021, after Alois P. Heinz *)

Formula

G.f. of column k: Product_{j>=1} (1+x^j)^(j^k).

A261053 Expansion of Product_{k>=1} (1+x^k)^(k^k).

Original entry on oeis.org

1, 1, 4, 31, 289, 3495, 51268, 891152, 17926913, 409907600, 10499834497, 297793199060, 9262502810645, 313457634240463, 11464902463397642, 450646709610954343, 18943070964019019671, 847932498252050293971, 40266255926484893366914, 2021845081107882645459639
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 08 2015

Keywords

Crossrefs

Programs

  • Magma
    m:=20; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[(1+x^k)^(k^k): k in [1..(m+2)]]))); // G. C. Greubel, Nov 08 2018
  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(i^i, j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..25);  # Alois P. Heinz, Aug 08 2015
  • Mathematica
    nmax=20; CoefficientList[Series[Product[(1+x^k)^(k^k),{k,1,nmax}],{x,0,nmax}],x]
  • PARI
    m=20; x='x+O('x^m); Vec(prod(k=1,m, (1+x^k)^(k^k))) \\ G. C. Greubel, Nov 08 2018
    

Formula

a(n) ~ n^n * (1 + exp(-1)/n + (exp(-1)/2 + 4*exp(-2))/n^2).
G.f.: exp(Sum_{k>=1} ( Sum_{d|k} (-1)^(k/d+1)*d^(d+1) ) * x^k/k). - Ilya Gutkovskiy, Nov 08 2018

A261052 Expansion of Product_{k>=1} (1+x^k)^(k!).

Original entry on oeis.org

1, 1, 2, 8, 31, 157, 915, 6213, 48240, 423398, 4147775, 44882107, 531564195, 6837784087, 94909482330, 1413561537884, 22482554909451, 380269771734265, 6815003300096013, 128992737080703803, 2571218642722865352, 53835084737513866662, 1181222084520177393143
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 08 2015

Keywords

Comments

Weigh transform of the factorial numbers. - Alois P. Heinz, Jun 11 2018

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(binomial(i!, j)*b(n-i*j,i-1), j=0..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..25);  # Alois P. Heinz, Aug 08 2015
  • Mathematica
    nmax=25; CoefficientList[Series[Product[(1+x^k)^(k!),{k,1,nmax}],{x,0,nmax}],x]
  • PARI
    seq(n)={Vec(exp(x*Ser(dirmul(vector(n, n, n!), -vector(n, n, (-1)^n/n)))))} \\ Andrew Howroyd, Jun 22 2018

Formula

a(n) ~ n! * (1 + 1/n + 2/n^2 + 10/n^3 + 57/n^4 + 401/n^5 + 3382/n^6 + 33183/n^7 + 371600/n^8 + 4685547/n^9 + 65792453/n^10).

A284897 Expansion of Product_{k>=1} 1/(1+x^k)^(k^3) in powers of x.

Original entry on oeis.org

1, -1, -7, -20, -8, 99, 455, 958, 715, -3606, -17450, -44157, -61852, 19546, 419786, 1442212, 3084950, 3756436, -2155907, -27112107, -88277693, -187777531, -251308697, -5153980, 1182558343, 4299818445, 9988792754, 16075200671, 12020651310, -29802956283
Offset: 0

Views

Author

Seiichi Manyama, Apr 05 2017

Keywords

Crossrefs

Cf. A248882.
Product_{k>=1} 1/(1+x^k)^(k^m): A081362 (m=0), A255528 (m=1), A284896 (m=2), this sequence (m=3), A284898 (m=4), A284899 (m=5).

Programs

  • Mathematica
    CoefficientList[Series[Product[1/(1 + x^k)^(k^3) , {k, 40}], {x, 0, 40}], x] (* Indranil Ghosh, Apr 05 2017 *)
  • PARI
    x= 'x + O('x^40); Vec(prod(k=1, 40, 1/(1 + x^k)^(k^3))) \\ Indranil Ghosh, Apr 05 2017

Formula

a(0) = 1, a(n) = -(1/n)*Sum_{k=1..n} A284900(k)*a(n-k) for n > 0. - Seiichi Manyama, Apr 06 2017

A261051 Expansion of Product_{k>=1} (1+x^k)^(Lucas(k)).

Original entry on oeis.org

1, 1, 3, 7, 14, 33, 69, 148, 307, 642, 1314, 2684, 5432, 10924, 21841, 43431, 85913, 169170, 331675, 647601, 1259737, 2441706, 4716874, 9083215, 17439308, 33387589, 63749174, 121409236, 230658963, 437198116, 826838637, 1560410267, 2938808875, 5524005110
Offset: 0

Views

Author

Vaclav Kotesovec, Aug 08 2015

Keywords

Crossrefs

Programs

  • Maple
    L:= n-> (<<0|1>, <1|1>>^n. <<2, 1>>)[1, 1]:
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
           add(binomial(L(i), j)*b(n-i*j, i-1), j=0..n/i)))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..50);  # Alois P. Heinz, Aug 08 2015
  • Mathematica
    nmax=40; CoefficientList[Series[Product[(1+x^k)^LucasL[k],{k,1,nmax}],{x,0,nmax}],x]

Formula

a(n) ~ phi^n / (2*sqrt(Pi)*n^(3/4)) * exp(-1 + 1/(2*sqrt(5)) + 2*sqrt(n) + s), where s = Sum_{k>=2} (-1)^(k+1) * (2 + phi^k)/((phi^(2*k) - phi^k - 1)*k) = -0.590290697526802161885355317939144642488927381134222996704542... and phi = A001622 = (1+sqrt(5))/2 is the golden ratio.
G.f.: exp(Sum_{k>=1} (-1)^(k+1)*x^k*(1 + 2*x^k)/(k*(1 - x^k - x^(2*k)))). - Ilya Gutkovskiy, May 30 2018

A309335 a(n) = n^3 if n odd, 7*n^3/8 if n even.

Original entry on oeis.org

0, 1, 7, 27, 56, 125, 189, 343, 448, 729, 875, 1331, 1512, 2197, 2401, 3375, 3584, 4913, 5103, 6859, 7000, 9261, 9317, 12167, 12096, 15625, 15379, 19683, 19208, 24389, 23625, 29791, 28672, 35937, 34391, 42875, 40824, 50653, 48013, 59319, 56000, 68921, 64827, 79507, 74536, 91125
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 24 2019

Keywords

Comments

Moebius transform of A007331.

Crossrefs

Programs

  • Mathematica
    a[n_] := If[OddQ[n], n^3, 7 n^3/8]; Table[a[n], {n, 0, 45}]
    nmax = 45; CoefficientList[Series[x (1 + 7 x + 23 x^2 + 28 x^3 + 23 x^4 + 7 x^5 + x^6)/(1 - x^2)^4, {x, 0, nmax}], x]
    LinearRecurrence[{0, 4, 0, -6, 0, 4, 0, -1}, {0, 1, 7, 27, 56, 125, 189, 343}, 46]
    Table[n^3 (15 - (-1)^n)/16, {n, 0, 45}]

Formula

G.f.: x * (1 + 7*x + 23*x^2 + 28*x^3 + 23*x^4 + 7*x^5 + x^6)/(1 - x^2)^4.
G.f.: Sum_{k>=1} J_3(k) * x^k/(1 - x^(2*k)), where J_3() is the Jordan function (A059376).
Dirichlet g.f.: zeta(s-3) * (1 - 1/2^s).
a(n) = n^3 * (15 - (-1)^n)/16.
a(n) = Sum_{d|n, n/d odd} J_3(d).
Sum_{n>=1} 1/a(n) = 57*zeta(3)/56 = 1.223522205001729897639...
Multiplicative with a(2^e) = 7*2^(3*e-3), and a(p^e) = p^(3*e) for odd primes p. - Amiram Eldar, Oct 26 2020
Euler transform is A248882. - Georg Fischer, Nov 10 2020
Previous Showing 11-20 of 23 results. Next