OFFSET
0,6
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Alois P. Heinz)
FORMULA
a(n) ~ exp(Pi*sqrt(2*n/3)) / (4*sqrt(3)*n). - Vaclav Kotesovec, Jan 28 2022
EXAMPLE
The qualifying partitions of 8 are 521, 431, 332, 421, 3221, 32111, 22211, 221111, 2111111, so that a(8) = 9.
MAPLE
g:= proc(n, i) option remember; `if`(n=0, 1,
`if`(i<1, 0, add(g(n-i*j, i-1), j=0..n/i)))
end:
b:= proc(n, i, l) option remember; `if`(n=0 or i<1, 0,
b(n, i-1, 0) +add(`if`(i+1=l, g(n-i*j, i-1),
b(n-i*j, i-1, i)), j=1..n/i))
end:
a:= n-> b(n$2, 0):
seq(a(n), n=0..60); # Alois P. Heinz, Feb 14 2014
MATHEMATICA
Map[Length[Cases[Map[Differences[DeleteDuplicates[#]] &, IntegerPartitions[#]], {___, -1, ___}]] &, Range[50]] (* Peter J. C. Moses, Feb 09 2014 *)
g[n_, i_] := g[n, i] = If[n==0, 1, If[i<1, 0, Sum[g[n-i*j, i-1], {j, 0, n/i}]]]; b[n_, i_, l_] := b[n, i, l] = If[n==0 || i<1, 0, b[n, i-1, 0] + Sum[If[i+1 == l, g[n-i*j, i-1], b[n-i*j, i-1, i]], {j, 1, n/i}]]; a[n_] := b[n, n, 0]; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Sep 01 2016, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Feb 11 2014
STATUS
approved