Bertram Felgenhauer,<a href=“/2013年3月/b321713号_1.txt“>n表,n=1..100000时a(n)(Gheorghe Coserea的前3023项)
检验过的
经核准的
OEIS服务器:已将新的b文件安装为b321713.txt。旧的b文件现在为b321713_1.txt。
提出
编辑
乔格·阿恩特:太好了,谢谢!
Gheorghe公司 科塞雷亚, 伯特伦 费尔根豪尔, <a href=“/A321713飞机/b321713号_1.txt“>n表,n=1时为a(n)。。10000</一> (第一 3023</一> 条款 通过 Gheorghe公司 科塞雷亚)
伯特伦·费尔根豪尔:是的,确实如此;只需要几秒钟。我已经按照建议上传了一个扩展的b文件。
乔格·阿恩特:该代码是否可以计算前10k项?如果是这样,请随意上传一个更大的b文件(这是第一个链接)。
(PARI)
b(n)={满足lambda(k)|n的k个数
my(R=1);
fordiv(n,d,if(i素数(d+1),
我的(e=1);而(n%(d+1)==0,n/=d+1;e++);
如果(d==1&e>1,e++);R*=e+1));
R(右)
};
a(n)=如果(n<=0,0,n==1,2,n%2,0,sumdiv(n,d,moebius(n/d)*b(d)));
向量(100,n,a(n))\\伯特伦·费尔根豪尔2022年3月27日
伯特伦·费尔根豪尔:添加更快的PARI代码,避免枚举所有解决方案。
invlambda(n)={\\A321712飞机A270562型
囊性纤维变性。A002322号,A321712飞机, A270562型, A321714飞机.