A118199 Number of partitions of n having no parts equal to the size of their Durfee squares.
1, 0, 1, 1, 1, 1, 2, 3, 5, 7, 10, 13, 18, 23, 31, 40, 53, 68, 89, 113, 146, 184, 234, 293, 369, 458, 572, 706, 874, 1073, 1320, 1611, 1970, 2393, 2909, 3518, 4255, 5122, 6167, 7394, 8862, 10585, 12637, 15038, 17886, 21213, 25141, 29723, 35112, 41383, 48737, 57278
Offset: 0
Keywords
Examples
a(7) = 3 because we have [7] with size of Durfee square 1, [4,3] with size of Durfee square 2 and [3,3,1] with size of Durfee square 2.
Links
- Vaclav Kotesovec, Table of n, a(n) for n = 0..5000 (terms 0..1000 from Alois P. Heinz)
Programs
-
Maple
g:=1+sum(x^(k^2+k)/(1-x^k)/product((1-x^i)^2,i=1..k-1),k=1..20): gser:=series(g,x=0,60): seq(coeff(gser,x,n),n=0..54); # second Maple program:: b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0, b(n, i-1)+`if`(i>n, 0, b(n-i, i)))) end: a:= n-> add(add(b(k, d) *b(n-d*(d+1)-k, d-1), k=0..n-d*(d+1)), d=0..floor(sqrt(n))): seq(a(n), n=0..70); # Alois P. Heinz, Apr 09 2012
-
Mathematica
b[n_, i_] := b[n, i] = If[n == 0, 1, If[i<1, 0, b[n, i-1] + If[i>n, 0, b[n-i, i]]]]; a[n_] := Sum[Sum[b[k, d]*b[n-d*(d+1)-k, d-1], {k, 0, n-d*(d+1)}], {d, 0, Floor[Sqrt[n]]}]; Table[a[n], {n, 0, 70}] (* Jean-François Alcover, May 22 2015, after Alois P. Heinz *) pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]]; conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]]; Table[Length[Select[IntegerPartitions[n],pq[#]>0&&pq[conj[#]]==0&]],{n,0,30}] (* a(0) = 0, Gus Wiseman, May 21 2022 *)
Formula
G.f.: 1+sum(x^(k^2+k)/[(1-x^k)*product((1-x^i)^2, i=1..k-1)], k=1..infinity).
a(n) ~ exp(Pi*sqrt(2*n/3)) / (16*n*sqrt(3)). - Vaclav Kotesovec, Jun 12 2025
Comments