%I#14 2022年1月12日03:21:19

%S 1,2,1,3,1,4,1,5,2,6,1,7,1,8,2,9,10,10,11,1,3,12,1,13,4,14,3,1,15,1，

%电话：16,5,17,6,18,1,19,7,20,1,21,1,22,1,23,1,24,2,25,8,2,1,26,9,27,10,28，

%U 1,29,1,30,11,31,2,1,1,32,12,12,3,1,33,1,34,13,4,14,35,1,36,4,37,1,38,39,15,40,1,41,2,42,16

%N k的值的数目，1<=k<=N，其中A063994（k）=A063995（N），其中A063 994（N）=Product_{primes p除以N}gcd（p-1，N-1）。

%C A063994的顺序转换。

%H Antti Karttunen，n的表格，n=1..65537的a（n）</a>

%t A063994[n_]：=如果[n==1，1，Times@@GCD[n-1，First/@FactorInteger[n]-1]；

%t模[{b}，b[_]=0；

%t a[n_]：=使用[{t=A063994[n]}，b[t]=b[t]+1]]；

%t数组[a，105]（*_Jean-François Alcover_，2022年1月12日*）

%o（PARI）

%o up_to=65537；

%o A063994（n）={my（f=系数（n）[，1]）；prod（i=1，#f，gcd（f[i]-1，n-1））；}；\\来自A063994

%o ordinal_transform（invec）={my（om=Map（），outvec=vector（length（invec）），pt）；对于（i=1，length（invec），if（mapisdefined（om，invec[i]），pt=mapget（om，invec[i]），pt=0）；outvec[i]=（1+pt）；mapput（om，invec[i]，（1+pt））；outvec；}；

%o v330756=序号变换（向量（up_to，n，A063994（n）））；

%o A330756（n）=v330756【n】；

%Y参考A063994，A209211。

%Y另请参阅A081373、A303756、A330747。

%K nonn公司

%O 1,2号机组

%A _Antti Karttunen，2019年12月30日