結果

問題 No.2896 Monotonic Prime Factors
ユーザー ねしん
提出日時 2024-09-20 22:21:13
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 619 ms / 2,000 ms
コード長 1,445 bytes
コンパイル時間 1,143 ms
コンパイル使用メモリ 85,332 KB
最終ジャッジ日時 2025-02-24 10:23:40
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <cmath>

using namespace std;

const int MAXN = 17 * 100000 + 1;
const int MOD = 998244353;

vector<long long> fact(MAXN, 1), invfact(MAXN, 1);

// 素因数分解して因数の数を返す関数
int factorize(int A) {
    int c = 0;
    int M = A;
    for (int i = 2; i <= sqrt(A) + 2; i++) {
        if (M < i) break;
        while (M % i == 0) {
            M /= i;
            c++;
        }
    }
    if (M != 1) c++;
    return c;
}

// 組み合わせC(n, r)を計算する関数
long long C(int n, int r) {
    return fact[n] * invfact[r] % MOD * invfact[n - r] % MOD;
}

// 繰り返し二乗法による累乗計算 (mod MOD)
long long mod_pow(long long base, long long exp, long long mod) {
    long long result = 1;
    while (exp > 0) {
        if (exp % 2 == 1) result = result * base % mod;
        base = base * base % mod;
        exp /= 2;
    }
    return result;
}

int main() {
    // 階乗と階乗の逆元を事前に計算
    for (int i = 2; i < MAXN; i++) {
        fact[i] = fact[i - 1] * i % MOD;
        invfact[i] = invfact[i - 1] * mod_pow(i, MOD - 2, MOD) % MOD;
    }

    int Q;
    cin >> Q;
    int cnt;
    cnt=0;
    while (Q--) {
        int a, b;
        cin >> a >> b;
        cnt = cnt+factorize(a);
        if (cnt < b) {
            cout << 0 << endl;
            continue;
        }
        cout << C(cnt - 1, cnt - b) << endl;
    }

    return 0;
}
0