[bzoj 4424]Cf19E Fairy

给定 n 个点，m 条边的无向图，可以从图中删除一条边，问删除哪些边可以使图变成 一个二分图。

这道题还是不错的。
看到无向图的题目，应该一开始要往dfs树上面想。其实是因为最近吃够教训了
那么根据上面套路，我们先弄出一棵生成树。对于那些非树边，如果x与y之间的点数为偶数，称之为偶边，如果为奇数，就为奇边。
二分图的充要条件是不存在奇环，且点数要在2个或以上（这个出题人好像没卡），所以我们的目标就是让所有奇环都被破坏掉。而我们通过手玩就可以发现，加进非树边后，奇边本身自成一个奇环，要想破坏它，要么删除奇边，要么删除x到y的树上路径上的任意一条边。
但非树边之间也会有影响，偶边跟偶边之间只会形成偶环，奇边跟奇边之间也只会形成偶环，但是奇边和偶边就会形成新的奇环（边与边之间的意思是两条边有交点），所以偶边的x到y的树上路径上的边删除了也没用，就要排除掉。而题目要求我们只能删除一条边，所以所有的奇边的x到y的树上路径都要经过这条边。满足了这两个条件就是合法边了。但是如何快速判断呢？用离线lca加上树上差分就可以做到O(n)了。
还有一些坑点，注意重边和自环，还有原图不一定联通，要好好关注细节。那这道题就做完了。
代码写的贼丑，关键是一开始想错了，发现之后就在上面修修补补，其实应该重构代码，这样也可以更简洁。但我懒啊，AC就完事了

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
inline void write(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
}
inline void pr1(int x){write(x),putchar(' ');}
inline void pr2(int x){write(x),puts("");}
vector<int>q[1000010],qq[1000010],he[1000010];
struct node
{
    int x,y,next,id;
}a[2000010];int len,last[1000010];
inline void ins(int x,int y,int id)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].id=id;
    a[len].next=last[x];last[x]=len;
}
struct edge
{
    int x,y,id;
}e[1000010];
int fa[1000010],f[1000010];
inline int findfa(int x)
{
    if(fa[x]!=x)fa[x]=findfa(fa[x]);
    return fa[x];
}
int rt,id,v[1000010],s[2][1000010],d[1000010],cnt[1000010],p[1000010];
int top,ans[1000010];
bool wy[1000010],jd;
inline void dfs(int x)
{
    v[x]=1;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(v[y])continue;
        d[y]=d[x]+1,f[y]=a[k].id;
        dfs(y);
        fa[y]=x;
    }
    for(int i=0;i<q[x].size();i++)
    {
        int y=q[x][i];
        if(x==y)cnt[rt]++,p[rt]=qq[x][i];
        if(v[y]==2)
        {
            int lca=findfa(y),w=(d[x]+d[y]-2*d[lca]+1)&1;
            if(w==1)cnt[rt]++,p[rt]=qq[x][i],jd=true;
            else he[rt].push_back(qq[x][i]);
            s[w][x]++,s[w][y]++;s[w][lca]-=2;
        }
    }v[x]=2;
}
inline void solve(int x)
{
    v[x]=1;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(v[y])continue;
        solve(y);
        s[0][x]+=s[0][y],s[1][x]+=s[1][y];
    }
}
int main()
{
    //freopen("4424.in","r",stdin);
    //freopen("4424.out","w",stdout);
    int n=read(),m=read(),tot=0;
    for(int i=1;i<=n;i++)fa[i]=i;
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        int fx=findfa(x),fy=findfa(y);
        if(fx!=fy){fa[fx]=fy;ins(x,y,i),ins(y,x,i);}
        else {tot++;e[tot].x=x,e[tot].y=y,e[tot].id=i;}
    }
    for(int i=1;i<=n;i++)fa[i]=i;
    for(int i=1;i<=tot;i++)
    {
        q[e[i].x].push_back(e[i].y);
        if(e[i].x!=e[i].y)q[e[i].y].push_back(e[i].x);
        qq[e[i].x].push_back(e[i].id);
        if(e[i].x!=e[i].y)qq[e[i].y].push_back(e[i].id);
    }
    for(int i=1;i<=n;i++)if(v[i]==0)rt=i,dfs(i),wy[i]=true;
    memset(v,0,sizeof(v));
    for(int i=1;i<=n;i++)if(v[i]==0)solve(i);
    memset(v,0,sizeof(v));
    if(jd==false)
    {
        for(int i=1;i<=n;i++)
        {
            if(wy[i]==true)
            {
                for(int j=0;j<he[i].size();j++)ans[++top]=he[i][j];
            }
        }
    }
    for(int i=1;i<=n;i++)if(cnt[i]==1)ans[++top]=p[i];
    for(int i=1;i<=n;i++)
    {
        int uf=findfa(i);
        if(cnt[uf]==0)
        {
            if(f[i] && jd==false)ans[++top]=f[i];
        }
        else if(!s[0][i] && s[1][i]==cnt[uf])ans[++top]=f[i],v[uf]=1;
    }
    bool bk=true;
    for(int i=1;i<=n;i++)if(cnt[i]>1 && !v[i])bk=false;
    int num=0;
    for(int i=1;i<=n;i++)if(cnt[i])num++;
    if(num>1)bk=false;
    if(bk==false)pr2(0);
    else
    {
        pr2(top);sort(ans+1,ans+top+1);
        for(int i=1;i<=top;i++)pr1(ans[i]);
    }
    return 0;
}