[BZOJ1475]方格取数 网络流 最小割

1475: 方格取数

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1025  Solved: 512
[Submit][Status][Discuss]

Description

在一个n*n的方格里，每个格子里都有一个正整数。从中取出若干数，使得任意两个取出的数所在格子没有公共边，且取出的数的总和尽量大。

Input

第一行一个数n；（n<=30） 接下来n行每行n个数描述一个方阵

Output

仅一个数，即最大和

Sample Input

2
1 2
3 5

Sample Output

6

 

黑白染色求最小割

 #include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 int tx[]={,,,-};
 int ty[]={,,-,};
 struct data {
     int to,next,f;
 }e[];
 int head[],cnt;
 void add(int u,int v,int f){e[cnt].to=v;e[cnt].next=head[u];e[cnt].f=f;head[u]=cnt++;}
 int n,m,s,t;
 bool vis[];
 int q[],dis[];
 bool bfs() {
     memset(dis,-,sizeof(dis));
     int h=,tail=;
     q[h]=t;
     dis[t]=;
     while(h!=tail) {
         int now=q[h++];if(h==) h=;
         for(int i=head[now];i>=;i=e[i].next) {
             if(dis[e[i].to]>-||!e[i^].f) continue;
             dis[e[i].to]=dis[now]-;
             q[tail++]=e[i].to;if(tail==) tail=;
         }
     }
     return dis[s]>=-;
 }
 int dfs(int now,int a) {
     int f=,flow=;
     if(now==t) return a;
     for(int i=head[now];i>=;i=e[i].next) {
         int to=e[i].to;
         if(dis[to]==dis[now]+&&e[i].f>) {
             f=dfs(to,min(a,e[i].f));
             flow+=f;
             e[i].f-=f;
             e[i^].f+=f;
             a-=f;
             if(a==) break;
         }
     }
     return flow;
 }
 int num=;
 int ans=;
 int sum=;
 void dinic() {
     while(bfs()) {
         sum+=dfs(s,);
     }
     printf("%d",num-sum);
 }
 int tt=;
 int a[][];
 int main() {
     memset(head,-,sizeof(head));
     scanf("%d%d",&m,&n);
     s=,t=;
     for(int i=;i<=m;i++)
         for(int j=;j<=n;j++) {
             scanf("%d",&a[i][j]);
             num+=a[i][j];
         }
     for(int i=;i<=m;i++) {
         for(int j=;j<=n;j++) {
             if((i+j)&) {
                 add(s,(i-)*n+j,a[i][j]),add((i-)*n+j,s,);
                 for(int k=;k<=;k++) {
                     int tox=i+tx[k],toy=j+ty[k];
                     if(tox==||tox>m||toy==||toy>n) continue;
                     add((i-)*n+j,(tox-)*n+toy,1e8);
                     add((tox-)*n+toy,(i-)*n+j,);
                 }
             }
             else add((i-)*n+j,t,a[i][j]),add(t,(i-)*n+j,);

         }
     }
     dinic();
 }