%I#25 2019年1月28日17:14:19
%S 0,1,0,2,2,2,1,0,3,2,2,0,2,5,4,6,1,4,2,3,4,3,0,2,6,2,2,12,6,
%T 10,1,0,4,4,8,4,9,16,2,4,0,3,6,2,6,4,37,6,14,2,4,12,3,12,18,4,10,3,18,
%U 4,5,0,2,12,5,14,6,42,2,28,26,16,3,0
%N A(N,k)是模m的个数,使得k mod m的乘法阶等于N;正方形数组A(n,k),n>=1,k>=1。
%H Alois P.Heinz,反对角线n=1..60</a>
%H维基百科,<a href=“https://en.wikipedia.org/wiki/Miplicative_order(英文)“>乘法顺序</a>
%F A(n,k)=|{m:k模m=n}|的乘法阶。
%F A(n,k)=Sum_{d|n}mu(n/d)*tau(k^d-1),mu=A008683,tau=A000005。
%对于{5,10,16,20,40,80}中的m,e A(4,3)=6:3^4=81=1(mod m)。
%e方阵A(n,k)开始:
%e 0、1、2、2、3、2、4、2。。。
%e 0、1、2、2、5、2、6、4。。。
%e 0、1、2、4、3、2、8、2。。。
%e 0、2、6、4、12、4、26、18。。。
%e 0、1、4、6、9、4、4、六。。。
%e 0、3、10、16、37、10、42、24。。。
%e 0、1、2、6、3、2、12、10。。。
%e 0、4、14、8、28、8、48、72。。。
%p(数字理论):
%pA:=(n,k)->加法(mobius(n/d)*tau(k^d-1),d=除数(n)):
%p序列(序列(A(n,1+d-n),n=1..d),d=1..15);
%t a[n_,k_]:=总和[MoebiusMu[n/d]*除数Sigma[0,k^d-1],{d,除数[n]}];a[1,1]=0;表[a[n-k+1,k],{n,1,12},{k,n,1,-1}]//Flatten(*_Jean-François Alcover_,2012年12月12日*)
%Y列k=1-10给出:A000004、A059499、A059885、A059896、A05988、A0591887、A0598.89、A05.9890、A0591891、A05989。
%Y行n=1-10表示:A000005、A059907、A059808、A05990%、A059910、A059911、A218256、A218257、A218258、A218259。
%Y主对角线表示A252760。
%Y参考A000005、A008683。
%K nonn,表格
%O 1,4型
%A _Alois P.Heinz,2012年6月1日
|