http://acm.hdu.edu.cn/showproblem.php?pid=1569?

題解：網(wǎng)絡(luò)流+最大流+最小割

把格子染成黑白兩色，形如下圖： ?有點(diǎn)丑。。好吧，不止一點(diǎn)

其實(shí)黑白兩色也可以反過來，不過這不是重點(diǎn)。

我們把每個(gè)格標(biāo)個(gè)號(hào)：

然后我們可以把格子抽象成一個(gè)點(diǎn)，然后分開：

最后就是連邊：

左邊那個(gè)藍(lán)藍(lán)的是源點(diǎn)，然后那個(gè)白色的大圈是匯點(diǎn)。

從源點(diǎn)連向黑點(diǎn)的邊的權(quán)值（容量）是那個(gè)黑點(diǎn)的值，然后從白點(diǎn)連向匯點(diǎn)的邊也是權(quán)值。

兩邊的點(diǎn)時(shí)候不能同時(shí)選的（因?yàn)橛泄策?#xff09;，邊權(quán)是INF（盡量大）。

題目要求最大點(diǎn)權(quán)獨(dú)立集（注意：是點(diǎn)權(quán)），就等于所有點(diǎn)的值減去最小點(diǎn)權(quán)覆蓋集，就是減去最小割（最大流）。

/* *@Author: STZG *@Language: C++ */ #include <bits/stdc++.h> #include<iostream> #include<algorithm> #include<cstdlib> #include<cstring> #include<cstdio> #include<string> #include<vector> #include<bitset> #include<queue> #include<deque> #include<stack> #include<cmath> #include<list> #include<map> #include<set> //#define DEBUG #define RI register int #define endl "\n" using namespace std; typedef long long ll; //typedef __int128 lll; const int N=10000+10; const int M=100000+10; const int MOD=1e9+7; const double PI = acos(-1.0); const double EXP = 1E-8; const ll INF = 0x3f3f3f3f3f3f3f; int s,t,n,m,k,p,l,r,u,v; ll w; ll ans,cnt,flag,temp,sum; int dis[N]; struct node{int u,v;ll c;node(){};node(int form,int to,ll cap):u(form),v(to),c(cap){} }; vector<node>edge; vector<int> G[N]; void Addedge(int u,int v,ll cap){edge.push_back({u,v,cap});edge.push_back({v,u,0});int sz=edge.size();//cout<<sz<<endl;G[u].push_back(sz-2);G[v].push_back(sz-1); } bool bfs(int u){memset(dis,-1,sizeof(dis));dis[u]=0;queue<int>q;q.push(u);while(!q.empty()){int u=q.front();q.pop();for(int i=0;i<G[u].size();i++){node e=edge[G[u][i]];//cout<<u<<" "<<e.v<<endl;if(dis[e.v]<0&&e.c>0){dis[e.v]=dis[u]+1;q.push(e.v);}}}return dis[t]>0; } ll dfs(int u,ll flow){if(u==t)return flow;ll now;for(int i=0;i<G[u].size();i++){node e=edge[G[u][i]];if(e.c>0&&dis[u]+1==dis[e.v]&&(now=dfs(e.v,min(flow,e.c)))){edge[G[u][i]].c-=now;edge[G[u][i]^1].c+=now;return now;}}return 0; } void dinic(){while(bfs(s)){ll res=0;while((res=dfs(s,INF))){ans+=res;}} } void init(){s=0;t=m*n+1;for(int i=0;i<=t;i++)G[i].clear();edge.clear();ans=0;sum=0; } int main() { #ifdef DEBUGfreopen("input.in", "r", stdin);//freopen("output.out", "w", stdout); #endif//ios::sync_with_stdio(false);//cin.tie(0);//cout.tie(0);//scanf("%d",&t);//int T=0;while(~scanf("%d%d",&m,&n)){init();for(int i=1;i<=m;i++){for(int j=1;j<=n;j++){scanf("%lld",&w);sum+=w;int id=(i-1)*n+j;if((i+j)%2){Addedge(s,id,w);if(j<n)Addedge(id,id+1,INF);if(j>1)Addedge(id,id-1,INF);if(i<m)Addedge(id,id+n,INF);if(i>1)Addedge(id,id-n,INF);}else{Addedge(id,t,w);}}}dinic();printf("%lld\n",sum-ans);}#ifdef DEBUGprintf("Time cost : %lf s\n",(double)clock()/CLOCKS_PER_SEC); #endif//cout << "Hello world!" << endl;return 0; }

?

總結(jié)

以上是生活随笔為你收集整理的方格取数(2)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。