HDU 1565格子取数(1)(最小割-Dinic)

8513 ワード

Descriptionはあなたにn*nの格子の碁盤をあげて、各格子の中に非負の数があります.任意の2つの数が存在する格子に共通のエッジがないように、いくつかの数を取り出します.つまり、取得した数が存在する2つの格子が隣接することができず、取得した数の和が最大になります.Inputは複数の試験例を含み、各試験例は1つの整数nとn*n個の負でない数(n<=20)Outputを含み、各試験例について、取得可能な最大値とSample Input 3 75 15 21 75 15 28,725 Sample Output 188 Solutionの最大獲得数=総数-最小不要値、すなわちsum-最小割を出力する.奇数格子の点連源点では容量はその上の値、偶数格子の点向匯点連容量はその上の値、奇数点向の到達可能な点(全偶数点)連容量はINFの点である.(INFのエッジは絶対に最小カットとして削除されません).1つの奇数点が当を保持し、彼と連係している偶数点が削除された場合(そうでなければS、T連通)の同理偶数点のみが削除される.Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
#define maxn 444
#define maxm 33333
#define INF 0x3f3f3f3f
int head[maxn],cur[maxn],d[maxn],st[maxm],s,e,no,n;
struct point
{
    int u,v,flow,next;
    point(){};
    point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)
{
    p[no]=point(x,y,head[x],z); 
    head[x]=no++;
    p[no]=point(y,x,head[y],0); 
    head[y]=no++;
}
void init()
{
    memset(head,-1,sizeof(head));
    no=0;
}
bool bfs()
{
    int i,x,y;
    queue<int>q;
    memset(d,-1,sizeof(d));
    d[s]=0; 
    q.push(s);
    while(!q.empty())
    {
        x=q.front();    
        q.pop();
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&& d[y = p[i].v]<0)
            {
                d[y]=d[x]+1;
                if(y==e)    
                    return true;
                q.push(y);
            }
        }
    }
    return false;
}
int dinic()
{
    int i,loc,top,x=s,nowflow,maxflow=0;
    while(bfs()){
        for(i=s;i<=e;i++)   
            cur[i]=head[i];
        top=0;
        while(true)
        {
            if(x==e)
            {
                nowflow=INF;
                for(i=0;i<top;i++)
                {
                    if(nowflow>p[st[i]].flow)
                    {
                        nowflow=p[st[i]].flow;
                        loc=i;
                    }
                }
                for(i=0;i<top;i++)
                {
                    p[st[i]].flow-=nowflow;
                    p[st[i]^1].flow+=nowflow;
                }
                maxflow+=nowflow;
                top=loc;    
                x=p[st[top]].u;
            }
            for(i=cur[x];i!=-1;i=p[i].next)
                if(p[i].flow&&d[p[i].v]==d[x]+1) 
                    break;
            cur[x]=i;
            if(i!=-1)
            {
                st[top++]=i;
                x=p[i].v;
            }
            else 
            {
                if(!top)    
                    break;
                d[x]=-1;
                x=p[st[--top]].u;
            }
        }
    }
    return maxflow;
}
int N,map[22][22];
int main()
{
    while(~scanf("%d",&N))
    {
        init();//  
        for(int i=1;i<=N;i++)
            for(int j=1;j<=N;j++)
                scanf("%d",&map[i][j]); 
        s=0;// 0 
        e=N*N+1;// N*N+1 
        n=N*N;
        int sum=0;
        for(int i=1;i<=N;i++)
            for(int j=1;j<=N;j++)
            {
                if((i+j)%2)//  
                {
                    add(s,(i-1)*N+j,map[i][j]);// 
                    // INF  
                    if(j>1)
                        add((i-1)*N+j,(i-1)*N+j-1,INF);
                    if(i>1)
                        add((i-1)*N+j,(i-2)*N+j,INF);
                    if(j<N)
                        add((i-1)*N+j,(i-1)*N+j+1,INF);
                    if(i<N)
                        add((i-1)*N+j,i*N+j,INF);
                }
                else//  
                    add((i-1)*N+j,e,map[i][j]);//  
                sum+=map[i][j];
            }
        int ans=sum-dinic();//sum-  
        printf("%d
"
,ans); } return 0; }