南邮 OJ 1356 方格取数问题

方格取数问题

时间限制(普通/Java) :  1000 MS/ 3000 MS          运行内存限制 : 65536 KByte
总提交 : 43            测试通过 : 18 

比赛描述

 

在一个有m*n 个方格的棋盘中，每个方格中有一个正整数。现要从方格中取数，使任意2 个数所在方格没有公共边，且取出的数的总和最大。试设计一个满足要求的取数算法。

    对于给定的方格棋盘，按照取数要求编程找出总和最大的数。

 

输入

        第1行有2个正整数m和n，分别表示棋盘的行数和列数。接下来的m行，每行有n个正整数，表示棋盘方格中的数

输出

    程序运行结束时，将取数的最大总和输出。

样例输入

3 3
1 2 3
3 2 3
2 3 1

样例输出

11

提示

 

题目来源

算法设计与实验题解

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>

using namespace std;

const int VM=2520;
const int EM=5010;
const int INF=0x3f3f3f3f;

struct Edge{
    int u,v,nxt;
    int flow;
}edge[EM<<1];

int n,m,cnt,head[VM];
int src,des,dep[VM];

void addedge(int cu,int cv,int cf){
    edge[cnt].u=cu;  edge[cnt].v=cv;  edge[cnt].flow=cf;
    edge[cnt].nxt=head[cu];  head[cu]=cnt++;

    edge[cnt].u=cv;  edge[cnt].v=cu;  edge[cnt].flow=0;
    edge[cnt].nxt=head[cv];  head[cv]=cnt++;
}

int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};

int legal(int i,int j,int k){
    int x=i+dir[k][0];
    int y=j+dir[k][1];
    return x>=1 && x<=n && y>=1 && y<=m;
}

int BFS(){  // 重新 建 图 （按 层数 建图）
    queue<int> q;
    while(!q.empty())
        q.pop();
    memset(dep,-1,sizeof(dep));
    dep[src]=0;
    q.push(src);
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=edge[i].nxt){
            int v=edge[i].v;
            if(edge[i].flow>0 && dep[v]==-1){  // 如果 可以  可以到达 但 还没有 访问
                dep[v]=dep[u]+1;
                q.push(v);
            }
        }
    }
    return dep[des]!=-1;
}
/*
int DFS(int u,int minx){  // 查找  路径上的 最小 的 流量
     if(u==des)
         return minx;
     int tmp;
     for(int i=head[u];i!=-1;i=edge[i].nxt){
         int v=edge[i].v;
         if(edge[i].flow>0 && dep[v]==dep[u]+1 && (tmp=DFS(v,min(minx,edge[i].flow)))){
             edge[i].flow-=tmp;
             edge[i^1].flow+=tmp;
             return tmp;
         }
         //printf("!!!!!!\n");
     }
     return 0;
 }
*/

int DFS(int u,int minx){
    int ans=0;
    if(u==des)
        return minx;
    for(int i=head[u];i!=-1 && ans<minx;i=edge[i].nxt){
        int v=edge[i].v;
        if(edge[i].flow>0 && dep[v]==dep[u]+1){
            int tmp=min(edge[i].flow,minx-ans);
            tmp=DFS(v,tmp);
            ans+=tmp;
            edge[i].flow-=tmp;
            edge[i^1].flow+=tmp;
        }
    }
    if(!ans)
        dep[u]=-2;
    return ans;
}

int Dinic(){
    int ans=0,tmp;
    while(BFS()){
        while(1){
            tmp=DFS(src,INF);
            if(tmp==0)
                break;
            ans+=tmp;
        }
    }
    return ans;
}

int main(){

    //freopen("input.txt","r",stdin);

    while(~scanf("%d%d",&n,&m)){
        cnt=0;
        memset(head,-1,sizeof(head));
        int x,sum=0;
        src=0;  des=n*m+1;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++){
                scanf("%d",&x);
                sum+=x;
                if((i+j)%2==0){
                    addedge(src,(i-1)*m+j,x);
                    for(int k=0;k<4;k++){
                        if(legal(i,j,k))
                            addedge((i-1)*m+j,(i+dir[k][0]-1)*m+(j+dir[k][1]),INF);
                    }
                }else{
                    addedge((i-1)*m+j,des,x);
                    for(int k=0;k<4;k++){
                        if(legal(i,j,k))
                            addedge((i+dir[k][0]-1)*m+(j+dir[k][1]),(i-1)*m+j,INF);
                    }
                }
            }
        int maxflow=Dinic();
        printf("%d\n",sum-maxflow);
    }
    return 0;
}




