結果
問題 | No.2658 L-R Nim |
ユーザー | ecottea |
提出日時 | 2024-03-06 16:20:36 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
CE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 7,990 bytes |
コンパイル時間 | 51 ms |
最終ジャッジ日時 | 2024-09-29 18:16:39 |
合計ジャッジ時間 | 1,304 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
43f36c4c4bd3 [/j_bin/judge_tool judge 0 30000 ../CompileMemory.txt /dev/null sud /dev/null _ g++-12 -I /boost_git/ -O2 -lm -Wuninitialized -DONLINE_JUDGE -o a.out main.cpp] open /home/yuki2006/gopath/src/yukicoder/judge/lang/lang.csv: no such file or directory goroutine 1 [running]: runtime/debug.Stack() /usr/local/go/src/runtime/debug/stack.go:24 +0x5e main.main.func1() /home/yuki2006/gopath/src/yukicoder/judge/main.go:20 +0x57 panic({0x7661e0?, 0xc000070b70?}) /usr/local/go/src/runtime/panic.go:770 +0x132 yukicoder/env.InitLangCommands(0x0) /home/yuki2006/gopath/src/yukicoder/env/lang.go:57 +0x3a5 main.main() /home/yuki2006/gopath/src/yukicoder/judge/main.go:42 +0x65
ソースコード
#ifndef HIDDEN_IN_VS // 折りたたみ用 // 警告の抑制 #define _CRT_SECURE_NO_WARNINGS // ライブラリの読み込み #include <bits/stdc++.h> using namespace std; // 型名の短縮 using ll = long long; using ull = unsigned long long; // -2^63 ~ 2^63 = 9 * 10^18(int は -2^31 ~ 2^31 = 2 * 10^9) using pii = pair<int, int>; using pll = pair<ll, ll>; using pil = pair<int, ll>; using pli = pair<ll, int>; using vi = vector<int>; using vvi = vector<vi>; using vvvi = vector<vvi>; using vvvvi = vector<vvvi>; using vl = vector<ll>; using vvl = vector<vl>; using vvvl = vector<vvl>; using vvvvl = vector<vvvl>; using vb = vector<bool>; using vvb = vector<vb>; using vvvb = vector<vvb>; using vc = vector<char>; using vvc = vector<vc>; using vvvc = vector<vvc>; using vd = vector<double>; using vvd = vector<vd>; using vvvd = vector<vvd>; template <class T> using priority_queue_rev = priority_queue<T, vector<T>, greater<T>>; using Graph = vvi; // 定数の定義 const double PI = acos(-1); const vi DX = { 1, 0, -1, 0 }; // 4 近傍(下,右,上,左) const vi DY = { 0, 1, 0, -1 }; int INF = 1001001001; ll INFL = 4004004003104004004LL; // (int)INFL = 1010931620; // 入出力高速化 struct fast_io { fast_io() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(18); } } fastIOtmp; // 汎用マクロの定義 #define all(a) (a).begin(), (a).end() #define sz(x) ((int)(x).size()) #define lbpos(a, x) (int)distance((a).begin(), std::lower_bound(all(a), x)) #define ubpos(a, x) (int)distance((a).begin(), std::upper_bound(all(a), x)) #define Yes(b) {cout << ((b) ? "Yes\n" : "No\n");} #define rep(i, n) for(int i = 0, i##_len = int(n); i < i##_len; ++i) // 0 から n-1 まで昇順 #define repi(i, s, t) for(int i = int(s), i##_end = int(t); i <= i##_end; ++i) // s から t まで昇順 #define repir(i, s, t) for(int i = int(s), i##_end = int(t); i >= i##_end; --i) // s から t まで降順 #define repe(v, a) for(const auto& v : (a)) // a の全要素(変更不可能) #define repea(v, a) for(auto& v : (a)) // a の全要素(変更可能) #define repb(set, d) for(int set = 0, set##_ub = 1 << int(d); set < set##_ub; ++set) // d ビット全探索(昇順) #define repis(i, set) for(int i = lsb(set), bset##i = set; i >= 0; bset##i -= 1 << i, i = lsb(bset##i)) // set の全要素(昇順) #define repp(a) sort(all(a)); for(bool a##_perm = true; a##_perm; a##_perm = next_permutation(all(a))) // a の順列全て(昇順) #define smod(n, m) ((((n) % (m)) + (m)) % (m)) // 非負mod #define uniq(a) {sort(all(a)); (a).erase(unique(all(a)), (a).end());} // 重複除去 #define EXIT(a) {cout << (a) << endl; exit(0);} // 強制終了 #define inQ(x, y, u, l, d, r) ((u) <= (x) && (l) <= (y) && (x) < (d) && (y) < (r)) // 半開矩形内判定 // 汎用関数の定義 template <class T> inline ll powi(T n, int k) { ll v = 1; rep(i, k) v *= n; return v; } template <class T> inline bool chmax(T& M, const T& x) { if (M < x) { M = x; return true; } return false; } // 最大値を更新(更新されたら true を返す) template <class T> inline bool chmin(T& m, const T& x) { if (m > x) { m = x; return true; } return false; } // 最小値を更新(更新されたら true を返す) template <class T> inline T get(T set, int i) { return (set >> i) & T(1); } // 演算子オーバーロード template <class T, class U> inline istream& operator>>(istream& is, pair<T, U>& p) { is >> p.first >> p.second; return is; } template <class T> inline istream& operator>>(istream& is, vector<T>& v) { repea(x, v) is >> x; return is; } template <class T> inline vector<T>& operator--(vector<T>& v) { repea(x, v) --x; return v; } template <class T> inline vector<T>& operator++(vector<T>& v) { repea(x, v) ++x; return v; } #endif // 折りたたみ用 #if __has_include(<atcoder/all>) #include <atcoder/all> using namespace atcoder; #ifdef _MSC_VER #include "localACL.hpp" #endif //using mint = modint1000000007; using mint = modint998244353; //using mint = modint; // mint::set_mod(m); namespace atcoder { inline istream& operator>>(istream& is, mint& x) { ll x_; is >> x_; x = x_; return is; } inline ostream& operator<<(ostream& os, const mint& x) { os << x.val(); return os; } } using vm = vector<mint>; using vvm = vector<vm>; using vvvm = vector<vvm>; using vvvvm = vector<vvvm>; #endif #ifdef _MSC_VER // 手元環境(Visual Studio) #include "local.hpp" #else // 提出用(gcc) inline int popcount(int n) { return __builtin_popcount(n); } inline int popcount(ll n) { return __builtin_popcountll(n); } inline int lsb(int n) { return n != 0 ? __builtin_ctz(n) : -1; } inline int lsb(ll n) { return n != 0 ? __builtin_ctzll(n) : -1; } inline int msb(int n) { return n != 0 ? (31 - __builtin_clz(n)) : -1; } inline int msb(ll n) { return n != 0 ? (63 - __builtin_clzll(n)) : -1; } #define dump(...) #define dumpel(v) #define dump_list(v) #define dump_mat(v) #define input_from_file(f) #define output_to_file(f) #define Assert(b) { if (!(b)) while (1) cout << "OLE"; } #endif bool naive(int n, int k, vi a) { rep(i, n) { repi(na, max(1, a[i] - k), a[i] + k) { vi a2(a); a2[i] = na; vi acc(n + 1); rep(i, n) acc[i + 1] = acc[i] ^ a2[i]; uniq(acc); if (sz(acc) == n + 1) return true; } } return false; } /* * 組の XOR の辞書が Binary Trie の魔改造で作れるのか??? → 無理でした. * もしやオンラインアダマール変換とかできる??? → 無理でした. * *【解説 AC】 * 典型「値域が狭い変数を決め打ち」をすればよかったが, * そもそも a XOR (a+x) (x∈[-100..100]) の値域が狭いことに気づいていなかった. * * 典型「XOR は繰り上がり[下がり]のない足し[引き]算」は頭にあったものの, * じゃあその差ってどれくらい?を考えたことがなかったのが敗因.ぐぬぬ. */ bool solve(int n, int k, vi a) { int A = 1 << (msb(*max_element(all(a))) + 1); vi acc(n + 1); rep(i, n) acc[i + 1] = acc[i] ^ a[i]; dump("acc:", acc); vi imos(n + 1); int v = 0; vi l_min(A, -1); repi(r, 0, n) { int l = l_min[acc[r]]; if (l != -1) { imos[l]++; imos[r]--; v++; } l_min[acc[r]] = r; } if (v == 0) return true; rep(i, n) imos[i + 1] += imos[i]; dump("imos:", imos); int L = INF, R = -INF; rep(i, n) { if (imos[i] == v) { chmin(L, i); chmax(R, i); } } dump("L, R:", L, R); if (L == INF) return false; unordered_map<int, vi> fb; repi(i, L, R) { repi(x, max(a[i] - k, 1), a[i] + k) { int y = a[i] ^ x; if (!fb.count(y)) { vi imos(n + 1); vi l_min(A, -1); repi(r, 0, n) { if ((acc[r] ^ y) < A) { int l = l_min[acc[r] ^ y]; if (l != -1) { imos[l]++; imos[r]--; } } if (l_min[acc[r]] == -1) l_min[acc[r]] = r; } rep(i, n) imos[i + 1] += imos[i]; dump(y); dump(imos); fb[y] = move(imos); } if (fb[y][i] == 0) return true; } } return false; } void bug_find() { #ifdef _MSC_VER // 合わない入力例を見つける. mute_dump = true; mt19937_64 mt; mt.seed((int)time(NULL)); uniform_int_distribution<ll> rnd(0LL, 1LL << 60); rep(hoge, 10000) { int n = rnd(mt) % 10 + 1; int k = rnd(mt) % 10 + 1; vi a(n); rep(i, n) a[i] = rnd(mt) % 18 + 1; auto res_naive = naive(n, k, a); auto res_solve = solve(n, k, a); if (res_naive != res_solve) { cout << "----------error!----------" << endl; cout << "input:" << endl; cout << n << " " << k << endl; cout << a << endl; cout << "results:" << endl; cout << res_naive << endl; cout << res_solve << endl; cout << "--------------------------" << endl; } } mute_dump = false; exit(0); #endif } int main() { // input_from_file("input.txt"); // output_to_file("output.txt"); bug_find(); int n, k; cin >> n >> k; vi a(n); rep(i, n) cin >> a[i]; Yes(solve(n, k, a)); }