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

• 問題: 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

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#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;
}
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX