[POJ1637][网络流]Sightseeing tour

POJ1637

首先我们把无向边随便定向，得出每个点的入度in和出度out
在欧拉回路中，每个点的in和out必须相等
如果有一个点的$|in-out|$为奇数，那么不可能存在欧拉回路
现在问题转化为可以改变无向边的方向，使得每个点的in和out相等，有向边不能改变，不用考虑
无向边就直接建成一开始定向的方向，流量为1
对于一个点，如果出度大于入度，就从s到它连一条流量为差值/2的边，反之就从它到t连一条流量为差值/2的边
那么只要满流就行了
证明：对于每个向t连边的点，有差值/2条边有流，就把这些边反向就好了
因为入度和出度是对称的，所以对于s向它连边的点也是这样
然后既没有往s也没有往t的点是中间的点，不能保留流量，所以是平衡的

构造方案：看哪些边有流，有流的边就反向，就做完了

Code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
inline int read(){
    int res=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-f;ch=getchar();}
    while(isdigit(ch)) {res=(res<<1)+(res<<3)+(ch^48);ch=getchar();}
    return res*f;
}
const int N=205,M=5005,INF=1e9;
int vis[M],head[N],nxt[M],c[M],tot;
int n,m,s,t;
int in[N],out[N];
inline void add(int x,int y,int z){
    vis[++tot]=y;nxt[tot]=head[x];head[x]=tot;c[tot]=z;
    vis[++tot]=x;nxt[tot]=head[y];head[y]=tot;c[tot]=0;
}
int st[N];
inline bool bfs(){
    memset(st,-1,sizeof(st));
    queue<int>q;
    q.push(s);st[s]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=nxt[i]){
            if(c[i] && st[vis[i]]==-1){
                q.push(vis[i]);
                st[vis[i]]=st[x]+1;
                if(vis[i]==t) return 1;
            }
        }
    }
    return 0;
}
int dfs(int v,int flow){
    if(v==t) return flow;
    int rest=flow,k;
    for(int i=head[v];i && rest;i=nxt[i])
        if(c[i] && st[vis[i]]==st[v]+1){
             k=dfs(vis[i],min(rest,c[i]));
             if(!k) st[vis[i]]=-1;
             rest-=k;c[i]-=k;c[i^1]+=k; 
        }
    return flow-rest;
}
inline int dinic(){
    int res,mxflow=0;
    while(bfs())
        while(res=dfs(s,INF)) mxflow+=res;
    return mxflow;
}
int main(){
    int T=read();
    while(T--){
        memset(head,0,sizeof(head));
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        tot=1;
        n=read();m=read();
        s=0;t=n+1;
        for(int x,y,z,i=1;i<=m;i++){
            x=read();y=read();z=read(); 
            ++in[y],++out[x];
            if(!z) add(x,y,1);
        }
        int ans=0,pt=1,all=0;
        for(int i=1;i<=n && pt;i++){
            if(in[i]!=out[i]){
                if(in[i]>out[i]){
                    int del=in[i]-out[i];
                    if(del&1) pt=0;
                    else add(i,t,del/2),all+=del; 
                }
                else{
                    int del=out[i]-in[i];
                    if(del&1) pt=0;
                    else add(s,i,del/2); 
                }
            }
        }
        all/=2;
        if(!pt) puts("impossible");
        else{
            ans=dinic();
            if(ans==all) puts("possible");
            else puts("impossible");
        }
    }
    return 0;
}