[BZOJ3514]Codechef MARCH14 GERALD07加强版 题解

[BZOJ3514]Codechef MARCH14 GERALD07加强版 题解

题目地址:BZOJ:Problem 3514. — Codechef MARCH14 GERALD07加强版

题目描述

N个点M条边的无向图,询问保留图中编号在[l,r]的边的时候图中的联通块个数。

输入输出格式

输入格式:
第一行四个整数N、M、K、type,代表点数、边数、询问数以及询问是否加密。
接下来M行,代表图中的每条边。
接下来K行,每行两个整数L、R代表一组询问。对于type=0的测试点,读入的L和R即为询问的L、R;对于type=1的测试点,每组询问的L、R应为L xor lastans和R xor lastans。

输出格式:
K行每行一个整数代表该组询问的联通块个数。

输入输出样例

输入样例#1:

3 5 4 0
1 3
1 2
2 1
3 2
2 2
2 3
1 5
5 5
1 2

输出样例#1:

2
1
3
1

说明

对于100%的数据,1≤N、M、K≤200,000。
2016.2.26提高时限至60s

题解

首先这题强制在线了。考虑原来的GERALD07是怎么做的,用LCT来维护生成树,LCT记录下当前连通块内编号最小的边,离线处理答案,权值线段树来查符合条件的边数,答案即为n-边数。该题链接见:[CC-GERALD07]Chef and Graph Queries 题解 | KSkun’s Blog
现在搞在线了怎么搞,那我们考虑把权值线段树改成主席树,然后查另外一个答案。我们在加边的过程中检查是否成环,成环了就把环中编号最小的边给删掉,然后把这条边的编号记下来。对于编号为[l, r]的边,问题的答案即n-这一段编号的边插入的时候删掉的编号小于l的边的数量。
下面我们证明它的正确性。对于一个加边时删边编号不小于l的边,其实删不删这一块本来就是连通的,并不能对答案产生贡献;反过来,删掉了前面那条边,其实这里并不连通了,此时把这条边加上以后变联通,其实对答案来说有-1的贡献。
至于n-的原因,如果原始有n个点没有边,那就是n个连通块,此时加1条边,就会合并掉某两个连通块,答案变为n-1,以此类推。
这个代码的常数巨大,复杂度没有用,就算BZOJ开了60s我也是蹭着限制AC的,太吓人了。

代码

// Code by KSkun, 2018/3
#include <cstdio>

inline void swap(int &a, int &b) {
    register int t = a;
    a = b;
    b = t;
}

inline char fgc() {
    static char buf[100000], *p1 = buf, *p2 = buf;
    return p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1++;
}

inline int readint() {
    register int res = 0, neg = 1;
    char c = fgc();
    while(c < '0' || c > '9') {
        if(c == '-') neg = -1;
        c = fgc();
    }
    while(c >= '0' && c <= '9') {
        res = res * 10 + c - '0';
        c = fgc();
    }
    return res * neg;
}

const int MAXN = 400005, INF = 1e9;

struct LCTNode {
    int ch[2], fa, val, mn;
    bool rev;
} lct[MAXN];

inline bool isleft(int p) {
    return lct[lct[p].fa].ch[0] == p;
}

inline bool isroot(int p) {
    register int fa = lct[p].fa;
    return lct[fa].ch[0] != p && lct[fa].ch[1] != p;
}

inline void update(int p) {
    register int ls = lct[p].ch[0], rs = lct[p].ch[1];
    lct[p].mn = p;
    if(lct[lct[ls].mn].val < lct[lct[p].mn].val) lct[p].mn = lct[ls].mn;
    if(lct[lct[rs].mn].val < lct[lct[p].mn].val) lct[p].mn = lct[rs].mn;
}

inline void reverse(int p) {
    swap(lct[p].ch[0], lct[p].ch[1]);
    lct[p].rev ^= 1;
}

inline void pushdown(int p) {
    register int ls = lct[p].ch[0], rs = lct[p].ch[1];
    if(lct[p].rev) {
        if(ls) reverse(ls);
        if(rs) reverse(rs);
        lct[p].rev ^= 1;
    }
}

int sta[MAXN], stop;

