登录
A155717号
一些正整数a,b的形式为N=a^2+7b^2的数。
11
8, 11, 16, 23, 29, 32, 37, 43, 44, 53, 56, 64, 67, 71, 72, 77, 79, 88, 92, 99, 107, 109, 112, 113, 116, 121, 127, 128, 137, 144, 148, 149, 151, 161, 163, 172, 176, 179, 184, 191, 193, 197, 200, 203, 207, 211, 212, 224, 232, 233, 239, 253, 256, 259, 261, 263, 268
抵消
1,1
评论
的后续A020670型(允许a和b为零)。
如果N=a^2+7*b^2是一个术语,那么7*N=(7*b)^2+7*a^2也是一个术语。相反,如果7*N是一个项,那么N就是一个项。示例:N=56;N/7=8是项,N*7=7^2+7*7^2是项。序列A154777号,A092572号A154778号具有相同的性质,7分别被素数2、3和5取代-Jerzy R Borysowicz公司,2020年5月22日
数学
选择[Range[300],Reduce[a>0&&b>0&#==a^2+7b^2,{a,b},Integers]=!=错误&](*Jean-François Alcover公司2016年11月17日*)
黄体脂酮素
(PARI)isA155717(n,/*可选的第二个参数允许我们获取其他序列*/c=7)={for(b=1,sqrtint((n-1)\c),issquare(n-c*b^2)&return(1))}
对于(n=1300,isA155717(n)&打印1(n“,”))
(Python)
定义缺陷(极限):
cands=范围(1,int(极限**.5)+2)
nums=[a**2+7*b**2表示a in cands表示b in cands]
返回排序(设置(如果k<=限制,则k为nums))
打印(aupto(268))#迈克尔·布拉尼基2021年8月11日
关键词
容易的,非n
作者
M.F.哈斯勒2009年1月25日
状态
经核准的