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-10 of 10 results.

A175804 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the n-th term of the k-th differences of partition numbers A000041.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, -1, 0, 1, 3, 2, 1, 1, 2, 5, -4, -2, -1, 0, 2, 7, 9, 5, 3, 2, 2, 4, 11, -21, -12, -7, -4, -2, 0, 4, 15, 49, 28, 16, 9, 5, 3, 3, 7, 22, -112, -63, -35, -19, -10, -5, -2, 1, 8, 30, 249, 137, 74, 39, 20, 10, 5, 3, 4, 12, 42, -539, -290, -153, -79, -40, -20, -10, -5, -2, 2, 14, 56
Offset: 0

Views

Author

Alois P. Heinz, Dec 04 2010

Keywords

Comments

Odlyzko showed that the k-th differences of A000041(n) alternate in sign with increasing n up to a certain index n_0(k) and then stay positive.
Are there any zeros after the first four, which all lie in columns k = 1, 2? - Gus Wiseman, Dec 15 2024

Examples

			Square array A(n,k) begins:
   1,  0,  1, -1,  2,  -4,   9,  ...
   1,  1,  0,  1, -2,   5, -12,  ...
   2,  1,  1, -1,  3,  -7,  16,  ...
   3,  2,  0,  2, -4,   9, -19,  ...
   5,  2,  2, -2,  5, -10,  20,  ...
   7,  4,  0,  3, -5,  10, -20,  ...
  11,  4,  3, -2,  5, -10,  22,  ...
		

Crossrefs

Columns k=0-5 give: A000041, A002865, A053445, A072380, A081094, A081095.
Main diagonal gives A379378.
For primes we have A095195 or A376682.
Row n = 0 is A281425.
Row n = 1 is A320590 except first term.
For composites we have A377033.
For squarefree numbers we have A377038.
For nonsquarefree numbers we have A377046.
For prime powers we have A377051.
Antidiagonal sums are A377056, absolute value version A378621.
The version for strict partitions is A378622, first column A293467.
A000009 counts strict integer partitions, differences A087897, A378972.

Programs

  • Maple
    A41:= combinat[numbpart]:
    DD:= proc(p) proc(n) option remember; p(n+1) -p(n) end end:
    A:= (n,k)-> (DD@@k)(A41)(n):
    seq(seq(A(n, d-n), n=0..d), d=0..11);
  • Mathematica
    max = 11; a41 = Array[PartitionsP, max+1, 0]; a[n_, k_] := Differences[a41, k][[n+1]]; Table[a[n, k-n], {k, 0, max}, {n, 0, k}] // Flatten (* Jean-François Alcover, Aug 29 2014 *)
    nn=5;Table[Table[Sum[(-1)^(k-i)*Binomial[k,i]*PartitionsP[n+i],{i,0,k}],{k,0,nn}],{n,0,nn}] (* Gus Wiseman, Dec 15 2024 *)

Formula

A(n,k) = (Delta^(k) A000041)(n).
A(n,k) = Sum_{i=0..k} (-1)^(k-i) * binomial(k,i) * A000041(n+i). In words, row x is the inverse zero-based binomial transform of A000041 shifted left x times. - Gus Wiseman, Dec 15 2024

A281425 a(n) = [q^n] (1 - q)^n / Product_{j=1..n} (1 - q^j).

Original entry on oeis.org

1, 0, 1, -1, 2, -4, 9, -21, 49, -112, 249, -539, 1143, -2396, 5013, -10550, 22420, -48086, 103703, -223806, 481388, -1029507, 2187944, -4625058, 9742223, -20490753, 43111808, -90840465, 191773014, -405523635, 858378825, -1817304609, 3845492204, -8129023694, 17162802918, -36191083386
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 05 2017

Keywords

Comments

a(n) is n-th term of the Euler transform of -n + 1, 1, 1, 1, ...
Inverse zero-based binomial transform of A000041. The version for strict partitions is A380412, or A293467 up to sign. - Gus Wiseman, Feb 06 2025

