結果

問題 No.391 CODING WAR
ユーザー konishikonishi
提出日時 2023-04-23 18:42:44
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 52 ms / 2,000 ms
コード長 1,922 bytes
コンパイル時間 3,255 ms
コンパイル使用メモリ 275,256 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-25 13:58:06
合計ジャッジ時間 4,165 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
6,812 KB
testcase_01 AC 3 ms
6,940 KB
testcase_02 AC 4 ms
6,940 KB
testcase_03 AC 4 ms
6,944 KB
testcase_04 AC 4 ms
6,944 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 3 ms
6,940 KB
testcase_07 AC 4 ms
6,940 KB
testcase_08 AC 3 ms
6,940 KB
testcase_09 AC 52 ms
6,944 KB
testcase_10 AC 52 ms
6,944 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 4 ms
6,944 KB
testcase_13 AC 51 ms
6,940 KB
testcase_14 AC 42 ms
6,940 KB
testcase_15 AC 46 ms
6,944 KB
testcase_16 AC 28 ms
6,940 KB
testcase_17 AC 33 ms
6,940 KB
testcase_18 AC 22 ms
6,944 KB
testcase_19 AC 23 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 extGcd(ll a,ll b,ll &x,ll &y){
    if(b==0){
        x=1;y=0;
        return a;
    }
    ll d=extGcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
ll modinv(ll a){
    ll x,y;
    extGcd(a,MOD,x,y);
    return (x%MOD+MOD)%MOD;
}
ll moddiv(ll a,ll b){//除算
    ll ap_2=modinv (b);
    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;
}
0