題解:可以看出每個格子有一些前驅,只有前驅都被消滅了才能走到這里。因為要求最大的權值和,所以我們可以用最大權閉合子圖來求解這題。最大權閉合子圖點這里!
然后讓蒟蒻講一講自己掉的坑。
首先,根據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;
}
轉載于:https://www.cnblogs.com/2016gdgzoi471/p/9476860.html
總結
以上是生活随笔為你收集整理的【bzoj1565】[NOI2009]植物大战僵尸 【网络流】【最大权闭合子图】的全部內容,希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。