%I#28 2023年6月19日12:13:20
%编号:1,246241560039060097649762441400006103499376152587874400,
%电话:381469685940095367431234400238418578084440059604644765235024,
%电话:14901161191304700003725290298436486000093132257460912111062423283064365380605500600582076609134515127375600
%N GL中的共轭类的数量(N,25)。
%H Alois P.Heinz,n表,n=0..250时的a(n)</a>
%F G.F.:产品{k>=1}(1-x^k)/(1-25*x^k_Alois P.Heinz,2012年11月3日
%p(数字理论):
%p b:=程序(n)b(n):=加法(φ(d)*25^(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]*25^(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,25)):n;
%o(PARI)
%o N=66;x='x+O('x^N);
%o gf=产品(n=1,n,(1-x^n)/(1-25*x^n;
%o v=Vec(gf)
%o/*_Joerg Arndt_,2013年1月24日*/
%Y请参阅A006951、A006952、A049314、A04931.5、A049316、A182603、A182601、A1182605、A182066、A182670、A182680、A182 609、A182610、A182 612。
%K nonn公司
%0、2
%A _Klaus Brockhaus,2010年11月23日
%E更多条款,来自_Alois P.Heinz,2012年11月3日
%E MAGMA代码由_Vincozo Librandi_编辑,2013年1月23日
|