(Python)
从sympy导入isprime,prime
从itertools导入链,组合为C
def功率集:#跳过空集和singleton
返回链.from_iterable(范围(2,len(s)+1)中r的C(s,r))
定义a(n):
alst,next_set=[范围(1,n)中i的素数(i+1)],{素数(n+1)}
while len(next_set):
alst.append(最小值(next_set));next_set=设置()
对于powerset中的s(alst[-n:]):
ss=总和
如果len(next_set):
如果ss>min(next_set):继续
如果ss>alst[-1]:
if isprime(ss):next_set.add(ss)
return len(alst)#return alst on a(11)forA338134飞机
对于范围(1,12)中的n: