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

A197126 Triangle T(n,k), n>=1, 1<=k<=n, read by rows: T(n,k) is the number of cliques of size k in all partitions of n.

Original entry on oeis.org

1, 1, 1, 3, 0, 1, 4, 2, 0, 1, 8, 2, 1, 0, 1, 11, 4, 2, 1, 0, 1, 19, 5, 3, 1, 1, 0, 1, 26, 10, 3, 3, 1, 1, 0, 1, 41, 11, 7, 3, 2, 1, 1, 0, 1, 56, 20, 8, 5, 3, 2, 1, 1, 0, 1, 83, 25, 13, 6, 5, 2, 2, 1, 1, 0, 1, 112, 38, 17, 11, 5, 5, 2, 2, 1, 1, 0, 1, 160, 49, 25, 13, 9, 5, 4, 2, 2, 1, 1, 0, 1
Offset: 1

Views

Author

Alois P. Heinz, Oct 10 2011

Keywords

Comments

All parts of a number partition with the same value form a clique. The size of a clique is the number of elements in the clique.

Examples

			T(4,1) = 4: [1,1,(2)], [(1),(3)], [(4)].
T(8,3) = 3: [1,1,(2,2,2)], [(1,1,1),2,3], [(1,1,1),5].
T(12,4) = 11: [(1,1,1,1),(2,2,2,2)], [1,(2,2,2,2),3], [(1,1,1,1),2,3,3], [(3,3,3,3)], [(1,1,1,1),2,2,4], [(2,2,2,2),4], [(1,1,1,1),4,4], [(1,1,1,1),3,5], [(1,1,1,1),2,6], [(1,1,1,1),8].  Here the first partition contains 2 cliques.
Triangle begins:
   1;
   1,  1;
   3,  0, 1;
   4,  2, 0, 1;
   8,  2, 1, 0, 1;
  11,  4, 2, 1, 0, 1;
  19,  5, 3, 1, 1, 0, 1;
  26, 10, 3, 3, 1, 1, 0, 1;
  ...

Links

Crossrefs

Columns k=1-10 give: A024786(n+1), A116646, A117524, A222704, A222705, A222706, A222707, A222708, A222709, A222710.
Row sums give: A000070(n-1). Diagonal gives: A000012. Limit of reversed rows: T(2*n+1,n+1) = A002865(n).
Cf. A213180.

Programs

  • Maple
    b:= proc(n, p, k) option remember; `if`(n=0, [1, 0], `if`(p<1, [0, 0],
      add((l->`if`(m=k, l+[0, l[1]], l))(b(n-p*m, p-1, k)), m=0..n/p)))
    end:
T:= (n, k)-> b(n, n, k)[2]:
seq(seq(T(n, k), k=1..n), n=1..20);
  • Mathematica
    Table[CoefficientList[ 1/q* Tr[Flatten[q^Map[Length, Split /@ IntegerPartitions[n], {2}]]], q], {n, 24}] (* Wouter Meeussen, Apr 21 2012 *)
b[n_, p_, k_] := b[n, p, k] = If[n == 0, {1, 0}, If[p < 1, {0, 0}, Sum[ Function[l, If[m == k, l + {0, l[[1]]}, l]][b[n - p*m, p - 1, k]], {m, 0, n/p}]]]; T[n_, k_] := b[n, n, k][[2]]; Table[Table[T[n, k], {k, 1, n}], {n, 1, 20}] // Flatten (* Jean-François Alcover, Aug 29 2016, after Alois P. Heinz *)

Formula

G.f. of column k: (x^k/(1-x^k)-x^(k+1)/(1-x^(k+1)))/Product_{j>0}(1-x^j).
Column k is asymptotic to exp(Pi*sqrt(2*n/3)) / (k*(k+1)*Pi*2^(3/2)*sqrt(n)). - Vaclav Kotesovec, May 24 2018

A118806 Triangle read by rows: T(n,k) is the number of partitions of n having k parts of multiplicity 3 (n,k>=0).

Original entry on oeis.org

