結果

問題 No.2365 Present of good number
ユーザー 鴨志田卓
提出日時 2023-07-16 22:06:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,137 bytes
コンパイル時間 1,533 ms
コンパイル使用メモリ 168,636 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-17 17:10:19
合計ジャッジ時間 2,482 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge6
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

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;

    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