結果
問題 | No.391 CODING WAR |
ユーザー | shogier1 |
提出日時 | 2019-10-09 21:37:09 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 145 ms / 2,000 ms |
コード長 | 3,774 bytes |
コンパイル時間 | 972 ms |
コンパイル使用メモリ | 101,044 KB |
実行使用メモリ | 79,192 KB |
最終ジャッジ日時 | 2024-11-17 17:46:40 |
合計ジャッジ時間 | 3,792 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 69 ms
79,068 KB |
testcase_01 | AC | 69 ms
79,064 KB |
testcase_02 | AC | 69 ms
79,064 KB |
testcase_03 | AC | 67 ms
78,936 KB |
testcase_04 | AC | 67 ms
79,064 KB |
testcase_05 | AC | 68 ms
78,940 KB |
testcase_06 | AC | 69 ms
78,936 KB |
testcase_07 | AC | 69 ms
79,192 KB |
testcase_08 | AC | 68 ms
79,064 KB |
testcase_09 | AC | 145 ms
79,068 KB |
testcase_10 | AC | 140 ms
78,936 KB |
testcase_11 | AC | 108 ms
79,064 KB |
testcase_12 | AC | 69 ms
79,064 KB |
testcase_13 | AC | 134 ms
79,064 KB |
testcase_14 | AC | 125 ms
79,064 KB |
testcase_15 | AC | 132 ms
79,064 KB |
testcase_16 | AC | 107 ms
78,940 KB |
testcase_17 | AC | 115 ms
79,192 KB |
testcase_18 | AC | 102 ms
79,064 KB |
testcase_19 | AC | 99 ms
79,064 KB |
ソースコード
#include <algorithm> #include <bitset> #include <complex> #include <deque> #include <iostream> #include <istream> #include <iterator> #include <map> #include <queue> #include <set> #include <stack> #include <string> #include <vector> #include <tuple> #include <cstdint> #include <iomanip> using namespace std; typedef long long ll; #define rep(i, n) for(ll i = 0; i < (n); i++) #define revrep(i, n) for(ll i = (n)-1; i >= 0; i--) #define pb push_back #define f first #define s second const ll INFL = 1LL << 60;//10^18 = 2^60 ll MOD = 1000000007; //ll MOD = 998244353; //vector<mint> dp(N, Mint); //vector<vector<mint>> dp2(N, vector<mint>(N, Mint)); //vector<vector<vector<mint>>> dp3(N, vector<vector<mint>>(N, vector<mint>(N, Mint))); struct mint{ ll x; mint(ll x):x(x % MOD){} mint& operator+=(const mint a){ (x += a.x) %= MOD; return *this; } mint& operator-=(const mint a){ (x += MOD-a.x) %= MOD; return *this; } mint& operator*=(const mint a){ (x *= a.x) %= MOD; return *this; } mint& operator%=(const mint a){ (x %= a.x); return *this; } mint& operator++ (int){ (x += 1) %= MOD; return *this; } mint& operator-- (int){ (x += MOD-1) %= 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 operator%(const mint a) const{ mint res(*this); return res%=a; } mint po(ll t) const{ if(!t) return 1; mint a = po(t>>1); a *= a; if(t&1) a *= *this; return a; } mint inverse() const{ return po(MOD-2); } mint& operator/=(const mint a){ return (*this) *= a.inverse(); } mint operator/(const mint a) const{ mint res(*this); return res/=a; } bool operator == (const mint a){ return this->x == a.x; } bool operator != (const mint a){ return this->x != a.x; } void get(){ cout << x << endl; } };mint Mint = 0; mint pow_mod(ll x, ll k){ mint res = 1; mint a = x; while(k > 0){ if(k % 2){ res *= a; } a *= a; k /= 2; } return res; } //二項演算 const int MAXcomb = 200010; ll fac[MAXcomb], finv[MAXcomb], inv[MAXcomb]; //facはn!,finvは1/n! //invは逆元 void COMinit(){ fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; for(int i = 2; i < MAXcomb; i++){ fac[i] = fac[i-1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD/i) % MOD; finv[i] = finv[i-1] * inv[i] % MOD; } } mint comb(int n, int k){ if(n < k) return 0; if(n < 0 || k < 0) return 0; mint res = fac[n]; res *= finv[k] * finv[n-k]; return res; } //第二種スターリング数 const ll MAXStir2 = 3010; vector<vector<mint>> Stir2memo(MAXStir2, vector<mint>(MAXStir2, Mint)); vector<mint> Bellmemo(MAXStir2, Mint); void Stir2init(){ Stir2memo[0][0] = 1; rep(i, MAXStir2-1)rep(j, i+1)Stir2memo[i+1][j+1] = Stir2memo[i][j] + Stir2memo[i][j+1] * (j+1); rep(i, MAXStir2){ Bellmemo[i] = 0; rep(j,i+1) Bellmemo[i] += Stir2memo[i][j]; } } mint Stir2(ll i, ll j){//区別できるi個をjグループに分ける場合の数 if(i < 0 || j < 0 || i < j) return 0; return Stir2memo[i][j]; } mint Bell(ll x){//区別できるx個をグループ分けする方法全ての場合の数 if(x < 0) return 0; return Bellmemo[x]; } mint kai_mod(ll K){ if(K < 0) return 0; if(K == 0) return 1; return kai_mod(K-1) * K; } int main(){ ll N, M; cin >> N >> M; COMinit(); mint ans = 0; rep(i, M+1){ if(i & 1) ans -= pow_mod(M-i, N) * comb(M, i); else ans += pow_mod(M-i, N) * comb(M, i); } ans.get(); }