#https://oeis.org/a338161 #(Python 3.8) S=[1] nomb=5000 terme=1 multipl=10**4 cpt=1 derterm=1 bon=0 fini=0 somm1=0 somm2=0 import copy while len(S)multipl: print(but, " plus de ",terme, " termes") fini=1 break termlst.append(terme) binaires = [bin(x) for x in range(2**(cpt),2**(cpt+1)-1,1) if bin(x).count("1")==int(cpt/2)+1 or bin(x).count("1")==int((cpt+1)/2)+1] bon=0 for ind in range(terme+1,multipl): if bon==1: break if ind in S: continue wlist=copy.deepcopy(termlst) wlist.append(ind) for ind2 in range(len(binaires)): if bon==1: break wbinaire=binaires[ind2] liste1=[wlist[ind3-3] for ind3 in range(3,len(wbinaire),1) if wbinaire[ind3]=="1"] liste2= [x for x in wlist if x not in liste1] somm1 =sum(liste1) somm2=sum(liste2) if somm2-somm1==but: bon=1 liste=liste1+liste2 liste.sort() for i in range(len(liste)): if (somm2-somm1==but and liste[i] in liste1): liste[i]=-liste[i] S=S+liste derterm=min(max(liste1),max(liste2)) break continue continue if bon !=1: print(ind, termlst, "pas trouve") continue print(S)