HDU 6311 Cover（欧拉路径）

Description

给出一个$n$个点$m$条边的无向图，要求用最少的不交的路径覆盖该图，输出这些路径

Input

多组用例，每组用例首先输入两个整数$n,m$表示点数和边数，之后$m$行每行两个整数$u,v$表示一条无向边

$(1\le n,m\le 10^5)$

Output

输出所用最少路径数以及每条路径所经过边的编号和方向（若和输入方向相反则在编号前加负号）

Sample Input

3 3
1 2
1 3
2 3

Sample Output

1
3 1 3 -2

Solution

对于一个连通分支，假设有$odd$个奇度点，由于一条路径至多解决两个奇度点，故至少需要$min(\frac{odd}{2},1)$条路径覆盖，$odd=0,2$时显然该连通块存在欧拉路径，此时一条路径即可覆盖，$odd>2$时，奇度点显然有偶数个，留下两个奇度点，其余奇度点两两连虚边，那么加完边后的图必然存在欧拉路径，求出欧拉路径后沿所连$\frac{odd-2}{2}$条虚边将路径分割即得到$\frac{odd}{2}$条路径

Code

#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
#define maxn 100005
int n,m,head[maxn],tot,du[maxn],vis[maxn];
vector<int>odd,sta,ans[maxn];
struct Edge
{
    int to,next,id,flag;
}e[3*maxn];
void init()
{
    tot=0;
    memset(head,-1,sizeof(head));
    memset(du,0,sizeof(du));
    memset(vis,0,sizeof(vis));
}
void add(int u,int v,int id)
{
    e[tot].to=v,e[tot].flag=0,e[tot].id=id;
    e[tot].next=head[u],head[u]=tot++;
}
void dfs1(int u)
{
    vis[u]=1;
    if(du[u]&1)odd.push_back(u);
    for(int i=head[u];~i;i=e[i].next)
    {
        int v=e[i].to;
        if(!vis[v])dfs1(v);
    }
}
void dfs2(int u)
{
    for(int i=head[u];~i;i=e[i].next)
        if(!e[i].flag)
        {
            e[i].flag=e[i^1].flag=1;
            dfs2(e[i].to);
            sta.push_back(e[i].id);
        }
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        init();
        for(int i=1;i<=m;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            add(u,v,i),add(v,u,-i);
            du[u]++,du[v]++;
        }
        int res=0;
        for(int i=1;i<=n;i++)
            if(!vis[i])
            {
                odd.clear();
                dfs1(i);
                if(odd.size()>2)
                    for(int j=2;j<odd.size();j+=2)
                        add(odd[j],odd[j+1],0),add(odd[j+1],odd[j],0);
                sta.clear();
                int s=odd.size()?odd[0]:i;
                dfs2(s);
                ans[res].clear();
                if(sta.size())
                {
                    for(int j=sta.size()-1;j>=0;j--)
                        if(sta[j])
                        {
                            if(j<sta.size()-1&&!sta[j+1])ans[++res].clear();
                            ans[res].push_back(sta[j]);
                        }
                    res++;
                }
            }
        printf("%d\n",res);
        for(int i=0;i<res;i++)
        {
            printf("%d",ans[i].size());
            for(int j=0;j<ans[i].size();j++)printf(" %d",ans[i][j]);
            printf("\n");
        }
    }
    return 0;
}