結果

問題 No.2440 Accuracy of Integer Division Approximate Functions
ユーザー 👑 rin204rin204
提出日時 2023-09-18 20:56:32
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 304 ms / 2,000 ms
コード長 846 bytes
コンパイル時間 216 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 79,236 KB
最終ジャッジ日時 2024-07-05 10:14:40
合計ジャッジ時間 6,767 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
52,224 KB
testcase_01 AC 287 ms
78,356 KB
testcase_02 AC 300 ms
78,196 KB
testcase_03 AC 270 ms
77,704 KB
testcase_04 AC 280 ms
78,376 KB
testcase_05 AC 283 ms
78,204 KB
testcase_06 AC 194 ms
77,312 KB
testcase_07 AC 195 ms
77,184 KB
testcase_08 AC 204 ms
77,372 KB
testcase_09 AC 198 ms
77,316 KB
testcase_10 AC 200 ms
77,456 KB
testcase_11 AC 271 ms
78,076 KB
testcase_12 AC 296 ms
79,236 KB
testcase_13 AC 304 ms
78,880 KB
testcase_14 AC 275 ms
77,964 KB
testcase_15 AC 293 ms
79,236 KB
testcase_16 AC 208 ms
77,496 KB
testcase_17 AC 209 ms
77,828 KB
testcase_18 AC 196 ms
78,128 KB
testcase_19 AC 187 ms
78,080 KB
testcase_20 AC 207 ms
77,328 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

def floor_sum(n, m, a, b):
    """
    return \\sum_{i=0}^{n-1} ((a*i+b)//m)
    """
    ret = 0
    while True:
        if a >= m:
            ret += (n - 1) * n // 2 * (a // m)
            a %= m
        if b >= m:
            ret += n * (b // m)
            b %= m
        y_max = (a * n + b) // m
        if y_max == 0:
            return ret
        x_max = y_max * m - b
        ret += (n - (x_max + a - 1) // a) * y_max
        n, m, a, b = y_max, a, m, -x_max % a


def solve(N, D, M, S):
    S = 1 << S
    if S == M * D:
        return N

    if S < M * D:
        th = -S * D // (S - M * D)
    else:
        th = S * D // (S - M * D)
    th = min(th, N)
    l = floor_sum(th, D, 1, 1)
    r = floor_sum(th, S, M, M)
    return th - abs(l - r)


for _ in range(int(input())):
    ans = solve(*map(int, input().split()))
    print(ans)
0