提出
经核准的
编辑
定义P(n):
f=因子(n)
定义 P(P)(n个): 返回已排序([保理商(n个)(f)[i] 对于我来说保理商(n个)(f)])
(Python)
来自sympy import factorint,totiten,gcd
定义T(n,m):返回((n+m)**2-n-3*m+2)/2
def P(n):返回排序的([factorint(n)[i]for i in factorent(n)])
定义a046523(n):
x=1
为True时:
如果P(n)==P(x):返回x
其他:x+=1
定义a(n):返回T(a046523(n),n/gcd(n,totient(n))#因德拉尼尔·戈什,2017年5月5日
复合过滤器:a(n)=T(A046523号(n) ,A109395号(n) ),其中T(n,k)是序列A000027号用作配对函数.
表[(2+(#1+#2)^2-#1-3#2)/2&@@{Times@@MapIndexed[Prime[First@#2]^#1&,Sort[FactorInteger[n][[All,-1]],Greater]]-Boole[n==1],分母[EulerPhi[n]/n]},{n,73}](*迈克尔·德弗利格2017年5月4日*)