bzoj 1565 [NOI2009]植物大战僵尸
拓扑排序+最小割。
每个植物有它保护的一些植物，等价于选一些点的先决条件是其它一些点，最大权闭合子图问题。
发现图里有环，所以环里的点都不能选，先决条件是环里的点的点也不能选，所以把所有边反向拓扑排序，拓扑不到的点不能选。
无视掉不能选的点后裸最大权闭合子图建图就好了。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<algorithm> 5 #include<queue> 6 using namespace std; 7 const int dian=605; 8 const int bian=1000005; 9 const int INF=0x3f3f3f3f; 10 int h[dian],nxt[bian],ver[bian],val[bian],ch[dian],cr[dian]; 11 int qz[dian],num[dian],dbh[dian][dian],in[dian],v[dian]; 12 int n,m,tot,aa,bb,cc,ans; 13 int S,T; 14 int bh(int a,int b){ 15 return (a-1)*m+b; 16 } 17 void add1(int a,int b){ 18 tot++;ver[tot]=b;nxt[tot]=h[a];h[a]=tot;in[b]++; 19 } 20 void toposort(){ 21 queue<int>q; 22 for(int i=1;i<=n*m;i++) 23 if(!in[i]) 24 q.push(i); 25 while(!q.empty()){ 26 int x=q.front(); 27 q.pop(); 28 v[x]=1; 29 for(int i=h[x];i;i=nxt[i]){ 30 int y=ver[i]; 31 in[y]--; 32 if(!in[y]) 33 q.push(y); 34 } 35 } 36 } 37 void add(int a,int b,int c){ 38 tot++;ver[tot]=b;val[tot]=c;nxt[tot]=h[a];h[a]=tot; 39 tot++;ver[tot]=a;val[tot]=0;nxt[tot]=h[b];h[b]=tot; 40 } 41 bool tell(){ 42 memset(ch,-1,sizeof(ch)); 43 queue<int>q; 44 q.push(S); 45 ch[S]=0; 46 while(!q.empty()){ 47 int t=q.front(); 48 q.pop(); 49 for(int i=h[t];i;i=nxt[i]) 50 if(ch[ver[i]]==-1&&val[i]){ 51 ch[ver[i]]=ch[t]+1; 52 q.push(ver[i]); 53 } 54 } 55 return ch[T]!=-1; 56 } 57 int zeng(int a,int b){ 58 if(a==T) 59 return b; 60 int r=0; 61 for(int i=cr[a];i&&b>r;i=nxt[i]) 62 if(ch[ver[i]]==ch[a]+1&&val[i]){ 63 int t=zeng(ver[i],min(b-r,val[i])); 64 val[i]-=t,r+=t,val[i^1]+=t; 65 if(val[i]) 66 cr[a]=i; 67 } 68 if(!r) 69 ch[a]=-1; 70 return r; 71 } 72 int dinic(){ 73 int r=0,t; 74 while(tell()){ 75 for(int i=1;i<=n*m+2;i++) 76 cr[i]=h[i]; 77 while(t=zeng(S,INF)) 78 r+=t; 79 } 80 return r; 81 } 82 int main(){ 83 memset(h,0,sizeof(h)); 84 memset(nxt,0,sizeof(nxt)); 85 memset(num,0,sizeof(num)); 86 memset(in,0,sizeof(in)); 87 memset(v,0,sizeof(v)); 88 tot=0,ans=0; 89 scanf("%d%d",&n,&m); 90 S=n*m+1,T=n*m+2; 91 for(int i=1;i<=n;i++) 92 for(int j=1;j<=m;j++){ 93 scanf("%d",&aa); 94 qz[bh(i,j)]=aa; 95 scanf("%d",&aa); 96 for(int k=1;k<=aa;k++){ 97 scanf("%d%d",&bb,&cc); 98 num[bh(i,j)]++; 99 dbh[bh(i,j)][num[bh(i,j)]]=bh(bb+1,cc+1); 100 add1(bh(i,j),bh(bb+1,cc+1)); 101 } 102 } 103 for(int i=1;i<=n;i++) 104 for(int j=2;j<=m;j++) 105 add1(bh(i,j),bh(i,j-1)); 106 toposort(); 107 memset(h,0,sizeof(h)); 108 memset(nxt,0,sizeof(nxt)); 109 tot=1; 110 for(int i=1;i<=n;i++) 111 for(int j=1;j<=m;j++) 112 if(v[bh(i,j)]){ 113 if(qz[bh(i,j)]>0){ 114 ans+=qz[bh(i,j)]; 115 add(S,bh(i,j),qz[bh(i,j)]); 116 } 117 else if(qz[bh(i,j)]<0) 118 add(bh(i,j),T,-qz[bh(i,j)]); 119 for(int k=1;k<=num[bh(i,j)];k++) 120 if(v[dbh[bh(i,j)][k]]) 121 add(dbh[bh(i,j)][k],bh(i,j),INF); 122 } 123 for(int i=1;i<=n;i++) 124 for(int j=1;j<m;j++) 125 if(v[bh(i,j)]&&v[bh(i,j+1)]) 126 add(bh(i,j),bh(i,j+1),INF); 127 printf("%d",ans-dinic()); 128 return 0; 129 }
