#来自在线整数序列百科全书的问候!本次搜索:id:a09912 显示1-1的1 ;%I a089912;%S a089912 1,3,4,4,5,5,6,11,8,8,9,10,16,12,12,18,14,21,23,17,18,18,26,20,28,30,30,31,31,24,33,26,36,28,37,30,57,32,33,45,45,46,47,45,45,45,47,45,45,45,47,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,45,38,51,,58,61,48,62, %U a089912 50,66,67,64,54,71,70,68,74,76,60100,62,81,77,65,82113,68,82,89118,72 %N a089912 a(N)给出了n和m的公酉因子,其中m从1到n。 %H A089912 Amiram Eldar,n=1..10000的n,a(n)表%A089912 udiv[n[n]以“:=选择[除数[n]、GCD[#、n/#]==1&];a[n Uu]:=模块[{d=udiv[n]},总和[长度[交叉点[d,udiv[k]]],{k,1,n}]]];阵列[a,100](*阿米拉姆Eldar_,2019年8月10日*);%o a0889912(PARI)a(n)={sdivn=Set();fordiv(n,d,如果n,d,如果他们的话,n,d,如果{k,d,d,如果他们的话,d,如果他们的话,d,如果他们的(GCD(d,n/d)==1,sdivn=setunion(sdivn,Set(d)));s=0;对于(m=1,n,sdivm=Set();fordiv(m,d,如果(gcd(d,m/d)==1,sdivm=setunion(sdivm,SETIVM(d))));s+=长度(setintersect(sdivn,sdivm)););返回(s);\\\\\米歇尔马库斯年7月15日,2013年7月15日,;%Y A089912 Cf.A034444,A034448.2013年7月15日;%Y A089912比较A034444,A034448.;%K A089912 easy,Non;%O A089912 1,2;%A A089912 %A A089912 %A A089912 %A A089912 Naoho Naoho 10;#内容根据OEIS最终用户许可协议提供:http://OEIS.org/License