【P2944】地震损失（最大流，洛谷）

绝对难度虚高的一题

看到题目，至少损坏几个房子，开始考虑最小割，建的是双向边，所以拆点，边权除了自己与自己的之外都连inf。然后把所有求救的点都连到超级源上，跑一遍最大流就可以了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define ll long long
#define inf 50000000
#define re register
using namespace std;
struct po
{
    int to,dis,nxt;
}edge[400001];
int head[100001],cur[100001],dep[10001],n,m,s,t,u,num=-1,x,y,l,tot,sum,k;
inline int read()
{
    int x=0,c=1;
    char ch=' ';
    while((ch>'9'||ch<'0')&&ch!='-')ch=getchar();
    while(ch=='-')c*=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-'0',ch=getchar();
    return x*c;
}
inline void add_edge(int from,int to,int dis)
{
    edge[++num].nxt=head[from];
    edge[num].to=to;
    edge[num].dis=dis;
    head[from]=num;
}
inline void add(int from,int to,int dis)
{
    add_edge(from,to,dis);
    add_edge(to,from,0);
}
inline bool bfs()
{
    memset(dep,0,sizeof(dep));
    queue<int> q;
    while(!q.empty())
    q.pop();
    dep[s]=1;
    q.push(s);
    while(!q.empty())
    {
        int now=q.front();
        q.pop();
        for(re int i=head[now];i!=-1;i=edge[i].nxt)
        {
            int v=edge[i].to;
            if(dep[v]==0&&edge[i].dis>0)
            {
                dep[v]=dep[now]+1;
                if(v==t)
                return 1;
                q.push(v); 
            }
        }
    }
    return 0;
}
inline int dfs(int u,int dis)
{
    if(u==t)
    return dis;
    int diss=0;
    for(re int& i=cur[u];i!=-1;i=edge[i].nxt)
    {
        int v=edge[i].to;
        if(dep[v]==dep[u]+1&&edge[i].dis!=0)
        {
            int check=dfs(v,min(dis,edge[i].dis));
            if(check>0)
            {
                diss+=check;
                dis-=check;
                edge[i].dis-=check;
                edge[i^1].dis+=check;
                if(dis==0) break;
            }
        }
    }
    return diss;
}
inline int dinic()
{
    int ans=0;
    while(bfs())
    {
        for(re int i=0;i<=2*n;i++)
        cur[i]=head[i];
        while(int d=dfs(s,inf))
        ans+=d;
    }
    return ans;
}
int main()
{
    memset(head,-1,sizeof(head));
    n=read();m=read();k=read();
    s=0;t=1;
    for(re int i=1;i<=m;i++)
    {
        x=read();y=read();
        if(x!=y)
        {
            add(x+n,y,inf);
            add(y+n,x,inf);
        }
    }
    for(re int i=1;i<=k;i++)
    {
        x=read();
        add(s,x+n,inf);
    }
    for(re int i=1;i<=n;i++)
    add(i,i+n,1);
    cout<<dinic();
}