显示找到的2个结果中的1-2个。
第页1
a(n)=n*(n+1)/2 mod(n的素因子之和(具有多重性))。
+10 4
1, 0, 2, 0, 1, 0, 0, 3, 6, 0, 1, 0, 6, 0, 0, 0, 3, 0, 3, 1, 6, 0, 3, 5, 6, 0, 10, 0, 5, 0, 8, 1, 6, 6, 6, 0, 6, 12, 6, 0, 3, 0, 0, 1, 6, 0, 10, 7, 3, 6, 1, 0, 0, 4, 10, 3, 6, 0, 6, 0, 6, 1, 4, 3, 3, 0, 15, 23, 7, 0, 0, 0, 6, 3, 5, 15, 3, 0, 3, 9, 6, 0, 0, 3, 6, 20, 6, 0, 0, 6, 12, 19, 6, 0, 2, 0
评论
如果n是奇数素数,a(n)=0。
如果n是素数,a(n^2)=n。
MAPLE公司
seq((n*(n+1)/2)mod加法(t[1]*t[2],t=ifactors(n)[2]),n=2..100);
数学
a[n_]:=Mod[n*(n+1)/2,加上@@Times@@FactorInteger[n]];数组[a,100,2](*阿米拉姆·埃尔达尔2022年4月15日*)
黄体脂酮素
(Python)
来自sympy导入因子
定义a(n):返回(n*(n+1)//2)%sum(p*e代表因子(n).items()中的p,e)
打印([a(n)代表范围(2,98)中的n])#迈克尔·布拉尼基,2022年4月24日
对k进行编号,使第k个三角形数==1修改k的整数对数。
+10 2
2, 6, 12, 21, 33, 45, 52, 63, 136, 162, 201, 205, 208, 225, 245, 253, 301, 304, 344, 441, 494, 531, 550, 637, 697, 720, 742, 806, 901, 910, 918, 1078, 1233, 1242, 1274, 1333, 1376, 1540, 1566, 1573, 1625, 1680, 1695, 1792, 1863, 1909, 2025, 2041, 2107, 2295, 2466, 2497, 2774, 2896, 2926, 2965
例子
a(3)=12=2*2*3是一个项,因为12*13/2=78==1(mod 2+2+3=7)。
MAPLE公司
过滤器:=proc(n)局部t;(n*(n+1)/2)mod add(t[1]*t[2],t=ifactors(n)[2])=1结束进程:
选择(过滤器,[2..3000]美元);
数学
选择[Range[3000],Mod[#*(#+1)/2,Plus@@Times@@FactorInteger[#]]==1&](*阿米拉姆·埃尔达尔2022年4月14日*)
搜索在0.004秒内完成
|