結果

問題 No.2262 Fractions
ユーザー Carpenters-CatCarpenters-Cat
提出日時 2023-04-07 23:18:08
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,542 bytes
コンパイル時間 2,016 ms
コンパイル使用メモリ 204,764 KB
実行使用メモリ 50,788 KB
最終ジャッジ日時 2024-10-02 20:29:18
合計ジャッジ時間 8,393 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 -- -
testcase_02 -- -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
int main () {
    int T;
    cin >> T;
    vector<vector<int>> X(300030);
    for (int i = 1; i <= 300000; i ++) {
        for (int j = 1; i * j <= 300000; j ++) {
            X[i * j].push_back(i);
        }
    }
    int meb[300030];
    fill(meb, meb + 300030, -2);
    meb[1] = 1;
    for (int i = 2; i <= 300000; i ++) {
        if (meb[i] == -2) {
            meb[i] = -1;
        }
        for (int j = 2; i * j <= 300000; j ++) {
            meb[i * j] *= -1;
            if (j % i == 0) {
                meb[i * j] = 0;
            }
        }
    }
    while (T --) {
        ll N, K;
        cin >> N >> K;
        if (N == 1) {
            puts("1/1");
            continue;
        }
        ll lb = 1, ls = 0, rb = 0, rs = 1;
        bool ok = false;
        while(true) {
            ll b = lb + rb, s = ls + rs;
            if (b > N || s > N) {
                break;
            }
            ll sum = 0;
            for (ll m = 1; m <= N; m ++) {
                ll n = min(N, (m * s) / b);
                for (auto a : X[m]) {
                    sum += (n / a) * meb[a];
                }
            }
            if (sum <= K) {
                lb = b; ls = s;
            } else {
                rb = b; rs = s;
            }
            if (sum == K) {
                ok = true;
                break;
            }
        }
        if (ok) {
            printf("%lld/%lld\n", ls, lb);
        } else {
            puts("-1");
        }
    }
}
0