#413766 (C++14) No.876 Range Compress Query

提出ソース
結果

問題	No.876 Range Compress Query
ユーザー	naoya_t
提出日時	2020-01-03 04:12:20
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	WA
実行時間	-
コード長	12,485 bytes
コンパイル時間	2,280 ms
コンパイル使用メモリ	183,560 KB
実行使用メモリ	14,336 KB
最終ジャッジ日時	2024-05-02 07:04:19
合計ジャッジ時間	4,356 ms
ジャッジサーバーID （参考情報）	judge1 / judge4
このコードへのチャレンジ
（要ログイン）
テストケース

テストケース表示
入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	WA	-
testcase_02	AC	2 ms 5,376 KB
testcase_03	WA	-
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	WA	-
testcase_08	AC	2 ms 5,376 KB
testcase_09	WA	-
testcase_10	AC	2 ms 5,376 KB
testcase_11	AC	137 ms 13,696 KB
testcase_12	AC	114 ms 13,184 KB
testcase_13	AC	114 ms 13,440 KB
testcase_14	AC	138 ms 13,696 KB
testcase_15	AC	97 ms 13,824 KB
testcase_16	AC	131 ms 14,336 KB
testcase_17	AC	132 ms 14,336 KB
testcase_18	AC	147 ms 14,336 KB
権限があれば一括ダウンロードができます
ソースコード

raw source code
#include <bits/stdc++.h>
using namespace std;

#define NDEBUG
#include <cassert>


typedef long long ll;
typedef long double Double;
typedef unsigned long long ull;
typedef pair<int,int> ii;
typedef pair<ll,ll> llll;
typedef pair<double,double> dd;

typedef vector<int> vi;
typedef vector<vector<int>> vvi;
typedef vector<ii> vii;
typedef vector<vector<ii>> vvii;
typedef vector<ll> vll;
typedef vector<vector<ll>> vvll;
typedef vector<llll> vllll;
typedef vector<bool> vb;
typedef vector<string> vs;
typedef vector<double> vd;
typedef vector<long double> vD;

#define sz(a)  int((a).size())
#define pb  push_back
#define eb  emplace_back
#define FOR(var,from,to) for(int var=(from);var<=(to);++var)
#define rep(var,n)  for(int var=0;var<(n);++var)
#define rep1(var,n)  for(int var=1;var<=(n);++var)
#define repC2(vari,varj,n)  for(int vari=0;vari<(n)-1;++vari)for(int varj=vari+1;varj<(n);++varj)
#define repC3(vari,varj,vark,n)  for(int vari=0;vari<(n)-2;++vari)for(int varj=vari+1;varj<(n)-1;++varj)for(int vark=varj+1;vark<(n);++vark)
#define ALL(c)  (c).begin(),(c).end()
#define RALL(c)  (c).rbegin(),(c).rend()
#define tr(i,c)  for(auto i=(c).begin(); i!=(c).end(); ++i)
#define found(s,e)  ((s).find(e)!=(s).end())
#define mset(arr,val)  memset(arr,val,sizeof(arr))
#define mid(x,y) ((x)+((y)-(x))/2)
#define IN(x,a,b) ((a)<=(x)&&(x)<=(b))
#define cons make_pair
#define clamp(v,lo,hi) min(max(v,lo),hi)

template<typename T1, typename T2> inline void amin(T1 & a, T2 const & b) { if (a>b) a=b; }
template<typename T1, typename T2> inline void amax(T1 & a, T2 const & b) { if (a<b) a=b; }
template<typename X, typename T> auto vectors(X x, T a) { return vector<T>(x, a); }
template<typename X, typename Y, typename Z, typename... Zs> auto vectors(X x, Y y, Z z, Zs... zs) { auto cont = vectors(y, z, zs...); return vector<decltype(cont)>(x, cont); }

inline ll square(ll x) { return x * x; }
inline ll gcd(ll a, ll b) { while(a) swap(a, b%=a); return b; }
template <typename T>
inline T mod(T a, T b) { return ((a % b) + b) % b; }

template <typename T>
int find_left(vector<T>& v, T elem) {
    return (int)(upper_bound(v.begin(), v.end(), elem) - v.begin()) - 1;
}
template <typename T>
int find_right(vector<T>& v, T elem) {
    return (int)(lower_bound(v.begin(), v.end(), elem) - v.begin());
}

