Exca王者之劍
阿新 • • 發佈:2019-01-10
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王者之劍