其实这题前前后后的思考时间加起来应该有两天之久了,dp状态,转移方式等等都还是比较好想,然而左看右看觉得spfa复杂度未免太爆炸……然后选择看了一篇题解,发现在多重优化之下,其实是可以过的……

  首先建立状态,这个应该比较明显:\(f[l][r][x][y]\) 代表合并完区间 \(l\) ~\(r\) 之后,机器人停在 \(x,y\) 处所需要的最少移动次数。转移状态即为:

  \(f[l][r][x][y] = f[l][k][x][y] + f[k + 1][r][x][y] \left ( l <= k <= r \right )\)

\(f[l][r][x][y] = f[l][r][x'][y'] + 1 \)

  其中第二个转移发生的条件是 \(x',y'\) 可以一步到达 \(x,y\)。第二个转移状态就是在之前的博客中所提及的那样:1.满足三角形不等式;2.不满足拓扑序;针对这样的转移,我们用 spfa 来优化 dp 的转移。注意在这张图中,边权均为1。在单源的最短路中,这样的图spfa可以优化为bfs, 在多源最短路中我们可以使用两个队列来进行优化。这两个队列分别存储新增的节点 & 被松弛所以要去松弛其余节点的节点。这样将节点分类之后,每一次取出队首元素权值更小的进行松弛操作。我们会发现第二个队列中节点的权值是单调的(在边权为1的图中,先访问到的节点权值更小),而第一个队列中的元素我们使用基数排序来排。(并不知道为什么要用基数排序,或许就是比较快吧?)

  然后这份代码是我抄的大佬的代码,非常感谢了。其中有一个小小的技巧:memset的时候默认赋给节点当前数据类型的最大值,相加会溢出。但对于这种没有正负要求的,我们可以利用 unsigned 自然溢出使得结果依然是最大值。(・ω<)☆ 感觉这题还是挺毒的,差点就被毒死了……

