%I#73 2021年6月13日13:35:11
%S 1841832693832046137831916552903831908477106883,
%电话83190743068779769831907383078281024371831907373418800027750413,
%电话:83190737272244910014741448783190735258907312489948783173583190732581823023465031521920149831907372 580768386561159866725731911
%N gcd(x,y,z)=1且1<={x,y、z}<=10^N的有序三元组(x,y,z)的数目。
%D Joachim von zur Gathern和Jürgen Gerhard,《现代计算机代数》,剑桥大学出版社,2003年第二版,第53-54页。
%H Chai Wah Wu,<a href=“/A342841/b342841.txt”>n表,n=0..15的a(n)</a>
%H Karl-Heinz Hofmann,n=1立方体的动画。
%F Lim_{n->无穷大}a(n)/10^(3*n)=1/zeta(3)=1/Apéry常数。
%F a(n)=A071778(10^n)。
%e为了可视化,集合(x,y,z)内接在立方体矩阵中。
%e“o”代表gcd=1。
%e“.”代表gcd>1。
%e、。
%e对于n=1,立方体矩阵的大小为10 X 10 X 10:
%e、。
%电子/::::
%e/100总和(z=1)
%e z=7 |/1 2 3 4 5 6 7 8 9 10|
%e--+--------------------75总和(z=2)
%e 1/|o o o o o o o o o o o o o 10|
%e 2/|o o o o oo o o o o o o 10 91总和(z=3)
%2010年5月|
%e z=8|/1 2 3 4 5 6 7 8 9 10 10 75总和(z=4)
%e--+--------------------10/
%e 1/|o o o o oo o o o o o o 10 10 96总和(z=5)
%e 2/月。o、。o、。o、。o、。5 9 /
%e/10 10 67总和(z=6)
%e z=9 |/1 2 3 4 5 6 7 8 9 10 5 10/
%e--+--------------------10 10/
%e 1/|o o o o oo o o o o o o 10 5--/
%e 2/|o o o o oo o o o o o o 10 10 99总和(z=7)
%e/7 5/
%e z=10 |/1 2 3 4 5 6 7 8 9 10 10/
%e--+--------------------10 5/
%e 1|o o o o o o o o o o o o 10 7--/
%e 2 | o。o、。o、。o、。o、。5 10 75总和(z=8)
%e 3 | o o o o o o o o 0 o o o 10 10/
%e 4 | o。o、。o、。o、。o、。5 7 /
%e 5 | o o o o。哦哦。8 10 /
%e 6 | o。o、。o、。o、。o、。5 --/
%e 7 | o o o o o o o o 0 o o o 10 91总和(z=9)
%e 8 | o。o、。o、。o、。o、。5 /
%e 9 | o o o o o o o o 0 o o o 10/
%e 10 | o。o。o、。o、。4 /
%e(电子)--/
%e 72总和(z=10)
%电子/
%e(电子)|
%e(电子)------
%e 841立方和(z=1..10)
%o(Python)
%o导入数学
%o表示范围(0,10)内的n:
%o计数器=0
%o对于范围(1,pow(10,n)+1)内的x:
%o对于范围(1,pow(10,n)+1)中的y:
%o对于范围(1,pow(10,n)+1)中的z:
%o如果math.gcd(math.gccd(y,x),z)==1:
%o计数器+=1
%o打印(n,计数器)
%Y参考A342586(用于10^n X 10^n)、A018805、A002117(zeta(3))、A071778。
%Y相关的k元组计数:
%Y对:A018805、A342632、A342586;
%Y三联件:A071778、A342935、A342841;
%Y四线:A082540、A343527、A343193;
%Y 5元组:A343282;
%Y 6元组:A343978、A344038.-_N.J.A.Sloane,2021年6月13日
%K nonn,硬
%0、2
%A _Karl-Heinz Hofmann,2021年3月24日
%E a(5)-a(10)摘自Hhugo Pfoertner_,2021年3月25日
%E a(11)摘自Hhugo Pfoertner_,2021年3月26日
%E a(12),来自布鲁斯·加内尔,2021年3月29日