bzoj3562 瞎搞+并查集+dfs

题目分析

听说的做法是这样的：
1.离线处理询问，那么有些边就是始终没有被处理的，只看这些边，并查集缩点。
2.加上有被处理的边。
3.离线处理询问，dfs暴力搞答案。
嗯嗯，是比直接暴力时间复杂度不知道低到哪里去了，可是——
极限数据还是过不了啊！！！！
还有那个“SHTSC 在变化过程中始终保持着一种特殊的性质:即不存在这样的原子序列 a1,a2,…,an(n>3)满足 a1 与 a2、a2 与a3、……、an-1 与 an 以及 an 与 a1 都通过化学键相连,但它们之间却没有其他化学键相连的情况。”这句话有什么用啊啊啊啊？！！！
注意两个连通块之间可能有多条边。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<climits>
using namespace std;
#define LL long long
const int N=5005,Q=10005,M=200005;
int n,m,T,tot;
int f[N];
int find(int x){
    if(x==f[x])return x;
    f[x]=find(f[x]);return f[x];
}
int s[M],t[M],ma[N][N],vis[N];
int h[N],ne[M],to[M],w[M];
struct node{int u,v,bj;}a[Q];
void dfs(int x){
    vis[x]=1;
    for(int i=h[x];i;i=ne[i])
        if(!vis[to[i]]&&w[i]>0)dfs(to[i]);
}
void add(int x,int y){
    ++tot,to[tot]=y,ne[tot]=h[x],h[x]=tot,w[tot]=1;
    ma[x][y]=tot;
}
void work(){
    int i,j,r1,r2;
    for(i=1;i<=m;++i){
        r1=f[s[i]],r2=f[t[i]];
        if(r1==r2)continue;
        if(!ma[r1][r2])add(r1,r2),add(r2,r1);
        else ++w[ma[r1][r2]],++w[ma[r2][r1]];
    }
    for(i=1;i<=T;++i){
        if(a[i].bj==1){
            for(j=1;j<=n;++j)vis[j]=0;
            for(j=1;j<=n;++j)
                if(f[j]==j&&!vis[j])++a[i].u,dfs(j);
        }
        else if(a[i].bj==2){
            r1=f[a[i].u],r2=f[a[i].v];
            if(r1==r2)continue;
            if(!ma[r1][r2])add(r1,r2),add(r2,r1);
            else ++w[ma[r1][r2]],++w[ma[r2][r1]];
        }
        else {
            r1=f[a[i].u],r2=f[a[i].v];
            --w[ma[r1][r2]],--w[ma[r2][r1]];
        }
    }
}
int main(){
    int i,j,x,y;char ch[10];
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;++i)f[i]=i;
    for(i=1;i<=m;++i){
        scanf("%d%d",&s[i],&t[i]);
        ma[s[i]][t[i]]=ma[t[i]][s[i]]=1;
    }
    scanf("%d",&T);
    for(i=1;i<=T;++i){
        scanf("%s",ch);
        if(ch[0]=='Q')a[i].bj=1;
        else {
            scanf("%d%d",&a[i].u,&a[i].v);
            if(ch[0]=='A')a[i].bj=2;
            else a[i].bj=3,ma[a[i].u][a[i].v]=ma[a[i].v][a[i].u]=0;
        }
    }
    for(i=1;i<=n;++i)
        for(j=i+1;j<=n;++j)
        if(ma[i][j])f[find(i)]=find(j);
    memset(ma,0,sizeof(ma));
    for(i=1;i<=n;++i)find(i);
    work();
    for(i=1;i<=T;++i)
        if(a[i].bj==1)printf("%d\n",a[i].u);
    return 0;
}