inline void pushto(int p) {
    stop = 0;
    while(!isroot(p)) {
        sta[stop++] = p;
        p = lct[p].fa;
    }
    pushdown(p);
    while(stop) {
        pushdown(sta[--stop]);
    }
}

inline void rotate(int p) {
    register bool t = !isleft(p); register int fa = lct[p].fa, ffa = lct[fa].fa;
    lct[p].fa = ffa; if(!isroot(fa)) lct[ffa].ch[!isleft(fa)] = p;
    lct[fa].ch[t] = lct[p].ch[!t]; lct[lct[fa].ch[t]].fa = fa;
    lct[p].ch[!t] = fa; lct[fa].fa = p;
    update(fa);
}

inline void splay(int p) {
    pushto(p);
    for(register int fa = lct[p].fa; !isroot(p); rotate(p), fa = lct[p].fa) {
        if(!isroot(fa)) rotate(isleft(fa) == isleft(p) ? fa : p);
    }
    update(p);
}

inline void access(int p) {
    for(register int q = 0; p; q = p, p = lct[p].fa) {
        splay(p);
        lct[p].ch[1] = q;
        update(p);
    }
}

inline void makert(int p) {
    access(p);
    splay(p);
    reverse(p);
}

inline int findrt(int p) {
    access(p);
    splay(p);
    while(lct[p].ch[0]) p = lct[p].ch[0];
    return p;
}

inline void split(int u, int v) {
    makert(u);
    access(v);
    splay(v);
}

inline void link(int u, int v) {
    split(u, v);
    lct[u].fa = v;
}

inline void cut(int u, int v) {
    split(u, v);
    if(lct[v].ch[0] != u || lct[lct[v].ch[0]].ch[1]) return;
    lct[u].fa = lct[v].ch[0] = 0;
}

inline int query(int u, int v) {
    split(u, v);
    return lct[v].mn;
}

struct SGTNode {
    int val, lch, rch;
} sgt[MAXN * 20];
int stot = 0, rt[MAXN];

inline void insert(int &o, int l, int r, int x) {
    sgt[++stot] = sgt[o]; o = stot;
    sgt[o].val++;
    if(l == r) return;
    register int mid = (l + r) >> 1;
    if(x <= mid) insert(sgt[o].lch, l, mid, x);
    else insert(sgt[o].rch, mid + 1, r, x);
}

inline int query(int o1, int o2, int l, int r, int rr) {
    if(r == rr) return sgt[o2].val - sgt[o1].val;
    register int mid = (l + r) >> 1;
    if(rr <= mid) return query(sgt[o1].lch, sgt[o2].lch, l, mid, rr);
    else return sgt[sgt[o2].lch].val - sgt[sgt[o1].lch].val + query(sgt[o1].rch, sgt[o2].rch, mid + 1, r, rr);
}

int n, m, q, type, ut, vt, lst[MAXN], u[MAXN], v[MAXN];

int main() {
    n = readint(); m = readint(); q = readint(); type = readint();
    lct[0].val = INF;
    for(register int i = 1; i <= n; i++) {
        lct[i].mn = i;
        lct[i].val = INF;
    }
    register int tot = n;
    for(register int i = 1; i <= m; i++) {
        ut = u[i] = readint(); vt = v[i] = readint();
        if(ut == vt) {
            lst[i] = i;
            continue;
        }
        if(findrt(ut) == findrt(vt)) {
            int t = query(ut, vt), tv = lct[t].val;
            lst[i] = tv;
            cut(u[tv], t);
            cut(v[tv], t);
        }
        tot++;
        lct[tot].mn = tot;
        lct[tot].val = i;
        link(ut, tot);
        link(vt, tot);
    }
    for(register int i = 1; i <= m; i++) {
        rt[i] = rt[i - 1];
        insert(rt[i], 0, m, lst[i]);
    }
    register int lastans = 0;
    while(q--) {
        ut = readint(); vt = readint();
        if(type) {
            ut ^= lastans;
            vt ^= lastans;
        }
        printf("%d\n", lastans = n - query(rt[ut - 1], rt[vt], 0, m, ut - 1));
    }
    return 0;
}


发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注

This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据