loj #121. 「離線可過」動態圖連通性

嘟嘟嘟


我用的方法是線段樹+並差集，簡單易懂，常數略大。


咱們按操做時間創建線段樹，而且每個節點開一個vector，記錄這個區間內有哪些邊。而後算出每一條邊的出現時間（這個時間必須是一段連續的區間，若是他加入後刪除又加入，那就算兩條邊），像區間更新打標記同樣把這條邊放進對應節點的vector。


對於詢問，咱們只用在這個時間點對應的節點上標記一下便可。顯然每個葉子結點最多隻能有一個詢問。


接下來，咱們dfs整棵線段樹，走到一個節點，就把這個節點的vector裏面的全部邊加進並差集中。到了葉子結點的時候，若是有詢問，就並差集查一下連通性便可。而後回溯的時候，咱們要撤銷全部的合併操做，即維護一個可撤銷的並差集（撤銷最後一條加進來的邊）。這個也不難，咱們每一次按秩合併的時候用啓發式合併，而且不要路徑壓縮，而後撤銷的時候恢復兩個集合原來的表明元和大小便可。
寫起來很簡單：

int p[maxn], siz[maxn], top = 0;
pr st[maxm << 2];
In int Find(int x) {return x == p[x] ? x : Find(p[x]);}
In int merge(int x, int y)
{
    int px = Find(x), py = Find(y);
    if(px == py) return 0;
    if(siz[px] > siz[py]) swap(px, py);
    p[px] = py, siz[py] += siz[px];
    st[++top] = mp(px, py);
    return 1;
}
In void cancel()
{
    int x = st[top].first, y = st[top].second; --top;
    p[x] = x, siz[y] -= siz[x];
}

分析一下複雜度，對於每個區間，在線段樹上最多隻會被分紅\(log\)個，所以區間總數量\(mlogm\)個，而後並差集複雜度\(logn\)，而遍歷線段樹的複雜度和並差集操做是獨立的，因此不影響，所以總複雜度\(O(mlogmlogn)\)。不過並差集的\(logn\)很小，實際上跑的仍是蠻快的。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cctype>
#include<map>
#include<queue>
#include<vector>
#include<assert.h>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define In inline
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 5e3 + 5;
const int maxm = 5e5 + 5;
In ll read()
{
    ll ans = 0;
    char ch = getchar(), las = ' ';
    while(!isdigit(ch)) las = ch, ch = getchar();
    while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();
    if(las == '-') ans = -ans;
    return ans;
}
In void write(ll x)
{
    if(x < 0) putchar('-'), x = -x;
    if(x >= 10) write(x / 10);
    putchar(x % 10 + '0');
}
In void MYFILE()
{
#ifndef mrclr
    freopen("ha.in", "r", stdin);
    freopen("ha.out", "w", stdout);
#endif
}

#define pr pair<int, int>
#define mp make_pair
int n, m, tim[maxn][maxn];
pr q[maxm];

int l[maxm << 2], r[maxm << 2], pos[maxm];
bool vis[maxm << 2];
vector<pr> v[maxm << 2];
In void build(int L, int R, int now)
{
    l[now] = L, r[now] = R;
    if(L == R) {pos[L] = now; return;}
    int mid = (L + R) >> 1;
    build(L, mid, now << 1), build(mid + 1, R, now << 1 | 1);
}
In void update_E(int L, int R, int now, pr E)
{
    if(l[now] == L && r[now] == R) {v[now].push_back(E); return;}
    int mid = (l[now] + r[now]) >> 1;
    if(R <= mid) update_E(L, R, now << 1, E);
    else if(L > mid) update_E(L, R, now << 1 | 1, E);
    else update_E(L, mid, now << 1, E), update_E(mid + 1, R, now << 1 | 1, E);
}


int p[maxn], siz[maxn], top = 0;
pr st[maxm << 2];
In int Find(int x) {return x == p[x] ? x : Find(p[x]);}
In int merge(int x, int y)
{
    int px = Find(x), py = Find(y);
    if(px == py) return 0;
    if(siz[px] > siz[py]) swap(px, py);
    p[px] = py, siz[py] += siz[px];
    st[++top] = mp(px, py);
    return 1;
}
In void cancel()
{
    int x = st[top].first, y = st[top].second; --top;
    p[x] = x, siz[y] -= siz[x];
}

int ans[maxm];
In void dfs(int now)
{
    int cnt = 0;
    for(auto i : v[now]) cnt += merge(i.first, i.second);
    if(l[now] == r[now])
    {
        if(vis[l[now]]) ans[l[now]] = Find(q[l[now]].first) == Find(q[l[now]].second);
        for(int i = 1; i <= cnt; ++i) cancel();
        return;
    }
    dfs(now << 1), dfs(now << 1 | 1);
    for(int i = 1; i <= cnt; ++i) cancel();
}

int main()
{
//  MYFILE();
    n = read(), m = read();
    build(1, m, 1);
    int FLG = 0;
    for(int i = 1; i <= m; ++i)
    {
        int op = read(), x = read(), y = read();
        if(x > y) swap(x, y);
        if(!op) tim[x][y] = i;
        else if(op == 1) update_E(tim[x][y], i, 1, mp(x, y)), tim[x][y] = 0;
        else q[i] = mp(x, y), FLG = vis[i] = 1;
    }
    if(!FLG) return 0;
    for(int i = 1; i <= n; ++i)
        for(int j = 1; j <= n; ++j) if(tim[i][j]) update_E(tim[i][j], m, 1, mp(i, j));
    for(int i = 1; i <= n; ++i) p[i] = i, siz[i] = 1;
    dfs(1);
    for(int i = 1; i <= m; ++i) if(vis[i]) puts(ans[i] ? "Y" : "N");
    return 0;
}