1. 程式人生 > >Exca王者之劍

Exca王者之劍

include des while void image cstring pri har 最小

Description
技術分享圖片
技術分享圖片

Input
第一行給出數字N,M代表行列數.N,M均小於等於100 下面N行M列用於描述數字矩陣

Output
輸出最多可以拿到多少塊寶石

Sample Input
2 2
1 2
2 1

Sample Output
4

可以發現,所有能取得寶石的點必定是不相鄰的,如果本題點權都為1,我們只需要建立二分圖然後求最大獨立點集即可。但是這題有點權,也就是變成了求最大帶權獨立點集,那麽我們考慮不用二分圖,使用網絡流,從源點向所有白點連一條流量為點權的邊,所有黑點向匯點連一條流量為點權的邊,所有白點向其相鄰的黑點連邊,答案即為總點權減去最小割

/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
#define lowbit(x) ((x)&-(x))
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline char gc(){
    static char buf[1000000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
inline int frd(){
    int x=0,f=1; char ch=gc();
    for (;ch<'0'||ch>'9';ch=gc())   if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline int read(){
    int x=0,f=1; char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar())  if (ch=='-')    f=-1;
    for (;ch>='0'&&ch<='9';ch=getchar())    x=(x<<3)+(x<<1)+ch-'0';
    return x*f;
}
inline void print(int x){
    if (x<0)    putchar('-'),x=-x;
    if (x>9)    print(x/10);
    putchar(x%10+'0');
}
const int N=1e4,M=6e4;
const int dx[4]={0,0,-1,1};
const int dy[4]={-1,1,0,0};
int pre[M+10],now[N+10],child[M+10],val[M+10];
int h[N+10],deep[N+10];
int n,m,S,T,Ans,tot=1;
void join(int x,int y,int z){pre[++tot]=now[x],now[x]=tot,child[tot]=y,val[tot]=z;}
void insert(int x,int y,int z){join(x,y,z),join(y,x,0);}
bool in_map(int x,int y){return x>0&&x<=n&&y>0&&y<=m;}
int G(int x,int y){return (x-1)*m+y;}
bool bfs(){
    int head=1,tail=1;
    memset(deep,255,sizeof(deep));
    deep[h[head]=S]=0;
    for (;head<=tail;head++){
        int Now=h[head];
        for (int p=now[Now],son=child[p];p;p=pre[p],son=child[p]){
            if (!~deep[son]&&val[p]){
                deep[h[++tail]=son]=deep[Now]+1;
                if (son==T) return 1;
            }
        }
    }
    return 0;
}
int dfs(int x,int v){
    if (x==T)   return v;
    int Ans=0;
    for (int p=now[x],son=child[p];p;p=pre[p],son=child[p]){
        if (deep[son]>deep[x]&&val[p]){
            int t=dfs(son,min(val[p],v));
            val[p]-=t,v-=t;
            val[p^1]+=t,Ans+=t;
            if (!v) break;
        }
    }
    if (!Ans)   deep[x]=-1;
    return Ans;
}
int main(){
    n=read(),m=read(),S=n*m+1,T=S+1;
    for (int i=1;i<=n;i++){
        for (int j=1;j<=m;j++){
            int v=read(); Ans+=v;
            if ((i+j)&1){
                insert(S,G(i,j),v);
                for (int k=0;k<4;k++){
                    int tx=i+dx[k],ty=j+dy[k];
                    if (!in_map(tx,ty)) continue;
                    insert(G(i,j),G(tx,ty),inf);
                }
            }else   insert(G(i,j),T,v);
        }
    }
    while (bfs())   Ans-=dfs(S,inf);
    printf("%d\n",Ans);
    return 0;
}

Exca王者之劍