首页 > 代码库 > HDU 4859 海岸线(最大流最小割)
HDU 4859 海岸线(最大流最小割)
难得的中文题,就不翻译了。
输入第一行为T,表示有T组测试数据。
每组数据以两个整数N和M开始,表示地图的规模。接下来的N行,每一行包含一个长度为M的字符串,表示地图,‘.’表示陆地,’E’表示浅海域,’D’表示深海域。
[Technical Specification]
1. 1 <= T <= 100
2. 1 <= N, M <= 47
题意:假设地图为一个N*M的格子,其中一些是陆地,一些是可以填充的浅海域,一些是不可填充的深海域。这里定义海岸线的长度为一个联通块陆地(可能包含浅海域填充变为的陆地)的边缘长度,两个格子至少有一个公共边,则视为联通。
值得注意的是,这里Z市的陆地区域可以是不联通的,并且整个地图都处在海洋之中,也就是说,Z市是由一些孤岛组成的,比如像,夏威夷?
你的任务是,填充某些浅海域,使得所有岛屿的海岸线之和最长。
输出最长海岸线和。
Solution:
一开始以为是贪心,后来发现这个规划问题好像不可以。算法是最小割,然后最小割==最大流,算法是最大流。
由于海岸线一定是海与岸的边界(废话哈哈哈~~~),不妨在给出的网格图外层加多层深海(的孤独~~~)
显然(i+j)%2==0的格子只可能和(i+j)%2==1的格子有海岸线。不妨二部图。
我们需要求的是最可能多的相邻不同对(<‘D‘, ‘.‘>),也就是求尽可能少的相邻相同对(<‘D‘, ‘D‘>或<‘.‘, ‘.‘>)。
建模如图,感觉这样画图比较直观。只画出部分,还有EE相连之类的。
左右流量均为inf,中间流量为1。
显然最大流跑出来的是最少的相同对数。
答案就是 sum - mf
#pragma comment (linker,"/STACK:102400000,102400000")#include <iostream>#include <cstdio>#include <cstring>#include <string>#include <vector>#include <algorithm>#include <cmath>#include <queue>#include <set>#include <map>#include <stack>using namespace std;#define mxn 2600#define mxe 26000#define inf 0x3f3f3f3fstruct SAP{ int dis[mxn],pre[mxn],gap[mxn],arc[mxn],f[mxe],cap[mxe]; int head[mxn],nxt[mxe],vv[mxe],e; void init(){e=0;memset(head,-1,sizeof(head));} void addedge(int u,int v,int c){ vv[e]=v,cap[e]=c,nxt[e]=head[u],head[u]=e++; vv[e]=u,cap[e]=0,nxt[e]=head[v],head[v]=e++; } int max_flow(int s,int t,int n){ int q[mxn],j,mindis,ans=0,ht=0,tl=1; int u,v,low; bool found,vis[mxn]; memset(dis,0,sizeof(dis)); memset(gap,0,sizeof(gap)); memset(vis,0,sizeof(vis)); memset(arc,0,sizeof(arc)); memset(f,0,sizeof(f)); q[0]=t,vis[t]=true,dis[t]=0,gap[0]=1; while(ht<tl){ int u = q[ht++]; for(int i=head[u];i!=-1;i=nxt[i]){ v = vv[i]; if(!vis[v]){ vis[v]=true; dis[v]=dis[u]+1; q[tl++]=v; gap[dis[v]]++; arc[v]=head[v]; } } } u=s;low=inf;pre[s]=s; while(dis[s]<n){ found = false; for(int &i = arc[u];i!=-1;i=nxt[i]){ if(dis[vv[i]]==dis[u]-1 && cap[i]>f[i]){ found = true; v=vv[i]; low = min(low, cap[i]-f[i]); pre[v]=u;u=v; if(u==t){ while(u!=s){ u=pre[u]; f[arc[u]]+=low; f[arc[u]^1]-=low; } ans+=low;low=inf; } break; } } if(found) continue; mindis = n; for(int i=head[u];i!=-1;i=nxt[i]){ if(mindis>dis[vv[i]] && cap[i]>f[i]){ mindis = dis[vv[j=i]]; arc[u]=i; } } if(--gap[dis[u]]==0) return ans; dis[u] = mindis+1; gap[dis[u]]++; u=pre[u]; } return ans; }}sap;char maze[55][55];int dx[]={0,1,0,-1};int dy[]={1,0,-1,0};int main(){ int t,n,m,ca=0; scanf("%d",&t); while(t--){ scanf("%d%d",&n,&m); for(int i=1;i<=n;++i) scanf("%s",maze[i]+1); for(int i=0;i<=n+1;++i) maze[i][0]=maze[i][m+1]=‘D‘; for(int j=0;j<=m+1;++j) maze[0][j]=maze[n+1][j]=‘D‘; sap.init(); int src = http://www.mamicode.com/(n+2)*(m+2);"Case %d: %d\n",++ca,sum-mf); } return 0;}
HDU 4859 海岸线(最大流最小割)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。