#include <bits/stdc++.h>
using namespace std;
#define maxn 505
#define maxk 400000
#define uns unsigned short
#define INF 32639
int n, W, H, ans = INF;
int mark[maxn][maxn][];
int L, R;
uns f[][][maxn][maxn];
int cnt, top, tank[maxk], S[maxk];
char Map[maxn][maxn];
bool vis[maxn][maxn];
int dxy[][] = {{,}, {,}, {,-}, {-,}}; struct node
{
int x, y;
node(int xx = , int yy = ) { x = xx, y = yy; }
}g[maxn][maxn][], pos[], q[maxk];
queue <node> q1, q2; int read()
{
int x = , k = ;
char c;
c = getchar();
while(c < '' || c > '') { if(c == '-') k = -; c = getchar(); }
while(c >= '' && c <= '') x = x * + c - '', c = getchar();
return x * k;
} void gmin(uns &x, uns y) { x = x > y ? y : x; } node dfs(int x, int y, int k)
{
if(mark[x][y][k] == cnt) return g[x][y][k] = node (-, -);
mark[x][y][k] = cnt;
if(g[x][y][k].x != && g[x][y][k].y != ) return g[x][y][k];
int pre = k;
if(Map[x][y] == 'A') k = (k + ) % ;
else if(Map[x][y] == 'C') k = (k + ) % ;
int xx = x + dxy[k][], yy = y + dxy[k][];
if(xx < || yy < || xx > H || yy > W || Map[xx][yy] == 'x') return g[x][y][pre] = node(x, y);
return g[x][y][pre] = dfs(xx, yy, k);
} void spfa()
{
memset(tank, , sizeof(tank));
for(int i = ; i <= top; i ++) tank[f[L][R][q[i].x][q[i].y]] ++;
for(int i = ; i <= INF; i ++) tank[i] += tank[i - ];
for(int i = ; i <= top; i ++) S[tank[f[L][R][q[i].x][q[i].y]] --] = i;
for(int i = ; i <= top; i ++) q1.push(q[S[i]]);
top = ;
while(!q1.empty() || !q2.empty())
{
node now;
if(q1.empty()) now = q2.front(), q2.pop();
else if(q2.empty()) now = q1.front(), q1.pop();
else
{
int x1 = q1.front().x, y1 = q1.front().y;
int x2 = q2.front().x, y2 = q2.front().y;
if(f[L][R][x1][y1] <= f[L][R][x2][y2]) now = q1.front(), q1.pop();
else now = q2.front(), q2.pop();
}
vis[now.x][now.y] = ;
for(int i = ; i < ; i ++)
{
node v = g[now.x][now.y][i];
if(v.x == - || v.y == -) continue;
if(f[L][R][v.x][v.y] > f[L][R][now.x][now.y] + )
{
f[L][R][v.x][v.y] = f[L][R][now.x][now.y] + ;
if(!vis[v.x][v.y]) vis[v.x][v.y] = , q2.push(v);
}
}
}
} int main()
{
n = read(), W = read(), H = read();
memset(f, , sizeof(f));
for(int i = ; i <= H; i ++)
{
scanf("%s", Map[i] + );
for(int j = ; j <= W; j ++)
if(Map[i][j] > '' && Map[i][j] <= '')
pos[Map[i][j] - ''] = node(i, j);
}
for(int i = ; i <= H; i ++)
for(int j = ; j <= W; j ++)
if(Map[i][j] != 'x')
for(int k = ; k < ; k ++)
++ cnt, dfs(i, j, k);
for(int i = ; i <= n; i ++)
{
vis[pos[i].x][pos[i].y] = ;
q[++ top] = pos[i];
L = R = i; f[i][i][pos[i].x][pos[i].y] = ;
spfa();
}
for(int l = , j; l <= n; l ++)
for(int i = ; (j = i + l - ) <= n; i ++)
{
for(int x = ; x <= H; x ++)
for(int y = ; y <= W; y ++)
{
for(int k = i; k < j; k ++)
gmin(f[i][j][x][y], f[i][k][x][y] + f[k + ][j][x][y]);
if(f[i][j][x][y] < INF) q[++ top] = node(x, y), vis[x][y] = ;
}
L = i, R = j; spfa();
}
unsigned short ans = INF;
for(int i = ; i <= H; i ++)
for(int j = ; j <= W; j ++)
ans = min(ans, f[][n][i][j]);
if(ans < INF) printf("%u\n", ans);
else printf("-1\n");
return ;
}

【题解】APIO2013机器人的更多相关文章

  1. bzoj3205 [Apio2013]机器人

    3205: [Apio2013]机器人 Time Limit: 15 Sec  Memory Limit: 128 MBSubmit: 953  Solved: 227[Submit][Status] ...

  2. [Bzoj3205][Apio2013]机器人(斯坦纳树)(bfs)

    3205: [Apio2013]机器人 Time Limit: 15 Sec  Memory Limit: 128 MBSubmit: 977  Solved: 230[Submit][Status] ...

  3. [APIO2013]机器人(斯坦纳树)

    题目描述 VRI(Voltron 机器人学会)的工程师建造了 n 个机器人.任意两个兼容的机 器人站在同一个格子时可以合并为一个复合机器人. 我们把机器人用 1 至 n 编号(n ≤ 9).如果两个机 ...

  4. bzoj千题计划230:bzoj3205: [Apio2013]机器人

    http://www.lydsy.com/JudgeOnline/problem.php?id=3205 历时一天,老子终于把它A了 哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈哈 因为不懂spfa ...

  5. [APIO2013]机器人[搜索、斯坦纳树]

    题意 题目链接 分析 记 g(d,x,y) 表示从 (x,y) 出发,方向为 d 到达的点,这个可以通过记忆化搜索求出,注意如果转移成环(此时向这个方向走没有意义)要特判. 记 f(l,r,x,y) ...

  6. BZOJ3205/UOJ107 [Apio2013]机器人

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...

  7. [APIO2013]机器人

    题目描述 VRI(Voltron 机器人学会)的工程师建造了 n 个机器人.任意两个兼容的机 器人站在同一个格子时可以合并为一个复合机器人. 我们把机器人用 1 至 n 编号(n ≤ 9).如果两个机 ...

  8. bzoj 3205: [Apio2013]机器人【dfs+斯坦纳树+spfa】

    第一次听说斯坦纳树这种东西 先dfs预处理出来dis[i][j][k]表示格子(i,j)向k方向转移能到哪,记忆话搜索预处理,注意如果有环的话特判一下 设f[i][j][x][y]表示复合机器人i-j ...

  9. 【BZOJ3205_洛谷3638】[APIO2013]机器人(动态规划)

    题目: 洛谷3638 分析: 卡了一天的神题--(OrzJumpmelon) 首先预处理出从点\(p\)向\(d\)方向出发最终能到达的点\(nxt[p][d]\).这个可以直接记忆化搜索解决.如果出 ...

