bzoj 1565 [NOI2009]植物大战僵尸解题报告

1565: [NOI2009]植物大战僵尸

Time Limit: 10 Sec Memory Limit: 64 MB
Submit: 2161 Solved: 1000
[Submit][Status][Discuss]

Description

Input

Output

仅包含一个整数，表示可以获得的最大能源收入。注意，你也可以选择不进行任何攻击，这样能源收入为0。

Sample Input

3 2
10 0
20 0
-10 0
-5 1 0 0
100 1 2 1
100 0

Sample Output

HINT

在样例中, 植物P1,1可以攻击位置(0,0), P2, 0可以攻击位置(2,1)。
一个方案为，首先进攻P1,1, P0,1，此时可以攻击P0,0 。共得到能源收益为(-5)+20+10 = 25。注意, 位置(2,1)被植物P2,0保护，所以无法攻击第2行中的任何植物。
【大致数据规模】
约20%的数据满足1 ≤ N, M ≤ 5；
约40%的数据满足1 ≤ N, M ≤ 10；
约100%的数据满足1 ≤ N ≤ 20，1 ≤ M ≤ 30，-10000 ≤ Score ≤ 10000 。

Source

[Submit][Status][Discuss]

HOME Back

————————————————我是分割线——————————————————————————

好题啊好题

这题折腾了三天，差点没吐出血来。

先拓扑排序求环，把环删掉。

之后重建图求最大权闭合子图，把一个植物向保护它的植物连边，跑一遍网络流即可。

第一次只过了4个点。不明所以。

换方法重写了拓扑，在网络流中加限制乱搞。