const ll MOD=1000000007LL;

inline ll ADD(ll x, ll y) { return (x+y) % MOD; }
inline ll SUB(ll x, ll y) { return (x-y+MOD) % MOD; }
inline ll MUL(ll x, ll y) { return x*y % MOD; }
inline ll POW(ll x, ll e) { ll v=1; for(; e; x=MUL(x,x), e>>=1) if (e&1) v = MUL(v,x); return v; }
inline ll INV(ll y) { /*assert(y%MOD!=0);*/ return POW(y, MOD-2); }
inline ll DIV(ll x, ll y) { return MUL(x, INV(y)); }

#define INTSPACE 12
char _buf[INTSPACE*1000000 + 3];

int loadint() {
    if (fgets(_buf, INTSPACE+3, stdin)==NULL) return 0;
    return atoi(_buf);
}

int loadvec(vector<int>& v, int N=-1) {
    if (N == 0) {
        v.clear();
        return 0;
    }
    if (N == -1) {
        N = loadint();
        if (N==0) return 0;
    }
    int bufsize = INTSPACE*N + 3;
    if (fgets(_buf, bufsize, stdin)==NULL) return 0;
    v.resize(N);

    int i=0;
    bool last = false;
    for (char *p=&_buf[0]; ;) {
        char *q = p;
        while (*q > ' ') ++q;
        if (*q == 0x0D || *q == 0x0A) last = true;
        *q = 0;
        v[i++] = atoi(p);
        if (last || i == N) break;
        p = q+1;
    }
    return i;
}
void read_cr() {
    fgets(_buf, 256, stdin);
}

void horizontal(vector<int>& v) {
    int L = v.size();
    for (int i=0; i<L; ++i) {
        printf("%d%c", v[i], (i==L-1)?'\n':' ');
    }
}
void horizontall(vector<long long>& v) {
    int L = v.size();
    for (int i=0; i<L; ++i) {
        printf("%lld%c", v[i], (i==L-1)?'\n':' ');
    }
}


template <typename Elem, typename LazyOperand, typename MERGE_DATA_PROC, typename APPLY_LAZY_PROC, typename MERGE_LAZY_PROC>
class LazySegmentTree {
 public:
    MERGE_DATA_PROC f;
    APPLY_LAZY_PROC g;
    MERGE_LAZY_PROC h;
    Elem elem_ident;
    LazyOperand lazy_operand_ident;
    vector<Elem> elems;
    vector<LazyOperand> lazy_operands;
    int n, height;

    inline void assign_merged_lazy(LazyOperand& dest, LazyOperand x) {
        dest = h(dest, x);
    }

    LazySegmentTree(MERGE_DATA_PROC f, APPLY_LAZY_PROC g, MERGE_LAZY_PROC h,
                    Elem elem_ident, LazyOperand lazy_operand_ident)
        : f(f), g(g), h(h), elem_ident(elem_ident), lazy_operand_ident(lazy_operand_ident) {
    }

    void init(int n_temp) {
        n = 1; height = 0;
        while (n < n_temp) { n <<= 1; ++height; }
        elems.assign(2*n, elem_ident);
        lazy_operands.assign(2*n, lazy_operand_ident);
    }

    void build(const vector<Elem>& v){
        int n_temp = v.size();
        init(n_temp);
        rep(i,n_temp) elems[n+i] = v[i];
        for (int i=n-1; i>0; --i) {
            elems[i] = f(elems[i*2], elems[i*2+1]);
        }
    }

    inline Elem reflect(int k){
        return (lazy_operands[k] == lazy_operand_ident) ? elems[k] : g(elems[k], lazy_operands[k]);
    }

    inline void _eval(int k){
        if (lazy_operands[k] != lazy_operand_ident) {
            assign_merged_lazy(lazy_operands[k*2], lazy_operands[k]);
            assign_merged_lazy(lazy_operands[k*2+1], lazy_operands[k]);

            elems[k] = reflect(k);
            lazy_operands[k] = lazy_operand_ident;
        }
    }

    inline void eval_down(int k) {
        for (int i=height; i>0; --i) {
            _eval(k >> i);
        }
    }

    inline void recalc(int k) {
        while (k >>= 1) {
            elems[k] = f(reflect(k*2), reflect(k*2+1));
        }
    }