随机推荐

  1. 富文本编辑器 summernote.js

    1.引用js  可在 https://summernote.org/ 官网下载 ,并查看详细的API  引入:summernote.js 和 summernote-zh-CN.js 以及样式文件:su ...

  2. django 面试题

    面试题1:migrate怎么判断哪些迁移脚本需要执行: 他会将代码中的迁移脚本和数据库中django_migrations中的迁移脚本进行对比,如果发现数据库中,没有这个迁移脚本,那么就会执行这个迁移 ...

  3. 关于使用array_rand随机取出数组的值

    代码如下 <?php echo "<meta charset='utf-8'/>";//选择解码方式,防止乱码现象 $a = array("abc&qu ...

  4. elasticsearch 5.x 系列之六 文档索引,更新,查询,删除流程

    一.elasticsearch index 索引流程 步骤: 客户端向Node1 发送索引文档请求 Node1 根据文档ID(_id字段)计算出该文档应该属于shard0,然后请求路由到Node3的P ...

  5. 基于pyecharts的IT各行业薪资展示

    我们的项目是一个信息采集系统,采集的是51job招聘网站,我爬取了Python,Java,C++,PHP还有北京各地区的职位数量,以及经验要求,和学历要求等等. 网页头; <!DOCTYPE h ...

  6. JSOI2018 R1 & 九省联考2018 滚粗记

    在NOIP与PKUWC相继滚粗后,rp守恒定律似乎终于开始起作用了…… (尽管Day2依然滚粗?) Day1: 本着前40min不写代码的准则,先把三道题大致过了一遍,似乎都比较喜闻乐见? T1:对抗 ...

  7. centos7下安装elasticSearch错误总结(单节点模式)

    1.首先确定你安装了jdk,版本需要1.8以上 2.上传elasticsearchjar包,只需配置一个文件即可 修改配置文件config/elasticsearch.yml    network.h ...

  8. 利尔达NB-IOT的PSM和eDRX低功耗模式笔记

    1. NB-IOT的技术优势,广覆盖,NB-IOT与GPRS和LTE相比较,最大链路预算提升了20dB,相当于提升了100倍,即使在地车车库.地下室.地下管道等普通无线网络信号难以到达的地方也容易覆盖 ...

  9. https 通信流程和Charles 抓包原理

    1. https 通信流程 ①客户端的浏览器向服务器传送客户端SSL 协议的版本号,加密算法的种类,产生的随机数,以及其他服务器和客户端之间通讯所需要的各种信息.②服务器向客户端传送SSL 协议的版本 ...

  10. html5特效库

    Swiper是纯javascript打造的滑动特效插件,面向手机.平板电脑等移动终端.能实现触屏焦点图.触屏Tab切换.触屏多图切换等常用效果. delaunay.js是一款能在图片中形成无数个三角形 ...