結果

問題 No.2326 Factorial to the Power of Factorial to the...
コンテスト
ユーザー Naoto Fujiwara
提出日時 2023-05-28 19:43:32
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 1,469 bytes
コンパイル時間 982 ms
コンパイル使用メモリ 82,844 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-27 16:31:34
合計ジャッジ時間 2,018 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

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)%mod;
        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%mod;
        res%=mod;
    }
    return res;
}

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

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

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

    ll e=factorial(N,MOD);
    
   

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

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


    cout<<ans<<endl;


}
0