    void update(int a, int b, LazyOperand x) {
        a += n; b += n;
        eval_down(a);
        eval_down(b-1);
        for (int l=a,r=b; l<r; l>>=1,r>>=1) {
            if (l & 1) assign_merged_lazy(lazy_operands[l++], x);
            if (r & 1) assign_merged_lazy(lazy_operands[--r], x);
        }
        recalc(a);
        recalc(b-1);
    }

    void set_val(int a, Elem e) {
        a += n;
        eval_down(a);
        elems[a]         = e;
        lazy_operands[a] = lazy_operand_ident;
        recalc(a);
    }

    Elem query(int a, int b) {
        a += n; b += n;
        eval_down(a);
        eval_down(b-1);
        Elem vl = elem_ident, vr = elem_ident;
        for (int l=a,r=b; l<r; l>>=1,r>>=1) {
            if (l & 1) vl = f(vl, reflect(l++));
            if (r & 1) vr = f(reflect(--r), vr);
        }
        Elem merged = f(vl, vr);
        return merged;
    }

    void desc() {
    }
};

/***
    using ll2 = pair<ll, int>;
    auto f_sum = [](ll2 a, ll2 b){ return ll2(a.first + b.first, a.second + b.second); };
    auto g = [](ll2 a, ll b){
        ll p = a.second * b;
        return ll2(a.first + p, a.second); };
    auto h = [](ll a, ll b){ return a+b; };
    LazySegmentTree<ll2,ll,decltype(f_sum),decltype(g),decltype(h)> st(f_sum, g, h, ll2(0,0), 0);
    st.build(vector<ll2>(N, ll2(0, 1)));

    auto g = [](ll2 a, ll b){
        return ll2((a.second % 2) ? a.first^b : a.first, a.second); }
    auto h = [](ll a, ll b){ return a^b; };
    LazySegmentTree<ll,ll,decltype(f_sum),decltype(g),decltype(h)> st(f_sum, g, h, ll2(0,0), 0);
    st.build(vector<ll>(N, ll2(0,1));

    auto g = [](ll2 a, ll b){
         ll p = a.second * b;
         return ll2(p ? p-a.first : a.first, a.second); }
    auto h = [](ll a, ll b){ return a^b; };
    LazySegmentTree<ll,ll,decltype(f_sum),decltype(g),decltype(h)> st(f_sum, g, h, ll2(0,0), 0);
    st.build(vector<ll>(N, ll2(0,1));

    auto f = [](ll a, ll b){ return min(a,b); }
    auto g = [](ll a, ll b){ return b; }
    auto h = [](ll a, ll b){ return min(a,b); };
    LazySegmentTree<ll,ll,decltype(f),decltype(g),decltype(h)> st(f, g, h, LLONG_MAX, LLONG_MAX);
    st.build(vector<ll>(N, LLONG_MAX));

    auto f = [](ll a, ll b){ return min(a,b); };
    auto g = [](ll a, ll b){ return min(a,b); };
    auto h = [](ll a, ll b){ return min(a,b); };
    LazySegmentTree<ll,ll,decltype(f),decltype(g),decltype(h)> st(f, g, h, LLONG_MAX, LLONG_MAX);
    st.build(vector<ll>(N, LLONG_MAX));
***/

template <typename T>
class SegmentTree {
public:
    int N;
    vector<T> buf_;
    using MERGER = function<T(T,T)>;
    MERGER merge;
    T ident;

    int ceil2(int size) {
        int n = 1;
        while (n < size) n *= 2;
        return n;
    }

    SegmentTree(int size, MERGER merge, T ident) : merge(merge), ident(ident) {
        N = ceil2(size);
        buf_.assign(N*2, ident);
    }
    SegmentTree(vector<T> ar, MERGER merge, T ident) : merge(merge), ident(ident) {
        int size = ar.size();
        N = ceil2(size);
        buf_.assign(N*2, ident);
        for (int i=0; i<size; ++i) update(i, ar[i]);
    }

    void update(int i, T x) {
        i = N + i - 1;

        if (merge(buf_[i], x) != x) return;

        buf_[i] = x;
        while (i > 0) {
			i = (i - 1) / 2;
			buf_[i] = merge(buf_[i*2+1], buf_[i*2+2]);
        }
    }

