(Python)
来自sympy import factorint,nextprime
从运算符导入mul
定义P(n):
f=因子(n)
返回排序([f[i]代表f中的i)
定义a046523(n):
x=1
为True时:
如果P(n)==P(x):返回x
其他:x+=1
定义a048673(n):
f=因子(n)
如果n==1,则返回1(1+reduce(mul,[nextprime(i)**f[i]代表f中的i))//2
定义a(n):返回a046523(a048673(n))
打印([a(n)代表范围(1101)中的n)]#因德拉尼尔·戈什,2017年6月12日
|