（PARI）
isok（n，k）=bigomega（k）==n&&bigomeka（k+1）==n+1；
a（n）=对于（k=1，oo，if（isok（n，k），return（k）））\\丹尼尔·苏图2022年5月5日
（PARI）
生成（A，B，n，k）=A=最大值（A，2^n）；（f（m，p，n）=my（list=list（））；如果（n==1，对于素数（q=max（p，ceil（A/m）），B\m，如果（bigomega（m*q-1）==k，listput（list，m*q-1-）），对于素音（q=p，sqrtnint（B\m、n），list=concat（list、f（m*q，q，n-1））））；列表）；向量（Vec（f（1，2，n）））；
a（n）=我的（x=2^n，y=2*x）；而（1，my（v=生成（x，y，n+1，n））；如果（#v>=1，返回（v[1]））；x=y+1；y=2*x）\\丹尼尔·苏图2023年7月9日
