BZOJ1475方格取数

1475: 方格取数
Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 710 Solved: 361
Description
在一个n*n的方格里，每个格子里都有一个正整数。从中取出若干数，使得任意两个取出的数所在格子没有公共边，且取出的数的总和尽量大。
Input
第一行一个数n；（n<=30） 接下来n行每行n个数描述一个方阵
Output
仅一个数，即最大和
Sample Input
2
1 2
3 5
Sample Output
6
如果x+y为奇数，(x,y)->T，容量v[x,y]
如果x+y为偶数，S->(x,y)，容量v[x,y]
并向四个方向的点连边，容量inf
答案为总收益减去最小割
附上本蒟蒻的代码：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define inf 0x7fffffff
int n,xx[4]={0,0,1,-1},yy[4]={1,-1,0,0},map[31][31],mark[31][31],h[10001],dis[10001],q[10001],head,tail,cnt=1,ans=0,sum;
struct kx
{
    int to,next,v;
}edge[10001];

bool pd(int x,int y)
{
    if (x<1 || y<1 || x>n || y>n) return false;
    else return true;
}

int read()
{
    int w=0,c=1; char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

void add(int u,int v,int w)
{
    cnt++,edge[cnt].next=h[u],h[u]=cnt,edge[cnt].to=v,edge[cnt].v=w;
    cnt++,edge[cnt].next=h[v],h[v]=cnt,edge[cnt].to=u,edge[cnt].v=0;
}

bool bfs()
{
    int j,p;
    memset(dis,-1,sizeof(dis));
    q[0]=0;
    dis[0]=0;
    head=0;
    tail=1;
    while (head<tail)
      {
        head++;
        j=q[head];
        p=h[j];
        while (p)
          {
            if (dis[edge[p].to]<0 && edge[p].v>0)
              {
                  dis[edge[p].to]=dis[j]+1;
                  tail++;
                  q[tail]=edge[p].to;
              }
            p=edge[p].next;
          }
      }
    if (dis[n*n+1]>0)
      return true;
    else
      return false;
}

int dfs(int x,int f)
{
    int w,used=0,i=h[x];
    if (x==n*n+1)
      return f;
    while (i)
      {
        if (edge[i].v && dis[edge[i].to]==dis[x]+1)
          {
            w=f-used;
            w=dfs(edge[i].to,min(w,edge[i].v));
            edge[i].v-=w;
            edge[i^1].v+=w;
            used+=w;
            if (used==f)
              return f;
          }
        i=edge[i].next;
      }
    if (!used)
      dis[x]=-1;
    return used;
}

int main()
{
    int i,j,b=0,w,k;
    n=read();
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
        map[i][j]=read(),ans+=map[i][j];
    w=(n*n+1)/2;
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
        if ((i+j)%2==0)
          b++,mark[i][j]=b;
        else
          w++,mark[i][j]=w;
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
        if ((i+j)%2==0)
          {
            add(0,mark[i][j],map[i][j]);
            for (k=0;k<4;k++)
              if (pd(i+xx[k],j+yy[k]))
                add(mark[i][j],mark[i+xx[k]][j+yy[k]],inf);
          }
        else
          add(mark[i][j],n*n+1,map[i][j]);
    while (bfs())
      while (sum=dfs(0,inf))
        ans-=sum;
    printf("%d",ans);
    return 0;
}