結果

問題 No.2066 Simple Math !
ユーザー GOTKAKOGOTKAKO
提出日時 2024-04-26 00:29:02
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 385 ms / 2,000 ms
コード長 1,550 bytes
コンパイル時間 2,138 ms
コンパイル使用メモリ 196,144 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-04-26 00:29:17
合計ジャッジ時間 11,922 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 82 ms
6,816 KB
testcase_02 AC 112 ms
6,944 KB
testcase_03 AC 55 ms
6,940 KB
testcase_04 AC 363 ms
6,944 KB
testcase_05 AC 385 ms
6,944 KB
testcase_06 AC 359 ms
6,940 KB
testcase_07 AC 360 ms
6,944 KB
testcase_08 AC 370 ms
6,940 KB
testcase_09 AC 361 ms
6,940 KB
testcase_10 AC 362 ms
6,944 KB
testcase_11 AC 360 ms
6,944 KB
testcase_12 AC 373 ms
6,940 KB
testcase_13 AC 359 ms
6,944 KB
testcase_14 AC 324 ms
6,940 KB
testcase_15 AC 331 ms
6,944 KB
testcase_16 AC 319 ms
6,944 KB
testcase_17 AC 324 ms
6,940 KB
testcase_18 AC 322 ms
6,940 KB
testcase_19 AC 225 ms
6,940 KB
testcase_20 AC 217 ms
6,940 KB
testcase_21 AC 212 ms
6,948 KB
testcase_22 AC 215 ms
6,944 KB
testcase_23 AC 217 ms
6,940 KB
testcase_24 AC 326 ms
6,940 KB
testcase_25 AC 314 ms
6,940 KB
testcase_26 AC 316 ms
6,944 KB
testcase_27 AC 313 ms
6,944 KB
testcase_28 AC 327 ms
6,944 KB
testcase_29 AC 62 ms
6,940 KB
testcase_30 AC 54 ms
6,940 KB
testcase_31 AC 96 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

__int128_t safemod(__int128_t a,__int128_t m){
    a %= m; a += m;
    return a>=m?a-m:a;
}
__int128_t floor(__int128_t a,__int128_t b){
    if(b < 0) a *= -1,b *= -1;
    return a<0?(a+1)/b-1:a/b; 
}
__int128_t ceil(__int128_t a,__int128_t b){
    if(b < 0) a *= -1,b *= -1;
    return a>0?(a-1)/b+1:a/b; 
}
__int128_t floorsum(__int128_t N,__int128_t M,__int128_t a,__int128_t b){
    //return Σ(i=0~N-1)floor((a*i+b)/M).
    __int128_t ret = 0;
    if(a >= M) ret += a/M*(N*(N-1)/2),a %= M;
    if(b >= M) ret += b/M*N,b %= M;
    __int128_t maxy = a*N+b;
    if(maxy >= M) ret += floorsum(maxy/M,a,M,maxy%M);
    return ret;
}
long long floorsum2(__int128_t N,__int128_t M,__int128_t a,__int128_t b){
    if(N <= 0) return 0;
    __int128_t ret = 0;
    if(a < 0 || a >= M) ret += floor(a,M)*(N*(N-1)/2),a = safemod(a,M);
    if(b < 0 || b >= M) ret += floor(b,M)*N,b = safemod(b,M);
    __int128_t maxy = a*N+b;
    if(maxy >= M) ret += floorsum(maxy/M,a,M,maxy%M);
    return ret;
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int T; cin >> T;
    while(T--){
        long long P,Q,K; cin >> P >> Q >> K;
        long long g = gcd(P,Q); P /= g,Q /= g;
        long long low = 0,high = 1e18;
        while(high-low > 1){
            long long mid = (high+low)/2;
            long long now = floorsum2(min(Q,mid/P+1),Q,-P,mid)+min(Q,mid/P+1);
            if(now >= K+1) high = mid;
            else low = mid;
        }
        cout << high*g << "\n";
    }    
}
0