結果

問題 No.1659 Product of Divisors
コンテスト
ユーザー ぷら
提出日時 2021-08-27 21:50:42
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 183 ms / 2,000 ms
コード長 1,139 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,136 ms
コンパイル使用メモリ 198,892 KB
最終ジャッジ日時 2025-01-24 02:53:04
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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

long long N,K;
constexpr int mod = 1000000007;
long long f1[2000001],f2[2000001],sum[1000001];

long long modpow(long long a,long long b) {
    long long ans = 1;
    while(b) {
        if(b & 1) {
            (ans *= a) %= mod;
        }
        (a *= a) %= mod;
        b /= 2;
    }
    return ans;
}

long long com(long long n,long long k) {
    return f2[n-K+1]*modpow(f1[k],mod-2)%mod;
}

int main() {
    cin >> N >> K;
    map<long long,int>tmp;
    for(long long i = 2; i*i <= N; i++) {
        while (N%i == 0) {
            N /= i;
            tmp[i]++;
        }
    }
    if(N > 1) {
        tmp[N]++;
    }
    f1[1] = 1;
    for(long long i = 2; i <= 2000000; i++) {
        f1[i] = f1[i-1]*i%mod;
    }
    f2[1] = K%mod;
    for(long long i = 2; i <= 2000000; i++) {
        f2[i] = f2[i-1]*((K+i-1)%mod)%mod;
    }
    sum[0] = 1;
    for(long long i = 1; i <= 1000000; i++) {
        sum[i] = sum[i-1]+com(K+i-1,i);
        sum[i] %= mod;
    }
    long long ans = 1;
    for(auto x:tmp) {
        ans *= sum[x.second];
        ans %= mod;
    }
    cout << ans << endl;
}
0