%我#16 2022年2月28日01:53:13

%S 1,9,9999099909989199990099989019999891099999890109999989010，

%电话99999 8891109999998901099999999889020099999998891100，

%电话：9999999 88901199999999988902090999998889120909999999998890119099999988990

%N分区总和（q=10），见第一条注释。

%C集q=10和f（m）=q^（m-1）*（q-1），则a（n）是所有乘积Product_{k=1..L}f（m_k）上n的所有分区P的和，其中L是分区P=[P_1^m_1，P_2^m_2，…，P_L^m_L]中不同部分的数量。

%C将q设置为素数幂，得出序列“GL（n，q）中的共轭类数”：

%Cq=3:A006952，q=4:A049314，q=5:A049315，q=7:A049316，q=8:A182603，

%Cq=9:A182604，q=11:A182605，q=13:A182606，q=16:A182607，q=17:A182608，

%Cq=19:A182609，q=23:A182610，q=25:A182611，q=27:A182612。

%C序列，其中q不是主幂：

%Cq=6:A221578，q=10:A221579，q=12:A221580，

%Cq=14:A221581，q=15:A221582，q=18:A221583，q=20:A221584。

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

%p（数字理论）：

%p b：=程序（n）b（n）：=加法（φ（d）*10^（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，2013年1月25日

%t b[n_]：=总和[EulerPhi[d]*10^（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（PARI）

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

%o gf=产品（n=1，n，（1-x^n）/（1-10*x^n；

%o v=Vec（gf）

%K nonn公司

%0、2

%A _Joerg Arndt_，2013年1月20日