login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A116861 Triangle read by rows: T(n,k) is the number of partitions of n such that the sum of the parts, counted without multiplicities, is equal to k (n>=1, k>=1). 47

%I #51 Jan 11 2024 14:46:04

%S 1,1,1,1,0,2,1,1,1,2,1,0,2,1,3,1,1,3,1,1,4,1,0,3,2,2,2,5,1,1,3,3,2,4,

%T 2,6,1,0,5,2,3,4,4,3,8,1,1,4,3,4,7,4,5,3,10,1,0,5,3,4,7,7,6,6,5,12,1,

%U 1,6,4,3,12,6,8,7,9,5,15,1,0,6,4,5,10,10,9,10,11,10,7,18,1,1,6,4,5,15,11,13,9,16,11,13,8,22

%N Triangle read by rows: T(n,k) is the number of partitions of n such that the sum of the parts, counted without multiplicities, is equal to k (n>=1, k>=1).

%C Conjecture: Reverse the rows of the table to get an infinite lower-triangular matrix b with 1's on the main diagonal. The third diagonal of the inverse of b is minus A137719. - _George Beck_, Oct 26 2019

%C Proof: The reversed rows yield the matrix I+N where N is strictly lower triangular, N[i,j] = 0 for j >= i, having its 2nd diagonal equal to the 2nd column (1, 0, 1, 0, 1, ...): N[i+1,i] = A000035(i), i >= 1, and 3rd diagonal equal to the 3rd column of this triangle, (2, 1, 2, 3, 3, 3, ...): N[i+2,i] = A137719(i), i >= 1. It is known that (I+N)^-1 = 1 - N + N^2 - N^3 +- .... Here N^2 has not only the second but also the 3rd diagonal zero, because N^2[i+2,i] = N[i+2,i+1]*N[i+1,i] = A000035(i+1)*A000035(i) = 0. Therefore the 3rd diagonal of (I+N)^-1 is equal to -A137719 without leading 0. - _M. F. Hasler_, Oct 27 2019

%C From _Gus Wiseman_, Aug 27 2023: (Start)

%C Also the number of ways to write n-k as a nonnegative linear combination of a strict integer partition of k. Also the number of ways to write n as a (strictly) positive linear combination of a strict integer partition of k. Row n=7 counts the following:

%C 7*1 . 1*2+5*1 1*3+4*1 1*3+2*2 1*5+2*1 1*7

%C 2*2+3*1 2*3+1*1 1*4+3*1 1*3+1*2+2*1 1*4+1*3

%C 3*2+1*1 1*5+1*2

%C 1*6+1*1

%C 1*4+1*2+1*1

%C (End)

%H Alois P. Heinz, <a href="/A116861/b116861.txt">Rows n = 1..141, flattened</a>

%H P. J. Rossky, M. Karplus, <a href="https://doi.org/10.1063/1.432387">The enumeration of Goldstone diagrams in many-body perturbation theory</a>, J. Chem. Phys. 64 (1976) 1569, equation (16)(1).

%F G.f.: -1 + Product_{j>=1} (1 + t^j*x^j/(1-x^j)).

%F Sum_{k=1..n} T(n,k) = A000041(n).

%F T(n,n) = A000009(n).

%F Sum_{k=1..n} k*T(n,k) = A014153(n-1).

%F T(n,1) = 1. T(n,2) = A000035(n+1). T(n,3) = A137719(n-2). - _R. J. Mathar_, Oct 27 2019

%F T(n,4) = A002264(n-1) + A121262(n). - _R. J. Mathar_, Oct 28 2019

%e T(10,7) = 4 because we have [6,1,1,1,1], [4,3,3], [4,2,2,1,1] and [4,2,1,1,1,1] (6+1=4+3=4+2+1=7).

%e Triangle starts:

%e 1;

%e 1, 1;

%e 1, 0, 2;

%e 1, 1, 1, 2;

%e 1, 0, 2, 1, 3;

%e 1, 1, 3, 1, 1, 4;

%e 1, 0, 3, 2, 2, 2, 5;

%e 1, 1, 3, 3, 2, 4, 2, 6;

%e 1, 0, 5, 2, 3, 4, 4, 3, 8;

%e 1, 1, 4, 3, 4, 7, 4, 5, 3, 10;

%e 1, 0, 5, 3, 4, 7, 7, 6, 6, 5, 12;

%e 1, 1, 6, 4, 3, 12, 6, 8, 7, 9, 5, 15;

%e ...

%p g:= -1+product(1+t^j*x^j/(1-x^j), j=1..40): gser:= simplify(series(g,x=0,18)): for n from 1 to 14 do P[n]:=sort(coeff(gser,x^n)) od: for n from 1 to 14 do seq(coeff(P[n],t^j),j=1..n) od; # yields sequence in triangular form

%p # second Maple program:

%p b:= proc(n, i) option remember; local f, g, j;

%p if n=0 then [1] elif i<1 then [ ] else f:= b(n, i-1);

%p for j to n/i do

%p f:= zip((x, y)->x+y, f, [0$i, b(n-i*j, i-1)[]], 0)

%p od; f

%p fi

%p end:

%p T:= n-> subsop(1=NULL, b(n, n))[]:

%p seq(T(n), n=1..20); # _Alois P. Heinz_, Feb 27 2013

%t max = 14; s = Series[-1+Product[1+t^j*x^j/(1-x^j), {j, 1, max}], {x, 0, max}, {t, 0, max}] // Normal; t[n_, k_] := SeriesCoefficient[s, {x, 0, n}, {t, 0, k}]; Table[t[n, k], {n, 1, max}, {k, 1, n}] // Flatten (* _Jean-François Alcover_, Jan 17 2014 *)

%t Table[Length[Select[IntegerPartitions[n],Total[Union[#]]==k&]],{n,0,10},{k,0,n}] (* _Gus Wiseman_, Aug 29 2023 *)

%o (PARI) A116861(n,k,s=0)={forpart(X=n,vecsum(Set(X))==k&&s++,k);s} \\ _M. F. Hasler_, Oct 27 2019

%Y Cf. A000041 (row sums), A000009 (diagonal), A014153.

%Y Cf. A114638 (count partitions with #parts = sum(distinct parts)).

%Y Column 1: A000012, column 2: A000035(1..), column 3: A137719(1..).

%Y For subsets instead of partitions we have A026820.

%Y This statistic is ranked by A066328.

%Y The central diagonal is T(2n,n) = A364910(n), non-strict A364907.

%Y Partial sums of columns are columns of A364911.

%Y Same as A364916 (offset 0) with rows reversed.

%Y A008284 counts partitions by length, strict A008289.

%Y A364912 counts linear combinations of partitions.

%Y A364913 counts combination-full partitions, strict A364839.

%Y Cf. A237667, A364272, A364915, A365002, A365004, A365006.

%K nonn,tabl

%O 1,6

%A _Emeric Deutsch_, Feb 27 2006

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 10:22 EDT 2024. Contains 371967 sequences. (Running on oeis4.)