結果

問題 No.2365 Present of good number
ユーザー 鴨志田卓鴨志田卓
提出日時 2023-07-17 06:07:55
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 1,164 bytes
コンパイル時間 1,540 ms
コンパイル使用メモリ 168,708 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-18 00:17:15
合計ジャッジ時間 7,517 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

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

using i64 = long long;
const int P = 1e9 + 7;

int qpow(int b, i64 k, int p = P) {
    int ret = 1;
    while(k > 0) {
        if(k & 1) ret = (i64)ret * b % p;
        b = (i64)b * b % p;
        k >>= 1;
    }
    return ret;
}

int main() {
    int n;
    i64 k;
    cin >> n >> k;
	assert(n <= 10 && k > 10);
    int c2 = 0, c3 = 0;
    while(n % 2 == 0) {
        n /= 2;
        c2 ++;
    }
    while(n % 3 == 0) {
        n /= 3;
        c3 ++;
    }

    while(k > 0 && n > 1) {
        swap(c2, c3);
        (c2 *= 2) %= P;
        int nxtn = 1;
        for(int j = 2; j <= n; j ++) {
            int cnt = 0;
            while(n % j == 0) {
                nxtn *= j + 1;
                n /= j;
            }
        }
        n = nxtn;
        while(n % 2 == 0) {
            n /= 2;
            c2 ++;
        }
        while(n % 3 == 0) {
            n /= 3;
            c3 ++;
        }
        k --;
    }

    c3 = c3 * (i64)qpow(2, k - k / 2, P - 1) % (P - 1);
    c2 = c2 * (i64)qpow(2, k / 2, P - 1) % (P - 1);
    if(k & 1) swap(c2, c3);
    cout << (i64)qpow(2, c2) * qpow(3, c3) % P * n % P;
}
0