OFFSET
0,3
COMMENTS
a(n) is the expected value of |tr(U)|^{2n} where U is drawn uniformly at random from the 120-element binary icosahedral group, viewed as a subgroup of SU(2) (or equivalently, the unit quaternion group). Note that |tr(U)| takes values in 0, phi^{-1}, 1, phi, 2 (with phi the golden ratio) with probabilities 1/4, 1/5, 1/3, 1/5, 1/60 respectively.
Is a reasonably good match with A000108 which corresponds to the case where U is drawn from all of SU(2) with the Haar distribution.
LINKS
Colin Barker, Table of n, a(n) for n = 0..1000
Terence Tao, A function field analogue of Riemann zeta statistics
Index entries for linear recurrences with constant coefficients, signature (8,-20,17,-4).
FORMULA
a(n) = 0^n/4 + (phi^{2n} + phi^{-2n})/5 + 1/3 + 4^n/60.
From Colin Barker, May 20 2019: (Start)
G.f.: (1 - 7*x + 14*x^2 - 8*x^3 + x^4) / ((1 - x)*(1 - 4*x)*(1 - 3*x + x^2)).
a(n) = 8*a(n-1) - 20*a(n-2) + 17*a(n-3) - 4*a(n-4) for n>4.
a(n) = (20 + 4^n + 3*2^(2-n)*((3-sqrt(5))^n + (3+sqrt(5))^n)) / 60 for n>0.
(End)
EXAMPLE
a(1) = (phi^2 + phi^{-2})/5 + 1/3 + 4/60 = 3/5 + 1/3 + 1/15 = 1.
PROG
(PARI) Vec((1 - 7*x + 14*x^2 - 8*x^3 + x^4) / ((1 - x)*(1 - 4*x)*(1 - 3*x + x^2)) + O(x^30)) \\ Colin Barker, May 20 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Terence Tao, May 20 2019
STATUS
approved