（Python）
导入系统
导入数学
来自症状理论输入因子
从sympy.theory导入素数范围
定义交集（lst1、lst2）：
lst3=[lst1中的值的值，如果lst2中的值]
返回长度（lst3）
n_primes=1000000
factors=[范围（0，n_primes）中n的素数（n）]
素数=列表（素数范围（0，n素数））
序列=[4]
总和=[序列[0]]
prime_factors=[f代表因子[sequence[0]]]中的f
大n=8
而len（序列）<bign：
new_a=错误
a=序列[-1]+1
while intersection（factors[a]，prime_factors）=0:
a+=1
n=长度（序列）
虽然不是新的a：
new_sum=[a+总和的总和]
prime_sum=假
对于new_sum中的总和：
如果素数和：
prime_sum=真
如果不是质数sum和a不在质数中：
序列追加（a）
打印（a，end=“，”）
sys.stdout.flush（）
总和=总和+新总和+[a]
sums=列表（dict.fromkeys（sums））
主因子=主因子+因子[a]
new_a=正确
其他：
a+=1
而a在素数或交集中（因子[a]，素数因子）=0:
a+=1
打印（）
（Python）
从数学导入gcd
从sympy导入isprime
从itertools导入islice
def-agen（start=4）：术语的生成器
alst，k，sums=[start]，start+1，｛0｝|｛start｝
为True时：
产量也[-1]
while any（gcd（k，an）！=1代表阿尔斯特）或\
任意（k+s不在和中，s的isprime（k+s）在和中）：
k+=1
alst.附加（k）
sums.update（[k+s代表总和中的s）
k+=1
打印（列表（islice（agen（），8））#迈克尔·布拉尼基2022年12月16日
