結果

問題 No.2326 Factorial to the Power of Factorial to the...
ユーザー Be11T_Be11T_
提出日時 2023-05-29 23:55:12
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 8 ms / 2,000 ms
コード長 2,901 bytes
コンパイル時間 5,013 ms
コンパイル使用メモリ 263,860 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-28 00:12:22
合計ジャッジ時間 5,769 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
4,380 KB
testcase_01 AC 8 ms
4,376 KB
testcase_02 AC 4 ms
4,376 KB
testcase_03 AC 7 ms
4,380 KB
testcase_04 AC 5 ms
4,376 KB
testcase_05 AC 4 ms
4,380 KB
testcase_06 AC 5 ms
4,380 KB
testcase_07 AC 8 ms
4,380 KB
testcase_08 AC 8 ms
4,380 KB
testcase_09 AC 7 ms
4,380 KB
testcase_10 AC 3 ms
4,376 KB
testcase_11 AC 4 ms
4,380 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 4 ms
4,380 KB
testcase_14 AC 7 ms
4,380 KB
testcase_15 AC 3 ms
4,376 KB
testcase_16 AC 2 ms
4,376 KB
testcase_17 AC 7 ms
4,376 KB
testcase_18 AC 4 ms
4,380 KB
testcase_19 AC 3 ms
4,376 KB
testcase_20 AC 2 ms
4,380 KB
testcase_21 AC 2 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#define rep(i, n) for(ll i = 0; i < n; i++)
#define rep2(i, m, n) for(ll i = m; i <= n; i++)
#define rrep(i, m, n) for(ll i = m; i >= n; i--)
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define MAX(x) *max_element(all(x))
#define MIN(x) *min_element(all(x))
#define SZ(a) ((ll)(a).size())
#define bit(n, k) ((n >> k) & 1)
using namespace std;
using ll = long long;
using P = pair<ll, ll>;
const int dx[4] = { 1, 0, -1, 0};
const int dy[4] = { 0, -1, 0, 1};
const int inf = 1e9 + 1;
const ll infll = 1e18;
//const double pi = acos(-1);
     
using Graph = vector<vector<ll>>;
using REV_PQ = priority_queue<ll, vector<ll>, greater<ll>>;
using PQ = priority_queue<ll>;

const int SIZE = 400010;
ll powll(ll n, ll x) {
    // return n ^ x;
    ll ret = 1;
    rep(i, x) ret *= n;
    return ret;
}
void OutputYesNo(bool val) {
    if (val) cout << "Yes" << endl;
    else cout << "No" << endl;
}
void OutputTakahashiAoki(bool val) {
    if (val) cout << "Takahashi" << endl;
    else cout << "Aoki" << endl;
}
void OutPutInteger(ll x) { cout << x << endl; }
void OutPutString(string x) { cout << x << endl; }

int pop_count(ll n){
    int ret = 0;
    while(n > 0){
        if(n % 2 == 1) ret++;
        n /= 2; 
    }
    return ret; 
}

//using PAIR_REV_PQ = priority_queue<pair<ll, ll>, vector<pair<ll, ll>>, greater<pair<ll, ll>>>;
//typedef atcoder::modint998244353 mint;
//typedef atcoder::modint1000000007 mint;
//typedef modint mint;
//typedef vector<mint> ModintVec;
//typedef vector<ModintVec> ModintMat; 

const int MAX_N = 100100;
vector<int> prime; // i.th prime
vector<int> first_prime(MAX_N + 1); // その数を最初に割り切る素数を格納　i = first_prime[i]ならiは素数(i≧2)

void sieve() {
    for (int i = 0; i <= MAX_N; i++) first_prime[i] = 1;
    first_prime[0] = first_prime[1] = 1;
    for (int i = 2; i <= MAX_N; i++) {
        if (first_prime[i] == 1) {
            prime.push_back(i);
            for (int j = i; j <= MAX_N; j += i) first_prime[j] = i;
        }
    }
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    int n, p; cin >> n >> p;
    sieve();

    int cnt = 0;
    rep2(i, 1, n){
        int j = i;
        while(j > 1){
            int x = first_prime[j];
            while(first_prime[j] == x){
                if(x == p) cnt++;
                j /= x;
            }
        }
    }


    using mint1 = dynamic_modint<0>;
    mint1::set_mod(1000000007);
    using mint2 = dynamic_modint<1>;
    mint2::set_mod(1000000006);

    mint1 res = cnt;
    mint2 x_m_1 = 1;
    mint1 x_m = 1;
    rep2(i, 1, n){
        x_m *= i;
        x_m_1 *= i;
    }

    res *= mint1(x_m.val()).pow(x_m_1.val());

    cout << res.val() << endl;


}
0