登录
OEIS由OEIS基金会的许多慷慨捐赠者.

 


A362961型
a(n)=和{b=0..层(sqrt(n)),n-b^2是平方}b。
2
1, 1, 0, 2, 3, 0, 0, 2, 3, 4, 0, 0, 5, 0, 0, 4, 5, 3, 0, 6, 0, 0, 0, 0, 12, 6, 0, 0, 7, 0, 0, 4, 0, 8, 0, 6, 7, 0, 0, 8, 9, 0, 0, 0, 9, 0, 0, 0, 7, 13, 0, 10, 9, 0, 0, 0, 0, 10, 0, 0, 11, 0, 0, 8, 20, 0, 0, 10, 0, 0, 0, 6, 11, 12, 0, 0, 0, 0, 0, 12, 9, 10, 0
抵消
1,4
评论
如果n in,a(n)=0A022544号.
如果n in,a(n)>0A001481号.
链接
斯特凡诺·斯佩齐亚,n=1..10000时的n,a(n)表
数学
a[n_]:=总和[b Boole[IntegerQ[Sqrt[n-b^2]],{b,0,Floor[Sqrt[n]]}];数组[a,83](*斯特凡诺·斯佩齐亚2023年5月15日*)
黄体脂酮素
(Python)
从gmpy2导入*
a=λn:总和([b对于范围(0,isqrt(n)+1)中的b,如果是平方(n-(b*b))])
打印([a(n)代表范围(1,84)中的n)]
(Python)
从sympy导入除数
来自sympy.solvers.diophantine.diophatine进口玉米
定义362961美元(n) :
c=0
对于除数(n)中的d:
如果(k:=d**2)>n:
打破
q、 r=divmod(n,k)
如果不是r:
c+=cornacchia(1,1,q)或[]中a的总和(d*(a[0]+(a[1]if a[0]!=a[1]else 0))
返回c#柴华武2023年5月15日
(PARI)a(n)=总和(b=0,平方(n),如果(发行方(n-b^2),b))\\米歇尔·马库斯2023年5月16日
交叉参考
囊性纤维变性。A022544号,A001481号.
囊性纤维变性。A143574号(b^2之和),A000925号.
关键词
非n,
作者
达里奥·克拉维乔,2023年5月10日
状态
经核准的

查找|欢迎光临|维基|注册|音乐|地块2|演示|索引|浏览|网络摄像头
贡献新序列。或评论|格式|样式表|变换|超级搜索|最近
OEIS社区|维护人OEIS基金会。

许可协议、使用条款、隐私政策。.

上次修改时间:2024年9月23日12:43 EDT。包含376164个序列。(在oeis4上运行。)