混合图是指既有有向边又有无向边的一类图
我们一般处理无向边是将它拆成两条有向边,而在求欧拉回路时不能这么做,因为一条边只能经过一次
因此用网络流解决这个问题
黑书p324讲了两种做法
我写的第二种
首先将无向边任意定向,将其当做一条有向边,判断这个图是否可能存在欧拉回路(对于每个点,abs(入度-出度)是偶数)
接下来构建网络:
保留原来的无向边,方向为你定的方向,流量为1
对于点i
如果出度>入度,从st连一条边到i,流量为(出度-入度)/2
如果入度>出度,从i连一条边到ed,流量为(入度-出度)/2
对于这样的一个网络,跑网络流
最后判断是否满流(st连出的边是否都增广完毕)
是的话就存在欧拉回路 将那些没有被增广的弧反向,就得到了一个欧拉回路
否则不存在欧拉回路
POJ1637 裸欧拉回路
//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i)
#define rrep(i,a,b) for(int i=a,tt=b;i>=tt;--i)
#define erep(i,e,x) for(int i=x;i;i=e[i].next)
#define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define read() (strtol(ipos,&ipos,10))
#define sqr(x) ((x)*(x))
#define pb push_back
#define PS system("pause");
typedef long long ll;
typedef pair<int,int> pii;
const int oo=~0U>>1;
const double inf=1e20;
const double eps=1e-6;
string name="",in=".in",out=".out";
//Var
struct E
{
int next,node,cap;
}e[5000];
struct ED
{
int a,b,d;
}edge[1008];
queue<int> q;
int h[208],tot,dis[208],gap[208];
int n,m,st,ed,T,din[208],dout[208];
void add(int a,int b,int c)
{
e[++tot].next=h[a];e[tot].node=b;e[tot].cap=c;h[a]=tot;
e[++tot].next=h[b];e[tot].node=a;e[tot].cap=0;h[b]=tot;
}
void Init()
{
memset(h,0,sizeof h);tot=1;
memset(din,0,sizeof din);
memset(dout,0,sizeof dout);
scanf("%d%d",&n,&m);int a,b,c;
rep(i,1,m)
{
scanf("%d%d%d",&a,&b,&c);
dout[a]++,din[b]++;
edge[i].a=a;edge[i].b=b;edge[i].d=c;
}
}
bool BFS()
{
memset(dis,-1,sizeof dis);
rep(i,1,n+2)gap[i]=h[i];
int x,y;q.push(ed);dis[ed]=0;bool flag=false;
while(!q.empty())
{
x=q.front();q.pop();
erep(i,e,h[x])
{
y=e[i].node;
if(dis[y]!=-1||!e[i^1].cap)continue;
if(y==st)flag=true;
dis[y]=dis[x]+1;
q.push(y);
}
}
return flag;
}
int DFS(int now,int low)
{
if(now==ed)return low;
int ret=low,tmp,y;
for(int i=gap[now];i;i=e[i].next)
{
y=e[i].node;
if(!e[i].cap||dis[y]!=dis[now]-1)continue;
tmp=DFS(y,min(low,e[i].cap));
e[i].cap-=tmp;
e[i^1].cap+=tmp;
low-=tmp;
if(!low)break;
gap[now]=e[i].next;
}
if(low==ret)dis[now]=-1;
return ret-low;
}
bool Work()
{
st=n+1,ed=n+2;
rep(i,1,n)if(abs(dout[i]-din[i])&1)return false;
rep(i,1,n)
if(dout[i]>din[i])add(st,i,dout[i]-din[i]>>1);
else if(din[i]>dout[i])add(i,ed,din[i]-dout[i]>>1);
rep(i,1,m)if(!edge[i].d)add(edge[i].a,edge[i].b,1);
while(BFS())
while(DFS(st,oo));
erep(i,e,h[st])if(e[i].cap)return false;
return true;
}
int main()
{
//freopen((name+in).c_str(),"r",stdin);
//freopen((name+out).c_str(),"w",stdout);
for(scanf("%d",&T);T;--T)
{
Init();
puts(Work()?"possible":"impossible");
}
return 0;
}