結果

問題 No.2066 Simple Math !
ユーザー 👑 rin204rin204
提出日時 2022-09-02 23:29:12
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 825 ms / 2,000 ms
コード長 925 bytes
コンパイル時間 547 ms
コンパイル使用メモリ 82,272 KB
実行使用メモリ 78,524 KB
最終ジャッジ日時 2024-11-16 06:31:14
合計ジャッジ時間 21,525 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import gcd

# return (\sum_{i=0}^{n-1} ((a*i+b)//m))
def floor_sum(n, m, a, b):
    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(p, q, k):
    l = 0
    r = 1 << 60

    def cnt(x):
        n = min(q, x // p + 1)
        return floor_sum(n, q, p, q + x - p * (n - 1))

    while r - l > 1:
        mid = (l + r) // 2        
        if cnt(mid) > k:
            r = mid
        else:
            l = mid
    return r


for _ in range(int(input())):
    p, q, k = map(int, input().split())
    g = gcd(p, q)
    p //= g
    q //= g
    print(g * solve(p, q, k))
0