|
%I #28 Jun 19 2023 09:05:46
%S 1,10,120,1320,14630,160920,1771440,19485720,214357440,2357931730,
%T 25937408640,285311493720,3138428201160,34522710196920,
%U 379749831637440,4177248147997440,45949729842155150,505447028263532520,5559917313256631160,61159090445821012920
%N Number of conjugacy classes in GL(n,11).
%H Alois P. Heinz, <a href="/A182605/b182605.txt">Table of n, a(n) for n = 0..300</a>
%F G.f.: Product_{k>=1} (1-x^k)/(1-11*x^k). - _Alois P. Heinz_, Nov 03 2012
%p with(numtheory):
%p b:= proc(n) b(n):= add(phi(d)*11^(n/d), d=divisors(n))/n-1 end:
%p a:= proc(n) a(n):= `if`(n=0, 1,
%p add(add(d*b(d), d=divisors(j)) *a(n-j), j=1..n)/n)
%p end:
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Nov 03 2012
%t b[n_] := Sum[EulerPhi[d]*11^(n/d), {d, Divisors[n]}]/n-1; a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d*b[d], {d, Divisors[j]}]*a[n-j], {j, 1, n}]/n]; Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Feb 17 2014, after _Alois P. Heinz_ *)
%o (Magma) N := 300; R<x> := PowerSeriesRing(Integers(), N);
%o Eltseq( &*[ (1-x^k)/(1-11*x^k) : k in [1..N] ] ); // _Volker Gebhardt_, Dec 07 2020
%o (PARI)
%o N=66; x='x+O('x^N);
%o gf=prod(n=1,N, (1-x^n)/(1-11*x^n) );
%o v=Vec(gf)
%o /* _Joerg Arndt_, Jan 24 2013 */
%Y Cf. A006951, A006952, A049314, A049315, A049316, A182603, A182604, A182606, A182607, A182608, A182609, A182610, A182611, A182612.
%K nonn
%O 0,2
%A _Klaus Brockhaus_, Nov 23 2010
%E More terms from _Alois P. Heinz_, Nov 03 2012
|
