显示找到的4个结果中的1-4个。
第页1
让我来=A344005型(n) =最小m,使得n除以m*(m+1);a(n)=gcd(n,m)。
+10 26
1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 5, 1, 1, 2, 1, 4, 3, 2, 1, 8, 1, 2, 1, 7, 1, 5, 1, 1, 11, 2, 7, 4, 1, 2, 3, 5, 1, 6, 1, 11, 9, 2, 1, 3, 1, 2, 17, 4, 1, 2, 5, 7, 3, 2, 1, 15, 1, 2, 9, 1, 5, 11, 1, 4, 23, 14, 1, 8, 1, 2, 3, 19, 7, 6, 1, 5, 1, 2, 1, 4, 17, 2, 29, 8, 1, 9, 13, 23, 3, 2, 19, 32, 1, 2, 11, 4, 1, 17
评论
当n/2在A344000型。这是真的,但本质上是微不足道的,无法深入了解这两个序列。
经验:对于n>=3,a(n)<=n/3,a(n)=n/3当n是3*{2^奇数,素数==-1模6}。
如果n是素数幂,则a(n)=1-柴华武,2022年6月1日
(结束)
MAPLE公司
ans:=[];
对于从1到40的n,做t1:=Findm(n)[1];ans:=[op(ans),igcd(n,t1)];日期:
ans;
数学
smd[n_]:=模块[{m=1},而[Mod[m(m+1),n]=0,m++];GCD[n,m]];数组[smd,110](*哈维·P·戴尔2022年1月7日*)
黄体脂酮素
(PARI)f(n)=我的(m=1);while((m*(m+1))%n,m++);米\\A344005型
(Python 3.8+)
从数学导入gcd,prod
从itertools导入组合
来自sympy导入因子
从sympy.theory.模块导入crt
如果n==1:
返回1
plist=元组(p**q代表p,q在因子(n).items()中)
如果len(plist)==1,则返回1,否则gcd(n,int(min(min(crt((m,n//m),(0,-1#柴华武,2022年6月1日
4, 7, 13, 16, 19, 25, 31, 37, 43, 49, 61, 64, 67, 73, 79, 97, 103, 109, 121, 127, 139, 151, 157, 163, 169, 181, 193, 199, 211, 223, 229, 241, 256, 271, 277, 283, 289, 307, 313, 331, 337, 343, 349, 361, 367, 373, 379, 397, 409, 421, 433, 439, 457, 463, 487, 499
评论
数字k不是3的幂,因此x^2+x+1(或x^2-x+1)在GF(k)上是可约的-宋嘉宁2019年9月24日
黄体脂酮素
(PARI)isok(n)=isprimepower(n)&&((n%3)==1)\\米歇尔·马库斯2019年9月24日
正方形数组,其中第n行列出了所有非负数kA345992型(k) =n,通过下降的反对偶来读取。
+10 5
1, 2, 6, 3, 10, 12, 4, 14, 21, 20, 5, 18, 39, 36, 15, 7, 22, 48, 52, 30, 42, 8, 26, 57, 68, 40, 78, 28, 9, 34, 75, 84, 55, 114, 35, 24, 11, 38, 93, 100, 65, 150, 56, 72, 45, 13, 46, 111, 116, 80, 186, 77, 88, 63, 110, 16, 50, 129, 148, 115, 222, 105, 136, 90, 310, 33, 17, 54, 147, 164, 130, 258, 133, 152, 126, 410, 44, 156
评论
数组由(n,k)=(1,1),(1,2),(2,1)。。。其中A(n,k)是第k个解x到A345992型(x) =个。
例子
阵列的左上角15x9:
否| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
----+--------------------------------------------------------------------------
1 | 1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25,
2 | 6, 10, 14, 18, 22, 26, 34, 38, 46, 50, 54, 58, 62, 74, 82,
3 | 12, 21, 39, 48, 57, 75, 93, 111, 129, 147, 183, 192, 201, 219, 237,
4 | 20, 36, 52, 68, 84, 100, 116, 148, 164, 196, 212, 228, 244, 292, 324,
5 | 15, 30, 40, 55, 65, 80, 115, 130, 155, 180, 205, 215, 230, 255, 265,
6 | 42, 78, 114, 150, 186, 222, 258, 294, 366, 402, 438, 474, 582, 618, 654,
7 | 28, 35, 56, 77, 105, 133, 154, 175, 203, 224, 273, 301, 329, 350, 371,
8 | 24, 72, 88, 136, 152, 200, 216, 264, 328, 344, 392, 456, 472, 536, 584,
9 | 45, 63, 90, 126, 144, 171, 207, 225, 252, 288, 333, 369, 387, 414, 450,
黄体脂酮素
(PARI)
up_to=105;
A345992型(n) =对于(m=1,oo,如果((m*(m+1))%n==0,返回(gcd(n,m)));
memoA354930sq=地图();
A354930sq(n,k)={my(v=0);如果(!mapisdefined(memoA354930 sq,[n,k-1],&v),如果(1==k,v=0,v=A354930msq(n,k-1)));对于(i=1+v,oo,如果(A345992型(i) ==n,地图输入(备忘录354930sq,[n,k],i);返回(i));};
A354930list(up_to)={my(v=向量(up_to),i=0);对于(a=1,oo,对于(col=1,a,i++;如果(i>up_to,返回(v));v[i]=A354930平方码(col,(a-(col-1))));(v);};
v354930=A354930列表(up_to);
a(n)=1,如果n是与1(mod 3)同余的素数幂的三倍,否则为0。
+10 三
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1
数学
a[n_]:=如果[Mod[n,9]==3&&PrimePowerQ[n/3],1,0];数组[a,100](*阿米拉姆·埃尔达尔,2022年6月15日*)
黄体脂酮素
(PARI)A354983型(n) =((3==(n%9))&&初始功率(n/3));
搜索在0.006秒内完成
|