bzoj 2132: 圈地计划

2132: 圈地计划

Time Limit: 2 Sec  Memory Limit: 256 MB
Submit: 918  Solved: 417
[Submit][Status][Discuss]

Description

最近房地产商GDOI(Group of Dumbbells Or Idiots)从NOI(Nuts Old Idiots)手中得到了一块开发土地。据了解，这块土地是一块矩形的区域，可以纵横划分为N×M块小区域。GDOI要求将这些区域分为商业区和工业区来开发。根据不同的地形环境，每块小区域建造商业区和工业区能取得不同的经济价值。更具体点，对于第i行第j列的区域，建造商业区将得到Aij收益，建造工业区将得到Bij收益。另外不同的区域连在一起可以得到额外的收益，即如果区域(I,j)相邻（相邻是指两个格子有公共边）有K块（显然K不超过4）类型不同于(I,j)的区域，则这块区域能增加k×Cij收益。经过Tiger.S教授的勘察，收益矩阵A,B,C都已经知道了。你能帮GDOI求出一个收益最大的方案么？

Input

输入第一行为两个整数，分别为正整数N和M，分别表示区域的行数和列数；第2到N+1列，每行M个整数，表示商业区收益矩阵A；第N+2到2N+1列，每行M个整数，表示工业区收益矩阵B；第2N+2到3N+1行，每行M个整数，表示相邻额外收益矩阵C。第一行，两个整数，分别是n和m（1≤n，m≤100）；

任何数字不超过1000”的限制

Output

输出只有一行，包含一个整数，为最大收益值。

Sample Input

3 3
1 2 3
4 5 6
7 8 9
9 8 7
6 5 4
3 2 1
1 1 1
1 3 1
1 1 1

Sample Output

81
【数据规模】
对于100%的数据有N,M≤100

HINT

 数据已加强，并重测--2015.5.15

在下面对拍，最大数据都不会错，交上去很快就WA了，不知道为什么，很奇怪。

#include<iostream>    #include<cstdio>    #include<cstring>    #include<algorithm>    #include<queue>    #include<cmath>    #define N 120    #define M 1000003    using namespace std;    int n,m,tot=-1,sum;    int next[M],point[M],remain[M],v[M];    int a[N][N],b[N][N],c[N][N],p[N][N];    int deep[M],cur[M];  int x[10]={0,1,0,-1},y[10]={1,0,-1,0};  const int inf=1e9;  void add(int x,int y,int z)    {      tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;      tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0;     }   void addedge(int x,int y,int z){  tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;  tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=z;} bool bfs(int s,int t)    {      memset(deep,0x7f,sizeof(deep));      for (int i=s;i<=t;i++) cur[i]=point[i];      deep[s]=0;       queue<int> p; p.push(s);      while(!p.empty())      {       int now=p.front(); p.pop();       for (int i=point[now];i!=-1;i=next[i])       if (deep[v[i]]>inf&&remain[i])        deep[v[i]]=deep[now]+1,p.push(v[i]);      }      if (deep[t]>inf) return false;      return true;    }    int dfs(int now,int t,int limit)  {    if (now==t||!limit) return limit;    int flow=0,f;    for(int i=cur[now];i!=-1;i=next[i])    {      cur[now]=i;      if (deep[v[i]]==deep[now]+1&&(f=dfs(v[i],t,min(limit,remain[i]))))      {        flow+=f;        limit-=f;      remain[i]-=f;        remain[i^1]+=f;        if (!limit) break;        }     }     return flow;  }  int dinic(int s,int t)    {      int ans=0;      while (bfs(s,t))       ans+=dfs(s,t,inf);      return ans;    }  int main()    {    freopen("a.in","r",stdin); freopen("my.out","w",stdout);  memset(point,-1,sizeof(point)); memset(next,-1,sizeof(next)); sum=0; scanf("%d%d",&n,&m);     for (int i=1;i<=n;i++)      for (int j=1;j<=m;j++) scanf("%d",&a[i][j]),sum+=a[i][j];    for (int i=1;i<=n;i++)      for (int j=1;j<=m;j++) scanf("%d",&b[i][j]),sum+=b[i][j];     for (int i=1;i<=n;i++)      for (int j=1;j<=m;j++)   {   scanf("%d",&c[i][j]);   if (i==1&&j==1||i==1&&j==m||i==n&&j==1||i==n&&j==m)    sum+=c[i][j]*2;   else    if (i==1||j==1||i==n||j==m)   sum+=c[i][j]*3;    else sum+=c[i][j]*4;   } for (int i=1;i<=n;i++)//黑白染色   for (int j=1;j<=m;j++)  {   if (i%2==1&&j%2==1)  p[i][j]=1;   if (i%2==0&&j%2==0)  p[i][j]=1;  } /*for (int i=1;i<=n;i++)  {    for (int j=1;j<=m;j++)     cout<<p[i][j]<<" ";    cout<<endl;  }*/ int num=n*m;     for (int i=1;i<=n;i++)      for (int j=1;j<=m;j++)      if (p[i][j])  add(0,(i-1)*m+j,a[i][j]),add((i-1)*m+j,num+1,b[i][j]);   else add(0,(i-1)*m+j,b[i][j]),add((i-1)*m+j,num+1,a[i][j]); for (int i=1;i<=n;i++)      for (int j=1;j<=m;j++)     if (p[i][j])     {        int x=(i-1)*m+j;        if (i-1>0)   addedge(x,(i-2)*m+j,c[i][j]+c[i-1][j]);       if (i+1<=n)   addedge(x,i*m+j,c[i][j]+c[i+1][j]);        if (j-1>0)   addedge(x,x-1,c[i][j]+c[i][j-1]);        if (j+1<=m)   addedge(x,x+1,c[i][j]+c[i][j+1]);      }      int t=dinic(0,num+1);   // cout<<t<<endl;   //cout<<sum<<endl;   printf("%d\n",sum-t);}