結果
問題 | No.868 ハイパー部分和問題 |
ユーザー | firiexp |
提出日時 | 2019-08-17 11:17:22 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,202 ms / 7,000 ms |
コード長 | 3,142 bytes |
コンパイル時間 | 1,194 ms |
コンパイル使用メモリ | 102,340 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-14 21:35:28 |
合計ジャッジ時間 | 17,168 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,816 KB |
testcase_02 | AC | 2 ms
6,816 KB |
testcase_03 | AC | 2 ms
6,816 KB |
testcase_04 | AC | 2 ms
6,816 KB |
testcase_05 | AC | 2 ms
6,820 KB |
testcase_06 | AC | 2 ms
6,820 KB |
testcase_07 | AC | 2 ms
6,816 KB |
testcase_08 | AC | 2 ms
6,820 KB |
testcase_09 | AC | 2 ms
6,816 KB |
testcase_10 | AC | 2 ms
6,816 KB |
testcase_11 | AC | 2 ms
6,816 KB |
testcase_12 | AC | 2 ms
6,820 KB |
testcase_13 | AC | 2 ms
6,824 KB |
testcase_14 | AC | 688 ms
6,816 KB |
testcase_15 | AC | 689 ms
6,824 KB |
testcase_16 | AC | 690 ms
6,820 KB |
testcase_17 | AC | 686 ms
6,816 KB |
testcase_18 | AC | 689 ms
6,820 KB |
testcase_19 | AC | 459 ms
6,820 KB |
testcase_20 | AC | 459 ms
6,816 KB |
testcase_21 | AC | 462 ms
6,816 KB |
testcase_22 | AC | 458 ms
6,820 KB |
testcase_23 | AC | 457 ms
6,820 KB |
testcase_24 | AC | 457 ms
6,816 KB |
testcase_25 | AC | 457 ms
6,816 KB |
testcase_26 | AC | 458 ms
6,816 KB |
testcase_27 | AC | 456 ms
6,820 KB |
testcase_28 | AC | 457 ms
6,820 KB |
testcase_29 | AC | 589 ms
6,820 KB |
testcase_30 | AC | 590 ms
6,816 KB |
testcase_31 | AC | 592 ms
6,816 KB |
testcase_32 | AC | 684 ms
6,820 KB |
testcase_33 | AC | 685 ms
6,816 KB |
testcase_34 | AC | 1,199 ms
6,816 KB |
testcase_35 | AC | 1,202 ms
6,820 KB |
testcase_36 | AC | 704 ms
6,820 KB |
testcase_37 | AC | 701 ms
6,816 KB |
testcase_38 | AC | 2 ms
6,816 KB |
testcase_39 | AC | 2 ms
6,820 KB |
ソースコード
#include <iostream> #include <algorithm> #include <iomanip> #include <map> #include <set> #include <queue> #include <stack> #include <numeric> #include <bitset> #include <cmath> #include <limits> static const int MOD = 1012924417; using ll = long long; using u32 = uint32_t; using namespace std; template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208; template <class T, class U> vector<T> make_v(U size, const T& init){ return vector<T>(static_cast<size_t>(size), init); } template<class... Ts, class U> auto make_v(U size, Ts... rest) { return vector<decltype(make_v(rest...))>(static_cast<size_t>(size), make_v(rest...)); } template<class T> void chmin(T &a, const T &b){ a = (a < b ? a : b); } template<class T> void chmax(T &a, const T &b){ a = (a > b ? a : b); } template<ll M = 1000000007> struct modint{ ll val; modint(): val(0){} template<typename T> explicit modint(T t){val = t%M; if(val < 0) val += M;} modint pow(ll k){ modint res(1), x(val); while(k){ if(k&1) res *= x; x *= x; k >>= 1; } return res; } template<typename T> modint& operator=(T a){ val = a%M; if(val < 0) val += M; return *this; } modint inv() {return pow(M-2);} modint& operator+=(modint a){ val += a.val; if(val >= M) val -= M; return *this;} modint& operator-=(modint a){ val += M-a.val; if(val >= M) val -= M; return *this;} modint& operator*=(modint a){ val = 1LL*val*a.val%M; return *this;} modint& operator/=(modint a){ return (*this) *= a.inv();} modint operator+(modint a) const {return modint(val) +=a;} modint operator-(modint a) const {return modint(val) -=a;} modint operator*(modint a) const {return modint(val) *=a;} modint operator/(modint a) const {return modint(val) /=a;} modint operator-(){ return modint(-val);} bool operator==(const modint a) const {return val == a.val;} bool operator!=(const modint a) const {return val != a.val;} bool operator<(const modint a) const {return val < a.val;} }; using mint = modint<MOD>; auto half = mint(2).inv(); int main() { int n, k; cin >> n >> k; vector<int> v(n); for (auto &&i : v) scanf("%d", &i); auto dp = make_v(2, k+1, mint(0)); dp[0][0] = mint(1); for (int i = 1; i <= n; ++i) { int now = i&1, prv = now^1; for (int j = 0; j <= k; ++j) { dp[now][j] = dp[prv][j]; if(j >= v[i-1]) dp[now][j] += dp[prv][j-v[i-1]]; } } auto DP = dp[n%2]; int q; cin >> q; for (int i = 0; i < q; ++i) { int x, u; scanf("%d %d", &x, &u); vector<mint> newdp(DP); x--; if(v[x] == 0){ for (int j = 0; j <= k; ++j) { newdp[j] = DP[j] * half; } }else { for (int j = v[x]; j <= k; ++j) { newdp[j] -= newdp[j-v[x]]; } } for (int j = 0; j <= k; ++j) { DP[j] = newdp[j]; if(j >= u) DP[j] += newdp[j-u]; } v[x] = u; printf("%d\n", DP[k] != mint(0)); } return 0; }