题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2547

挺有意思的一道题,这道题可以划分成几个小题。。。。。。。

题目大意:

三个兵种在一个n*m的图上,图的每一个格子有一个高度,三个兵种一种走不下降的路,一种走不上升的路,一种随便走。数量分别是k,k,1

然后t个目的地,每个目的地必须有ri个兵,ri之和为2*k+1

然后你可以交换兵,让不能走的情况又变成能走,例如,第一种兵走到了山顶,四周都是下降的,这时候就可以和第二,三种兵交换位置,继续走下去。。。

【解题思路】

这道题分解下来就是先来一个spfa找每一个兵到每一个点的最小交换次数,然后就是二分这个一共的交换次数,接着就是二分图匹配。。。

我们跑spfa,以每个兵开始扩展,遇到不能走就换一次,记录每个兵到每个点的最小交换次数,然后把当前兵到指定点的值记录入一个新的数组(表示二分图的连接情况)

然后二分。。因为最少是不交换,最大交换2*k次,因为可以每个兵和天兵交换一次然后去终点,天兵是随便走,所以最多交换2*k次,然后二分次数。。二分当中的check就是用二分图匹配了。。。Check成立的判断时交换次数+最大匹配数,如果这个值大于了兵的数量,说明这个方案是可行的。。。

这是大致思路,但是还是要注意一些细节的地方。

Spfa要注意当前的兵种是在变化的,所以要结合到当前位置之前的交换次数和最初状态

最初为0,交换偶数次,当前还是0,交换奇数次,当前为1

最初为1,交换偶数次,当前还是1,交换奇数次,当前为0

假设开始是c,次数是t,这个最终结果就是c^(t&1)

另外一个就是要注意二分图匹配时,除了要求vis[i]==0之外,还有x兵到i点的次数小于我们二分出来的限制。。。

然后就还是结合代码理解吧,我看别人的思路都不是很懂,还是结合代码读懂的。。。

 #include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstdlib>
#include<queue>
#define maxn 105
using namespace std; const int dx[]={,-,,,};
const int dy[]={,,,-,}; int att[maxn][maxn],h[maxn][maxn];
int f[maxn][maxn];//表示到达这个位置的最小交换次数
int n,m,k,t,r;
pair<int,int>a[maxn];
pair<int,int>b[maxn];
queue<pair<int,int> >q; void init()
{
scanf("%d%d%d%d",&n,&m,&k,&t);
for(int i=;i<=*k+;i++)
{
scanf("%d%d",&a[i].first,&a[i].second);
}
int tot=;
for(int i=;i<=t;i++)
{
scanf("%d%d%d",&b[tot].first,&b[tot].second,&r);
r--;tot++;
for(r;r>=;r--,tot++){
b[tot]=b[tot-];
}
}
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
scanf("%d",&h[i][j]);
} void bfs(int x,int y,int cnt){
memset(f,0x3f3f3f,sizeof(f));
f[x][y]=;
q.push(make_pair(x,y));
while(!q.empty()){
int nx=q.front().first,ny=q.front().second;
for(int i=;i<=;i++)
{
int nowx=nx+dx[i],nowy=ny+dy[i];
if(nowx<||nowx>n||nowy<||nowy>m)continue;
int tt=;
if(cnt^(f[nx][ny]&)){
if(h[nowx][nowy]>h[nx][ny])tt=;
}else if(h[nowx][nowy]<h[nx][ny])tt=;
if(f[nowx][nowy]>f[nx][ny]+tt){
f[nowx][nowy]=f[nx][ny]+tt;
q.push(make_pair(nowx,nowy));
} }
q.pop();
} } int vis[maxn*],ck[maxn*]; int ok(int x,int limit){
for(int i=;i<=*k+;i++)
{
if(vis[i]==&&att[x][i]<=limit){
vis[i]=;
if(ck[i]==||ok(ck[i],limit)){
ck[i]=x;return ;
}
}
}
return ;
} int check(int x){
memset(ck,,sizeof(ck));
int tot=;
for(int i=;i<=*k;i++){
memset(vis,,sizeof(vis));
if(ok(i,x))tot++;
}
if(tot+x>=*k)return ;
else return ;
} int ans(){
int l=,r=*k;//最坏情况就是天兵挨着和每一个交换,最多换2*k次
while(l<r){
int mid=(l+r)>>;
if(check(mid))r=mid;else l=mid+;
}
return l;
} int main()
{
init();
for(int i=;i<=*k;i++)
{
if(i<=k)bfs(a[i].first,a[i].second,);
else bfs(a[i].first,a[i].second,);
for(int j=;j<=*k+;j++){
att[i][j]=f[b[j].first][b[j].second];
}
}
printf("%d\n",ans());
}

总结:

把一个大题打乱成多个小题是一个很有效的方式

注意分析状态的变化,比如题中的兵种在spfa的时候是变化的

