CF #325 div1 C + D + hdu 5458

hdu 5458

给定一个连通图，动态删边，并询问任意两点u,v之间的桥的数量(i.e删掉这条边后，u，v不连通)。保证任意时刻图是联通的。

把原图叫G0，第i个操作结束后的图叫Gi，所有询问结束后的图叫Gq。

把删边改为反向加边，然后询问。那么就变成初始状态为Gq，最终状态为G0的加边过程。

当时想的是先把删完边后的图缩点，构成树。添加边(u,v)时又要缩点，但此时只需要将路径(u,v)边标记一下即可，可以用树链剖分。

缩点比较复杂的地方是要重新构图。要处理好原图和新图的映射关系，比较繁琐。

好一点的做法是，对Gq做个一个生成树T，把Gq - T 中的边放到操作里面去。就变成从T->Gq->G0的过程。

Gq->T只要并查集便可以搞定，相对于tarjan缩点更易实现。

并查集+树链剖分+线段树

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <vector>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include <map>
#include <complex>
#include <limits>

#define PI acos(-1)
#define gcd(A,B) __gcd(A,B)
#define lcm(A,B) (A/gcd(A,B)*B)

#define clr(A,x) memset(A,x,sizeof(A))
#define eps (1e-6)

using namespace std;

typedef long long LL;
typedef pair<int,int> P;
typedef pair<int,P> pp;
typedef complex<double> cpx;
const int mm = 200005;
const int inf = 1e9+7;

int n,m,q;
int op[mm],a[mm],b[mm];
vector<int> G[mm];
multiset<P> ES;
int Sz[mm],hson[mm],idx[mm],top[mm],dp[mm],f[mm];
int id;
int F[mm];
vector<int> ans;
struct segment{
    int col[mm<<2];
    void init(int n)
    {
        int N = 1;
        while(N>>1 < n) N <<= 1;
        fill(col,col+N+1,1);
    }
    void pushdown(int rt)
    {
        if(col[rt] == 0) col[rt<<1] = col[rt<<1|1] = 0;
    }

    void color(int rt,int L,int R,int l,int r)
    {
        if(r < L || R < l) return ;
        if(l <= L && R <= r)
        {
            col[rt] = 0;
            return ;
        }
        pushdown(rt);
        int mid = (L+R)>>1;
        color(rt<<1, L,mid,l,r);
        color(rt<<1|1,mid+1,R,l,r);
        if(col[rt<<1] == col[rt<<1|1]) col[rt] = col[rt<<1];
        else col[rt] = -1;
    }
    int query(int rt,int L,int R,int l,int r)
    {
        if( r < L || R < l) return 0;
        if( l <= L && R <= r && col[rt] != -1)
        {
            if(col[rt] == 0) return 0;
            return R-L+1;
        }
        pushdown(rt);
        int mid = (L+R)>>1;
        return query(rt<<1,L,mid,l,r) + query(rt<<1|1,mid+1,R,l,r);
    }
}seg;

void input()
{
    ES.clear();
    cin >> n >> m >> q;
    for(int i = 0; i < m; i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        if(x < y) swap(x,y);
        ES.insert(P(x,y));
    }
    for(int i = 0; i < q; i++)
    {
        scanf("%d%d%d",op+i,a+i,b+i);
        if(a[i] < b[i]) swap(a[i],b[i]);
        if(op[i] == 1){
            ES.erase(ES.find(P(a[i],b[i])));
        }
    }
}

int ba(int k)
{
    if(k != F[k]) F[k] = ba(F[k]);
    return F[k];
}

void insert(int x,int y)
{
    F[x] = y;
}

void clear()
{
//    cout << ES.size() << endl;
    for(int i = 0; i <= n; i++) F[i] = i;
    for(int i = 0; i <= n ; i++) G[i].clear();
    for(multiset<P>::iterator it = ES.begin(); it != ES.end(); it++)
    {
        int x = it->first, y = it->second;
        if(ba(x) == ba(y)){
            op[q] = 1;a[q] = x;b[q] = y;
            q++;
        }else{
            insert(ba(x),ba(y));
            G[x].push_back(y);
            G[y].push_back(x);
        }
    }
}

void dfs_Sz(int rt)
{
    Sz[rt] = 1;dp[rt] = dp[f[rt]] + 1;
    int &son = hson[rt];
    son = 0;
    for(int i = 0; i < G[rt].size();i++)
    if(G[rt][i] != f[rt]){
        int &u = G[rt][i];
        f[u] = rt;
        dfs_Sz(u);
        Sz[rt] += Sz[u];
        if(son == 0 || Sz[u] > Sz[son]) son = u;
    }
}

void dfs_Chain(int rt,int ancestor)
{
    idx[rt] = ++id;top[rt] = ancestor;
    if(hson[rt]) dfs_Chain(hson[rt],ancestor);
    for(int i = 0; i < G[rt].size(); i++)
    if(G[rt][i] != f[rt] && G[rt][i] != hson[rt]){
        dfs_Chain(G[rt][i],G[rt][i]);
    }
}

void color(int u,int v)
{
    while(top[u] != top[v])
    {
        if(dp[top[u]] < dp[top[v]]) swap(u,v);
        seg.color(1,1,n,idx[top[u]],idx[u]);
        u = f[top[u]];
    }
    if(dp[u] < dp[v]) swap(u,v);
    seg.color(1,1,n,idx[v]+1,idx[u]);
}

int query(int u,int v)
{
    int res = 0;
    while(top[u] != top[v])
    {
        if(dp[top[u]] < dp[top[v]]) swap(u,v);
        res += seg.query(1,1,n,idx[top[u]],idx[u]);
        u = f[top[u]];
    }
    if(dp[u] < dp[v]) swap(u,v);
    return res + seg.query(1,1,n,idx[v]+1,idx[u]);
}

void solve()
{
    seg.init(n);
    int root = 1;
    dp[root] = -1;f[root] = root;
    id = 0;
    dfs_Sz(root);
    dfs_Chain(root,root);

    ans.clear();
    for(int i = q-1; i >= 0; i--)
    {
        if(op[i] == 1){
            color(a[i],b[i]);
        }else{
            ans.push_back(query(a[i],b[i]));
        }
    }
    for(int i = ans.size()-1; i >= 0; i--)
        printf("%d\n",ans[i]);
}


int main()
{
    freopen("wcbao.in","r",stdin);
    int T,cas = 0;
    cin >> T;
    while(T--)
    {
        input();
        clear();
        printf("Case #%d:\n",++cas);
//        cout << q << endl;
        solve();
    }
    return 0;
}

#325 div1 C

x了狗了

打了表发现x，y必须互质才有解。自然想到模拟辗转相除。

当时想的是输出看看，正着反着都试一下。

反序输出后样例不对蒙了忘了反过来试一遍。

#325 div1 D

1<=n <= 25

meet in the middle

sb了