Crossrefs

Programs

  • Maple
    b:= proc(n, k) option remember; `if`(k=0,
          combinat[numbpart](n), b(n, k-1)-b(n-1, k-1))
        end:
    a:= n-> b(n$2):
    seq(a(n), n=0..35);  # Alois P. Heinz, Dec 21 2024
  • Mathematica
    Table[SeriesCoefficient[(1 - q)^n / Product[(1 - q^j), {j, 1, n}], {q, 0, n}], {n, 0, 35}]
    Table[SeriesCoefficient[(1 - q)^n QPochhammer[q^(1 + n), q]/QPochhammer[q, q], {q, 0, n}], {n, 0, 35}]
    Table[SeriesCoefficient[1/QFactorial[n, q], {q, 0, n}], {n, 0, 35}]
    Table[Differences[PartitionsP[Range[0, n]], n], {n, 0, 35}] // Flatten
    Table[Sum[(-1)^j*Binomial[n, j]*PartitionsP[n-j], {j, 0, n}], {n, 0, 30}] (* Vaclav Kotesovec, Oct 06 2017 *)

Formula

a(n) = [q^n] 1/((1 + q)*(1 + q + q^2)*...*(1 + q + ... + q^(n-1))).
a(n) = Sum_{j=0..n} (-1)^j * binomial(n, j) * A000041(n-j). - Vaclav Kotesovec, Oct 06 2017
a(n) ~ (-1)^n * 2^(n - 3/2) * exp(Pi*sqrt(n/12) + Pi^2/96) / (sqrt(3)*n). - Vaclav Kotesovec, May 07 2018

A320590 Expansion of Product_{k>=1} 1/(1 - x^k/(1 + x)^k).

Original entry on oeis.org

1, 1, 1, 0, 1, -2, 5, -12, 28, -63, 137, -290, 604, -1253, 2617, -5537, 11870, -25666, 55617, -120103, 257582, -548119, 1158437, -2437114, 5117165, -10748530, 22621055, -47728657, 100932549, -213750621, 452855190, -958925784, 2028187595, -4283531490, 9033779224
Offset: 0

Views

Author

Ilya Gutkovskiy, Oct 16 2018

Keywords

Comments

The zero-based binomial transform of this sequence is A000070, and if we remove first terms it becomes A000041.

Crossrefs

Row n=1 of A175804 (except first term). Row n=0 is A281425.
The version for strict partitions is A320591, row n=1 of A378622, first column A293467.
A000009 counts strict integer partitions, differences A087897, A378972.
A000041 counts integer partitions, differences A002865.

Programs

  • Magma
    m:=50; R:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1 - x^k/(1 + x)^k): k in [1..(m+2)]]) )); // G. C. Greubel, Oct 29 2018
  • Maple
    seq(coeff(series(mul(1/(1-x^k/(1+x)^k),k=1..n),x,n+1), x, n), n = 0 .. 35); # Muniru A Asiru, Oct 16 2018
  • Mathematica
    nmax = 34; CoefficientList[Series[Product[1/(1 - x^k/(1 + x)^k), {k, 1, nmax}], {x, 0, nmax}], x]
    nmax = 34; CoefficientList[Series[Exp[Sum[DivisorSigma[1, k] x^k/(k (1 + x)^k), {k, 1, nmax}]], {x, 0, nmax}], x]
  • PARI
    m=50; x='x+O('x^m); Vec(prod(k=1, m+2, 1/(1 - x^k/(1 + x)^k))) \\ G. C. Greubel, Oct 29 2018
    

Formula

G.f.: exp(Sum_{k>=1} x^k/(k*((1 + x)^k - x^k))).
G.f.: exp(Sum_{k>=1} sigma(k)*x^k/(k*(1 + x)^k)).

