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 0x7fffffffint 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;}