BZOJ1585 Earthquake Damage 2 地震伤害

题目链接：传送门

题目大意：自行参考原题

我承认这个题如果不告诉我是最小割我是想不到的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;}