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了