登录
OEIS由
OEIS基金会的许多慷慨捐赠者
.
提示
(来自的问候
整数序列在线百科全书
!)
A121083号
对于k=-2且0<c<=10^n的本原勾股线样三元组a^2+b^2=c^2+k的数目。
三
3, 19, 182, 1779, 17697, 176794, 1768021, 17676780, 176776851, 1767763756
(
列表
;
图表
;
参考
;
听
;
历史
;
文本
;
内部格式
)
抵消
1,1
评论
第一作者推测a(n)/10^n作为n->inf是1/(4*sqrt(2))=0.17677。。。
链接
n=1..10时的n,a(n)表。
例子
a(1)=3,因为有3个解(a,b,c)as(1,1,2),(3,5,6),(7,7,10),其中0<c<=10^1。
数学
countTriples[m_,k_]:=模块[{c2,c2add,total=0,fax,g},Do[c2=c^2+k;如果[c2<2,继续[]];
c2添加=c2;
而[EvenQ[c2od],c2odd/=2];
如果[c2od==1,如果[OddQ[Log[2,c2],total++];
继续[]];
如果[Mod[c2add,4]==3,继续[]];
g=GCD【c2添加,100947】;
如果[g!=1&g^2!=GCD[c2add,10190296809],继续[]];
fax=映射[{Mod[#[[1]],4],#[[2]]}&,FactorInteger[c2add]];
如果[Apply[Or,Map[#[[1]]==3&&OddQ[#[2]]&,fax]],继续[]];
fax=案例[fax,{1,aa_}:>aa+1];
fax=上限[应用[次数,传真]/2];
总计+=传真;
,{c,m}];
总计](*由Wolfram Research的Daniel Lichtblau提供*)
交叉参考
囊性纤维变性。
A101931号
.
上下文中的序列:
A343672
A161630型
A371320型
*
A343685型
A213533型
A203133型
相邻序列:
A121080型
A121081号
A121082号
*
A121084号
A121085号
A121086号
关键词
更多
,
非n
作者
蒂托·皮耶扎斯三世
2006年8月11日
扩展
发现的前几个术语
蒂托·皮耶扎斯三世
詹姆斯·沃尔德比(j-Waldby(AT)pat7.com)。
随后的术语由Andrzej Kozlowski(akoz(AT)mimuw.edu.pl)、Daniel Lichtblau(danl(AT)wolfram.com)发现。
a(8)-a(10)来自
山口Hiroaki Yamanouchi
2015年10月17日
状态
经核准的