【bzoj1565】[NOI2009]植物大戰僵尸 【網絡流】【最大權閉合子圖】

題解：可以看出每個格子有一些前驅，只有前驅都被消滅了才能走到這里。因為要求最大的權值和，所以我們可以用最大權閉合子圖來求解這題。最大權閉合子圖點這里！
然后讓蒟蒻講一講自己掉的坑。
首先，根據WYC大佬的博客，我們要先進行一次拓撲排序來把因出現了環而無敵的格子排除出去。
然后我就掉了一個坑：拓撲排序建邊應該是與網絡流建邊反向的，因為按照最大權閉合子圖的建法，環上連出去的邊連到的是保護環上格子的，這些格子是可以被收集的，但是如果按照這種建法，拓撲排序時永遠都不會拓展到這些格子。所以要反向建邊。
我真得有點懷疑人生，我的網絡流模板是不是有大問題，常數大到上天！
UPD：網絡流模板真的有問題，work數組簡直就是害人的！去掉！
代碼：

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int inf=0x7fffffff;
int n,m,w,x,y,s,t,ans,sc[25][35],id[25][35];
int cnt,head[605],work[605],dep[605],in[605],to[400005],nxt[400005],dd[400005];
bool ck[605];
vector<int> v[605];
queue<int> q;
void addedge(int u,int v){in[u]++;to[++cnt]=u;nxt[cnt]=head[v];head[v]=cnt;
}
void adde(int u,int v,int d){to[++cnt]=v;nxt[cnt]=head[u];dd[cnt]=d;head[u]=cnt;to[++cnt]=u;nxt[cnt]=head[v];dd[cnt]=0;head[v]=cnt;
}
bool bfs(){memset(dep,0,sizeof(dep));dep[s]=1;while(!q.empty()){q.pop();}q.push(s);while(!q.empty()){int u=q.front(),v;q.pop();for(int i=head[u];i;i=nxt[i]){v=to[i];if(dd[i]&&!dep[v]){dep[v]=dep[u]+1;if(v==t){return true;}q.push(v);}}}return false;
}
int dfs(int u,int f){if(u==t){return f;}int tmp,res=0,v;for(int &i=work[u];i&&f;i=nxt[i]){v=to[i];if(dd[i]&&dep[v]==dep[u]+1&&(tmp=dfs(v,min(f,dd[i])))){dd[i]-=tmp;dd[i^1]+=tmp;f-=tmp;res+=tmp;}}if(!res){dep[u]=0;}return res;
}
int maxflow(){int res=0;while(bfs()){memcpy(work,head,sizeof(head));res+=dfs(s,0x7fffffff);}return res;
}
int main(){scanf("%d%d",&n,&m);for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){id[i][j]=(i-1)*m+j;}}for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){scanf("%d%d",&sc[i][j],&w);while(w--){scanf("%d%d",&x,&y);x++,y++;v[id[x][y]].push_back(id[i][j]);}}}for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){for(int k=0;k<(int)v[id[i][j]].size();k++){addedge(id[i][j],v[id[i][j]][k]);}if(j<m){addedge(id[i][j],id[i][j+1]);}}}for(int i=1;i<=n*m;i++){if(!in[i]){q.push(i);}}while(!q.empty()){int u=q.front(),v;q.pop();ck[u]=true;for(int i=head[u];i;i=nxt[i]){v=to[i];in[v]--;if(!in[v]){q.push(v);}}}s=0;t=n*m+1;cnt=1;memset(head,0,sizeof(head));for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){if(ck[id[i][j]]){if(sc[i][j]>0){ans+=sc[i][j];adde(s,id[i][j],sc[i][j]);}else if(sc[i][j]<0){adde(id[i][j],t,-sc[i][j]);}for(int k=0;k<(int)v[id[i][j]].size();k++){adde(id[i][j],v[id[i][j]][k],inf);}if(j<m){adde(id[i][j],id[i][j+1],inf);}}}}printf("%d\n",ans-maxflow());return 0;
}