poj2987 Firing[最小割]

题目

---

求选最少点个数的最大权闭合子图。（板子题）

---

最小割入门题，什么都不想说，丢个别人题解地址就跑。

附加几点个人理解：与s相通的S点集是闭合子图，剩下的与t相通的T点集是其他的。任意一个割都保证了有一个闭合子图是S集。

然后是选最少点个数，考虑加入最大权时存在多个方案，选的点个数不同，既然各方案都是闭合的，那我全选也是闭合的，就会创出更多权值，与之前最大权矛盾，所以最大权一定点数唯一，选法唯一。那就直接从s点找残量网络，把他的点集都找一下，计下数即可。←这里有个原理，做最小割（最大流），而最大流算好后因为最小割必然满流，所以最小割对应的边残量为0，则dfs时就不遍历残量为0的边。←刚才又提到一个结论，一定存在一个最大流使得最小割在其中满流，这里口胡证明一下，防止以后我忘了：最大流做完后s所能到的所有点S集，以及另一T集之间，应当是有一些边满流的也就是残量为0的，不然就可继续增广，就不是最大流了，那这些边对应了一个割CUT[S,T]，再由任意流小于等于任意割结论（打个比方，水流将流到割断的地方，而水流量小于等于容量，故得证）知，任意一个割大于等于这个最大流，即大于等于这个割，那这个割就是最小割，而她满流。证毕。

UPD:我现在看发现当时好像证的有点问题。。所以请不要信以为真，有空回来补详细的一些证明，当时只是个网络流初心者。

我还是太菜了啊。

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<queue>
 using namespace std;
 typedef long long ll;
 typedef pair<int,int> pii;
 template<typename T>inline char MIN(T&A,T B){return A>B?A=B,:;}
 template<typename T>inline char MAX(T&A,T B){return A<B?A=B,:;}
 template<typename T>inline T _min(T A,T B){return A<B?A:B;}
 template<typename T>inline T _max(T A,T B){return A>B?A:B;}
 template<typename T>inline void inc(T&A,T B){A+=B;}
 template<typename T>inline T read(T&x){
     x=;int f=;char c;while(!isdigit(c=getchar()))if(c=='-')f=;
     while(isdigit(c))x=x*10+c-'0',c=getchar();return f?x=-x:x;
 }
 const int N=+,M=+;ll INF=1ll<<;
 int Head[N],cur[N],Next[M<<],v[M<<],dis[N],vis[N],tot=,n,cnt,m,s,t;ll w[M<<];

 inline void Addedge(int x,int y,ll z){
     v[++tot]=y,Next[tot]=Head[x],Head[x]=tot,w[tot]=z;
     v[++tot]=x,Next[tot]=Head[y],Head[y]=tot,w[tot]=;
 }
 #define y v[j]
 inline int bfs(){
     queue<int> q;memset(dis,,sizeof dis),dis[s]=,q.push(s);
     for(register int i=;i<=n+;++i)cur[i]=Head[i];
     while(!q.empty()){
         int x=q.front();q.pop();
         for(register int j=Head[x];j;j=Next[j])if(w[j]&&!dis[y]){
             dis[y]=dis[x]+,q.push(y);
             if(y==t)return ;
         }
     }
     return ;
 }
 ll dinic(int x,ll flow){
     if(!flow||x==t)return flow;
     ll rest=flow,k;
     for(register int j=cur[x];j&&rest;cur[x]=j,j=Next[j])if(w[j]&&dis[y]==dis[x]+){
         if(!(k=dinic(y,_min(rest,w[j]))))dis[y]=;
         rest-=k,w[j]-=k,w[j^]+=k;
     }
     return flow-rest;
 }
 void dfs(int x){
     for(register int j=Head[x];j;j=Next[j])if(w[j]&&!vis[y])vis[y]=,++cnt,dfs(y);
 }
 #undef y
 int x,y;
 ll maxflow,ans,z;

 int main(){//freopen("tmp.in","r",stdin);freopen("tmp.out","w",stdout);
     read(n),read(m);s=n+,t=s+;
     for(register int i=;i<=n;++i)read(z)>?(inc(ans,z),Addedge(s,i,z)):Addedge(i,t,-z);
     for(register int i=;i<=m;++i)read(x),read(y),Addedge(x,y,INF);
     while(bfs())inc(maxflow,dinic(s,INF));
     ans-=maxflow;cnt=,vis[s]=,dfs(s);
     printf("%d %lld\n",--cnt,ans);
     return ;
 }

Upd:今天被巨佬嘲讽了说我快读用位运算太慢。orzorzorzTQL