HDU 1565 체크 수 (1)

체크 개수 (1)
Time Limit: 1000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2410 Accepted Submission(s): 919
Problem Description
n * n 칸 짜 리 바둑판 을 드 리 겠 습 니 다. 칸 마다 마이너스 가 있 습 니 다.
그 중에서 몇 개의 수 를 꺼 내 서 임의의 두 개의 수가 있 는 칸 이 공공 변 이 없 게 한다. 즉, 취 하 는 수가 있 는 두 개의 칸 이 서로 인접 하지 못 하고 취 하 는 수의 합 이 가장 크다 는 것 이다.
Input
여러 개의 테스트 인 스 턴 스 를 포함 하고 모든 테스트 인 스 턴 스 는 하나의 정수 n 과 n * n 개의 비 마이너스 (n < 20) 를 포함한다.
Output
모든 테스트 인 스 턴 스 에 대해 출력 이 얻 을 수 있 는 가장 큰 합
Sample Input

   
   
   
   

    
    
    
    3
75 15 21 
75 15 28 
34 70 5 
   
   
   
   

Sample Output

   
   
   
   

    
    
    
    188
   
   
   
   

Author
ailyanlu
Source
Happy 2007
Recommend
8600

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
#include <string>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <list>

#define INT_INF 0x3fffffff
#define LL_INF 0x3fffffffffffffff
#define EPS 1e-12
#define MOD 1000000007
#define PI 3.141592653579798
#define N 3000
#define E 100000

using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
typedef double DB;

struct Edge
{
    int en,cap,flow,next;
} edge[E];
int head[N] , tot , now[N];
int source,sink,tot_num;
int pre[N] , dis[N] , gap[N];

void add_edge(int st,int en,int cap)
{
    edge[tot].en=en;
    edge[tot].cap=cap;
    edge[tot].flow=0;
    edge[tot].next=head[st];
    head[st]=tot++;

    edge[tot].en=st;
    edge[tot].cap=0;
    edge[tot].flow=0;
    edge[tot].next=head[en];
    head[en]=tot++;
}

void augment(int flow)
{
    for(int i=source;i!=sink;i=edge[now[i]].en)
    {
        edge[now[i]].flow+=flow;
        edge[now[i]^1].flow-=flow;
    }
}

int sap()
{
    memset(dis,0,sizeof(dis));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(pre));
    for(int i=0;i<tot_num;i++)
        now[i]=head[i];
    gap[0]=tot_num;
    int point=source,flow=0,min_flow=INT_INF;
    while(dis[source]<tot_num)
    {
        bool fg=false;
        for(int i=now[point];i!=-1;i=edge[i].next)
            if(edge[i].cap-edge[i].flow>0 && dis[point]==dis[edge[i].en]+1)
            {
                min_flow=min(min_flow,edge[i].cap-edge[i].flow);
                now[point]=i;
                pre[edge[i].en]=point;
                point=edge[i].en;
                if(point==sink)
                {
                    flow+=min_flow;
                    augment(min_flow);
                    point=source;
                    min_flow=INT_INF;
                }
                fg=true;
                break;
            }
        if(fg) continue;
        if(--gap[dis[point]]==0) break;
        int Min=tot_num;
        for(int i=head[point];i!=-1;i=edge[i].next)
            if(edge[i].cap-edge[i].flow>0 && Min>dis[edge[i].en])
            {
                Min=dis[edge[i].en];
                now[point]=i;
            }
        gap[dis[point]=Min+1]++;
        if(point!=source) point=pre[point];
    }
    return flow;
}

int build(int n,int m)
{
    memset(head,-1,sizeof(head));
    tot=0;
    int ans=0;
    source=N-2; sink=N-1; tot_num=N;
    for(int i=0; i<n; i++)
        for(int j=0,pos,val; j<m; j++)
        {
            scanf("%d",&val);
            ans+=val;
            pos=i*m+j;
            if((i+j)%2==0) add_edge(pos,sink,val);
            else
            {
                add_edge(source,pos,val);
                if(i+1<n) add_edge(pos,(i+1)*m+j,INT_INF);
                if(i-1>=0) add_edge(pos,(i-1)*m+j,INT_INF);
                if(j+1<m) add_edge(pos,i*m+j+1,INT_INF);
                if(j-1>=0) add_edge(pos,i*m+j-1,INT_INF);
            }
        }
    return ans;
}

int main()
{
    int n,m;
    while(scanf("%d",&n)!=EOF)
    {
        m=n;
        int ans=build(n,m);
        ans-=sap();
        printf("%d
",ans);
    }
    return 0;
}