/*8MS
#include<iostream>
#include<cstring>
#include<cstdio>
#define INF 0x7fffffff
using namespace std;
int n,m,cnt=1,head[1001],h[1001],ans;
int xx[4]={0,0,1,-1},yy[4]={1,-1,0,0};
int mp[31][31],mark[31][31];
struct data{int to,next,v;}e[10001];
bool jud(int x,int y){if(x<1||y<1||x>n||y>m)return 0;return 1;}
void ins(int u,int v,int w)
{cnt++;e[cnt].to=v;e[cnt].v=w;e[cnt].next=head[u];head[u]=cnt;}
void insert(int u,int v,int w)
{ins(u,v,w);ins(v,u,0);}
void INS(){
     for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
           if((i+j)%2==0)
           {
              insert(0,mark[i][j],mp[i][j]);
              for(int k=0;k<4;k++)
              {
                  if(jud(i+xx[k],j+yy[k])){insert(mark[i][j],mark[i+xx[k]][j+yy[k]],INF);}
                   }
            }
           else insert(mark[i][j],n*m+1,mp[i][j]);
     }
bool bfs()
{
     int q[40002],t=0,w=1,i,now;
     memset(h,-1,sizeof(h));
     q[0]=h[0]=0;
     while(t<w)
     {
            now=q[t];t++;
            i=head[now];
            while(i)
            {
                  if(e[i].v&&h[e[i].to]<0){h[e[i].to]=h[now]+1;q[w++]=e[i].to;}
                  i=e[i].next;
             }
     }
     if(h[n*m+1]==-1)return 0;
     return 1;
     }
int dfs(int x,int f)
{
    if(x==n*m+1)return f;
    int i=head[x];
    int w,used=0;
    while(i)
    {
            if(e[i].v&&h[e[i].to]==h[x]+1)
            {
                w=f-used;
                w=dfs(e[i].to,min(w,e[i].v));   
                e[i].v-=w;
                e[i^1].v+=w;
                used+=w;
                if(used==f)return f;                      
                }
                i=e[i].next;
            }
    if(!used)h[x]=-1;
    return used;
    }
void dinic(){while(bfs())ans-=dfs(0,INF);}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
       for(int j=1;j<=m;j++)
       {
          scanf("%d",&mp[i][j]);
          ans+=mp[i][j];
          }
    int b=0,w=(n*m+1)/2; 
    for(int i=1;i<=n;i++)
       for(int j=1;j<=m;j++)
           if((i+j)%2==0){b++;mark[i][j]=b;}
               else{w++;mark[i][j]=w;}
    INS();
    dinic();
    printf("%d",ans);
    return 0;
    }
*/




/*Time Limit Exceed at Test 3
#include<stdio.h>
#define N 200
int m,n,max_val,current;
int a[N][N];
bool visited[N][N];

void handle(int i, int j){
	if( !visited[i][j] && 
		( (0==j) || !visited[i][j-1]) &&
		( (0==i) || !visited[i-1][j]) &&
		( (m-1==j) || !visited[i][j+1]) &&
		( (n-1==i) || !visited[i+1][j]) ){
		current += a[i][j];
		visited[i][j] = 1;
		if(n-1==j){
			if(m-1==i){
				if(max_val<current){
					max_val = current;
				}
			}else{
				handle(i+1,0);
			}
		}else{
			handle(i,j+1);
		}
		current -= a[i][j];
		visited[i][j]=0;
	}
	if(n-1==j){
		if(i<m-1){
			handle(i+1,0);
		}else{
			if(max_val<current){
				max_val = current;
			}
		}
	}else{
		handle(i,j+1);
	}
}


int main(){
	int i,j;
	scanf("%d%d",&m,&n);
	for(i=0; i<m; i++){
		for(j=0; j<n; j++){
			scanf("%d",&a[i][j]);
		}
	}
	handle(0,0);
	printf("%d\n",max_val);
}
*/