結果
問題 | No.391 CODING WAR |
ユーザー | Baplisca |
提出日時 | 2020-07-01 23:21:03 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 116 ms / 2,000 ms |
コード長 | 3,568 bytes |
コンパイル時間 | 2,273 ms |
コンパイル使用メモリ | 203,964 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-14 06:24:23 |
合計ジャッジ時間 | 3,786 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,812 KB |
testcase_01 | AC | 2 ms
6,940 KB |
testcase_02 | AC | 2 ms
6,944 KB |
testcase_03 | AC | 2 ms
6,940 KB |
testcase_04 | AC | 2 ms
6,940 KB |
testcase_05 | AC | 2 ms
6,940 KB |
testcase_06 | AC | 2 ms
6,944 KB |
testcase_07 | AC | 2 ms
6,940 KB |
testcase_08 | AC | 2 ms
6,940 KB |
testcase_09 | AC | 116 ms
6,944 KB |
testcase_10 | AC | 92 ms
6,940 KB |
testcase_11 | AC | 66 ms
6,940 KB |
testcase_12 | AC | 2 ms
6,944 KB |
testcase_13 | AC | 105 ms
6,944 KB |
testcase_14 | AC | 94 ms
6,944 KB |
testcase_15 | AC | 107 ms
6,944 KB |
testcase_16 | AC | 65 ms
6,944 KB |
testcase_17 | AC | 72 ms
6,940 KB |
testcase_18 | AC | 52 ms
6,940 KB |
testcase_19 | AC | 52 ms
6,940 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; //#define int long long #define all(v) (v).begin(), (v).end() #define rall(v) (v).rbegin(), (v).rend() #define rep(i,n) for(int i=0;i<n;++i) #define rep1(i,n) for(int i=1;i<n;++i) #define exrep(i, a, b) for(ll i = a; i < b; i++) #define out(x) cout << x << endl #define EPS (1e-7) #define gearup ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0); typedef long long ll; typedef long double ld; typedef unsigned long long ull; typedef vector<int> vi; typedef vector<char> vc; typedef vector<bool> vb; typedef vector<double> vd; typedef vector<string> vs; typedef vector<pair<int,int> > vpii; typedef vector<vector<int> > vvi; typedef vector<vector<char> > vvc; typedef vector<vector<bool> > vvb; typedef vector<vector<double> > vvd; typedef vector<vector<string> > vvs; typedef vector<ll> vl; typedef vector<vector<ll> > vvl; typedef vector<vector<vector<ll> > > vvvl; ll MOD = 1000000007; const long long L_INF = 1LL << 60; const int INF = 2147483647; // 2^31-1 const double PI = acos(-1); //cout<<fixed<<setprecision(10); template<class T> inline bool chmin(T& a, T b) {if (a > b) {a = b;return true;}return false;} template<class T> inline bool chmax(T& a, T b) {if (a < b) {a = b;return true;}return false;} template<class T> void debug(T v){rep(i,v.size()) cout<<v[i]<<" ";cout<<endl;} const ll dx[8] = {1, 1, 0, -1, -1, -1, 0, 1}; const ll dy[8] = {0, 1, 1, 1, 0, -1, -1, -1}; //以降 cin の入力元が 'input.txt' になる //std::ifstream in("input.txt"); //std::cin.rdbuf(in.rdbuf()); // a^b mod を計算する ll mod_pow(ll a, ll b, ll p) { if (b == 0) return 1; if (b & 1) return (a * mod_pow(a, b - 1, p)) % p; ll t = mod_pow(a, b / 2, p); return (t * t) % p; } struct mint { long long x; mint(long long x=0):x((x%MOD+MOD)%MOD){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= MOD) x -= MOD; return *this; } mint& operator-=(const mint a) { if ((x += MOD-a.x) >= MOD) x -= MOD; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= MOD; return *this; } mint operator+(const mint a) const { mint res(*this); return res+=a; } mint operator-(const mint a) const { mint res(*this); return res-=a; } mint operator*(const mint a) const { mint res(*this); return res*=a; } mint pow(long long t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime MOD mint inv() const { return pow(MOD-2); } mint& operator/=(const mint a) { return (*this) *= a.inv(); } mint operator/(const mint a) const { mint res(*this); return res/=a; } }; struct combination { vector<mint> fact, ifact; combination(long long n):fact(n+1),ifact(n+1) { assert(n < MOD); fact[0] = 1; for (int i = 1; i <= n; ++i) fact[i] = fact[i-1]*i; ifact[n] = fact[n].inv(); for (int i = n; i >= 1; --i) ifact[i-1] = ifact[i]*i; } mint operator()(long long n, long long k) { if (k < 0 || k > n) return 0; return fact[n]*ifact[k]*ifact[n-k]; } }; //ホウジョ原理 signed main() { gearup; ll n,m; cin >> n >> m; combination C(m+10); mint res = 0; rep(i,m+1){ mint now = C(m,i) * mod_pow(m-i,n,MOD); if(i & 1)res-=now; else res+=now; } out(res.x); }