Bzoj1324 Exca王者之剑

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 591  Solved: 296

Description

 

Input

第一行给出数字N,M代表行列数.N,M均小于等于100 下面N行M列用于描述数字矩阵

Output

输出最多可以拿到多少块宝石

Sample Input

2 2
1 2
2 1

Sample Output

4

HINT

 

Source

2007Amber国家队论文

图论 网络流

将图二分染色

求最大点权独立集

最大点权独立集就等于总权值和减去最小权值覆盖

最小权值覆盖就是建好二分图以后的最小割

也就是说这是道模板题

然而我居然没有1A，原因是反向弧的流量建成了w而不是0

真鸡儿丢人，我退群吧

/*by SilverN*/
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f;
const int mxn=100010;
const int mx[5]={0,1,0,-1,0};
const int my[5]={0,0,1,0,-1};
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct edge{
    int v,nxt,f;
}e[mxn<<1];
int hd[mxn],mct=1;
void add_edge(int u,int v,int w){
    e[++mct].v=v;e[mct].nxt=hd[u];e[mct].f=w; hd[u]=mct;return;
}
void insert(int u,int v,int w){
    add_edge(u,v,w);add_edge(v,u,0);return;
}
queue<int>q;
int d[mxn];
int n,m,S,T;
bool BFS(){
    memset(d,0,sizeof d);
    d[S]=1;
    q.push(S);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=hd[u];i;i=e[i].nxt){
            int v=e[i].v;
            if(e[i].f && !d[v]){
                d[v]=d[u]+1;
                q.push(v);
            }
        }
    }
    return d[T];
}
int DFS(int u,int lim){
    if(u==T)return lim;
    int f=0,tmp;
    for(int i=hd[u],v;i;i=e[i].nxt){
        v=e[i].v;
        if(d[v]==d[u]+1 && e[i].f && (tmp=DFS(v,min(lim,e[i].f)))){
            e[i].f-=tmp;  e[i^1].f+=tmp;
            lim-=tmp;   f+=tmp;
            if(!lim)return f;
        }
    }
    d[u]=0;
    return f;
}
int Dinic(){
    int res=0;
    while(BFS())res+=DFS(S,0x3f3f3f3f);
    return res;
}
int id[120][120];
int mp[120][120];
int smm=0,dt=0;
void solve(){
    int i,j;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
            id[i][j]=++dt;
    S=0;T=dt+1;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++){
            if((i+j)&1){
                insert(S,id[i][j],mp[i][j]);
                for(int k=1;k<=4;k++){
                    int nx=i+mx[k],ny=j+my[k];
                    if(nx>0 && nx<=n && ny>0 && ny<=m)
                        insert(id[i][j],id[nx][ny],INF);
                     
                }
            }
            else insert(id[i][j],T,mp[i][j]);
        }
    return;
}
int main(){
    int i,j;
    n=read();m=read();
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++){
            mp[i][j]=read();
            smm+=mp[i][j];
        }
    solve();
    int res=Dinic();
//  printf("Res:%d
",res);
    printf("%d
",smm-res);
    return 0;
}