|
（PARI）invecn（v，k，x）=for（i=1，k，if（v[i]==x，返回（i））；0
alist（n）=局部（v=向量（n），x，g）；v[1]=4；对于（k=2，n，x=4；while（invecn（v，k-1，x）|（g=gcd（v[k-1]，x））==1 | | isprime（g），x++）；v[k]=x）；v
（哈斯克尔）
导入数据。列表（删除）
a251756 n=a251756清单！！（n-1）
a251756_list=0:f 0 a002808_list其中
f x zs=g zs其中
g（y:ys）| d==1 | | a010051'd==1=g ys
|否则=y:f y（删除y zs）
式中d=gcd x y
--莱因哈德·祖姆凯勒2014年12月8日
蟒蛇
从gmpy2导入gcd，是
A251756号_列表，l，s，b=[0]，0，1，{}
对于范围内（10**3）：
……i=s
……虽然是真的：
……如果不是b中的i：
………m=gcd（i，l）
………如果不是（m==1或是_prime（m））：
................A251756号_list.append（一）
……l，b[i]=i，真实
……当s在b中时：
……….流行音乐
………s+=1
……休息
……i+=1#柴华武2014年12月8日
| 