codeforces 19E Fairy 树链剖分

题意:

给定一张无向图，求有多少条边保证删掉它之后能将该图转化成二分图。

---

解析:

首先二分图没有奇环。
我们可以dfs出一棵树。
这时候，非树边都是返祖边。
每一条非树边和树边都可以构成一个环。
我们只需要考虑所有只含一条返祖边的环即可。
如果有且只有一个奇环的话，那么显然删除环上所有的边都可行。
如果有多个奇环的话，显然是选取这些奇环共同经过的树边。
而同时我们需要注意，这些树边是否被任意一个偶环经过。
如果被任意一个偶环，且被奇环经过的话，那么该边是不可取的，因为如果取了它，那么经过它的奇环以及偶环中剩下的边可以组合成一个奇环。
所以我们可以用树剖来维护每个点的父边是否被偶环经过以及被多少奇环经过即可。
复杂度O(nlog^2n)。

---

当然可以用树形DP解决这个问题，请戳dashgua.coding.io

---

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 11000
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
using namespace std;
int n,m,cnt,tot;
int head[N],col[N],vis[N];
int deep[N],son[N],siz[N],top[N],fa[N],dfn[N];
struct node
{
    int from,to,next;
}edge[N<<1];
struct Line
{
    int u,v;
}L[N];
void init()
{
    memset(head,-1,sizeof(head));
    cnt=1;
}
void edgeadd(int from,int to)
{
    edge[cnt].from=from,edge[cnt].to=to,edge[cnt].next=head[from];
    head[from]=cnt++;
}
void dfs1(int now,int ff)
{
    fa[now]=ff,siz[now]=1,son[now]=-1;
    deep[now]=deep[ff]+1;
    for(int i=head[now];i!=-1;i=edge[i].next)
    {
        int to=edge[i].to;
        if(deep[to])continue;
        col[to]=(i/2)+(i%2),vis[(i/2)+(i%2)]=1;
        dfs1(to,now);
        siz[now]+=siz[to];
        if(son[now]==-1||siz[to]>siz[son[now]])son[now]=to;
    }
}
void dfs2(int now,int tp)
{
    top[now]=tp;
    dfn[now]=++tot;
    if(son[now]!=-1)dfs2(son[now],tp);
    for(int i=head[now];i!=-1;i=edge[i].next)
    {
        int to=edge[i].to;
        if(to==son[now]||to==fa[now]||fa[to]!=now)continue;
        dfs2(to,to);
    }
}
int sum[N<<2],cov[N<<2],flag[N<<2];
void pushup(int rt)
{
    sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void pushdown(int rt,int l,int r)
{
    if(cov[rt])
    {
        int mid=(l+r)>>1;
        cov[rt<<1]+=cov[rt];
        cov[rt<<1|1]+=cov[rt];
        sum[rt<<1]+=(mid-l+1)*cov[rt];
        sum[rt<<1|1]+=(r-mid)*cov[rt];
        cov[rt]=0;
    }
}
void update(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        cov[rt]++;
        sum[rt]+=r-l+1;
        return;
    }
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    if(L<=mid)update(L,R,lson);
    if(R>mid)update(L,R,rson);
    pushup(rt);
}
void update2(int L,int R,int l,int r,int rt)
{
    if(flag[rt])return;
    if(L<=l&&r<=R)
    {
        flag[rt]=1;
        return;
    }
    int mid=(l+r)>>1;
    if(L<=mid)update2(L,R,lson);
    if(R>mid)update2(L,R,rson);
    pushup(rt);
}
void Update(int x,int y,int tag)
{
    if(deep[x]<deep[y])swap(x,y);
    while(top[x]!=top[y])
    {
        if(deep[top[x]]<deep[top[y]])swap(x,y);
        int l=dfn[top[x]],r=dfn[x];
        if(tag)update(l,r,1,n,1);
        else update2(l,r,1,n,1);
        x=fa[top[x]];
    }
    if(x==y)return;
    if(deep[x]>deep[y])swap(x,y);
    int l=dfn[x]+1,r=dfn[y];
    if(tag)update(l,r,1,n,1);
    else update2(l,r,1,n,1);
}
int query(int p,int l,int r,int rt)
{
    if(l==r)return sum[rt];
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    if(p<=mid)return query(p,lson);
    else return query(p,rson);
}
bool query2(int p,int l,int r,int rt)
{
    if(flag[rt])return 1;
    if(l==r)return 0;
    int mid=(l+r)>>1;
    if(p<=mid)return query2(p,lson);
    else return query2(p,rson);
}
int sta[N];
int Top;
int main()
{
    init();
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        edgeadd(u,v),edgeadd(v,u);
        L[i].u=u,L[i].v=v;
    }
    for(int i=1;i<=n;i++)
    {
        if(!deep[i])
            dfs1(i,0),dfs2(i,i);
    }
    int cntodd=0,cnteven=0,no;
    for(int i=1;i<=m;i++)
    {
        if(!vis[i])
        {
            int u=L[i].u,v=L[i].v;
            if((deep[u]&1)^(deep[v]&1))
                Update(u,v,0),cnteven++;
            else Update(u,v,1),cntodd++,no=i;
        }
    }
    if(cntodd==0)
    {
        for(int i=1;i<=m;i++)
            sta[++Top]=i;
    }else
    {
        if(cntodd==1)sta[++Top]=no; 
        for(int i=2;i<=n;i++)
        {
            if(query(dfn[i],1,n,1)==cntodd&&!query2(dfn[i],1,n,1))
                sta[++Top]=col[i];
        }
        sort(sta+1,sta+Top+1);
    }
    printf("%d\n",Top);
    for(int i=1;i<=Top;i++)
        printf("%d ",sta[i]);
    puts("");
}