python3 - <<'PY' p=996979 g=379153 d=14641 B=list(range(509847,509896)) A1=list(range(221445,221491)) A2=list(range(884827,884872)) cliques=[] for a in A1: cliques.append([a]+B) for b in A2: cliques.append(B+[b]) C2=C3=0 bad=0 for X in cliques: orig=[(d*x)%p for x in X] for aa in orig: C2=(C2+aa*aa)%983 C3=(C3+aa*aa*aa)%983 for i in range(50): for j in range(i+1,50): lam=(g*(orig[i]+orig[j])+123801)%p if not(146517<=lam<=146611 or 521496<=lam<=521588 or 855093<=lam<=855186): bad+=1 print('cliques',len(cliques),'bad',bad,'C2,C3',C2,C3) # also verify all lambda registry B and C shift conditions raw bad2=0 for X in cliques: orig=[(d*x)%p for x in X] for i in range(50): for j in range(i+1,50): lam=(g*(orig[i]+orig[j])+123801)%p ok1=146517<=lam<=146611 ok2=793324<=(lam+271828)%p<=793416 ok3=172273<=(lam+314159)%p<=172366 if not(ok1 or ok2 or ok3): bad2+=1 print('bad2',bad2) PY