时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld
题目描述
随着海上运输石油泄漏的问题,一个新的有利可图的行业正在诞生,那就是撇油行业。如今,在墨西哥湾漂浮的大量石油,吸引了许多商人的目光。这些商人们有一种特殊的飞机,可以一瓢略过整个海面20米乘10米这么大的长方形。(上下相邻或者左右相邻的格子,不能斜着来)当然,这要求一瓢撇过去的全部是油,如果一瓢里面有油有水的话,那就毫无意义了,资源完全无法利用。现在,商人想要知道,在这片区域中,他可以最多得到多少瓢油。
地图是一个N×N的网络,每个格子表示10m×10m的正方形区域,每个区域都被标示上了是油还是水
输入描述:
测试输入包含多条测试数据
测试数据的第一行给出了测试数据的数目T(T<75)
每个测试样例都用数字N(N<50)来表示地图区域的大小,接下来N行,每行都有N个字符,其中符号’.’表示海面、符号’#’表示油面。
输出描述:
输出格式如下“Case X: M”(X从1开始),M是商人可以最多得到的油量。
示例1
输入
1 6 ...... .##... ...... .#..#. .#..## ......
输出
Case 1: 3分析: 本质是考察 二分图的最大匹配。 首先需要对油面进行 编号, 然后,用 邻接表建立二分图, 最后就是经典的 匈牙利算法 的运用
#include <cstdio>
#include <iostream>
#include <vector>
#include <cstring>
using namespace std;
const int N = 55;
char map[N][N];
bool vis[N*N];
int from[N*N]; //Not N!
vector<int> mp[N*N]; //邻接表
void Read(int &n){
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%s",map[i]+1);//从地址mp[i][1]开始接收字符串
}
void PreProcess(int n, int &cnt){
//对油量区域编号
int num[N][N];
for(int t=1;t<=n;t++)
for(int r=1;r<=n;r++){
if(map[t][r] == '#')
num[t][r] = ++cnt;
}
for(int i=1;i<=cnt;i++) mp[i].clear(); //清空元素(not n!)
//用邻接表构建二分图
for(int k=1;k<=n;k++)
for(int j=1;j<=n;j++){
if(map[k][j] == '#'){
if(k>1 && map[k-1][j]=='#') mp[num[k][j]].push_back(num[k-1][j]);
if(k<n && map[k+1][j]=='#') mp[num[k][j]].push_back(num[k+1][j]);
if(j>1 && map[k][j-1]=='#') mp[num[k][j]].push_back(num[k][j-1]);
if(j<n && map[k][j+1]=='#') mp[num[k][j]].push_back(num[k][j+1]);
}
}
}
bool AgumentPath(int i){
for(int k=0;k<mp[i].size();k++){
int tmp = mp[i].at(k);
if(!vis[tmp]){ //与i关联的点tmp不在增广路上
vis[tmp] = true;
//tmp是未盖点,或从tmp对应项出发有增广路
if(from[tmp]==-1 || AgumentPath(from[tmp])){
from[tmp] = i;
return true;
}
}
}
return false;
}
int Hungary(int cnt){
memset(from,-1,sizeof(from)); //初始值-1表示未盖点
int ans = 0;
for(int i=1;i<=cnt;i++){
memset(vis,false,sizeof(vis));//每次找增广路vis都要init为false
if(AgumentPath(i)) ans ++;
}
return ans;
}
int main(){
int t,n,ans,cnt;
scanf("%d",&t);
for(int i=1;i<=t;i++){
Read(n);
cnt = 0;
PreProcess(n,cnt);
//对每个区域的增广路会多计算一倍,且所有区域合并与分开计算增广路并无差别
ans = Hungary(cnt)/2;
printf("Case %d: %d\n",i,ans);
}
return 0;
}