#862400 (C++23) No.391 CODING WAR

• 問題: No.391 CODING WAR
• ユーザー: konishi
• 提出日時: 2023-04-23 18:29:17
• 言語: C++23 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 41 ms / 2,000 ms
• コード長: 1,698 bytes
• コンパイル時間: 3,042 ms
• コンパイル使用メモリ: 244,128 KB
• 実行使用メモリ: 5,248 KB
• 最終ジャッジ日時: 2024-11-08 00:50:40
• 合計ジャッジ時間: 3,924 ms
• ジャッジサーバーID （参考情報）: judge1 / judge4

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	4 ms; 5,248 KB
testcase_01	AC	4 ms; 5,248 KB
testcase_02	AC	4 ms; 5,248 KB
testcase_03	AC	4 ms; 5,248 KB
testcase_04	AC	4 ms; 5,248 KB
testcase_05	AC	2 ms; 5,248 KB
testcase_06	AC	4 ms; 5,248 KB
testcase_07	AC	4 ms; 5,248 KB
testcase_08	AC	4 ms; 5,248 KB
testcase_09	AC	41 ms; 5,248 KB
testcase_10	AC	39 ms; 5,248 KB
testcase_11	AC	2 ms; 5,248 KB
testcase_12	AC	4 ms; 5,248 KB
testcase_13	AC	38 ms; 5,248 KB
testcase_14	AC	31 ms; 5,248 KB
testcase_15	AC	33 ms; 5,248 KB
testcase_16	AC	22 ms; 5,248 KB
testcase_17	AC	26 ms; 5,248 KB
testcase_18	AC	18 ms; 5,248 KB
testcase_19	AC	19 ms; 5,248 KB

ソースコード

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef vector<int> VI;
typedef vector<ll> VL;
typedef vector<vector<int>> VII;
typedef vector<vector<ll>> VLL;
typedef vector<string> VS;

#define pii pair<int,int>
#define pll pair<ll,ll>
#define pis pair<int,string>
#define psi pair<string,int>
#define rep(i,num,n) for(int i=num;i<(int)(n);i++) //for_loop
#define REP(i,n) for(int i=0;i<(int)(n);i++)
#define rrep(i,num,n) for(int i=num-1;i>=(int)(n);i--) //reverse_for>
#define in(x,a,b) (a<=x && x<b)//範囲
#define reo return 0
#define all(v) v.begin(),v.end()
#define INF 2000000000
#define INFL 9000000000000000000
#define sz(c) ((int)(c).size())
#define ten(x) ((int)1e##x)
#define tenll(x) ((ll)1e##x)

const long double pi = 3.1415926535897932384626433832795028841971L;
const ll MOD=1000000007;
const int num_fac=200001;//階乗計算の最大数
ll _fact[200002];//階乗を保管する配列
ll modfact(ll x){//階乗
    if(_fact[x]!=0){
        return _fact[x];
    }
    if(_fact[0]==0){
        _fact[0]=1;
        for(int i=1;i<=num_fac;++i)_fact[i]=_fact[i-1]*i%MOD;
    }
    return _fact[x];
}
ll modpow(ll a,ll n){//累乗
    ll r=1;
    while(n)r=r*((n%2)?a:1)%MOD,a=a*a%MOD,n>>=1;
    return r;
}
ll moddiv(ll a,ll b){//除算
    ll ap_2=modpow(b,MOD-2);
    return (a*ap_2)%MOD;
}

ll aCb(ll a,ll b){
    return moddiv(modfact(a),modfact(a-b)*modfact((b))%MOD);
}

int main(){
    ll n,m;cin>>n>>m;
    if(n<m){
        cout<<0;reo;
    }
    ll ans=0;
    int m2=1;
    for(int i=m;i>=1;i--){
        ll t=(aCb(m,i)*modpow(i,n))%MOD;
        ans+=m2*t;
        ans=(ans+MOD)%MOD;
        m2*=-1;
    }
    cout<<ans<<endl;
}