結果
| 問題 |
No.1349 Subset Product Queries
|
| コンテスト | |
| ユーザー |
KoD
|
| 提出日時 | 2021-01-13 19:37:49 |
| 言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 1,725 bytes |
| コンパイル時間 | 5,925 ms |
| コンパイル使用メモリ | 107,644 KB |
| 最終ジャッジ日時 | 2025-01-17 17:33:07 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | RE * 1 |
| other | WA * 18 RE * 12 |
ソースコード
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include <cstdio>
#include <vector>
#include <bitset>
constexpr unsigned B = 10;
int main() {
unsigned N, Q;
std::scanf("%u %u", &N, &Q);
std::vector<unsigned> A(N);
for (auto &x: A) {
std::scanf("%u", &x);
}
const unsigned Size = N / B;
std::vector<std::vector<std::bitset<5001>>> memo(Size, std::vector<std::bitset<5001>>(Size));
for (unsigned i = 0; i < Size; ++i) {
memo[i][i].set(0);
for (unsigned j = i + 1; j < Size; ++j) {
memo[i][j] |= memo[i][j - 1];
for (unsigned k = B * (j - 1); k < B * j; ++k) {
memo[i][j] |= memo[i][j] << A[k];
}
}
}
while (Q--) {
unsigned l, r, k;
std::scanf("%u %u %u", &l, &r, &k);
l -= 1;
if (r - l <= B * 2) {
std::bitset<5001> dp;
dp.set(0);
for (unsigned i = l; i < r; ++i) {
dp |= dp << A[i];
}
if (dp.test(k)) {
std::puts("Yes");
}
else {
std::puts("No");
}
}
else {
r -= 1;
const unsigned lk = (l + B - 1) / B;
const unsigned rk = r / B;
auto dp = memo[lk][rk];
for (unsigned k = l; k < B * lk; ++k) {
dp |= dp << A[k];
}
for (unsigned k = B * rk; k <= r; ++k) {
dp |= dp << A[k];
}
if (dp.test(k)) {
std::puts("Yes");
}
else {
std::puts("No");
}
}
}
}
KoD