羅区[P2147]洞窟の調査

エフェクトタイトル:メンテナンスNポイント無向グラフには、接続性の2点を答えるために、ダイナミックかつ穿刺エッジの側面をサポートしています。

解決策:パーティション+セグメントツリーは失効をチェックしてするように設定することができ
、パーティションツリーラインの基本である、オフラインでお願いします。
動作セグメントツリーのタイムラインに基づいて、ツリーライン分割統治アルゴリズムと呼ばれます。
この問題ツリーラインのメンテナンス期間が発生し、現在の時間における辺の集合です。分析に示すように、一方の側のために、木まで線分が表示され\(O(関数logm)\)ノードの集合、最も\(M \)エッジ、したがって、セグメントツリーは、エッジの集合の合計サイズである\ (O(mlogm)\) サブツリーに再帰的に追加してピントを確認しながらノードを介して、それぞれが、現在の時間間隔で表示されます良いセグメントツリー、ルートDFSからツリー全体を、確立した後、現在のノードを離れる前に、あなたは、現在のノード側に集中する必要があります焦点は、側面とチェックから削除されます。リーフノードまでの時間間隔た場合、対応する時間が質問に答えます。ランクの併用によって相殺互いに素なセットバンド、単一の操作のための時間複雑\(O(logN個)\) そのため、全体的な時間の複雑さがある\(O(mlogmlogn)\)

コードは以下の通りです

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e4 + 10;
const int maxm = 2e5 + 10;
typedef pair<int, int> PI;

int n, m;
map<PI, int> mp;
map<int, PI> qry;
vector<int> ans;

struct node {
    int x, y, rx, ry;
};
stack<node> stk;
int fa[maxn], rk[maxn];
int find(int x) {
    while (x != fa[x]) {
        x = fa[x];
    }
    return x;
}
void merge(int x, int y) {
    x = find(x), y = find(y);
    if (x == y) {
        return;
    }
    if (rk[x] > rk[y]) {
        swap(x, y);
    }
    stk.push(node {x, y, rk[x], rk[y]});
    fa[x] = y;
    if (rk[x] == rk[y]) {
        ++rk[y];
    }
}
void undo() {
    node t = stk.top();
    fa[t.x] = t.x, fa[t.y] = t.y;
    rk[t.x] = t.rx, rk[t.y] = t.ry;
    stk.pop();
}
int ask(int x, int y) {
    return find(x) == find(y) ? 1 : 2;
}

struct SegTree {
#define ls(o) t[o].lc
#define rs(o) t[o].rc
    int lc, rc;
    vector<PI> s;
} t[maxm << 1];
int rt;
void build(int &o, int l, int r) {
    static int tot = 0;
    o = ++tot;
    if (l == r) {
        return;
    }
    int mid = l + r >> 1;
    build(ls(o), l, mid);
    build(rs(o), mid + 1, r);
}
void insert(int o, int l, int r, int x, int y, PI z) {
    if (l == x && r == y) {
        t[o].s.push_back(z);
        return;
    }
    int mid = l + r >> 1;
    if (y <= mid) {
        insert(ls(o), l, mid, x, y, z);
    } else if (x > mid) {
        insert(rs(o), mid + 1, r, x, y, z);
    } else {
        insert(ls(o), l, mid, x, mid, z);
        insert(rs(o), mid + 1, r, mid + 1, y, z);
    }
}
void dfs(int o, int l, int r) {
    for (auto e : t[o].s) {
        merge(e.first, e.second);
    }
    if (l == r) {
        if (qry.find(l) != qry.end()) {
            PI q = qry[l];
            ans.push_back(ask(q.first, q.second));
        }
        for (auto e : t[o].s) {
            undo();
        }
        return;
    }
    int mid = l + r >> 1;
    dfs(ls(o), l, mid);
    dfs(rs(o), mid + 1, r);
    for (auto e : t[o].s) {
        undo();
    }
}

int main() {
    static char s[10]; 
    scanf("%d %d", &n, &m);
    
    for (int i = 1; i <= n; i++) {
        fa[i] = i, rk[i] = 1;
    }
    build(rt, 1, m);
    
    for (int i = 1, x, y; i <= m; i++) {
        scanf("%s %d %d", s, &x, &y);
        if (x > y) {
            swap(x, y);
        }
        if (s[0] == 'C') {
            mp[PI(x, y)] = i;
        } else if (s[0] == 'D') {
            auto p = mp.find(PI(x, y));
            assert(p != mp.end());
            insert(rt, 1, m, p->second, i, p->first);
            mp.erase(p);
        } else {
            qry[i] = PI(x, y);
        }
    }
    for (auto e : mp) {
        insert(rt, 1, m, e.second, m, e.first);
    }
    dfs(rt, 1, m);
    for (int i = 0; i < ans.size(); i++) {
        puts(ans[i] == 1 ? "Yes" : "No");
    }
    return 0;
}

おすすめ

転載: www.cnblogs.com/wzj-xhjbk/p/11469655.html