login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A214808
Kostant partition function P(2*rho) for root system A_n, where 2*rho is the sum of the positive roots.
0
1, 1, 3, 26, 898, 128826, 82462230, 244868962698, 3481404811793685, 242641164238940173212
OFFSET
0,3
COMMENTS
This is the number of ways of writing the vector 2*rho = (n,n-2,...,-n+2,-n) as a sum of the vectors e_i - e_j for 1 <= i < j <= n+1 with nonnegative integer coefficients. - Alejandro H. Morales, May 24 2024
FORMULA
See MathOverflow link for a recurrence.
EXAMPLE
For n = 2, the a(2) = 3 solutions are
0*(1, -1, 0) + 2*(1, 0, -1) + 0*(0, 1, -1) = (2, 0, -2),
1*(1, -1, 0) + 1*(1, 0, -1) + 1*(0, 1, -1) = (2, 0, -2),
2*(1, -1, 0) + 0*(1, 0, -1) + 2*(0, 1, -1) = (2, 0, -2).
MAPLE
multcoeff:=proc(n, f, coeffv, k)
local i, currcoeff;
currcoeff:=f;
for i from 1 to n do
currcoeff:=`if`(coeffv[i]=0, coeff(series(currcoeff, x[i], k), x[i], 0), coeff(series(currcoeff, x[i], k), x[i]^coeffv[i]));
end do;
return currcoeff; end proc:F:=n->mul(mul((1-x[i]*x[j]^(-1))^(-1), j=i+1..n), i=1..n):
a := n -> multcoeff(n+1, F(n+1), [seq(n-2*i, i=0..n)], 2*n):
seq(a(i), i=0..5) # Alejandro H. Morales, May 24 2024
CROSSREFS
Sequence in context: A280222 A354652 A174811 * A378063 A219895 A088730
KEYWORD
nonn,more
AUTHOR
N. J. A. Sloane, Jul 30 2012
EXTENSIONS
a(0) and a(7)-a(9) from Alejandro H. Morales, May 24 2024
STATUS
approved