with（padic，ordp）：菱形：=proc（n）#奇数整数n局部i，j，s的色调菱形；s：={}；对于i从1乘2到n do，对于j从1乘2中到n dos：=s并集{r2d2（i/j）}odod；sort（convert（s，list））end:r2d2:=proc（q）#有理数q 2^的倍频程减缩（-floor（evalf（ln（q）/ln（2））））*q end:plim:=proc；r:=1；i：=0；当不是（r=q）时，做i:=i+1；p:=ithprime（i）；r:=r*p^ordp（q，p）od；i end:vai:=proc（n，i）#第i素数通过专利val对n轮（evalf（n*ln（ithprime（i））/ln（2））的映射end:via:=proch（n，l）#长度l局部i，v的n的专利val；对于i从1到l做v[i]：=vai（n，i）od；convert（convert（v，array），list）end:h:=proc（n，q）#如果q=1，则通过专利val n映射间隔q，然后返回（0）fi；dotprod（vec（q），via（n，plim（q）））end:consis:=proc（n，s）#edo n相对于调和集s局部i的一致性；对于i从1到nops，如果不是h（n，s[i]）=round（n*l2（s[i]]），则返回（false）fiod；RETURN（true）end:consl:=proc（n）#edo n local c的最高odd-limit一致性；c:=3；而consis（n，diamond（c））做c:=c+2od；c-2端：
