题意:在一个无向图中,一个人要从 A 点赶往 B 点,
之后再赶往 C 点,且要求中途不 能多次经过同一个点。问是否存在这样的路线。(3 <= N <= 30011, 1 <= M <= 50011)
解题:不错的无向图拆点最大流。
由于要求每个点只通过一次,可以把点约束转化为边约束。边流量为1就是了。
S是与2‘相连,而不是2相连。!在这道题目,我们需要从1到2再到3,不重复经过点。现在化成边了,也就是不重复经过边。很明显,从2进行两次增广,第一次到1,第二次到3.两次增广都造访了2号顶点,也就是如果S直接与2相连,由于2号点的约束,导致2号点只能造访一次,也就是2-2’这条边,所有只能增广一次。将S与2‘相连,就能进行多次增广了。
还有一个坑点是:存在一些不必要的边,需要自己判断一下
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=501110;
const int maxm=5101100;
const int inf=0x3f3f3f3f;
int head[maxn],cnt=0;
struct edge
{
int v,nxt,w;
}edge[maxm*4+100];
void add_edge(int u,int v,int w)
{
edge[cnt].v=v;
edge[cnt].w=w;
edge[cnt].nxt=head[u];
head[u]=cnt++;
edge[cnt].v=u;
edge[cnt].w=0;
edge[cnt].nxt=head[v];
head[v]=cnt++;
}
int numh[maxn],h[maxn],curedge[maxn],pre[maxn];
//最大流板子
int sap(int s,int t,int n)
{
memset(numh,0,sizeof(numh));
memset(h,0,sizeof(h));
memset(pre,-1,sizeof(pre));
int cur_flow,flow_ans=0,u,tmp,neck,i;
for(i=1;i<=n;i++)
{
curedge[i]=head[i];
}
numh[0]=n;
u=s;
while(h[s]<n)
{
if(u==t)
{
cur_flow=inf;
for(i=s;i!=t;i=edge[curedge[i]].v)
{
if(cur_flow>edge[curedge[i]].w)
{
neck=i;
cur_flow=edge[curedge[i]].w;
}
}
for(i=s;i!=t;i=edge[curedge[i]].v)
{
tmp=curedge[i];
edge[tmp].w-=cur_flow;
edge[tmp^1].w+=cur_flow;
}
flow_ans+=cur_flow;
u=neck;
}
for(i=curedge[u];i!=-1;i=edge[i].nxt)
{
if(edge[i].w&&h[u]==h[edge[i].v]+1)
{
break;
}
}
if(i!=-1)
{
curedge[u]=i;
pre[edge[i].v]=u;
u=edge[i].v;
}
else
{
if(0==--numh[h[u]])
break;
curedge[u]=head[u];
for(tmp=n,i=head[u];i!=-1;i=edge[i].nxt)
{
if(edge[i].w)
{
tmp=min(tmp,h[edge[i].v]);
}
}
h[u]=tmp+1;
++numh[h[u]];
if(u!=s)
u=pre[u];
}
}
return flow_ans;
}
int main ()
{
int t;
scanf("%d",&t);
while(t--)
{
int n,m;
memset(head,-1,sizeof(head));
cnt=0;
scanf("%d%d",&n,&m);
add_edge(2*n+1,2+n,2);
for(int i=1;i<=m;i++)
{
int xx,yy;
scanf("%d%d",&xx,&yy);
if(xx>n||xx<1||yy>n||yy<1)
continue;
add_edge(xx+n,yy,1);
add_edge(yy+n,xx,1);
}
for(int i=1;i<=n;i++)
{
add_edge(i,i+n,1);
add_edge(i+n,i,1);
}
add_edge(n+1,2*n+2,1);
add_edge(3+n,2*n+2,1);
int res=sap(2*n+1,2*n+2,2*n+2);
//printf("%d\n",res);
if(res>=2)
printf("YES\n");
else
printf("NO\n");
}
}