poj1637 Sightseeing tour（混合图欧拉回路，最大流）

判断混合图是否存在欧拉回路。

首先有关欧拉回路的一些定义：
如果图G中的一个路径包括每个边恰好一次，则该路径称为欧拉路径(Euler path)。
如果一个回路是欧拉路径，则称为欧拉回路(Euler circuit)。
具有欧拉回路的图称为欧拉图（简称E图）。具有欧拉路径但不具有欧拉回路的图称为半欧拉图。

有关欧拉回路的一些性质：
以下判断基于此图的基图连通。
无向图存在欧拉回路的充要条件：
一个无向图存在欧拉回路，当且仅当该图所有顶点度数都为偶数,且该图是连通图。
有向图存在欧拉回路的充要条件：
一个有向图存在欧拉回路，所有顶点的入度等于出度且该图是连通图。
混合图存在欧拉回路条件：
要判断一个混合图G（V,E）（既有有向边又有无向边）是欧拉图，方法如下：
假设有一张有向图G’，在不论方向的情况下它与G同构。并且G’包含了G的所有有向边。那么如果存在一个图G’使得G’存在欧拉回路，那么G就存在欧拉回路。

至于此题题解：orz网络流建模汇总

大概就是先把无向边随便定向，计算每个点的入度与出度之差，如果是奇数则一定不可行。如果全是偶数，用网络流来判断是否可行。首先原图中的有向边，方向是确定的，无用了，都删掉。然后对于无向边，就按你给他定的方向建边，容量为1。然后对于每个点，设du[i]表示入度-出度，如果入度>出度，则我们要把一些入边反过来变成出边才能使得入=出，需要翻几条边呢？du[i]/2.所以建边i->T,du[i]/2.同理，对于出度大于入度的，需要翻-du[i]/2条出边，建边S->i,-du[i]/2.看能否满流即可。如果满流了一定可行，一条可行的欧拉回路就是你把流量非0的双向边方向都反过来，就ok啦。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 210
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;
}
int n,m,h[N],num=1,cur[N],lev[N],T=201,du[N];
struct edge{
    int to,next,val;
}data[2500];
struct Edge{
    int x,y;bool op;
}e[1010];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;
    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=0;
}
inline bool bfs(){
    queue<int>q;memset(lev,0,sizeof(lev));
    q.push(0);lev[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(lev[y]||!data[i].val) continue;
            lev[y]=lev[x]+1;if(y==T) return 1;q.push(y);
        }
    }return 0;
}
inline int dinic(int x,int low){
    if(x==T) return low;int tmp=low;
    for(int &i=cur[x];i;i=data[i].next){
        int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue;
        int res=dinic(y,min(tmp,data[i].val));
        if(!res) lev[y]=0;else tmp-=res,data[i].val-=res,data[i^1].val+=res;
        if(!tmp) return low;
    }return low-tmp;
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();
    while(tst--){
        n=read();m=read();bool flag=1;int ans=0,tot=0;
        memset(h,0,sizeof(h));num=1;memset(du,0,sizeof(du));
        for(int i=1;i<=m;++i) e[i].x=read(),e[i].y=read(),e[i].op=read(),du[e[i].y]++,du[e[i].x]--;
        for(int i=1;i<=n;++i) if(du[i]%2) flag=0;
        if(!flag){puts("impossible");continue;}
        for(int i=1;i<=n;++i){
            if(du[i]>0) add(i,T,du[i]>>1),tot+=du[i]>>1;
            if(du[i]<0) add(0,i,(-du[i])>>1);
        }for(int i=1;i<=m;++i)
            if(!e[i].op) add(e[i].x,e[i].y,1);
        while(bfs()){memcpy(cur,h,sizeof(h));ans+=dinic(0,inf);}
        puts(ans==tot?"possible":"impossible");
    }return 0;
}