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

 

标志
提示
(来自的问候整数序列在线百科全书!)
A015616号 1<=i<j<k<=n且gcd(i,j,k)=1的三元组(i,j,k)的数量。 5
0, 0, 1, 4, 10, 19, 34, 52, 79, 109, 154, 196, 262, 325, 409, 493, 613, 712, 865, 997, 1171, 1336, 1567, 1747, 2017, 2251, 2548, 2818, 3196, 3472, 3907, 4267, 4717, 5125, 5665, 6079, 6709, 7222, 7858, 8410, 9190, 9748, 10609, 11299, 12127 (列表;图表;参考;;历史;文本;内部格式)
抵消
1,4个
链接
阿洛伊斯·海因茨,n=1..10000时的n,a(n)表
配方奶粉
a(n)=(A071778号(n) -3个*A018805型(n) +2)/6-弗拉德塔·乔沃维奇2004年12月1日
a(n)=和{i=1..n}A000741号(i) ●●●●-阿洛伊斯·海因茨2011年2月8日
对于n>1,a(n)=n(n-1)(n-2)/6-和{j=2..n}a(楼层(n/j))=A000292号(n-2)-总和{j=2..n}a(楼层(n/j))-柴华武2021年3月30日
例子
对于n=6,a(6)=19的解是二项式(6,3)=(6*5*4)/(1*2*3)=20个可能的三元组减去GCD=2的三元组(2,4,6)。
MAPLE公司
f: =程序(n)局部i,j,k,t1,t2,t3;t1:=0;对于i从1到n-2 do对于j从i+1到n-1 do t2:=gcd(i,j);对于从j+1到n的k,做t3:=gcd(t2,k);如果t3=1,则t1:=t1+1;fi;od:od:od:t1;结束;
#基于Moebius变换的程序A000741号:
带有(数字理论):
b: =proc(n)选项记忆;
加法(mobius(n/d)*(d-2)*(d-1)/2,d=除数(n))
结束时间:
a: =proc(n)选项记忆;
b(n)+`if`(n=1,0,a(n-1))
结束时间:
seq(a(n),n=1..100)#阿洛伊斯·海因茨2011年2月8日
数学
a[n_]:=(cnt=0;Do[cnt+=Boole[GCD[i,j,k]==1],{i,1,n-2},{j,i+1,n-1},},j+1,n}];cnt);表[a[n],{n,1,45}](*Jean-François Alcover公司2013年3月5日*)
黄体脂酮素
(PARI)打印1(c=0);对于(k=1,99,对于(j=1,k-1,gcd(j,k)==1&&(c+=j-1)&&next;对于(i=1,j-1,gcd([i,j,k])>1|c++));打印1(“,”c))
(Python)
从functools导入lru_cache
@lru_cache(最大大小=无)
定义A015616号(n) :
如果n<=1:
返回0
c、 j=n*(n-1)*(n-2)//6,2
k1=无
当k1>1时:
j2=无/无k1+1
c-=(j2-j)*A015616号(k1)
j、 k1=j2,n//j2
返回c#柴华武2021年3月30日
交叉参考
关键词
非n
作者
状态
经核准的

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

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

上次修改时间:美国东部夏令时2024年4月18日15:05。包含371780个序列。(在oeis4上运行。)