(PARI)
A003961号(n) =我的(f=系数(n));对于(i=1,#f~,f[i,1]=下一素数(f[i、1]+1));factorback(f);\\使用的代码米歇尔·马库斯
(Python)
来自sympy import factorint、nextprime、prod
定义a007814(n):返回1+bin(n-1)[2:].count(“1”)-bin(n)[2:].count
定义a048673(n):
f=因子(n)
如果n==1,则返回1(1+prod(nextprime(i)**f[i]代表f中的i)//2
def a(n):返回a007814(a048673(n))#因德拉尼尔·戈什2017年6月12日
|