結果

問題 No.1973 Divisor Sequence
ユーザー 🍮かんプリン🍮かんプリン
提出日時 2022-06-11 03:11:21
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,041 ms / 2,000 ms
コード長 3,981 bytes
コンパイル時間 1,643 ms
コンパイル使用メモリ 175,352 KB
実行使用メモリ 100,844 KB
最終ジャッジ日時 2024-09-21 10:32:29
合計ジャッジ時間 5,409 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 154 ms
29,920 KB
testcase_03 AC 9 ms
5,432 KB
testcase_04 AC 88 ms
18,648 KB
testcase_05 AC 30 ms
16,628 KB
testcase_06 AC 97 ms
25,284 KB
testcase_07 AC 18 ms
6,872 KB
testcase_08 AC 48 ms
14,728 KB
testcase_09 AC 75 ms
16,768 KB
testcase_10 AC 97 ms
20,348 KB
testcase_11 AC 13 ms
5,516 KB
testcase_12 AC 76 ms
20,840 KB
testcase_13 AC 23 ms
8,320 KB
testcase_14 AC 49 ms
11,904 KB
testcase_15 AC 133 ms
21,820 KB
testcase_16 AC 128 ms
25,980 KB
testcase_17 AC 83 ms
17,352 KB
testcase_18 AC 6 ms
5,376 KB
testcase_19 AC 110 ms
21,888 KB
testcase_20 AC 20 ms
6,388 KB
testcase_21 AC 120 ms
30,312 KB
testcase_22 AC 104 ms
21,888 KB
testcase_23 AC 1,041 ms
100,844 KB
testcase_24 AC 272 ms
38,492 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/**
 *   @FileName	a.cpp
 *   @Author	kanpurin
 *   @Created	2022.06.11 03:11:16
**/

#include "bits/stdc++.h" 
using namespace std; 
typedef long long ll;


template< int MOD >
struct mint {
public:
    unsigned int x;
    mint() : x(0) {}
    mint(long long v) {
        long long w = (long long)(v % (long long)(MOD));
        if (w < 0) w += MOD;
        x = (unsigned int)(w);
    }
    mint(std::string &s) {
        unsigned int z = 0;
        for (int i = 0; i < s.size(); i++) {
            z *= 10;
            z += s[i] - '0';
            z %= MOD;
        }
        x = z;
    }
    mint operator+() const { return *this; }
    mint operator-() const { return mint() - *this; }
    mint& operator+=(const mint &a) {
        if ((x += a.x) >= MOD) x -= MOD;
        return *this;
    }
    mint& operator-=(const mint &a) {
        if ((x -= a.x) >= MOD) x += MOD;
        return *this;
    }
    mint& operator*=(const mint &a) {
        unsigned long long z = x;
        z *= a.x;
        x = (unsigned int)(z % MOD);
        return *this;
    }
    mint& operator/=(const mint &a) {return *this = *this * a.inv(); }
    friend mint operator+(const mint& lhs, const mint& rhs) {
        return mint(lhs) += rhs;
    }
    friend mint operator-(const mint& lhs, const mint& rhs) {
        return mint(lhs) -= rhs;
    }
    friend mint operator*(const mint& lhs, const mint& rhs) {
        return mint(lhs) *= rhs;
    }
    friend mint operator/(const mint& lhs, const mint& rhs) {
        return mint(lhs) /= rhs;
    }
    friend bool operator==(const mint& lhs, const mint& rhs) {
        return lhs.x == rhs.x;
    }
    friend bool operator!=(const mint& lhs, const mint& rhs) {
        return lhs.x != rhs.x;
    }
    friend std::ostream& operator<<(std::ostream &os, const mint &n) {
        return os << n.x;
    }
    friend std::istream &operator>>(std::istream &is, mint &n) {
        unsigned int x;
        is >> x;
        n = mint(x);
        return is;
    }
    mint inv() const {
        assert(x);
        return pow(MOD-2);
    }
    mint pow(long long n) const {        
        assert(0 <= n);
        mint p = *this, r = 1;
        while (n) {
            if (n & 1) r *= p;
            p *= p;
            n >>= 1;
        }
        return r;
    }
    
    mint sqrt() const {
        if (this->x < 2) return *this;
        if (this->pow((MOD-1)>>1).x != 1) return mint(0);
        mint b = 1, one = 1;
        while (b.pow((MOD-1) >> 1) == 1) b += one;
        long long m = MOD-1, e = 0;
        while (m % 2 == 0) m >>= 1, e += 1;
        mint x = this->pow((m - 1) >> 1);
        mint y = (*this) * x * x;
        x *= (*this);
        mint z = b.pow(m);
        while (y.x != 1) {
            int j = 0;
            mint t = y;
            while (t != one) j += 1, t *= t;
            z = z.pow(1LL << (e-j-1));
            x *= z; z *= z; y *= z; e = j;
        }
        return x;
    }
};

constexpr int MOD = 1e9 + 7;

template<typename T>
vector<pair<T,int>> prime_factorization(T n) {
    vector<pair<T, int>> res;
    for (T i = 2; i*i <= n; i++) {
        int cnt = 0;
        while (n % i == 0) {
            n /= i;
            cnt++;
        }
        if (cnt > 0) res.push_back({i,cnt});
    }
    if (n > 1) res.push_back({n,1});
    return res;
}

int main() {
    int n;cin >> n;
    ll m;cin >> m;
    auto pf = prime_factorization(m);
    mint<MOD> ans = 1;
    for (int i = 0; i < pf.size(); i++) {
        vector<vector<mint<MOD>>> dp(n+1,vector<mint<MOD>>(pf[i].second+1));
        dp[0][0] = 1;
        for (int j = 0; j < n; j++) {
            for (int k = 0; k <= pf[i].second; k++) {
                for (int s = 0; s+k <= pf[i].second; s++) {
                    dp[j+1][k] += dp[j][s];
                }
            }
        }
        mint<MOD> sum = 0;
        for (int j = 0; j <= pf[i].second; j++) {
            sum += dp[n][j];
        }
        ans *= sum;
    }
    cout << ans << endl;
    return 0;
}
0