题目链接:传送门
题目大意:自行参考原题
我承认这个题如果不告诉我是最小割我是想不到的QwQ
这个题如果知道了是最小割就非常显然了:
S=1,T=2n+1,对每个点x拆成x和x2=x+n,
对于图上的每一条边(u,v),建边(u2,v,INF),(v2,u,INF)。
对于报告的点u,建边(u,u2,INF),(u2,T,INF)。
对于其它不是1的点u,建边(u,u2,1)
对于点1,建边(1,n+1,INF)。
最小割就是答案。
代码:
//BZOJ1585
#include<queue>
#include<cstdio>
#include<climits>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 7000
#define MAXM 200000
#define INF INT_MAX
#define build_edge(u,v,c) add_edge(u,v,c),add_edge(v,u,0)
#define P(x,y) (~-(x)*m+(y))
#define inrange(x,y) (1<=x&&x<=n&&1<=y&&y<=m)
using namespace std;
struct edges{
int to,pre,resf;
}e[MAXM];
int etop,h[MAXN],lev[MAXN],cur[MAXN];
bool vis[MAXN];queue<int> q;
int dx[]={0,-1,0,1,0},
dy[]={0,0,1,0,-1};
int add_edge(int u,int v,int c)
{
etop++;
e[etop].to=v;
e[etop].pre=h[u];
h[u]=etop;
e[etop].resf=c;
return 0;
}
bool bfs(int s,int t)
{
memset(vis,false,sizeof(vis));
memset(lev,0,sizeof(lev));
while(!q.empty()) q.pop();
q.push(s);lev[s]=0;vis[s]=true;
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=h[x];i;i=e[i].pre)
if(!vis[e[i].to]&&e[i].resf)
{
vis[e[i].to]=true;
lev[e[i].to]=lev[x]+1;
q.push(e[i].to);
}
}
return vis[t];
}
int dfs(int s,int t,int a)
{
if(s==t||!a) return a;
int flow=0,f;
for(int &i=cur[s];i;i=e[i].pre)
if(lev[e[i].to]==lev[s]+1&&(f=dfs(e[i].to,t,min(a,e[i].resf)))>0)
{
a-=f;flow+=f;e[i].resf-=f;
e[((i-1)^1)+1].resf+=f;
if(!a) break;
}
return flow;
}
bool rep[MAXN];
int main()
{
int n,m,p;scanf("%d%d%d",&n,&m,&p);
int s=1,t=2*n+1;
for(int i=1;i<=m;i++)
{
int x,y;scanf("%d%d",&x,&y);
if(x==y) continue;
int x2=x+n,y2=y+n;
build_edge(x2,y,INF);
build_edge(y2,x,INF);
}
while(p--)
{
int x,x2;scanf("%d",&x);x2=x+n;
build_edge(x,x2,INF);
build_edge(x2,t,INF);
rep[x]=true;
}
build_edge(1,n+1,INF);
for(int i=2;i<=n;i++)
if(!rep[i]) build_edge(i,i+n,1);
int flow=0;
while(bfs(s,t))
{
for(int i=s;i<=t;i++)
cur[i]=h[i];
flow+=dfs(s,t,INF);
}
printf("%d\n",flow);
return 0;
}