結果

問題 No.391 CODING WAR
コンテスト
ユーザー konishi
提出日時 2023-04-23 18:42:44
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 58 ms / 2,000 ms
コード長 1,922 bytes
コンパイル時間 3,039 ms
コンパイル使用メモリ 244,132 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-08 01:01:53
合計ジャッジ時間 4,298 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

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