题目链接:点这里!!!
题意:给你n(3<=n<=100)个点m(m<=n*(n-1)/2)t条无向边,让你找一条a->b的一条路线,且必须经过mid点,每个点最多走一遍,输出这条路径。
好久没写网络流的题了,遇到这个题就有点傻比了,比赛的时候没做出来,赛后补了一下知识点,就差不多啦!!
题解:
每个点最多只能经过一次,我们可以把i点拆成i和i+n两个点,且在i和i+n两个点间连一条容量为1(i=mid的时候容量为2)的有向边,假设x与y边相连,则建一条x+n->y容量为1的边。
然后我们设超级源为s,超级汇为t,建一条s->mid为2容量的边,a->t,b->t容量为1的边,题目就可以简化成mid到a,b两点路径不存在交叉。然后跑一边最大流。
我们得到了跑完最大流的图,我们知道如果边有经过的话,那么他是有流量的。我们可以先在图上dfs下,从s到t然后将有流量的边标记下或者减一下,然后再从s到t,就得到了我们想要的答案。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define LL long long
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
#pragma comment(linker, "/STACK:102400000000,102400000000")
const LL MOD = 1000000007;
const int N = 100+15;
const int maxn = 8e3+15;
const int letter = 130;
const int INF = 1000000;
const double pi=acos(-1.0);
const double eps=1e-10;
using namespace std;
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,head[2*N],tot;
int a,b,mid;
int dis[2*N],cur[2*N],vis[N*2];
int p[2*N];
vector<int>vp;
queue<int>q;
struct edges{
int to,next,ca,fl;
}e[8*N*N];
void addedges(int u,int v,int val){
e[tot].to=v,e[tot].ca=val,e[tot].fl=0,e[tot].next=head[u],head[u]=tot++;
e[tot].to=u,e[tot].ca=0, e[tot].fl=0,e[tot].next=head[v],head[v]=tot++;
}
bool bfs(int s,int t){
int x,to;
clr(dis,-1);
while(!q.empty()) q.pop();
dis[s]=0;
q.push(s);
while(!q.empty()){
x=q.front(),q.pop();
for(int i=head[x];i!=-1;i=e[i].next){
to=e[i].to;
if(dis[to]<0&&e[i].ca-e[i].fl>0){
dis[to]=dis[x]+1;
q.push(to);
}
}
}
return dis[t]!=-1;
}
int dfs(int t,int x,int f){
if(x==t) return f;
int w,used=0;
for(int i=cur[x];i!=-1;i=e[i].next){
int to=e[i].to;
if(dis[to]==dis[x]+1){
w=f-used;
w=dfs(t,to,min(w,e[i].ca-e[i].fl));
e[i].fl+=w;
e[i^1].fl-=w;
if(e[i].ca-e[i].fl>0) cur[x]=i;
used+=w;
if(used==f) return f;
}
}
if(!used) dis[x]=-1;
return used;
}
int maxflow(int s,int t){
int flow=0;
while(bfs(s,t)){
for(int i=0;i<=t;i++)cur[i]=head[i];
flow+=dfs(t,s,INF);
}
return flow;
}
void pp(int x){
if(x<=n&&x!=0) vp.pb(x);
if(x==2*n+2) return;
for(int i=head[x];i!=-1;i=e[i].next){
int to=e[i].to;
if(e[i].ca&&e[i].fl){
e[i].fl--;
pp(to);
return;
}
}
}
void solve(int s,int t){
vp.clear();
pp(s);
if(vp[vp.size()-1]==a){
for(int i=vp.size()-1;i>=0;i--){
if(i==vp.size()-1) printf("%d",vp[i]);
else printf(" %d",vp[i]);
}
vp.clear();
pp(s);
for(int i=1;i<vp.size();i++) printf(" %d",vp[i]);
puts("");
}
else {
pp(s);
reverse(vp.begin(),vp.end());
int flag=0;
for(int i=0;i<vp.size();i++){
if(vp[i]==b){
flag=i;
break;
}
else {
if(i==0) printf("%d",vp[i]);
else printf(" %d",vp[i]);
}
}
for(int i=vp.size()-2;i>=flag;i--){
printf(" %d",vp[i]);
}
puts("");
}
}
int main(){
int tc;
scanf("%d",&tc);
while(tc--){
tot=0,clr(head,-1),clr(p,0);
scanf("%d%d",&n,&m);
scanf("%d%d%d",&a,&b,&mid);
int s=0,t=2*n+2;
for(int i=1;i<=n;i++){
if(i==a||i==b){
addedges(i,i+n,1);
addedges(i+n,t,1);
}
else if(i==mid){
addedges(s,i,2);
addedges(i,i+n,2);
}
else {
addedges(i,i+n,1);
}
}
int x,y;
for(int i=0;i<m;i++){
scanf("%d%d",&x,&y);
addedges(x+n,y,1);
addedges(y+n,x,1);
}
maxflow(s,t);
solve(s,t);
}
return 0;
}