首页 > 代码库 > 网络流24题-方格取数
网络流24题-方格取数
题目
方格中取数若干,两两不相邻,求最大选数和。
样例:
3 3
1 2 3
3 2 3
2 3 1
输出:
11
ans=2+3+3+3=11
黑白染色
建成二分图:
中间这些边连成INF(即不限制流量)
其实就是求最大独立集
定理:|二分图最大独立集|=|顶点数|-|二分图最大匹配数|
这个,我不会证明。
所以:|二分图最大独立集|=|总价值|-|二分图最大匹配价值|
考虑构图:s=>黑点=>白点=>t
首先,对与第[i][j]格,有价值为val[i][j]
edge(s,黑点,val[i][j])
edge(白点,t,val[i][j])
然后,对所有黑点
edge(黑点,相邻白点,INF)
令,|总价值|=sum=∑∑val[i][j] (1<=i<=n,1<=j<=m)
|二分图最大匹配价值|<=>|最小割|<=>|最大流| 这个就不分析了。
ans=sum-MaxFlow(G)
于是,代码呼之欲出
1 #include <algorithm> 2 #include <iostream> 3 #include <cstring> 4 #include <cstdlib> 5 #include <cstdio> 6 #include <vector> 7 #include <cmath> 8 #include <queue> 9 #include <map> 10 #include <set> 11 using namespace std; 12 13 inline void read(int &ans){ 14 ans=0;char x=getchar();int f=1; 15 while(x<‘0‘||x>‘9‘){if(x==‘-‘)f=-1;x=getchar();} 16 while(x>=‘0‘&&x<=‘9‘)ans=ans*10+x-‘0‘,x=getchar(); 17 ans*=f; 18 } 19 const int MAXN=30+10; 20 const int MAXV=MAXN*MAXN; 21 const int MAXE=MAXV<<4; 22 const int INF=2e9; 23 struct net{int v,next,c;}E[MAXE]; 24 int head[MAXV],tot,n,m; 25 int vis[MAXN][MAXN],val[MAXN][MAXN],cnt; 26 void edge(int u,int v,int c){ 27 E[tot]=(net){v,head[u],c},head[u]=tot++; 28 E[tot]=(net){u,head[v],0},head[v]=tot++; 29 } 30 struct Net{ 31 int s,t; 32 int dep[MAXV],cur[MAXV]; 33 void init(int a,int b){ 34 s=a,t=b,memset(head,-1,sizeof(head)); 35 } 36 bool bfs(){ 37 memset(dep,0,sizeof(dep)); 38 queue<int>q; 39 dep[s]=1,q.push(s); 40 while(!q.empty()){ 41 int u=q.front();q.pop(); 42 for(int i=head[u];~i;i=E[i].next){ 43 int v=E[i].v; 44 if(!dep[v]&&E[i].c>0){ 45 dep[v]=dep[u]+1; 46 q.push(v); 47 } 48 } 49 if(dep[t])return true; 50 } 51 return false; 52 } 53 int dfs(int u,int f){ 54 if(u==t)return f; 55 int ans=0,cup; 56 for(int &i=cur[u];~i;i=E[i].next){ 57 int v=E[i].v; 58 if(E[i].c>0&&dep[v]==dep[u]+1){ 59 cup=dfs(v,min(f-ans,E[i].c)); 60 E[i].c-=cup; 61 E[i^1].c+=cup; 62 ans+=cup; 63 if(ans==f)return ans; 64 } 65 } 66 return ans; 67 } 68 int work(){ 69 int ans=0; 70 while(bfs()){ 71 for(int i=0;i<=t;i++)cur[i]=head[i]; 72 ans+=dfs(s,INF); 73 } 74 return ans; 75 } 76 }a; 77 78 int Tx[4]={1,-1,0,0}; 79 int Ty[4]={0,0,-1,1}; 80 int main(){ 81 read(n),read(m); 82 int s=n*m+1,t=s+1; 83 a.init(s,t); 84 for(int i=1;i<=n;i++) 85 for(int j=1;j<=m;j++) 86 read(val[i][j]),vis[i][j]=++cnt; 87 int SUM=0; 88 for(int i=1;i<=n;i++) 89 for(int j=1;j<=m;j++){ 90 SUM+=val[i][j]; 91 if((i+j)&1)edge(s,vis[i][j],val[i][j]); 92 else edge(vis[i][j],t,val[i][j]); 93 int I,J; 94 if((i+j)&1) 95 for(int k=0;k<4;k++){ 96 I=i+Tx[k],J=j+Ty[k]; 97 if(I<1||I>n||J<1||J>m)continue; 98 edge(vis[i][j],vis[I][J],INF); 99 } 100 } 101 cout<<SUM-a.work()<<endl; 102 return 0; 103 }
网络流24题-方格取数
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。