%I#25 2023年6月19日08:57:20
%S 1,16288489683504141955224137280410334726975752256,
%电话:11858778808020150599381203234271894799648582235726688,
%电话:990457800726556816837782653376518428624230507392518448661191875230982480827240261878925204256014850656
%N GL中共轭类的数量(N,17)。
%H Alois P.Heinz,n表,n=0..250时的a(n)</a>
%F G.F.:产品{k>=1}(1-x^k)/(1-17*x^k_Alois P.Heinz,2012年11月3日
%p(数字理论):
%p b:=proc(n)b(n):=加法(phi(d)*17^(n/d),d=除数(n))/n-1结束:
%p a:=进程(n)a(n):=`if`(n=0,1,
%p加(加(d*b(d),d=除数(j))*a(n-j),j=1..n)/n)
%p端:
%p序列(a(n),n=0..30);#_Alois P.Heinz,2012年11月3日
%t b[n_]:=总和[EulerPhi[d]*17^(n/d),{d,除数[n]}]/n-1;a[n_]:=a[n]=如果[n==0,1,Sum[Sum[d*b[d],{d,Divisors[j]}]*a[n-j],{j,1,n}]/n];表[a[n],{n,0,30}](*_Jean-François Alcover_,2014年2月17日,在_Alois P.Heinz_*之后)
%o(Magma)/*程序不适用于[1..4]中n>4:*/[1]cat[NumberOfClasses(GL(n,17)):n;
%o(PARI)
%o N=66;x='x+O('x^N);
%o gf=产品(n=1,n,(1-x^n)/(1-17*x^n;
%o v=Vec(gf)
%o/*_Joerg Arndt_,2013年1月24日*/
%Y参考A006951、A006952、A049314、A049315、A049316、A182603、A182604、A182605、A182606、A182607、A182609、A182610、A182611、A182612。
%K nonn公司
%0、2
%A _Klaus Brockhaus,2010年11月23日
%E更多条款,来自_Alois P.Heinz,2012年11月3日
%E MAGMA代码由_Winenzo Librandi_编辑,2013年1月24日
|