OFFSET
0,13
COMMENTS
Clearly a(p)=0, when p is prime.
Number of partitions p of n such that mean(p) = multiplicity(max(p)). For example, a(12) counts these two partitions: 441111, 332211. See the Mathematica program at A240200 for a count of partitions defined in this manner, along with related sequences. - Clark Kimberling, Apr 03 2014
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = coefficient of x^n in expansion of Sum_{d|n} x^(d+n/d)/Product(1-x^k, k=d..n/d). - Vladeta Jovovic, Nov 24 2008
EXAMPLE
a(21) = 4 since property holds for 4 partitions of 21: (7,7,4,3), (7,6,5,3), (7,5,3,3,3), (7,4,4,3,3).
MAPLE
b:= proc(n, i, m) option remember; `if`(n=0, 1,
`if`(i<m, 0, add(b(n-i*j, i-1, m), j=0..n/i)))
end:
a:= n-> add(b(n-d-n/d, d, n/d), d=select(x->
is(x>=sqrt(n)), numtheory[divisors](n))):
seq(a(n), n=0..80); # Alois P. Heinz, Apr 03 2014
MATHEMATICA
f[n_] := Length@ Select[ IntegerPartitions@n, (Length@ # > 1 && Last@# First@# == n) &]; Array[f, 72] (* Robert G. Wilson v, Mar 15 2006 *)
b[n_, i_, m_] := b[n, i, m] = If[n==0, 1, If[i<m, 0, Sum[b[n-i*j, i-1, m], {j, 0, n/i}]]]; a[n_] := Sum[b[n-d-n/d, d, n/d], {d, Select[Divisors[n], # >= Sqrt[n]&]}]; a[0] = 0; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Sep 12 2015, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Resta, Mar 14 2006
EXTENSIONS
More terms from Robert G. Wilson v, Mar 15 2006
STATUS
approved