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

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

25

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 。

我们先消环。然后就是最大权闭合子图了。

我的程序最后一个点TLE了= =谁能帮我看看哪里错了么。。

#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
int head[100001],exhead[100001];
struct map
{
	 int f;
	 int s,t;
     int next;
}a[800001],exa[800001];
int edge,exedge;
int p;
int q[400001],d[400001];
int sum;
inline void exadd(int s,int t)
{
     exa[exedge].next=exhead[s];
     exhead[s]=exedge;
     exa[exedge].s=s;
     exa[exedge].t=t;
}
inline void add(int s,int t,int f)
{
     a[edge].next=head[s];
     head[s]=edge;
     a[edge].s=s;
     a[edge].t=t;
     a[edge].f=f;
}
inline bool bfs()
{
     int l=0,r=0;
     memset(q,0,sizeof(q));
     r++;
     q[r]=0;
     memset(d,-1,sizeof(d));
     d[0]=0;
     int i,k;
     while(l<r)
     {
     	  l++;
          int k=q[l];
          for(i=head[k];i!=0;i=a[i].next)
          {
               if(a[i].f>0&&d[a[i].t]==-1)
               {
                    d[a[i].t]=d[k]+1;
                    r++;
                    q[r]=a[i].t;
               }
          }
     }
     if(d[p]>=0)
          return true;
     return false;
}
inline int dfs(int k,int s)
{
     if(k==p)
          return s;
     int t=s;
     int i;
     for(i=head[k];i!=0;i=a[i].next)
     {
          if(d[a[i].t]==d[k]+1&&a[i].f>0)
          {
               int xx=dfs(a[i].t,min(s,a[i].f));
               a[i].f-=xx;
               if(i%2==0)
                    a[i-1].f+=xx;
               else
                    a[i+1].f+=xx;
               s-=xx;
          }
     }
     return t-s;
}
inline int maxflow()
{
     int s=0;
     while(bfs())
     {
          s+=dfs(0,2100000000);
          if(s>sum)
          {
               s=sum;
               break;
          }
     }
     return s;
}
struct npc
{
     int s,a;
}tx[100001];
int indeg[100001];
bool v[100001];
int main()
{
	 freopen("pvz.in","r",stdin);
	 freopen("pvz.out","w",stdout);
     int n,m;
     scanf("%d%d",&n,&m);
     int i,j,x,y;
     int px=0;
     for(i=1;i<=n*m;i++)
     {
          scanf("%d%d",&tx[i].s,&tx[i].a);
          px=px%m+1;
          if(px!=1)
          {
               exedge++;
               exadd(i,i-1);
               indeg[i-1]++;
          }
          for(j=1;j<=tx[i].a;j++)
          {
          	   scanf("%d%d",&x,&y);
          	   x++;
          	   y++;
               exedge++;
               exadd(i,(x-1)*m+y);
               indeg[(x-1)*m+y]++;
          }
     }
     queue <int>Q;
     while(!Q.empty())
          Q.pop();
     for(i=1;i<=n*m;i++)
          if(indeg[i]==0)
               Q.push(i);
     memset(v,false,sizeof(v));
     while(!Q.empty())
     {
          int d=Q.front();
          Q.pop();
          v[d]=true;
          for(i=exhead[d];i!=0;i=exa[i].next)
          {
               int t=exa[i].t;
               indeg[t]--;
               if(indeg[t]==0)
                    Q.push(t);
          }
     }
     for(i=1;i<=n*m;i++)
     {
          for(j=exhead[i];j!=0;j=exa[j].next)
          {
               int t=exa[j].t;
               if(v[t])
               {
                    edge++;
                    add(t,i,2100000000);
                    edge++;
                    add(i,t,0);
               }
          }
     }
     p=n*m+1;
     for(i=1;i<=n*m;i++)
     {
          if(v[i])
          {
               if(tx[i].s>0)
               {
                    edge++;
                    add(0,i,tx[i].s);
                    edge++;
                    add(i,0,0);
                    sum+=tx[i].s;
               }
               else if(tx[i].s<0)
               {
                    edge++;
                    add(i,p,-tx[i].s);
                    edge++;
                    add(p,i,0);
               }
          }
     }
     printf("%d\n",sum-maxflow());
     return 0;
}