結果
問題 | No.2066 Simple Math ! |
ユーザー |
![]() |
提出日時 | 2024-05-09 12:26:15 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,082 bytes |
コンパイル時間 | 341 ms |
コンパイル使用メモリ | 81,996 KB |
実行使用メモリ | 164,668 KB |
最終ジャッジ日時 | 2024-12-16 01:29:35 |
合計ジャッジ時間 | 69,957 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 11 TLE * 20 |
ソースコード
import sysinput = sys.stdin.readlinefrom math import gcd# https://github.com/shakayami/ACL-for-python/blob/master/math.py# \sum_{i=0}^{n-1} (a * i + b)//mdef floor_sum(n,m,a,b):ans=0if a<0:a2=a%mreturn floor_sum(n,m,a2,b)-n*(n-1)*((a2-a)//m)//2if b<0:b2=b%mreturn floor_sum(n,m,a,b2)-n*((b2-b)//m)if a>=m:ans+=(n-1)*n*(a//m)//2a%=mif b>=m:ans+=n*(b//m)b%=my_max=(a*n+b)//mx_max=(y_max*m-b)if y_max==0:return ansans+=(n-(x_max+a-1)//a)*y_maxans+=floor_sum(y_max,a,m,(a-x_max%a)%a)return ansdef BinarySearch(yes, no, P, Q, K):while abs(yes - no) != 1:mid = (yes + no)//2num = min(P, mid//Q + 1)if floor_sum(num, P, -Q, mid) + num - 1 >= K:yes = midelse:no = midreturn yesT = int(input())for _ in range(T):P, Q, K = map(int, input().split())g = gcd(P, Q)P //= gQ //= gyes = 10 ** 20no = -1print(BinarySearch(yes, no, P, Q, K) * g)