|
|
A321233
|
|
a(n) is the number of reflectable bases of the root system of type D_n.
|
|
1
|
|
|
0, 4, 128, 4992, 241664, 14131200, 972521472, 77138231296, 6935178903552, 697359579217920, 77576992194560000, 9461629052252061696, 1255632936007234486272, 180144800985155488448512, 27786422394606966747955200, 4585649599904345055716966400, 806288164205933489807717040128
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The root systems of type D_n are only defined for n >= 4. See chapter 3 of the Humphreys reference. Sequence extended to n=1 using formula/recurrence.
|
|
REFERENCES
|
J. E. Humphreys, Introduction to Lie algebras and representation theory, 2nd ed, Springer-Verlag, New York, 1972.
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: Sum_{m>=2} (1/(4*m)) (Sum_{k>=1} k^(k-1)*(4*x)^k/k!)^m.
|
|
MATHEMATICA
|
Rest[With[{m = 25}, CoefficientList[Series[Sum[Sum[j^(j - 1)*(4*x)^j/j!, {j, 1, m + 1}]^k/(4*k), {k, 2, m}], {x, 0, m}], x]*Range[0, m]!]] (* G. C. Greubel, Dec 09 2018 *)
|
|
PROG
|
(PARI) a(n)={n!*polcoef(sum(m=2, n, (sum(k=1, n, k^(k-1)*(4*x)^k/k!) + O(x^(n-m+2)))^m/(4*m)), n)} \\ Andrew Howroyd, Nov 01 2018
(Magma) m:=25; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (&+[ (&+[ j^(j-1)*(4*x)^j/Factorial(j) :j in [1..m+3]])^k/(4*k) :k in [2..m+2]]) )); [0] cat [Factorial(n+1)*b[n]: n in [1..m-2]]; // G. C. Greubel, Dec 09 2018
(Python)
from math import comb
def A321233(n): return 0 if n<2 else ((sum(comb(n, k)*(n-k)**(n-k)*k**k for k in range(1, (n+1>>1)))<<1) + (0 if n&1 else comb(n, m:=n>>1)*m**n))//n<<(n-1<<1) # Chai Wah Wu, Apr 26 2023
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|