結果
問題 | No.1474 かさまJ |
ユーザー | SHIJOU |
提出日時 | 2021-04-09 23:35:33 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 143 ms / 2,500 ms |
コード長 | 3,921 bytes |
コンパイル時間 | 2,211 ms |
コンパイル使用メモリ | 204,164 KB |
実行使用メモリ | 141,656 KB |
最終ジャッジ日時 | 2023-09-07 18:32:57 |
合計ジャッジ時間 | 4,801 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 38 ms
141,612 KB |
testcase_01 | AC | 38 ms
141,456 KB |
testcase_02 | AC | 38 ms
141,532 KB |
testcase_03 | AC | 38 ms
141,464 KB |
testcase_04 | AC | 38 ms
141,464 KB |
testcase_05 | AC | 83 ms
141,620 KB |
testcase_06 | AC | 45 ms
141,468 KB |
testcase_07 | AC | 38 ms
141,636 KB |
testcase_08 | AC | 39 ms
141,568 KB |
testcase_09 | AC | 103 ms
141,656 KB |
testcase_10 | AC | 41 ms
141,612 KB |
testcase_11 | AC | 38 ms
141,468 KB |
testcase_12 | AC | 38 ms
141,644 KB |
testcase_13 | AC | 104 ms
141,452 KB |
testcase_14 | AC | 58 ms
141,652 KB |
testcase_15 | AC | 47 ms
141,564 KB |
testcase_16 | AC | 41 ms
141,580 KB |
testcase_17 | AC | 62 ms
141,572 KB |
testcase_18 | AC | 124 ms
141,524 KB |
testcase_19 | AC | 138 ms
141,560 KB |
testcase_20 | AC | 143 ms
141,568 KB |
ソースコード
//#define _GLIBCXX_DEBUG #include <bits/stdc++.h> using namespace std; #define rep(i, n) for(int i=0; i<n; ++i) #define all(v) v.begin(), v.end() #define rall(v) v.rbegin(), v.rend() using ll = int64_t; using ull = uint64_t; using ld = long double; using P = pair<int, int>; using vs = vector<string>; using vi = vector<int>; using vvi = vector<vi>; template<class T> using PQ = priority_queue<T>; template<class T> using PQG = priority_queue<T, vector<T>, greater<T>>; const int INF = 0xccccccc; const ll LINF = 0xcccccccccccccccLL; template<typename T1, typename T2> inline bool chmax(T1 &a, T2 b) {return a < b && (a = b, true);} template<typename T1, typename T2> inline bool chmin(T1 &a, T2 b) {return a > b && (a = b, true);} template<typename T1, typename T2> istream &operator>>(istream &is, pair<T1, T2> &p) { return is >> p.first >> p.second;} template<typename T1, typename T2> ostream &operator<<(ostream &os, const pair<T1, T2> &p) { return os << p.first << ' ' << p.second;} const unsigned int mod = 1000000007; //const unsigned int mod = 998244353; struct mint { unsigned int x; mint():x(0) {} mint(int64_t x_) { int64_t v = int64_t(x_ % mod); if(v < 0) v += mod; x = (unsigned int)v; } static mint row(int v) { mint v_; v_.x = v; return v_; } mint operator-() const { return row(x?mod-x:0);} 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) { uint64_t z = x; z *= a.x; x = (unsigned int)(z % mod); return *this; } mint operator+(const mint a) const { return mint(*this) += a;} mint operator-(const mint a) const { return mint(*this) -= a;} mint operator*(const mint a) const { return mint(*this) *= a;} friend bool operator==(const mint &a, const mint &b) {return a.x == b.x;} friend bool operator!=(const mint &a, const mint &b) {return a.x != b.x;} mint &operator++() { x++; if(x == mod) x = 0; return *this; } mint &operator--() { if(x == 0) x = mod; x--; return *this; } mint operator++(int) { mint result = *this; ++*this; return result; } mint operator--(int) { mint result = *this; --*this; return result; } mint pow(int64_t t) const { mint x_ = *this, r = 1; while(t) { if(t&1) r *= x_; x_ *= x_; t >>= 1; } return r; } //for prime mod mint inv() const { return pow(mod-2);} mint& operator/=(const mint a) { return *this *= a.inv();} mint operator/(const mint a) {return mint(*this) /= a;} }; istream& operator>>(istream& is, mint& a) { return is >> a.x;} ostream& operator<<(ostream& os, const mint& a) { return os << a.x;} string to_string(mint a) { return to_string(a.x); } struct combination { vector<mint> frac, ifrac; combination(int n):frac(n+1), ifrac(n+1) { assert((unsigned int)(n) < mod); frac[0] = 1; for (int i = 1; i <= n; ++i) frac[i] = frac[i-1]*i; ifrac[n] = frac[n].inv(); for (int i = n; i >= 1; --i) ifrac[i-1] = ifrac[i]*i; } mint operator()(int n, int k) { if (k < 0 || k > n) return 0; return frac[n]*ifrac[k]*ifrac[n-k]; } } c(100000); #define N 42 #define M 20010 //head mint dp[N][N][M]; int main() { ios::sync_with_stdio(false); cin.tie(0); int n, mp, mq, l; cin >> n >> mp >> mq >> l; vi s(n); rep(i, n) cin >> s[i]; dp[0][0][0] = 1; rep(i, n) { rep(j, i+2) { rep(k, mq+1) { if(j and k) { dp[i+1][j][k] += dp[i+1][j][k-1]+dp[i][j-1][k-1]; if(k > s[i]) dp[i+1][j][k] -= dp[i][j-1][k-s[i]-1]; } } rep(k, mq+1) { dp[i+1][j][k] += dp[i][j][k]; } } } mint ans; rep(j, n+1) { rep(k, mq+1) { ans += dp[n][j][k]*c(mp-l*j+k+n-1, n-1); } } cout << ans << endl; }