    T query(int lo, int hi) {
        return query(lo, hi, 0, 0, N);
    }
    T query(int a, int b, int k, int l, int r) {
        if (r <= a || b <= l) return ident;

        if (a <= l && r <= b) {
            return buf_[k];
        } else {
			T lv = query(a, b, 2*k+1, l, (l+r)/2);
			T rb = query(a, b, 2*k+2, (l+r)/2, r);
			return merge(lv, rb);
        }
    }

    void dump() {
    }
};

template <typename T>
class MinSegmentTree : public SegmentTree<T> {
 public:
    MinSegmentTree(int size) :
        SegmentTree<T>(size, [](T a,T b){return min(a,b);}, numeric_limits<T>::max()) {}
    MinSegmentTree(vector<T> ar) :
        SegmentTree<T>(ar, [](T a,T b){return min(a,b);}, numeric_limits<T>::max()) {}
};

template <typename T>
class MaxSegmentTree : public SegmentTree<T> {
 public:
    MaxSegmentTree(int size) :
        SegmentTree<T>(size, [](T a,T b){return max(a,b);}, numeric_limits<T>::min()) {}
    MaxSegmentTree(vector<T> ar) :
        SegmentTree<T>(ar, [](T a,T b){return max(a,b);}, numeric_limits<T>::min()) {}
};

template <typename T>
class XorSegmentTree : public SegmentTree<T> {
 public:
    XorSegmentTree(int size) :
        SegmentTree<T>(size, [](T a,T b){return a ^ b;}, 0) {}
    XorSegmentTree(vector<T> ar) :
        SegmentTree<T>(ar, [](T a,T b){return a ^ b;}, 0) {}
};

template <typename T>
class GcdSegmentTree : public SegmentTree<T> {
 public:
    GcdSegmentTree(int size) :
        SegmentTree<T>(size, [](T a,T b){return gcd(a,b);}, 0) {}
    GcdSegmentTree(vector<T> ar) :
        SegmentTree<T>(ar, [](T a,T b){return gcd(a,b);}, 0) {}
};

template <typename T>
class SumSegmentTree : public SegmentTree<T> {
 public:
    SumSegmentTree(int size) :
        SegmentTree<T>(size, [](T a,T b){return a+b;}, 0) {}
    SumSegmentTree(vector<T> ar) :
        SegmentTree<T>(ar, [](T a,T b){return a+b;}, 0) {}
};

int main() {
    using ll2 = pair<ll, int>;

    int N, Q; scanf("%d %d%*c", &N, &Q);
    vi a(N); loadvec(a,N);
    vector<ll2> ar(N); rep(i,N) ar[i] = ll2(a[i], 1);

    auto f_sum = [](ll2 a, ll2 b){ return ll2(a.first + b.first, a.second + b.second); };
    auto g = [](ll2 a, ll b){
        ll p = a.second * b;
        return ll2(a.first + p, a.second); };
    auto h = [](ll a, ll b){ return a+b; };
    LazySegmentTree<ll2,ll,decltype(f_sum),decltype(g),decltype(h)> st(f_sum, g, h, ll2(0,0), 0);
    st.build(ar);

    vi dif(N-1,0);
    rep(i,N-1) {
        dif[i] = (int)(a[i] != a[i+1]);
    }
    SumSegmentTree<int> sst(dif);

    int cnt = 0;
    rep(i,Q){
        int op; scanf("%d ", &op);
        vi b(3);
        switch (op){
            case 1:
                {
                    loadvec(b, 3);
                    int l = b[0]-1, r = b[1]-1, x = b[2];
                    st.update(l, r+1, (ll)x);
                    if (0 <= l-1 && l+1 <= N) {
                        ll u = st.query(l-1,l).first, v = st.query(l,l+1).first;
                        sst.update(l-1, (int)(u != v));
                    }
                    if (0 <= r && r+2 <= N) {
                        ll u = st.query(r,r+1).first, v = st.query(r+1,r+2).first;
                        sst.update(r, (int)(u != v));
                    }
                }
                break;
            case 2:
                {
                    loadvec(b, 2);
                    int l = b[0]-1, r = b[1]-1;
                    printf("%d\n", 1 + sst.query(l, r));
                }
                break;
        }

    }
    return 0;
}
yukicoder

結果

テストケース

ソースコード
