%I#22 2022年9月8日08:45:55

%S 1,26728196565314141434815238741976010460332792282429516096，

%电话76255969338902058911315435525559060551656248150094635282119528，

%电话：405255512616676888109418913111007878429543127065399715971847976643076861647780830

%N GL中的共轭类数（N，27）。

%H Alois P.Heinz，n的表格，n=0..250时的a（n）</a>

%F G.F.：产品{k>=1}（1-x^k）/（1-27*x^k_Alois P.Heinz，2012年11月3日

%p（数字理论）：

%p b：=程序（n）b（n）：=加法（φ（d）*27^（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..20）；#_Alois P.Heinz，2012年11月3日

%t b[n_]：=总和[EulerPhi[d]*27^（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，20}]（*_Jean-François Alcover_，2014年2月17日，在_Alois P.Heinz_*之后）

%o（Magma）/*程序不适用于[1..4]中n>4:*/[1]cat[NumberOfClasses（GL（n，27））：n；

%o（PARI）

%o N=66；x='x+O（'x^N）；

%o gf=产品（n=1，n，（1-x^n）/（1-27*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 611。

%K nonn公司

%0、2

%A _Klaus Brockhaus，2010年11月23日

%E更多条款，来自_Alois P.Heinz，2012年11月3日

%E MAGMA代码由_Winenzo Librandi_编辑，2013年1月24日