BZOJ1585: [Usaco2009 Mar]Earthquake Damage 2 地震伤害

n<=3000个点m<=20000条无向边的图，有p<=n个出发点，每个出发点都不可拆，现拆一些点使每个出发点都不能到达点1，求最小点数。

简单的最小割。每个点拆成两个x和y，无向边A--B即Ay->Bx,By->Ax，每个点拆成的x和y再连边容量1，然后建超级源向p个点连边，最大流，没了。

错误。没看题。p个点不可割。WA了一发。正确做法只要把p个点的x->y的边容量设inf即可。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<cstdlib>
#include<math.h>
//#include<iostream>
using namespace std;

int n,m,p;
#define maxn 6011
#define maxm 100011
const int inf=0x3f3f3f3f;
struct Network
{
    struct Edge{int to,next,cap,flow;}edge[maxm];
    int first[maxn],le,n;
    void clear(int n)
    {
        memset(first,0,sizeof(first));
        le=2;this->n=n;
    }
    void in(int x,int y,int cap)
    {
        Edge &e=edge[le];
        e.to=y;e.cap=cap;e.flow=0;
        e.next=first[x];first[x]=le++;
    }
    void insert(int x,int y,int cap)
    {
        in(x,y,cap);
        in(y,x,0);
    }
    int s,t,que[maxn],head,tail,dis[maxn],cur[maxn];
    bool bfs()
    {
        memset(dis,0,sizeof(dis));
        dis[s]=1;
        que[head=(tail=1)-1]=s;
        while (head!=tail)
        {
            const int now=que[head++];
            for (int i=first[now];i;i=edge[i].next)
            {
                Edge &e=edge[i];
                if (e.cap>e.flow && !dis[e.to])
                {
                    dis[e.to]=dis[now]+1;
                    que[tail++]=e.to;
                }
            }
        }
        return dis[t];
    }
    int dfs(int x,int a)
    {
        if (x==t || !a) return a;
        int flow=0,f;
        for (int &i=cur[x];i;i=edge[i].next)
        {
            Edge &e=edge[i];
            if (dis[e.to]==dis[x]+1 && (f=dfs(e.to,min(a,e.cap-e.flow)))>0)
            {
                e.flow+=f;
                edge[i^1].flow-=f;
                flow+=f;
                a-=f;
                if (!a) break;
            }
        }
        return flow;
    }
    int Dinic(int s,int t)
    {
        this->s=s,this->t=t;
        int flow=0;
        while (bfs())
        {
            for (int i=1;i<=n;i++) cur[i]=first[i];
            flow+=dfs(s,0x3f3f3f3f);
        }
        return flow;
    }
}g;
int x,y;int vis[maxn];
int main()
{
    scanf("%d%d%d",&n,&m,&p);
    g.n=(n<<1)^1;int s=g.n,t=1;
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        g.insert(x+n,y,inf);
        g.insert(y+n,x,inf);
    }
    memset(vis,0,sizeof(vis));
    for (int i=1;i<=p;i++)
    {
        scanf("%d",&x);vis[x]=1;
        g.insert(s,x,inf);
    }
    for (int i=1;i<=n;i++) if (!vis[i]) g.insert(i,i+n,1);else g.insert(i,i+n,inf);
    printf("%d
",g.Dinic(s,t));
    return 0;
}