| 数据
|
1, 2, 3, 5, 6, 7, 8, 9, 12, 15, 16, 18, 19, 22, 23, 24, 25, 26, 27, 31, 35, 36, 37, 41, 46, 53, 58, 72, 80, 87, 94, 111, 121, 130, 149, 159, 183, 217, 282, 311, 388, 422, 624, 935, 1171, 1323, 1578, 1600, 2554, 12348, 14842, 17461
|
| MAPLE公司
|
带有(padic,ordp):
diamond:=proc(n)#奇数整数n局部i,j,s的色调菱形;s:={};对于i从1乘2到n do,对于j从1乘2中到n dos:=s并集{r2d2(i/j)}odod;排序(转换,列表)结束:
r2d2:=proc(q)#有理数的倍频程减缩q 2^(-地板(evalf(ln(q)/ln(2)))*q结束:
plim:=proc(q)#有理数q的素数极限局部r,i,p;r:=1;i:=0;当不是(r=q)时,做i:=i+1;p:=ithprime(i);r:=r*p^ordp(q,p)od;i结束:
vai:=proc(n,i)#第i素数通过专利val对n轮(evalf(n*ln(ithprime(i))/ln(2))end的映射:
via:=proc(n,l)#长度为l的n的专利值局部i,v;对于i从1到l do v[i]:=vai(n,i)od;convert(转换(v,数组),列表)结束:
h:=proc(n,q)#区间q的专利val n映射,如果q=1,则返回(0)fi;dotprod(vec(q),via(n,plim(q)))端:
condit:=proc(n,s)#关于调和集s局部i,d的edon的不同一致性;对于i从1到nops,如果不是h(n,s[i])=round(n*l2(s[i]]),则返回(false)fiod;对于i从1到nops,做d[i]:=h(n,s[i])od;如果不是nops(convert(d,set))=nops(s),则返回(false)fi;返回(真)结束:
condl:=proc(n)#edo n局部c的最高独立奇数限制一致性;c:=3;而condi(n,diamond(c))做c:=c+2od;c-2端:
|