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 29 results. Next

A292137 G.f.: Im(1/(i*x; x)_inf), where (a; q)_inf is the q-Pochhammer symbol, i = sqrt(-1).

Original entry on oeis.org

0, 1, 1, 0, 0, 0, -1, -2, -2, -2, -2, -3, -3, -2, -2, -2, -1, 1, 2, 2, 4, 6, 7, 8, 10, 13, 14, 14, 15, 17, 17, 15, 15, 16, 14, 10, 8, 6, 1, -5, -10, -14, -21, -31, -38, -43, -53, -64, -71, -77, -86, -97, -104, -108, -115, -124, -127, -125, -127, -130, -125, -116
Offset: 0

Views

Author

Seiichi Manyama, Sep 09 2017

Keywords

Examples

			Product_{k>=1} 1/(1 - i*x^k) = 1 + (0+1i)*x + (-1+1i)*x^2 + (-1+0i)*x^3 + (-1+0i)*x^4 + (-1+0i)*x^5 + (-2-1i)*x^6 + (-1-2i)*x^7 + ...
		

Crossrefs

Programs

  • Maple
    N:= 100:
    S := convert(series( add( (-1)^n*x^(2*n+1)/(mul(1 - x^k,k = 1..2*n+1)), n = 0..N ), x, N+1 ), polynom):
    seq(coeff(S, x, n), n = 0..N); # Peter Bala, Jan 15 2021
  • Mathematica
    Im[CoefficientList[Series[1/QPochhammer[I*x, x], {x, 0, 100}], x]] (* Vaclav Kotesovec, Sep 17 2017 *)

Formula

1/(i*x; x)_inf is the g.f. for A292136(n) + i*a(n).
a(n) = Sum (-1)^((k - 1)/2) where the sum is over all integer partitions of n into an odd number of parts and k is the number of parts. - Gus Wiseman, Mar 08 2018
G.f.: Sum_{n >= 0} (-1)^n * x^(2*n+1)/Product_{k = 1..2*n+1} (1 - x^k). - Peter Bala, Jan 15 2021

A300354 Number of enriched p-trees of weight n with distinct leaves.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 8, 8, 13, 17, 54, 56, 98, 125, 195, 500, 606, 921, 1317, 1912, 2635, 6667, 7704, 12142, 16958, 24891, 33388, 47792, 106494, 126475, 195475, 268736, 393179, 523775, 750251, 979518, 2090669, 2457315, 3759380, 5066524, 7420874, 9726501, 13935546
Offset: 0

Views

Author

Gus Wiseman, Mar 03 2018

Keywords

Comments

An enriched p-tree of weight n > 0 is either a single node of weight n, or a sequence of two or more enriched p-trees with weakly decreasing weights summing to n.

Examples

			The a(6) = 8 enriched p-trees with distinct leaves: 6, (42), (51), ((31)2), ((32)1), (3(21)), ((21)3), (321).
		

