%I#46 2024年4月21日19:14:49

%S 0,0,0,1,0,0，0,0'，0,0，

%温度0,29,0,0,16,0,3,0,11,0,0,74,0,10,10,01,6,0,1,0,98,0,15,0，

%U 3,0,1,0,3,0181,0,0,1,1,03,0,74,1,0,10,0,98,0,,0,16,0,98，0,1,0，0,0220,0,11,29,0,0

%N N的除数（d_i，d_j，d_k，d_m）的有序四元组的个数，使得GCD（d_i,d_j、d_k、d_m）>1。

%C集合{（x，y，z，w）中的元素数：x|n，y|n，z|n，w|n，x<y<z<w，GCD（x，y，z，w）>1}。

%C序列的每一项都无限重复；例如：

%n=1、2、3、4、5、6、7、8、9、10、11、13、14、15、17、19、21、22、23、25、26…时，C a（n）=0。。。（数字k使得k的真除数的乘积小于等于k；即k的除数的积小于等于k^2；参见A007964）。

%当n=12、16、18、20、28、44、45、50、52、63、68、75、76、81、92、98、99。。。（素数的四次幂，或素数与不同素数的平方的乘积；参见A080258）。

%C a（n）=5，n=32，243，3125，16807。。。（素数的五次方；参见A050997）。

%对于n=64、729、15625、117649，C a（n）=15。。。（有7个除数的数字：素数的6次方；见A030516）。

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

%H Antti Karttunen，<a href=“/A336432/A336432.txt”>数据补充：n，a（n）计算n=1.65537</a>

%H<a href=“/index/Eu#epf”>根据n的因式分解中的指数计算序列的索引项</a>

%e a（30）=3，因为30的除数是{1，2，3，5，6，10，15，30}，对于以下三个四倍的除数：（2,6,10,30），（3,6,15,30）和（5,10,15,三十），GCD（d_i，d_j，d_k，d_m）>1。

%p（数字理论）：nn:=100：

%p代表从1到nn的n do：

%p it：=0：d：=除数（n）：n0：=nops（d）：

%i从1到n0-3的p为：

%p表示从i+1到n0-2的j do：

%p表示k从j+1到n0-1 do：

%p代表从k+1到n0的l do：

%如果igcd（d[i]，d[j]，d[C]，d[1]）>1，则为p

%那么p

%p it：=它+1：

%p其他

%p fi:

%日期：

%日期：

%日期：

%日期：

%p打印（`%d，`，it）：

%日期：

%t数组[Count[GCD@@#&/@Subsets[Divisors[#]，{4}]，_？（#>1&）]&，100]（*_Amiram Eldar_，2020年10月31日，在Michael De Vlieger_于A336530*）

%o（PARI）a（n）=my（d=除数（n））；总和（i=1，#d-3，总和（j=i+1，#d-2，总和（k=j+1，#d-1，总和（m=k+1，#d，gcd（[d[i]，d[j]，d[k]，d[m]）>1））；\\_米歇尔·马库斯，2020年10月31日

%o（PARI）a（n）={my（f=因子（n），vp=vecprod（f[，1]），d=除数（vp），res=0）；对于（i=2，#d，res-=二项式（numdiv（n/d[i]），4）*（-1）^omega（d[i]

%Y参见A275387，A336530。

%K nonn，已更改

%O 1,24号机组

%2020年10月5日，拉格瑙市

%E条款由_David A.Corneth修正，2020年10月31日