[bzoj2547]玩具兵<Spfa+二分+匈牙利算法>的更多相关文章

  1. [ACM_图论] The Perfect Stall 完美的牛栏(匈牙利算法、最大二分匹配)

    描述 农夫约翰上个星期刚刚建好了他的新牛棚,他使用了最新的挤奶技术.不幸的是,由于工程问题,每个牛栏都不一样.第一个星期,农夫约翰随便地让奶牛们进入牛栏,但是问题很快地显露出来:每头奶牛都只愿意在她们 ...

  2. hdu 2444 The Accomodation of Students(二分匹配 匈牙利算法 邻接表实现)

    The Accomodation of Students Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K ( ...

  3. 51 NOd 2006 飞行员配对(匈牙利算法二分匹配)

    题目来源: 网络流24题 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏  关注 第二次世界大战时期,英国皇家空军从沦陷国征募了大量外籍飞行员.由皇家空军派出的每一 ...

  4. 匈牙利算法模板 hdu 1150 Machine Schedule(二分匹配)

    二分图:https://blog.csdn.net/c20180630/article/details/70175814 https://blog.csdn.net/flynn_curry/artic ...

  5. hdu2063 最大二分匹配(匈牙利算法)

    过山车 Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...

  6. 【01染色法判断二分匹配+匈牙利算法求最大匹配】HDU The Accomodation of Students

    http://acm.hdu.edu.cn/showproblem.php?pid=2444 [DFS染色] #include<iostream> #include<cstdio&g ...

  7. [Bzoj 2547] [Ctsc2002] 玩具兵

    2547: [Ctsc2002]玩具兵 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 317  Solved: 152[Submit][Status] ...

  8. HDU5090——Game with Pearls(匈牙利算法|贪心)(2014上海邀请赛重现)

    Game with Pearls Problem DescriptionTom and Jerry are playing a game with tubes and pearls. The rule ...

  9. Ural1109_Conference(二分图最大匹配/匈牙利算法/网络最大流)

    解题报告 二分图第一题. 题目描写叙述: 为了參加即将召开的会议,A国派出M位代表,B国派出N位代表,(N,M<=1000) 会议召开前,选出K队代表,每对代表必须一个是A国的,一个是B国的; ...

随机推荐

  1. JZOJ 5246. 【NOIP2017模拟8.8A组】Trip(trip)

    5246. [NOIP2017模拟8.8A组]Trip(trip) (File IO): input:trip.in output:trip.out Time Limits: 1500 ms Memo ...

  2. 5W随想

    什么时候都要想5个W:what.where.when.why.who 比如:提bug的时候:这个问题是什么问题,在哪里发现的,什么时候发现的/什么时候修复好,怎么发现的/为什么是bug,谁发现的/指给 ...

  3. Redis07——Redis到底能用在什么地方(下)

    在前一篇文章中,我们已经介绍过Redis的一些实际应用.如KV缓存.分布式锁.消息队列,由于篇幅原因,并未介绍完全.接下来将继续为各位带来Redis的更多应用. bitmat(位图) 实现 位图的基本 ...

  4. 使用.Net Core编写命令行工具(CLI)

    命令行工具(CLI) 命令行工具(CLI)是在图形用户界面得到普及之前使用最为广泛的用户界面,它通常不支持鼠标,用户通过键盘输入指令,计算机接收到指令后,予以执行. 通常认为,命令行工具(CLI)没有 ...

  5. Lucene查询语法汇总

    目录 一.单词查询 二.通配符查询 三.模糊查询 四.近似查询 五.范围查询 六.优先级查询 七.逻辑操作 八.括号分组 九.转义特殊字符 Lucene是目前最为流行的开源全文搜索引擎工具包,提供了完 ...

  6. Simulink仿真入门到精通(十三) Simulink创建自定义库

    当用户自定义了一系列模块之后,可以自定义模块库将同类自定义模块显示到Simulink Browser中,作为库模块方便地拖曳到新建模型中. 建立这样的自定义库需要3个条件: 建立library的mdl ...

  7. 一个交互式可视化Python库——Bokeh

    本篇为<Python数据可视化实战>第十篇文章,我们一起学习一个交互式可视化Python库--Bokeh. Bokeh基础 Bokeh是一个专门针对Web浏览器的呈现功能的交互式可视化Py ...

  8. 分享几个 PHP 编码的最佳实践

    对于初学者而言,可能很难理解为什么某些做法更安全. 但是,以下一些技巧可能超出了 PHP 的范围. 始终使用大括号 让我们看下面的代码: if (isset($condition) && ...

  9. 小程序中内容审核功能的使用(后台使用thinkPHP5.1)

    本文包含文本和图片的检测 //接收要检测的文本内容并调用检测方法 public function textCheck(Request $request){ //内容安全识别 $data['conten ...

  10. MySQL InnoDB表的碎片量化和整理(data free能否用来衡量碎片?)

    网络上有很多MySQL表碎片整理的问题,大多数是通过demo一个表然后参考data free来进行碎片整理,这种方式对myisam引擎或者其他引擎可能有效(本人没有做详细的测试).对Innodb引擎是 ...