POJ 1637 Sightseeing tour混合欧拉图的判定

网络流简直爆炸了。

以下解释转自kuangbin博客，爆炸了。

【混合图】
混合图（既有有向边又有无向边的图）中欧拉环、欧拉路径的判定需要借助网络流！

（1）欧拉环的判定：
一开始当然是判断原图的基图是否连通，若不连通则一定不存在欧拉环或欧拉路径（不考虑度数为0的点）。

其实，难点在于图中的无向边，需要对所有的无向边定向（指定一个方向，使之变为有向边），使整个图变成一个有向欧拉图（或有向半欧拉图）。若存在一个定向满足此条件，则原图是欧拉图（或半欧拉图）否则不是。关键就是如何定向？

首先给原图中的每条无向边随便指定一个方向（称为初始定向），将原图改为有向图G'，然后的任务就是改变G'中某些边的方向（当然是无向边转化来的，原混合图中的有向边不能动）使其满足每个点的入度等于出度。
设D[i]为G'中(点i的出度 - 点i的入度）。可以发现，在改变G'中边的方向的过程中，任何点的D值的奇偶性都不会发生改变（设将边<i, j>改为<j, i>，则i入度加1出度减1，j入度减1出度加1，两者之差加2或减2，奇偶性不变）！而最终要求的是每个点的入度等于出度，即每个点的D值都为0，是偶数，故可得：若初始定向得到的G'中任意一个点的D值是奇数，那么原图中一定不存在欧拉环！

若初始D值都是偶数，则将G'改装成网络：设立源点S和汇点T，对于每个D[i]>0的点i，连边<S, i>，容量为D[i]/2；对于每个D[j]<0的点j，连边<j, T>，容量为-D[j]/2；G'中的每条边在网络中仍保留，容量为1（表示该边最多只能被改变方向一次）。求这个网络的最大流，若S引出的所有边均满流，则原混合图是欧拉图，将网络中所有流量为1的中间边（就是不与S或T关联的边）在G'中改变方向，形成的新图G''一定是有向欧拉图；若S引出的边中有的没有满流，则原混合图不是欧拉图。

为什么能这样建图？
考虑网络中的一条增广路径S-->i-->...-->j-->T，将这条从i到j的路径在G'中全部反向，则：i的入度加1出度减1，j的入度减1出度加1，路径中其它点的入度出度均不变。而i是和S相连的，因此初始D[i]>0，即i的出度大于入度，故这样反向之后D[i]减少2；同理，j是和T相连的，这样反向之后D[j]增加2。因此，若最大流中边<S, i>满流（流量为初始D[i]/2），此时D[i]值就变成了0，也就是i的入度等于出度。因此只要使所有S引出的边全部满流，所有初始D值>0的点的D值将等于0，又因为将边变向后所有点的D值之和不变，所有初始D值小于0的点的D值也将等于0，而初始D值等于0的D点既不与S相连也不与T相连，所以它们是网络中的中间点，而中间点的流入量等于流出量，故它们的入度和出度一直不变，即D值一直为0。因此，整个图G'成为欧拉图。

（2）欧拉路径的判定：
首先可以想到的是枚举欧拉路径的起点i和终点j，然后在图中添加边<j, i>，再求图中是否有欧拉回路即可。但是，该算法的时间复杂度达到了O(M * 最大流的时间)，需要优化。
前面已经说过，在将边变向的过程中任何点的D值的奇偶性都不会改变，而一个有向图有欧拉路径的充要条件是基图连通且有且只有一个点的入度比出度少1（作为欧拉路径的起点），有且只有一个点的入度比出度多1（作为终点），其余点的入度等于出度。这就说明，先把图中的无向边随便定向，然后求每个点的D值，若有且只有两个点的初始D值为奇数，其余的点初始D值都为偶数，则有可能存在欧拉路径（否则不可能存在）。进一步，检查这两个初始D值为奇数的点，设为点i和点j，若有D[i]>0且D[j]<0，则i作起点j作终点（否则若D[i]与D[j]同号则不存在欧拉路径），连边<j, i>，求是否存在欧拉环即可（将求出的欧拉环中删去边<j, i>即可）。这样只需求一次最大流。 

好像CSDN又啥比了

#include <iostream>
#include<time.h>
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<string>
#include<map>
#include<vector>
#include <algorithm>
#include <queue>


#define INF 0x3f3f3f3f


