POJ 1637 混合图的欧拉回路

       混合图欧拉回路用的是网络流。

　　把该图的无向边随便定向，计算每个点的入度和出度，包括那些定向边所关联的点。如果有某个点出入度之差为奇数，那么肯定不存在欧拉回路。因为欧拉回路要求每点入度= 出度，也就是总度数为偶数，存在奇数度点必不能有欧拉回路。

　　好了，现在每个点入度和出度之差均为偶数。那么将这个偶数除以2，得x。也就是说，对于每一个点，只要将x条边改变方向（入>出就是变入，出>入就是变出），就能保证出= 入。如果每个点都是出 = 入，那么很明显，该图就存在欧拉回路。

　 　现在的问题就变成了：我该改变哪些边，可以让每个点出 =入？构造网络流模型。首先，有向边是不能改变方向的，要之无用，删。一开始不是把无向边定向了吗？定的是什么向，就把网络构建成什么样，边长容量上限1。另新建s和t。对于入 > 出的点u，连接边(u, t)、容量为x，对于出> 入的点v，连接边(s,v)，容量为x（注意对不同的点x不同）。之后，察看是否有满流的分配。有就是能有欧拉回路，没有就是没有。欧拉回路是哪个？查看流值分配，将所有流量非0（上限是1，流值不是0就是1）的边反向，就能得到每点入度 = 出度的欧拉图。

　 　由于是满流，所以每个入> 出的点，都有x条边进来，将这些进来的边反向，OK，入 = 出了。对于出> 入的点亦然。那么，没和s、t连接的点怎么办？和s连接的条件是出 >入，和t连接的条件是入 > 出，那么这个既没和s也没和t连接的点，自然早在开始就已经满足入 =出了。那么在网络流过程中，这些点属于“中间点”。我们知道中间点流量不允许有累积的，这样，进去多少就出来多少，反向之后，自然仍保持平衡。

　　所以，就这样，混合图欧拉回路问题，解了。

代码：

#pragma comment(linker,"/STACK:102400000,102400000")
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
typedef long long ll;
const int maxn=5005;
const int maxm=1000005;
const int INF=1000000000;
struct EdgeNode
{
    int from;
    int to;
    int cost;
    int next;
}edge[maxm];
int head[maxn],cnt;
void add(int x,int y,int z)
{
    edge[cnt].from=x;edge[cnt].to=y;edge[cnt].cost=z;edge[cnt].next=head[x];head[x]=cnt++;
    edge[cnt].from=y;edge[cnt].to=x;edge[cnt].cost=0;edge[cnt].next=head[y];head[y]=cnt++;
}

void init()
{
    cnt=0;
    memset(head,-1,sizeof(head));
}

int S,T,n,m;
int d[maxn],gap[maxn],curedge[maxn],pre[maxn];
//curedge[]为当前弧数组，pre为前驱数组

int sap(int S,int T,int n)  //n为点数
{
    int cur_flow,flow_ans=0,u,tmp,neck,i;
    memset(d,0,sizeof(d));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(pre));
    for(i=0;i<=n;i++)curedge[i]=head[i]; //初始化当前弧为第一条邻接表
    gap[0]=n;
    u=S;
    while(d[S]<n)             //当d[S]>=n时，网络中肯定出现了断层
    {
        if(u==T)
        {
            cur_flow=INF;
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                           //增广成功，寻找瓶颈边
                if(cur_flow>edge[curedge[i]].cost)
                {
                    neck=i;
                    cur_flow=edge[curedge[i]].cost;
                }
            }
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                             //修改路径上的边容量
                tmp=curedge[i];
                edge[tmp].cost-=cur_flow;
                edge[tmp^1].cost+=cur_flow;
            }
            flow_ans+=cur_flow;
            u=neck;                     //下次增广从瓶颈边开始
        }
        for(i=curedge[u];i!=-1;i=edge[i].next)
            if(edge[i].cost&&d[u]==d[edge[i].to]+1)
               break;
        if(i!=-1)
        {
            curedge[u]=i;
            pre[edge[i].to]=u;
            u=edge[i].to;
        }
        else
        {
            if(0==--gap[d[u]])break;    //gap优化
            curedge[u]=head[u];
            for(tmp=n,i=head[u];i!=-1;i=edge[i].next)
                if(edge[i].cost)
                   tmp=min(tmp,d[edge[i].to]);
            d[u]=tmp+1;
            ++gap[d[u]];
            if(u!=S)u=pre[u];           //重标号并且从当前点前驱重新增广
        }
    }
    return flow_ans;
}

int dsc[maxn]; //出度入度差
int main()
{
    int t,x,y,k,i,flag,sum;
    scanf("%d",&t);
    while(t--){
        init();
        scanf("%d%d",&n,&m);
        memset(dsc,0,sizeof(dsc));
        for(i=1;i<=m;i++){
            scanf("%d%d%d",&x,&y,&k);
            dsc[x]--; dsc[y]++; //随意定向x->y
            if(!k)add(x,y,1); //怎么定向怎么建图
        }
        flag=1; S=0; T=n+1; sum=0;
        for(i=1;i<=n;i++)if(dsc[i]%2)flag=0;
        if(flag==0){printf("impossible\n");continue;}
        for(i=1;i<=n;i++){
            if(dsc[i]>0)add(i,T,dsc[i]/2);
            else if(dsc[i]<0)add(S,i,-dsc[i]/2),sum+=-dsc[i]/2;
        }
        if(sap(S,T,T+1)!=sum)printf("impossible\n");
        else printf("possible\n");
    }
    return 0;
}