登录
A015634号
没有全局因子的[1..n]中整数的有序四元组数。
7
1, 4, 13, 29, 63, 106, 189, 289, 444, 626, 911, 1203, 1657, 2130, 2766, 3462, 4430, 5359, 6688, 7992, 9670, 11405, 13704, 15840, 18730, 21548, 25037, 28521, 33015, 37067, 42522, 47690, 53940, 60108, 67760, 74748, 83886, 92433, 102629, 112469
抵消
1,2
链接
配方奶粉
通用公式:(1/(1-x))*Sum_{k>=1}mu(k)*x^k/(1-x^k)^4-伊利亚·古特科夫斯基2020年2月14日
a(n)=n*(n+1)*(n+2)*(n+3)/24-求和{j=2..n}a(楼层(n/j))=A000332号(n+3)-总和{j=2..n}a(楼层(n/j))-柴华武2021年4月18日
黄体脂酮素
(Python)
从functools导入lru_cache
@lru_cache(最大大小=无)
定义A015634号(n) :
如果n==0:
返回0
c、 j=1,2
k1=无
当k1>1时:
j2=无/无k1+1
c+=(j2-j)*A015634号(k1)
j、 k1=j2,n//j2
返回n*(n+1)*(n+2)*(n+3)//24-c+j-n#柴华武2021年4月18日
(PARI)a(n)=总和(k=1,n,sumdiv(k,d,moebius(k/d)*二项式(d+2,3))\\Seiichi Manyama先生2021年6月12日
(PARI)a(n)=二项式(n+3,4)-和(k=2,n,a(n\k))\\Seiichi Manyama先生2021年6月12日
(PARI)我的(N=66,x='x+O('x^N));Vec(总和(k=1,N,moebius(k)*x^k/(1-x^k)^4)/(1-x))\\Seiichi Manyama先生2021年6月12日
关键词
非n
状态
经核准的