A377054 First term of the n-th differences of the powers of primes. Inverse zero-based binomial transform of A000961.

Original entry on oeis.org

1, 1, 0, 0, 0, 1, -5, 15, -34, 63, -97, 115, -54, -251, 1184, -3536, 8736, -18993, 37009, -64545, 98442, -121393, 82008, 147432, -860818, 2710023, -7110594, 17077281, -38873146, 85085287, -179965647, 367885014, -725051280, 1372311999, -2481473550, 4257624252
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2024

Keywords

Examples

			The sixth differences of A000961 begin: -5, 10, -9, 1, 6, -10, 16, -18, ..., so a(6) = -5.
		

Crossrefs

The version for primes is A007442, noncomposites A030016, composites A377036.
For squarefree numbers we have A377041, nonsquarefree A377049.
This is the first column of the array A377051.
For antidiagonal-sums we have A377052, absolute A377053.
For positions of first zeros we have A377055.
A000040 lists the primes, differences A001223, seconds A036263.
A000961 lists the powers of primes, differences A057820.
A001597 lists perfect-powers, complement A007916.
A008578 lists the noncomposites, differences A075526.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.

Programs

  • Mathematica
    q=Select[Range[100],#==1||PrimePowerQ[#]&];
    Table[Sum[(-1)^(j-k)*Binomial[j,k]*q[[1+k]],{k,0,j}],{j,0,Length[q]/2}]

Formula

The inverse zero-based binomial transform of a sequence (q(0), q(1), q(2), ...) is the sequence p given by:
p(j) = sum_{k=0..j} (-1)^(j-k)*binomial(j,k)*q(k)

A377052 Antidiagonal-sums of the array A377051(n,k) = n-th term of k-th differences of powers of primes.

Original entry on oeis.org

1, 3, 4, 5, 6, 13, -6, 45, -50, 113, -98, 73, 274, -1159, 3563, -8707, 19024, -36977, 64582, -98401, 121436, -81961, -147383, 860871, -2709964, 7110655, -17077217, 38873213, -85085216, 179965720, -367884935, 725051361, -1372311916, 2481473639, -4257624155
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2024

Keywords

Comments

These are the row-sums of the triangle-version of A377051.

Examples

			The sixth antidiagonal of A377051 is (8, 1, -1, -2, -3, -4, -5), so a(6) = -6.
		

Crossrefs

The version for primes is A140119, noncomposites A376683, composites A377034.
For squarefree numbers we have A377039, nonsquarefree A377047.
These are the antidiagonal-sums of A377051.
The unsigned version is A377053.
For leaders we have A377054, for primes A007442 or A030016.
For first zero-positions we have A377055.
A version for partitions is A377056, cf. A175804, A053445, A281425, A320590.
A000040 lists the primes, differences A001223, seconds A036263.
A001597 lists perfect-powers, complement A007916.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.

Programs

  • Mathematica
    nn=20;
    t=Table[Differences[NestList[NestWhile[#+1&, #+1,!PrimePowerQ[#]&]&,1,2*nn],k],{k,0,nn}];
    Total/@Table[t[[j,i-j+1]],{i,nn},{j,i}]

A377053 Antidiagonal-sums of the absolute value of the array A377051(n,k) = n-th term of k-th differences of powers of primes.

Original entry on oeis.org

1, 3, 4, 5, 6, 13, 24, 45, 80, 123, 174, 229, 382, 1219, 3591, 8849, 19288, 37899, 67442, 108323, 156054, 206733, 311525, 860955, 2710374, 7111657, 17080759, 38884849, 85124764, 180097856, 368321633, 726482493, 1377039690, 2496856437, 4306569569, 7016267449
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2024

Keywords

Comments

These are the row-sums of the absolute value of the triangle-version of A377051.

Examples

			The sixth antidiagonal of A377051 is (8, 1, -1, -2, -3, -4, -5), so a(6) = 24.
		

Crossrefs

The version for primes is A376681, noncomposites A376684, composites A377035.
For squarefree numbers we have A377040, nonsquarefree A377048.
This is the antidiagonal-sums of the absolute value of A377051.
The signed version is A377052.
For leaders we have A377054, for primes A007442 or A030016.
For first zero-positions we have A377055.
A version for partitions is A377056, cf. A175804, A053445, A281425, A320590.
A000040 lists the primes, differences A001223, seconds A036263.
A008578 lists the noncomposites, differences A075526.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.

Programs

  • Mathematica
    nn=20;
    t=Table[Differences[NestList[NestWhile[#+1&, #+1,!PrimePowerQ[#]&]&,1,2*nn],k],{k,0,nn}];
    Total/@Abs[Table[t[[j,i-j+1]],{i,nn},{j,i}]]

A377055 Position of first appearance of zero in the n-th differences of the prime-powers (A246655), or 0 if it does not appear.

Original entry on oeis.org

0, 0, 1, 1, 4, 48, 61, 83, 29, 57, 290, 121, 7115, 14207, 68320, 14652, 149979, 122704, 481540, 980376, 632441, 29973, 25343678, 50577935, 7512418, 210836403, 67253056, 224083553, 910629561, 931524323, 452509699, 2880227533, 396690327, 57954538325, 77572935454, 35395016473
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2024

Keywords

Examples

			The fourth differences of A246655 begin: 1, -3, 3, 0, -2, 2, ... so a(4) = 4.
		

Crossrefs

The version for primes is A376678, noncomposites A376855, composites A377037.
For squarefree numbers we have A377042, nonsquarefree A377050.
These are the positions of first zeros in each row of A377051.
For antidiagonal-sums we have A377052, absolute A377053.
For leaders we have A377054, for primes A007442 or A030016.
A000040 lists the primes, differences A001223, seconds A036263.
A000961 lists the powers of primes, differences A057820.
A008578 lists the noncomposites, differences A075526.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.
A246655 lists the prime-powers, differences A057820 (except first term).

Programs

  • Mathematica
    nn=10000;
    u=Table[Differences[Select[Range[nn],PrimePowerQ],k],{k,2,16}];
    mnrm[s_]:=If[Min@@s==1,mnrm[DeleteCases[s-1,0]]+1,0];
    m=Table[Position[u[[k]],0][[1,1]], {k,mnrm[Union[First/@Position[u,0]]]}]

Extensions

a(12)-a(27) from Pontus von Brömssen, Oct 22 2024
a(28)-a(30) from Chai Wah Wu, Oct 23 2024
a(31)-a(35) from Lucas A. Brown, Nov 03 2024

A378970 Antidiagonal-sums of the array A378622(n,k) = n-th term of k-th differences of strict partition numbers (A000009).

Original entry on oeis.org

1, 1, 1, 5, -4, 18, -20, 47, -56, 110, -153, 309, -532, 1045, -1768, 2855, -3620, 2928, 2927, -20371, 62261, -148774, 314112, -613835, 1155936, -2175658, 4244218, -8753316, 19006746, -42471491, 95234915, -210395017, 453414314, -949507878, 1931940045
Offset: 0

Views

Author

Gus Wiseman, Dec 14 2024

Keywords

Examples

			Antidiagonal 4 of A378622 is (2, 0, -1, -2, -3), so a(4) = -4.
		

Crossrefs

For primes we have A140119 or A376683, absolute value A376681 or A376684.
For composites we have A377034, absolute value A377035.
For squarefree numbers we have A377039, absolute value A377040.
For nonsquarefree numbers we have A377047, absolute value A377048.
For prime powers we have A377052, absolute value A377053.
For partition numbers we have A377056, absolute value A378621.
Row-sums of the triangular form of A378622. See also:
- A175804 is the version for partitions.
- A293467 gives the first column (up to sign).
- A377285 gives position of first zero in each row.
The unsigned version is A378971.
A000009 counts strict integer partitions, differences A087897, A378972.
A000041 counts integer partitions, differences A002865, A053445.

Programs

  • Mathematica
    nn=30;
    t=Table[Take[Differences[PartitionsQ/@Range[0,2nn],k],nn],{k,0,nn}];
    Total/@Table[t[[j,i-j+1]],{i,nn/2},{j,i}]

A378971 Antidiagonal-sums of absolute value of the array A378622(n,k) = n-th term of k-th differences of strict partition numbers (A000009).

Original entry on oeis.org

1, 1, 1, 5, 8, 18, 30, 47, 70, 110, 177, 309, 574, 1063, 1892, 3107, 4598, 6166, 8737, 20603, 62457, 149132, 314116, 614093, 1155968, 2176048, 4244322, 8753864, 19006756, 42472117, 95235017, 210396059, 453414950, 949510166, 1931941261, 3826650257, 7400745917
Offset: 0

Views

Author

Gus Wiseman, Dec 14 2024

Keywords

Examples

			Antidiagonal 4 of A378622 is (2, 0, -1, -2, -3), so a(4) = 8.
		

Crossrefs

For primes we have A376681 or A376684, signed version A140119 or A376683.
For composites we have A377035, signed version A377034.
For squarefree numbers we have A377040, signed version A377039.
For nonsquarefree numbers we have A377048, signed version A377049.
For prime powers we have A377053, signed version A377052.
For partition numbers we have A378621, signed version A377056.
Row-sums of the triangular form of A378622. See also:
- A175804 is the version for partitions.
- A293467 gives the first column (up to sign).
- A377285 gives position of first zero in each row.
The signed version is A378970.
A000009 counts strict integer partitions, differences A087897, A378972.
A000041 counts integer partitions, differences A002865, A053445.

Programs

  • Mathematica
    nn=30;
    t=Table[Take[Differences[PartitionsQ/@Range[0,2nn],k],nn],{k,0,nn}];
    Total/@Abs/@Table[t[[j,i-j+1]],{i,nn/2},{j,i}]

A378621 Antidiagonal-sums of absolute value of the array A175804(n,k) = n-th term of k-th differences of partition numbers (A000041).

Original entry on oeis.org

1, 1, 4, 5, 11, 16, 36, 65, 142, 285, 595, 1210, 2497, 5134, 10726, 22637, 48383, 104066, 224296, 481985, 1030299, 2188912, 4626313, 9743750, 20492711, 43114180, 90843475, 191776658, 405528200, 858384333, 1817311451, 3845500427, 8129033837, 17162815092
Offset: 0

Views

Author

Gus Wiseman, Dec 14 2024

Keywords

Examples

			Antidiagonal i + j = 3 of A175804 is (3, 1, 0, -1), so a(3) = 5.
		

Crossrefs

These are the antidiagonal-sums of the absolute value of A175804.
First column of the same array is A281425.
For primes we have A376681 or A376684, signed A140119 or A376683.
For composites we have A377035, signed A377034.
For squarefree numbers we have A377040, signed A377039.
For nonsquarefree numbers we have A377048, signed A377049.
For prime powers we have A377053, signed A377052.
The signed version is A377056.
The corresponding array for strict partitions is A378622, see A293467, A377285, A378971, A378970.
A000009 counts strict integer partitions, differences A087897, A378972.
A000041 counts integer partitions, differences A002865, A053445.

Programs

  • Mathematica
    nn=30;
    q=Table[PartitionsP[n],{n,0,nn}];
    t=Table[Sum[(-1)^(j-k)*Binomial[j,k]*q[[i+k]],{k,0,j}],{j,0,Length[q]/2},{i,Length[q]/2}]
    Total/@Abs/@Table[t[[j,i-j+1]],{i,nn/2},{j,i}]
Showing 1-10 of 10 results.