結果
問題 |
No.2162 Copy and Paste 2
|
ユーザー |
|
提出日時 | 2022-12-25 20:51:38 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,044 ms / 7,000 ms |
コード長 | 3,609 bytes |
コンパイル時間 | 2,175 ms |
コンパイル使用メモリ | 180,964 KB |
実行使用メモリ | 24,376 KB |
最終ジャッジ日時 | 2024-11-19 11:41:21 |
合計ジャッジ時間 | 14,600 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; template <class T> std::vector<int> z_algorithm(const std::vector<T>& s) { int n = int(s.size()); if (n == 0) return {}; std::vector<int> z(n); z[0] = 0; for (int i = 1, j = 0; i < n; i++) { int& k = z[i]; k = (j + z[j] <= i) ? 0 : std::min(j + z[j] - i, z[i - j]); while (i + k < n && s[k] == s[i + k]) k++; if (j + z[j] < i + z[i]) j = i; } z[0] = n; return z; } std::vector<int> z_algorithm(const std::string& s) { int n = int(s.size()); std::vector<int> s2(n); for (int i = 0; i < n; i++) { s2[i] = s[i]; } return z_algorithm(s2); } //双対セグメント木 template <class S,S (*mapping)(S, S),S (*id)()> struct dual_segtree { public: dual_segtree() : dual_segtree(0) {} dual_segtree(int n) : dual_segtree(std::vector<S>(n, id())) {} dual_segtree(const std::vector<S>& v) : _n(int(v.size())) { log = ceil_pow2(_n); size = 1 << log; d = std::vector<S>(2 * size, id()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; } const S& operator[](int p) const { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); return d[p]; } S& operator[](int p) { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); return d[p]; } void apply(int p, S f) { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); d[p] = mapping(f, d[p]); } void apply(int l, int r, S f) { assert(0 <= l && l <= r && r <= _n); if (l == r) return; l += size; r += size; for (int i = log; i >= 1; i--) { if (((l >> i) << i) != l) push(l >> i); if (((r >> i) << i) != r) push((r - 1) >> i); } while (l < r) { if (l & 1) all_apply(l++, f); if (r & 1) all_apply(--r, f); l >>= 1; r >>= 1; } } private: int _n, size, log; std::vector<S> d; void all_apply(int k, S f) { d[k] = mapping(f, d[k]); } void push(int k) { all_apply(2 * k, d[k]); all_apply(2 * k + 1, d[k]); d[k] = id(); } int ceil_pow2(int n) { int x = 0; while ((1U << x) < (unsigned int)(n)) x++; return x; } }; //xにfを作用させた時の値の変化 int op(int f, int x){return max(x, f);} //作用させても変化させないもの int id(){return 0;} int main(){ ios::sync_with_stdio(false); cin.tie(0); string s; cin >> s; int n = s.size(); auto za = z_algorithm(s); vector<vector<int>> tb(n); set<int> S; for(int i = 0; i < n; i++){ za[i] = min(za[i], i); if(za[i] >= 2){ tb[za[i]].push_back(i); S.insert(i); } } S.insert(n + 1); dual_segtree<int, op, id> seg(n + 1); for(int i = 2; i < n; i++){ int cnt = seg[i] + i - 2; int pos = *S.begin() + i; while(pos <= n){ seg.apply(pos, n + 1, cnt); cnt += i - 1; pos = *S.lower_bound(pos) + i; } for(auto &&p:tb[i]) S.erase(p); } cout << n - seg[n] << '\n'; }