//** I/O Accelerator Interface .. **/
#define g (c=getchar())
#define d isdigit(g)
#define p x=x*10+c-'0'
#define n x=x*10+'0'-c
#define pp l/=10,p
#define nn l/=10,n
template<class T> inline T& RD(T &x)
{
    char c;
    while(!d);
    x=c-'0';
    while(d)p;
    return x;
}
template<class T> inline T& RDD(T &x)
{
    char c;
    while(g,c!='-'&&!isdigit(c));
    if (c=='-')
    {
        x='0'-g;
        while(d)n;
    }
    else
    {
        x=c-'0';
        while(d)p;
    }
    return x;
}
inline double& RF(double &x)      //scanf("%lf", &x);
{
    char c;
    while(g,c!='-'&&c!='.'&&!isdigit(c));
    if(c=='-')if(g=='.')
        {
            x=0;
            double l=1;
            while(d)nn;
            x*=l;
        }
        else
        {
            x='0'-c;
            while(d)n;
            if(c=='.')
            {
                double l=1;
                while(d)nn;
                x*=l;
            }
        }
    else if(c=='.')
    {
        x=0;
        double l=1;
        while(d)pp;
        x*=l;
    }
    else
    {
        x=c-'0';
        while(d)p;
        if(c=='.')
        {
            double l=1;
            while(d)pp;
            x*=l;
        }
    }
    return x;
}
#undef nn
#undef pp
#undef n
#undef p
#undef d
#undef g


using namespace std;


const int MAXN = 1010;//点数
const int MAXM = 400010;//正向弧和反向弧的数量


struct Edge
{
    int to,next,cap,flow;
} edge[MAXM];


int Top;
int head[MAXN];
int gap[MAXN],dep[MAXN],cur[MAXN];


void Init()
{
    Top = 0;
    memset(head,-1,sizeof(head));
}


void Link(int u,int v,int w)
{
    edge[Top].to = v,edge[Top].cap = w,edge[Top].flow = 0;
    edge[Top].next = head[u],head[u] = Top++;


    edge[Top].to = u,edge[Top].cap = 0,edge[Top].flow = 0;
    edge[Top].next = head[v],head[v] = Top++;
}


int Q[MAXN];


void BFS(int S,int E)
{
    memset(dep,-1,sizeof(dep));
    memset(gap,0,sizeof(gap));
    int u,v,i,Front = 0,Rear = 0;
    dep[E] = 0;
    Q[Rear++] = E;
    gap[0] = 1;


    while(Front != Rear)
    {
        u = Q[Front++];


        for(i = head[u]; i != -1; i = edge[i].next)
        {
            v = edge[i].to;
            if(dep[v] != -1) continue;
            Q[Rear++] = v;
            dep[v] = dep[u] + 1;
            gap[dep[v]]++;
        }
    }
}


int S[MAXN];


int ISAP(int Start,int End,int N)
{
    BFS(Start,End);
    if(dep[Start] == -1)
        return 0;
    memcpy(cur,head,sizeof(head));
    int Top = 0;
    int u = Start,v;
    int flow = 0;
    int Min;
    int inser;
    int i;
    while(dep[Start] < N)
    {
        if(u == End){
            Min = INF,inser = -1;


            for(i = 0;i < Top; ++i)
                if(Min > edge[S[i]].cap - edge[S[i]].flow)
                    inser = i,Min = edge[S[i]].cap - edge[S[i]].flow;


            for(i = 0 ; i < Top; ++i)
                edge[S[i]].flow += Min,edge[S[i]^1].flow -= Min;


            flow += Min;
            Top = inser;
            u = edge[S[Top]^1].to;
            continue;
        }


        bool flag = false;
        for(i =  cur[u];i != -1 ;i = edge[i].next)
        {
            v = edge[i].to;
            if(edge[i].cap-edge[i].flow && dep[v] + 1  == dep[u]){
                flag = true;
                cur[u] = i;
                break;
            }
        }


        if(flag){
            S[Top++] = cur[u];
            u = v;
            continue;
        }
        int Min = N;
        for(i = head[u] ; i != -1;i = edge[i].next)
        {
            if(edge[i].cap - edge[i].flow && dep[edge[i].to] < Min)
            {
                Min = dep[edge[i].to];
                cur[u] = i;
            }
        }
        gap[dep[u]]--;
        if(gap[dep[u]] == 0)
            return flow;
        dep[u] = Min+1;
        gap[dep[u]]++;
        if(u != Start) u = edge[S[--Top]^1].to;
    }
    return flow;
}


int deg[MAXN];


int main()
{
    int T;


    scanf("%d",&T);


    int u,v,w,n,m,i,sum;


    while(T--)
    {
        scanf("%d %d",&n,&m);


        int S = 0,E = n+1;


        memset(deg,0,sizeof(deg));


        Init();


        while(m--)
        {
            scanf("%d %d %d",&u,&v,&w);


            if(w == 1)
                deg[v]++,deg[u]--;
            else
                Link(u,v,1),deg[v]++,deg[u]--;
        }


        for(i = 1;i <= n && (abs(deg[i])&1) == 0; ++i)
            ;


        if(i <= n)
        {
            puts("impossible");
            continue;
        }
        for(i = 1,sum = 0;i <= n; ++i)
        {
            if(deg[i] < 0)
                Link(S,i,abs(deg[i])>>1);
            else
                sum += (deg[i]>>1),Link(i,E,deg[i]>>1);
        }


        if(ISAP(S,E,E+1) == sum)
            puts("possible");
        else
            puts("impossible");


    }


    return 0;
}