Crossrefs

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    ept[q_]:=ept[q]=If[Length[q]===1,1,Total[Times@@@Map[ept,Join@@Function[sptn,Join@@@Tuples[Permutations/@GatherBy[sptn,Total]]]/@Select[sps[q],Length[#]>1&],{2}]]];
    Table[Total[ept/@Select[IntegerPartitions[n],UnsameQ@@#&]],{n,1,30}]

Formula

a(n) = Sum_{i=1..A000009(n)} A299203(A246867(n,i)).

A220420 Express the Sum_{n>=0} p(n)*x^n, where p(n) is the partition function, as a product Product_{k>=1} (1 + a(k)*x^k).

Original entry on oeis.org

1, 2, 1, 4, 1, 0, 1, 14, 1, -4, 1, -8, 1, -16, 1, 196, 1, -54, 1, -92, 1, -184, 1, 144, 1, -628, 1, -1040, 1, -2160, 1, 41102, 1, -7708, 1, -12932, 1, -27592, 1, 54020, 1, -98496, 1, -173720, 1, -364720, 1, 853624, 1, -1341970, 1, -2383916, 1, -4918536, 1
Offset: 1

Views

Author

Michel Marcus, Dec 14 2012

Keywords

Comments

This is the PPE (power product expansion) of A000041.
When n is odd, a(n) = 1.
When n is even, a(n) = 2, 4, 0, 14, -4, -8, -16, 196, -54, -92, -184, 144, -628, -1040, -2160, 41102, ...
Alkauskas (2016, Problem 3, p. 3) conjectured that a(8*k+2), a(8*k+4), and a(8*k+6) are all negative, and a(8*k) is positive for k >= 1. [This statement is not wholly true for k = 0.] - Petros Hadjicostas, Oct 07 2019

Crossrefs

Programs

  • Mathematica
    terms = 55; sol[0] = {};
    sol[m_] := sol[m] = Join[sol[m - 1], If[OddQ[m], {a[m] -> 1}, First @ Solve[Thread[Table[PartitionsP[n], {n, 0, m}] == CoefficientList[ (Product[1 + a[n]*x^n, {n, 1, m}] /. sol[m - 1]) + O[x]^(m + 1), x]]]]];
    Array[a, terms] /. sol[terms] (* Jean-François Alcover, Dec 06 2018, corrected Oct 03 2019 *)
    (* Second program: *)
    A[m_, n_] := A[m, n] = Which[m == 1, PartitionsP[n], m > n >= 1, 0, True, A[m - 1, n] - A[m - 1, m - 1]*A[m, n - m + 1] ];
    a[n_] := A[n, n];
    a /@ Range[1, 55] (* Jean-François Alcover, Oct 03 2019, using the formula given by Petros Hadjicostas *)
  • PARI
    a(m) = {default(seriesprecision, m+1); ak = vector(m); pol = 1 / eta(x + x * O(x^m)); ak[1] = polcoeff(pol, 1); for (k=2, m, pol = taylor(pol / (1+ak[k-1]*x^(k-1)), x); ak[k] = polcoeff(pol, k, x);); for (k=1, m, print1(ak[k], ", "););}

Formula

From Petros Hadjicostas, Oct 04 2019: (Start)
Define (A(m,n): n,m >= 1) by A(m=1,n) = p(n) = A000041(n) for n >= 1, A(m,n) = 0 for m > n >= 1 (upper triangular), and A(m,n) = A(m-1,n) - A(m-1,m-1) * A(m,n-m+1) for n >= m >= 2. Then a(n) = A(n,n). [Theorem 3 in Gingold et al. (1988).]
a(n) = Sum_{s|n} s/n + Sum_{s|n, s > 1} (-a(n/s))^s/s. [Eq. (1) in Alkauskas (2008, 2009).]
(End)

A290262 Irregular triangle read by rows: rows give the (negated) nonzero coefficients of t in each term of the inverse power product expansion of 1 - t * x/(1-x).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 3, 4, 2, 1, 3, 5, 5, 3, 1, 1, 4, 9, 13, 13, 9, 4, 1, 1, 4, 9, 13, 13, 9, 4, 1, 1, 5, 14, 25, 30, 24, 12, 3, 1, 5, 15, 30, 42, 42, 30, 15, 5, 1, 1, 6, 21, 48, 75, 81, 60, 30, 10, 2
Offset: 1

Views

Author

Gus Wiseman, Jul 24 2017

Keywords

Comments

Row sums are A290261(n). A regular version is A290320.

Examples

			Triangle begins:
  1;
  1,  1;
  1,  1,
  1,  2,  2,  1;
  1,  2,  2,  1;
  1,  3,  4,  2;
  1,  3,  5,  5,  3,  1;
  1,  4,  9, 13, 13,  9,  4,  1;
  1,  4,  9, 13, 13,  9,  4,  1;
  1,  5, 14, 25, 30, 24, 12,  3;
  1,  5, 15, 30, 42, 42, 30, 15,  5,  1;
  1,  6, 21, 48, 75, 81, 60, 30, 10,  2;
		

Crossrefs

Programs

  • Mathematica
    eptrees[n_]:=Prepend[Join@@Table[Tuples[eptrees/@y],{y,Rest[IntegerPartitions[n]]}],n];
    eptrans[a_][n_]:=Sum[(-1)^Count[tree,_List,{0,Infinity}]*Product[a[i],{i,Flatten[{tree}]}],{tree,eptrees[n]}];
    Table[DeleteCases[CoefficientList[-eptrans[-t&][n],t],0],{n,12}]

A300862 Solution to 1 = Sum_y Product_{k in y} a(k) for each n > 0, where the sum is over all integer partitions of n with an odd number of parts.

Original entry on oeis.org

1, 1, 0, 0, -1, -1, 0, 1, 1, 0, -2, -3, -2, 2, 7, 6, -3, -15, -19, -2, 32, 54, 24, -64, -153, -123, 95, 389, 444, -43, -966, -1475, -516, 2066, 4414, 3092, -3874, -12480, -12936, 3847, 32445, 45494, 8950, -77282, -147663, -86313, 157456, 435623, 399041, -229616, -1211479, -1535700, -73132
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=1-Sum[Times@@a/@y,{y,Select[IntegerPartitions[n],Length[#]>1&&OddQ[Length[#]]&]}];
    Array[a,40]

A300863 Signed recurrence over enriched p-trees: a(n) = (-1)^(n - 1) + Sum_{y1 + ... + yk = n, y1 >= ... >= yk > 0, k > 1} a(y1) * ... * a(yk).

Original entry on oeis.org

1, 0, 2, 2, 6, 14, 34, 82, 214, 566, 1482, 4058, 10950, 30406, 83786, 235714, 658286, 1874254, 5293674, 15189810, 43312542, 125075238, 359185586, 1043712922, 3015569582, 8800146182, 25565402802, 74918274562, 218572345718, 642783954238, 1882606578002
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=(-1)^(n-1)+Sum[Times@@a/@y,{y,Select[IntegerPartitions[n],Length[#]>1&]}];
    Array[a,40]

Formula

O.g.f.: (-1/(1+x) + Product 1/(1-a(n)x^n))/2.

A290973 Write 2*x/(1-x) in the form Sum_{j>=1} ((1-x^j)^a(j) - 1).

Original entry on oeis.org

-2, 1, 2, 3, 4, 6, 6, 10, 8, 15, 10, 25, 12, 28, 10, 60, 16, 25, 18, 125, 0, 66, 22, 218, 24, 91, -30, 420, 28, -387, 30, 2011, -88, 153, 28, -1894, 36, 190, -182, 8902, 40, -3234, 42, 2398, -132, 276, 46, 2340, 48, -2678, -510, 4641, 52, -1754, -198, 108400
Offset: 1

Views

Author

Gus Wiseman, Aug 16 2017

Keywords

Examples

			2x/(1-x) = (1-x)^(-2) - 1 + (1-x^2)^1 - 1 + (1-x^3)^2 - 1 + (1-x^4)^3 - 1 + ...
		

Crossrefs

Programs

  • Maple
    a:= n-> add(binomial(n/d-1-a(d), n/d), d=
            numtheory[divisors](n) minus {n})-2:
    seq(a(n), n=1..60);  # Alois P. Heinz, Aug 27 2017
  • Mathematica
    nn=60;
    rus=SolveAlways[Normal[Series[2x/(1-x)==Sum[(1-x^n)^a[n]-1,{n,nn}],{x,0,nn}]],x];
    Array[a,nn]/.First[rus]

Formula

For all n > 0 we have: 2 = Sum_{d|n} binomial(-a(d) + n/d - 1, n/d).

A295632 Write 1/Product_{n > 1}(1 - 1/n^s) in the form Product_{n > 1}(1 + a(n)/n^s).

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1
Offset: 2

Views

Author

Gus Wiseman, Nov 24 2017

Keywords

Comments

First negative entry is a(1024) = -4.

Crossrefs

Programs

  • Mathematica
    nn=100;
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Solve[Table[Length[facs[n]]==Sum[Times@@a/@f,{f,Select[facs[n],UnsameQ@@#&]}],{n,2,nn}],Table[a[n],{n,2,nn}]][[1,All,2]]

A300864 Signed recurrence over strict trees: a(n) = -1 + Sum_{y1 + ... + yk = n, y1 > ... > yk > 0, k > 1} a(y1) * ... * a(yk).

Original entry on oeis.org

1, 1, 0, 1, 0, 0, 2, -1, 0, 4, -6, 6, 6, -24, 38, -17, -64, 188, -230, -6, 662, -1432, 1286, 1210, -6362, 10692, -5530, -18274, 57022, -74364, 174, 216703, -489544, 467860, 391258, -2256430, 3948206, -2234064, -6725362, 21920402, -29716570, 2095564, 84595798, -198418242, 197499846
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=-1+Sum[Times@@a/@y,{y,Select[IntegerPartitions[n],Length[#]>1&&UnsameQ@@#&]}];
    -Array[a,40]

A300865 Signed recurrence over binary enriched p-trees: a(n) = (-1)^(n-1) + Sum_{x + y = n, 0 < x <= y < n} a(x) * a(y).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 2, 2, 4, 6, 10, 16, 27, 46, 77, 131, 224, 391, 672, 1180, 2050, 3626, 6344, 11276, 19863, 35479, 62828, 112685, 200462, 360627, 644199, 1162296, 2083572, 3768866, 6777314, 12289160, 22158106, 40255496, 72765144, 132453122, 239936528, 437445448
Offset: 1

Views

Author

Gus Wiseman, Mar 13 2018

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_]:=a[n]=(-1)^(n-1)+Sum[a[k]*a[n-k],{k,1,n/2}];
    Array[a,50]
Previous Showing 11-20 of 29 results. Next