终于90，但第九个点一直RE，最后发现数组开小了。这题最多600个点，但边数组一直开到800005才A了。TAT~

 /*

     Problem: bzoj 1565 [NOI2009]植物大战僵尸

     OJ:         bzoj

     User:    S.B.S.

     Time:     264 MS

     Memory:     13884 KB

     Length:  6934 B

 */

 #include<iostream>

 #include<cstdio>

 #include<cstring>

 #include<cmath>

 #include<algorithm>

 #include<queue>

 #include<cstdlib>

 #include<iomanip>

 #include<cassert>

 #include<climits>

 #include<functional>

 #include<bitset>

 #include<vector>

 #include<list>

 #define F(i,j,k) for(int i=j;i<=k;++i)

 #define M(a,b) memset(a,b,sizeof(a))

 #define FF(i,j,k) for(int i=j;i>=k;i--)

 #define inf 0x3f3f3f3f

 #define maxm 1001

 #define mod 998244353

 //#define LOCAL

 using namespace std;

 int read(){

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}

     return x*f;

 }

 int n,m;

 int mp[];

 struct EDGE

 {

     int from;

     int to;

     int value;

     int next;

 }e[];

 int head[];

 int v[],in[],out[];

 int tot=;

 inline void addedge(int u,int v,int w)

 {

     e[++tot].from=u;

     e[tot].to=v;

     e[tot].value=w;

     e[tot].next=head[u];

     head[u]=tot;

 //    tot++;

 }

 //bool kan[2000][2000];

 int ans=;

 //inline int at(int x,int y)

 //{

 //    return m*x+y+1;

 //}

 //int stk[4001],top=0;

 //inline void topsort()

 //{

 //    int temp[4001];M(temp,0);

 //    F(i,1,m*n) temp[i]=in[i];

 //    F(i,1,m*n){

 //        if(in[i]==0){

 //            stk[top]=i;

 ////            cout<<i<<"********"<<endl;

 //            top++;

 //        }

 //    }

 ////    cout<<top<<": "<<endl;

 ////    F(i,0,top-1) cout<<stk[i]<<" ";cout<<endl;

 //    while(top!=0){

 //        top--;int u=stk[top];

 //        for(int i=head[u];i!=-1;i=e[i].next)

 //        {

 //            in[e[i].to]--;

 //            if(in[e[i].to]==0){

 //                stk[top]=e[i].to;

 //                top++;

 //            }

 //        }

 //    }

 //    F(i,1,m*n){

 //        if(in[i]!=0){

 //            stk[top]=i;

 //            top++;

 //        }

 //    }

 //    while(top!=0){

 //        top--;int u=stk[top];

 //        temp[u]=-1;out[u]=-1;

 //        v[u]=0;

 //    }

 //    F(i,1,m*n) in[i]=temp[i];

 //    return;

 //}

 inline void add(int u,int v,int w)

 {

     addedge(u,v,w);

     addedge(v,u,);

     ++in[u];

     return;

 }

 int map[];

 int S,T;

 bool ok()

 {

     F(i,S,T) if(map[i]!=-) map[i]=-;

     int que[];

     int hd=,tl=;

     que[tl++]=S;

     map[S]=;

     while (hd<tl){

         int u=que[hd++];

         for(int i=head[u];i;i=e[i].next)

         {

             if(map[e[i].to]==-&&e[i].value)

             {

                 map[e[i].to]=map[u]+;

                 que[tl++]=e[i].to;

             }

         }

     }

 //    if (map[T]!=-1) return true;

 //    else return false;

     return map[T]>;

 }

 int cur[];

 int zeng(int k,int now)

 {

     if (k==T) return now;

     int r=;

     for (int i=cur[k];i&&now>r;i=e[i].next)

     {

         if (map[k]+==map[e[i].to]&&e[i].value)

         {

             int t=zeng(e[i].to,min(now-r,e[i].value));

             e[i].value-=t;e[i^].value+=t;r+=t;

             if(e[i].value) cur[k]=i;//***

 //            cout<<e[i].from<<" ---> "<<e[i].to<<" $ "<<t<<endl;

         }

     }

     if (!r) map[k]=-;

     return r;

 }

 int dinic()

 {

     int r=,t;

     M(cur,);

     while (ok()){

         for(int i=S;i<=T;++i) cur[i]=head[i];

         r+=zeng(S,inf);

 //        while(t=zeng(S,inf)) r+=t;

     }

     return r;

 }

 //pair<int,int> kk[4001];

 //int cur;

 int main()

 {

     std::ios::sync_with_stdio(false);//cout<<setiosflags(ios::fixed)<<setprecision(1)<<y;

     #ifdef LOCAL

     freopen("111.in","r",stdin);

     freopen("data.out","w",stdout);

     #endif

     cin>>n>>m;

     S=;T=m*n+;

     M(mp,);M(head,);

 //    F(i,0,n-1)F(j,0,m-1){

 //        int x,y;

 //        cin>>x>>y;

 //        mp[i][j]=x;

 //        v[at(i,j)]=x;

 //        in[at(i,j)]++;

 //        if(j==m-1) in[at(i,j)]--;

 //        if(j!=0) add(at(i,j),at(i,j-1),v[at(i-1,j)]);

 //        if(y==0) continue;

 ////        **********************************************

 //        F(k,1,y){

 //            int a,b;

 //            cin>>a>>b;

 //            kk[cur].first=at(i,j);

 //            kk[cur].second=at(a,b);

 //            cur++;

 //            in[at(a,b)]++;

 //            addedge(at(i,j),at(a,b),0);

 ////            cout<<at(i,j)<<" "<<at(a,b)<<" &&&&&"<<endl;

 //        }

 //    }

 //    F(i,1,m*n) cout<<in[i]<<" ";cout<<endl;

 //    cout<<"/*****************************"<<endl;

 //    F(i,1,tot) cout<<e[i].from<<" ---> "<<e[i].to<<" $ "<<e[i].value<<endl;

 //    cout<<"*****************************/"<<endl;

     F(i,,n*m){

         int temp;

         cin>>mp[i];

         if(mp[i]>) add(S,i,mp[i]);

         else add(i,T,-mp[i]);

         cin>>temp;

         while(temp--){

             int x,y;

             cin>>x>>y;

             add(x*m+y+,i,inf);

         }

         if(i%m) add(i,i+,inf);

     }

 //    topsort();

 //    cout<<"))))))"<<endl;

 //    F(i,1,m*n) if(in[i]==-1) cout<<i<<" ";cout<<"((((("<<endl;

 //    M(e,0);M(head,-1);tot=0;

     int q[];M(q,);

     int cur1=,cur2=;

     for(int i=S;i<=T;++i){

         if(!in[i]) q[cur2++]=i;

         map[i]=-;

     }

     int sum=;

     while(cur1<cur2){

         int cnt=q[cur1++];map[cnt]=;

         if(mp[cnt]>) sum+=mp[cnt];

         for(int i=head[cnt];i;i=e[i].next)

         {

             if(i&) if(!--in[e[i].to]) q[cur2++]=e[i].to;

         }

     }

 //    F(i,1,m*n) cout<<in[i]<<" ";cout<<endl;

 //    cout<<"v: "<<endl;

 //    F(i,1,m*n) cout<<v[i]<<" ";cout<<endl;

 //    F(i,0,n-1)F(j,0,m-1){

 //        if(in[at(i,j)]==-1) continue;

 ////        cout<<at(i,j)<<" : "<<in[at(i,j)]<<"@ "<<endl;

 //        if(j!=0&&in[at(i,j-1)]!=-1) add(at(i,j-1),at(i,j),inf);

 //        if(v[at(i,j)]==0) continue;

 //        if(v[at(i,j)]>0) add(31,at(i,j),v[at(i,j)]),sum+=v[at(i,j)];

 //        else add(at(i,j),32,-v[at(i,j)]);

 //    }

 //    while(cur!=0){

 //        cur--;

 //        if(in[kk[cur].first]==-1||in[kk[cur].second]==-1) continue;

 //        int u=kk[cur].first;

 //        int v=kk[cur].second;

 //        add(v,u,inf);

 //    }

 //    cout<<"/*****************************"<<endl;

 //    F(i,1,tot) cout<<e[i].from<<" ---> "<<e[i].to<<" $ "<<e[i].value<<endl;

 //    cout<<"*****************************/"<<endl;

 //    cout<<sum<<endl;

     cout<<sum-dinic()<<endl;

     return ;

 }

 /*

 3 2

 10 0

 20 0

 -10 0

 -5 1 0 0

 100 1 2 1

 100 0

 */

植物大战僵尸