本文轉載至哲學與程序博客：http://zhexue.sinaapp.com/?p=105
題目來源：ZOJ1992
對于一個混合圖，即有有向邊又有無向邊的圖，判斷是否存在一條歐拉回路。
解法（轉）：混合圖歐拉回路用的是網絡流。把該圖的無向邊隨便定向，計算每個點的入度和出度。如果有某個點出入度之差為奇數，那么肯定不存在歐拉回路。因為歐拉回路要求每點入度 = 出度，也就是總度數為偶數，存在奇數度點必不能有歐拉回路。現在每個點入度和出度之差均為偶數。將這個偶數除以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連接的點，自然早在開始就已經滿足入 = 出了。那么在網絡流過程中，這些點屬于“中間點”。我們知道中間點流量不允許有累積的，這樣，進去多少就出來多少，反向之后，自然仍保持平衡。所以，就這樣，混合圖歐拉回路問題，解了。
#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<string.h>
using namespace std;
#define N 205
#define MAXN N
#define inf 100000000
int map[N][N];
int flow[N][N];
int max_flow(int n,int mat[][MAXN],int source,int sink,int flow[][MAXN]){
int pre[MAXN],que[MAXN],d[MAXN],p,q,t,i,j;
if (source==sink) return inf;
for (i=0;i<n;i++)
for (j=0;j<n;flow[i][j++]=0);
for (;;){
for (i=0;i<n;pre[i++]=0);
pre[t=source]=source+1,d[t]=inf;
for (p=q=0;p<=q&&!pre[sink];t=que[p++])
for (i=0;i<n;i++)
if (!pre[i]&&(j=mat[t][i]-flow[t][i]))
pre[que[q++]=i]=t+1,d[i]=d[t]<j?d[t]:j;
else if (!pre[i]&&(j=flow[i][t]))
pre[que[q++]=i]=-t-1,d[i]=d[t]<j?d[t]:j;
if (!pre[sink]) break;
for (i=sink;i!=source;)
if (pre[i]>0)
flow[pre[i]-1][i]+=d[sink],i=pre[i]-1;
else
flow[i][-pre[i]-1]-=d[sink],i=-pre[i]-1;
}
for (i=0;i<n;i++)
if(mat[source][i] > flow[source][i])
return 0;
return 1;
}
int main()
{
int T, degin[N],degout[N], n, m, x, y, z;
int flag;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&m);
memset(degin,0,sizeof(degin));
memset(degout,0,sizeof(degout));
memset(map,0,sizeof(map));
for(int i = 0; i < m; i++){
scanf("%d%d%d",&x,&y,&z);
degout[x]++;
degin[y]++;
if(!z)map[x][y]++;
}
flag = 0;
for(int i = 1; i <= n && !flag; i++){
if((degin[i]+degout[i])%2)flag=1;
if(degin[i] > degout[i])
map[i][n+1] = (degin[i]-degout[i])/2;
else
map[0][i] = (degout[i]-degin[i])/2;
}
if(flag)
printf("impossible\n");
else{
if(max_flow(n+2,map,0,n+1,flow))
printf("possible\n");
else
printf("impossible\n");
}
}
return 0;
}