您的位置:首页 > 其它

Codeforces 19E 树上差分

2018-07-08 22:55 330 查看

思路:

先随便建出来一棵搜索树(图可能不连通?)

每一条非树边(剩下的边)和树边都可以构成一个环。 

我们只看一个非树边和某些树边构成的这些环。

分成三种情况:

1.没有奇环  所有边都可以删

2.有一个奇环  奇环上的边可以删

3.有一堆奇环,一堆偶环

边化到点上

如果是奇环 非树边为(x,y) 在树上x,y这两个地方+1 lca(x,y)-2

偶环相反

做一遍树上递推上去

最后如果答案和奇环边的数量相等  这个边可以删。

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=20050;
int n,m,xx,yy,v
,next
,first
,tot,num,T,stk
,a
,odd,R
,remm,lastT,
rev
,fa
,vis
,size
,son
,dfn
,cnt,top
,deep
,ans,ra
;
struct Node{int x,y,wei;}rec
;
void add(int x,int y){v[tot]=y,next[tot]=first[x],first[x]=tot++;}
void dfs1(int x){
size[x]=1,vis[x]=1,rev[++T]=x,stk[T]=x;
for(int i=first[x];~i;i=next[i])if(v[i]!=fa[x]){
if(vis[v[i]]){rec[++num].x=x,rec[num].y=v[i],rec[num].wei=i/2+1;continue;}
R[v[i]]=i/2+1,deep[v[i]]=deep[x]+1,fa[v[i]]=x;
dfs1(v[i]),size[x]+=size[v[i]];
if(size[v[i]]>size[son[x]])son[x]=v[i];
}
}
void dfs2(int x,int tp){
vis[x]=1,dfn[x]=++cnt;top[x]=tp;
if(son[x])dfs2(son[x],tp);
for(int i=first[x];~i;i=next[i])if(v[i]!=fa[x]&&v[i]!=son[x]&&!vis[v[i]])dfs2(v[i],v[i]);
}
int lca(int x,int y){
int fx=top[x],fy=top[y];
while(fx!=fy){
if(deep[fx]<deep[fy])swap(fx,fy),swap(x,y);
x=fa[fx],fx=top[x];
}if(deep[x]<deep[y])return x;return y;
}
int main(){
memset(first,-1,sizeof(first));
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)scanf("%d%d",&xx,&yy),add(xx,yy),add(yy,xx);
for(int i=1;i<=n;i++)if(!vis[i]){
dfs1(i);for(int j=lastT;j<=T;j++)vis[rev[j]]=0;dfs2(i,i);lastT=T;
}for(int i=1;i<=num;i++){
if((deep[rec[i].x]%2)==(deep[rec[i].y]%2))
a[rec[i].x]++,a[rec[i].y]++,a[lca(rec[i].x,rec[i].y)]-=2,odd++,remm=rec[i].wei;
else a[rec[i].x]--,a[rec[i].y]--,a[lca(rec[i].x,rec[i].y)]+=2;
}if(!odd){printf("%d\n",m);for(int i=1;i<=m;i++)printf("%d ",i);return 0;}
for(int i=n;i;i--)a[fa[stk[i]]]+=a[stk[i]];
for(int i=1;i<=n;i++)if(a[i]==odd)ra[++ans]=R[i];
if(odd==2)ra[++ans]=remm;
printf("%d\n",ans),sort(ra+1,ra+1+ans);
for(int i=1;i<=ans;i++)printf("%d ",ra[i]);
}

 

内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: