結果

問題 No.2326 Factorial to the Power of Factorial to the...
ユーザー Naoto FujiwaraNaoto Fujiwara
提出日時 2023-05-28 19:33:00
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,396 bytes
コンパイル時間 584 ms
コンパイル使用メモリ 80,504 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-08-27 15:27:23
合計ジャッジ時間 1,669 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

//https://yukicoder.me/problems/no/2326
#include<iostream>
#include<set>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<map>
#include<numeric>
#include<queue>
#include<tuple>
#include<cmath>
using namespace std;
typedef long long ll;
const ll INF=1LL<<60;
typedef pair<int,int> P;
typedef pair<int,P> PP;
const ll MOD=1e9+7;

ll mod_pow(ll x,ll y, ll mod){
    ll res=1;
    while(y>0){
        if(y&1){
            res*=x;
            res%=mod;
        }
        x*=x;
        x%=mod;
        y/=2;
    }
    return res;
}

ll Legendre(ll n, ll p,ll mod){
    ll res=0;
    ll q=p;
    while (1)
    {
        ll d=n/q;
        if(d==0) break;
        res+=d;
        res%=mod;
        q=q*p;
    }
    return res;
}

//n!のmod
ll factorial(ll n,ll mod){
    ll res=1;
    for(ll t=1;t<=n;t++){
        res*=t;
        res%=mod;
    }
    return res;
}

int main(){
    ll N,P;
    cin>>N>>P;
    //N!の中のPで割れる個数
    ll cnt=Legendre(N,P,MOD);
    /*
    while(n>0){
        cnt+=(n/P)%MOD;
        cnt%=MOD;
        n/=P;
    }
    */

    //N!をMOD-1で割ったあまり
    ll p=factorial(N,MOD-1);


    ll k=p;
    for(ll t=N;t>=1;t--){
        k=mod_pow(k,t,MOD);
    }

    //cout<<"cnt="<<cnt<<endl;
    //cout<<"k="<<k<<endl;

    //これで kは　
    
    ll ans=cnt*k%MOD;
    //ll ans=mod_pow(cnt,k,MOD);


    cout<<ans<<endl;


}
0