Intergalactic Map SPOJ - IM

题意：在一个无向图中，一个人要从 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");
    }
}