1, 1, 2, 2, 1, 5, 6, 1, 9, 2, 12, 3, 19, 3, 24, 5, 1, 34, 8, 43, 13, 62, 13, 2, 77, 23, 1, 105, 28, 2, 132, 40, 4, 177, 49, 5, 220, 71, 6, 287, 89, 8, 1, 356, 123, 11, 462, 147, 18, 570, 198, 23, 1, 723, 249, 29, 1, 888, 329, 37, 1, 1121, 400, 50, 4, 1370, 518, 69, 1, 1705, 642, 85
Offset: 0

Views

Author

Emeric Deutsch, Apr 29 2006

Keywords

Comments

T(n,0)=A118807(n). T(n,1)=A118808(n). Row sums yield the partition numbers (A000041). Sum(k*T(n,k), k>=0)=A117524(n) (n>=1).

Examples

			T(12,2) = 2 because we have [3,3,3,1,1,1] and [3,2,2,2,1,1,1].
Triangle starts:
1;
1;
2;
2,  1;
5;
6,  1;
9,  2;
12, 3;

Links

Crossrefs

Cf. A000041, A118807, A118808, A117524, A116595, A116644.

Programs

  • Maple
    g:=product(1+x^j+x^(2*j)+t*x^(3*j)+x^(4*j)/(1-x^j),j=1..35): gser:=simplify(series(g,x=0,35)): P[0]:=1: for n from 1 to 30 do P[n]:=coeff(gser,x^n) od: for n from 0 to 30 do seq(coeff(P[n],t,j),j=0..degree(P[n])) od; # sequence given in triangular form

Formula

G.f.: product(1+x^j+x^(2j)+tx^(3j)+x^(4j)/(1-x^j), j=1..infinity).

A213180 Sum over all partitions lambda of n of Sum_{p:lambda} p^m(p,lambda), where m(p,lambda) is the multiplicity of part p in lambda.

Original entry on oeis.org

0, 1, 3, 7, 16, 28, 59, 91, 170, 269, 450, 655, 1162, 1602, 2527, 3793, 5805, 8034, 12660, 17131, 26484, 37384, 53738, 73504, 114683, 153613, 221225, 313339, 453769, 609179, 927968, 1223909, 1804710, 2522264, 3539835, 4855420, 7439870, 9765555, 14009545
Offset: 0

Views

Author

Alois P. Heinz, Feb 27 2013

Keywords

Examples

			a(6) = 59: (1^6) + (2+1^4) + (2^2+1^2) + (2^3) + (3+1^3) + (3+2+1) + (3^2) + (4+1^2) + (4+2) + (5+1) + (6) = 1+3+5+8+4+6+9+5+6+6+6 = 59.

Links

Crossrefs

Cf. A000070 (Sum 1), A006128 (Sum m), A014153 (Sum p), A024786 (Sum floor(1/m)), A066183 (Sum p^2*m), A066186 (Sum p*m), A073336 (Sum floor(m/p)), A116646 (Sum delta(m,2)), A117524 (Sum delta(m,3)), A103628 (Sum delta(m,1)*p), A117525 (Sum delta(m,2)*p), A197126, A213191.

Programs

  • Maple
    b:= proc(n, p) option remember; `if`(n=0, [1, 0], `if`(p<1, [0, 0],
      add((l->`if`(m=0, l, l+[0, l[1]*p^m]))(b(n-p*m, p-1)), m=0..n/p)))
    end:
a:= n-> b(n, n)[2]:
seq(a(n), n=0..40);
  • Mathematica
    b[n_, p_] := b[n, p] = If[n==0, {1, 0}, If[p<1, {0, 0}, Sum[Function[l, If[m==0, l, l+{0, l[[1]]*p^m}]][b[n-p*m, p-1]], {m, 0, n/p}]]]; a[n_] := b[n, n][[2]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Feb 15 2017, translated from Maple *)

Formula

From Vaclav Kotesovec, May 24 2018: (Start)
a(n) ~ c * 3^(n/3), where
c = 5.0144820680945600131204662934686439430547... if mod(n,3)=0
c = 4.6144523178014379613985400559486878971522... if mod(n,3)=1
c = 4.5237761454818383598444208605033385016299... if mod(n,3)=2
(End)
Showing 1-3 of 3 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.