結果
問題 | No.876 Range Compress Query |
ユーザー | hamray |
提出日時 | 2020-11-24 22:34:03 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 275 ms / 2,000 ms |
コード長 | 7,854 bytes |
コンパイル時間 | 3,196 ms |
コンパイル使用メモリ | 180,072 KB |
実行使用メモリ | 13,148 KB |
最終ジャッジ日時 | 2024-07-23 18:51:41 |
合計ジャッジ時間 | 5,034 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 3 ms
5,376 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 3 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 3 ms
5,376 KB |
testcase_07 | AC | 3 ms
5,376 KB |
testcase_08 | AC | 3 ms
5,376 KB |
testcase_09 | AC | 3 ms
5,376 KB |
testcase_10 | AC | 3 ms
5,376 KB |
testcase_11 | AC | 274 ms
12,732 KB |
testcase_12 | AC | 224 ms
12,520 KB |
testcase_13 | AC | 219 ms
12,716 KB |
testcase_14 | AC | 266 ms
12,728 KB |
testcase_15 | AC | 183 ms
12,916 KB |
testcase_16 | AC | 267 ms
13,000 KB |
testcase_17 | AC | 253 ms
13,148 KB |
testcase_18 | AC | 275 ms
12,992 KB |
ソースコード
#include <bits/stdc++.h> //#include <atcoder/all> //using namespace atcoder; #pragma GCC target ("avx2") #pragma GCC optimization ("O3") #pragma GCC optimization ("unroll-loops") const double pi = 3.141592653589793238462643383279; using namespace std; //typedef //------------------------------------------ typedef vector<int> VI; typedef vector<VI> VVI; typedef vector<string> VS; typedef pair<int, int> PII; typedef pair<long long, long long> PLL; typedef pair<int, PII> TIII; typedef long long LL; typedef unsigned long long ULL; typedef vector<LL> VLL; typedef vector<VLL> VVLL; //container util //------------------------------------------ #define ALL(a) (a).begin(), (a).endf() #define RALL(a) (a).rbegin(), (a).rend() #define PB push_back #define MP make_pair #define SZ(a) int((a).size()) #define SQ(a) ((a) * (a)) #define EACH(i, c) for (typeof((c).begin()) i = (c).begin(); i != (c).endf(); ++i) #define EXIST(s, e) ((s).find(e) != (s).endf()) #define SORT(c) sort((c).begin(), (c).endf()) //repetition //------------------------------------------ #define FOR(i, s, n) for (int i = s; i < (int)n; ++i) #define REP(i, n) FOR(i, 0, n) #define MOD 1000000007 #define rep(i, a, b) for (int i = a; i < (b); ++i) #define trav(a, x) for (auto &a : x) #define all(x) x.begin(), x.end() typedef long long ll; typedef pair<int, int> pii; typedef vector<int> vi; #define chmin(x, y) x = min(x, y) #define chmax(x, y) x = max(x, y) const double EPS = 1e-9, PI = acos(-1); //ここから編集 typedef string::const_iterator State; template< int mod > struct ModInt { int x; ModInt() : x(0) {} ModInt(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} ModInt &operator+=(const ModInt &p) { if((x += p.x) >= mod) x -= mod; return *this; } ModInt &operator-=(const ModInt &p) { if((x += mod - p.x) >= mod) x -= mod; return *this; } ModInt &operator*=(const ModInt &p) { x = (int) (1LL * x * p.x % mod); return *this; } ModInt &operator/=(const ModInt &p) { *this *= p.inverse(); return *this; } ModInt operator-() const { return ModInt(-x); } ModInt operator+(const ModInt &p) const { return ModInt(*this) += p; } ModInt operator-(const ModInt &p) const { return ModInt(*this) -= p; } ModInt operator*(const ModInt &p) const { return ModInt(*this) *= p; } ModInt operator/(const ModInt &p) const { return ModInt(*this) /= p; } bool operator==(const ModInt &p) const { return x == p.x; } bool operator!=(const ModInt &p) const { return x != p.x; } ModInt inverse() const { int a = x, b = mod, u = 1, v = 0, t; while(b > 0) { t = a / b; swap(a -= t * b, b); swap(u -= t * v, v); } return ModInt(u); } ModInt pow(int64_t n) const { ModInt ret(1), mul(x); while(n > 0) { if(n & 1) ret *= mul; mul *= mul; n >>= 1; } return ret; } friend ostream &operator<<(ostream &os, const ModInt &p) { return os << p.x; } friend istream &operator>>(istream &is, ModInt &a) { int64_t t; is >> t; a = ModInt< mod >(t); return (is); } static int get_mod() { return mod; } }; using modint = ModInt< 1000000007 >; template< typename T > struct Combination { vector< T > _fact, _rfact, _inv; Combination(int sz) : _fact(sz + 1), _rfact(sz + 1), _inv(sz + 1) { _fact[0] = _rfact[sz] = _inv[0] = 1; for(int i = 1; i <= sz; i++) _fact[i] = _fact[i - 1] * i; _rfact[sz] /= _fact[sz]; for(int i = sz - 1; i >= 0; i--) _rfact[i] = _rfact[i + 1] * (i + 1); for(int i = 1; i <= sz; i++) _inv[i] = _rfact[i] * _fact[i - 1]; } inline T fact(int k) const { return _fact[k]; } inline T rfact(int k) const { return _rfact[k]; } inline T inv(int k) const { return _inv[k]; } T P(int n, int r) const { if(r < 0 || n < r) return 0; return fact(n) * rfact(n - r); } T C(int p, int q) const { if(q < 0 || p < q) return 0; return fact(p) * rfact(q) * rfact(p - q); } T H(int n, int r) const { if(n < 0 || r < 0) return (0); return r == 0 ? 1 : C(n + r - 1, r); } }; /* 区間add, 区間min */ template<typename T> struct LazySegmentTree{ const int INF = numeric_limits<int>::max(); int n0; vector<T> node, lazy; LazySegmentTree(vector<T> &v){ int n = v.size(); n0 = 1; while(n0 < n) n0 <<= 1; node.resize(2*n0-1); lazy.assign(2*n0-1, 0); for(int i=0; i<n; i++) node[i+n0-1] = v[i]; for(int i=n0-2; i>=0; i--) node[i] = node[i*2+1] + node[i*2+2]; } void eval(int k){ if(lazy[k] == 0) return; node[k] += lazy[k]; if(k < n0-1){ lazy[k*2+1] += lazy[k]; lazy[k*2+2] += lazy[k]; } lazy[k] = 0; } void update(int a, int b, T x, int k, ll l, ll r){ eval(k); if(a <= l && r <= b){ lazy[k] += x; eval(k); }else if(a < r && l < b) { update(a, b, x, k*2+1, l, (l+r)/2); update(a, b, x, k*2+2, (l+r)/2, r); node[k] = node[k*2+1] + node[k*2+2]; } } void update(int l, int r, T x) { update(l, r, x, 0, 0, n0); } T query(int a, int b, int k, ll l, ll r){ eval(k); if(r <= a || b <= l) return 0; else if(a <= l && r <= b){ return node[k]; }else{ T vl = query(a, b, k*2+1, l, (l+r)/2); T vr = query(a, b, k*2+2, (l+r)/2, r); return vl+vr; } } T query(int l, int r) { return query(l, r, 0, 0, n0); } inline T operator[](int a) { return query(a, a + 1); } }; template<typename T> struct LazySegmentTree2{ const int INF = numeric_limits<int>::max(); int n0; vector<T> node, lazy; LazySegmentTree2(vector<T> &v){ int n = v.size(); n0 = 1; while(n0 < n) n0 <<= 1; node.resize(2*n0-1); lazy.assign(2*n0-1, INF); for(int i=0; i<n; i++) node[i+n0-1] = v[i]; for(int i=n0-2; i>=0; i--) node[i] = node[i*2+1] + node[i*2+2]; } void eval(int k, int p){ if(lazy[k] == INF) return; if(k < n0-1){ lazy[k*2+1] = lazy[k]; lazy[k*2+2] = lazy[k]; } node[k] = lazy[k] * p; lazy[k] = INF; } void update(int a, int b, T x, int k, ll l, ll r){ eval(k, r-l); if(a <= l && r <= b){ lazy[k] = x; eval(k, r-l); }else if(a < r && l < b) { update(a, b, x, k*2+1, l, (l+r)/2); update(a, b, x, k*2+2, (l+r)/2, r); node[k] = node[k*2+1] + node[k*2+2]; } } void update(int l, int r, T x) { update(l, r, x, 0, 0, n0); } T query(int a, int b, int k, ll l, ll r){ eval(k, r-l); if(r <= a || b <= l) return 0; else if(a <= l && r <= b){ return node[k]; }else{ T vl = query(a, b, k*2+1, l, (l+r)/2); T vr = query(a, b, k*2+2, (l+r)/2, r); return vl+vr; } } T query(int l, int r) { return query(l, r, 0, 0, n0); } inline T operator[](int a) { return query(a, a + 1); } }; int main(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(10); int N, Q; cin >> N >> Q; vector<ll> a(N); REP(i,N) cin >> a[i]; vector<ll> b; REP(i,N-1) { if(a[i] != a[i+1]) b.push_back(1); else b.push_back(0); } //cout << b.size() << endl; LazySegmentTree<ll> ori(a); LazySegmentTree2<ll> seg(b); while(Q--){ int t; cin >> t; if(t == 1){ int l, r, x; cin >> l >> r >> x; l--; r--; ori.update(l, r+1, x); if(l != 0){ int t1 = ori[l-1]; int t2 = ori[l]; if(t1 != t2){ seg.update(l-1, l, 1); }else{ seg.update(l-1, l, 0); } } if(r != N-1){ int t1 = ori[r+1]; int t2 = ori[r]; if(t1 != t2){ seg.update(r, r+1, 1); }else{ seg.update(r, r+1, 0); } } }else{ int l, r; cin >> l >> r; l--; r--; cout << seg.query(l, r) + 1 << endl; } } return 0; }