b2OJ_1565_[NOI2009]植物大战僵尸_拓扑排序+最大权闭合子
题意:n*m个植物,每个植物有分数(可正可负),和能保护植物的位置。只能从右往左吃,并且不能吃正被保护着的,可以一个不吃,求获得的最大分数。
分析:把每个植物向能保护它的植物连边。源点连正权点,负权点连汇点。
考虑在一个环上的植物是吃不到的,我们可以用拓扑排序确定哪些是能吃的。
然后求一遍最大权闭合子图就是答案。
代码:
1 #include <stdio.h> 2 #include <string.h> 3 #include <algorithm> 4 #include <queue> 5 using namespace std; 6 #define S (n*m+1) 7 #define T (n*m+2) 8 #define P(x,y) (m*(x-1)+y) 9 #define inf 100000000 10 int head[700],to[1000010],nxt[1000010],flow[1000010],c[1000010],cnt=1,dep[700]; 11 int n,m,can[700],sum,val[700]; 12 void add(int u,int v,int f) 13 { 14 to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;flow[cnt]=f; 15 to[++cnt]=u;nxt[cnt]=head[v];head[v]=cnt;flow[cnt]=0; 16 } 17 bool bfs() 18 { 19 queue <int> q; 20 memset(dep,0,sizeof(dep)); 21 q.push(S);dep[S]=1; 22 while(q.size()) 23 { 24 int x=q.front();q.pop(); 25 for(int i=head[x];i;i=nxt[i]) 26 { 27 if(!dep[to[i]]&&flow[i]&&can[to[i]]) 28 { 29 dep[to[i]]=dep[x]+1; 30 q.push(to[i]); 31 } 32 } 33 } 34 return dep[T]; 35 } 36 int dfs(int x,int mf) 37 { 38 if(x==T)return mf; 39 int nf=0; 40 for(int i=head[x];i;i=nxt[i]) 41 { 42 if(dep[to[i]]==dep[x]+1&&flow[i]&&can[to[i]]) 43 { 44 int tmp=dfs(to[i],min(flow[i],mf-nf)); 45 nf+=tmp; 46 flow[i]-=tmp; 47 flow[i^1]+=tmp; 48 if(nf==mf)break; 49 } 50 } 51 dep[x]=0; 52 return nf; 53 } 54 void dinic() 55 { 56 int ans=0,f; 57 while(bfs()) 58 { 59 while(f=dfs(S,inf)) 60 { 61 ans+=f; 62 } 63 } 64 printf("%d",sum-ans); 65 } 66 void topsort() 67 { 68 queue <int> q; 69 for(int i=1;i<=T;i++) 70 { 71 if(!c[i])q.push(i); 72 } 73 while(q.size()) 74 { 75 int x=q.front();q.pop(); 76 can[x]=1; 77 if(val[x]>0)sum+=val[x]; 78 for(int i=head[x];i;i=nxt[i]) 79 { 80 c[to[i]]--; 81 if(c[to[i]]==0)q.push(to[i]); 82 } 83 } 84 } 85 int main() 86 { 87 scanf("%d%d",&n,&m); 88 int x,y,z,w; 89 for(int i=1;i<=n*m;i++) 90 { 91 scanf("%d%d",&x,&y); 92 val[i]=x; 93 if(x>0)add(S,i,x),c[S]++; 94 else add(i,T,-x),c[i]++; 95 if(i%m)add(i,i+1,inf),c[i]++; 96 while(y--) 97 { 98 scanf("%d%d",&z,&w); 99 add(P(z+1,w+1),i,inf); 100 c[P(z+1,w+1)]++; 101 } 102 } 103 topsort(); 104 dinic(); 105 } 106 107 /*************************************************************** 108 Problem: 1580 109 User: 20170105 110 Language: C++ 111 Result: Accepted 112 Time:320 ms 113 Memory:16704 kb 114 ****************************************************************/