「CERC2017」Donut Drone
题目链接
题目翻译:
你正在模拟无人机探索一个不稳定的环状行星的过程。技术上说,无人机正在穿过一个环形网格———一个在两维上都首尾环绕在一起的矩形网格。格子的行号从上到下依次编号为\(1\)到\(r\),列号
从上到下依次编号为\(1\)到\(c\)。每个格子还有一个海拔——这是个正数。
无人机一开始位于第一行第一列的格子。每一步,无人机会考虑这样三个格子:右边、右上方、右下方(注意这个网格首尾相接)。无人机会飞到它们之中海拔最高的一个格子。模拟的过程中,共有两种可能的操作:
\(move\) \(k\)无人机移动k步
\(change\) \(a\) \(b\) \(e\)第\(a\)行第\(b\)列的格子海拔修改为\(e\)。
在每次\(move\)操作后,你都需要立刻找到无人机的位置。你可以认为,每次移动的三个目标位置海拔互不相同,因此每一步移动都是良定义的。
solution
每次移动的步数\(k\)很大,所以可以联想到倍增跳法,即每次跳\(2^k\)步。
但步数不太好维护,我们考虑将移动\(k\)步转化为先暴力跳到第一列,再从第一列出发跳若干圈,最后在暴力跳剩下的不到一圈的步数。显然2段暴力跳的复杂度仅为\(O(c)\),可以接受
那么问题转化为了要维护从第一列的每一个位置出发跳\(2^k\)圈后所处的位置。本题在\(y轴\)上的跳法不确定,但在\(x\)轴上一直是每次向右移动一格,于是可以在列上建立线段树,每个节点(设对应的列为\(l,r\))维护第\(l\)列上每个位置,跳到第\(r+1\)列后所处的位置,就可以
t[p][i]=t[p<<1|1][t[p<<1][i]]
实现转移(其中\(t[p][i]\)表示处在第\(i\)行的点,跳过\(p\)对应的这段区间后所处的位置,在代码中是\(T[p].t[i]\))
于是根节点维护的答案就是第一列中的每个位置跳一圈后所处的位置,接下来的只不过是倍增基本套路。
每次对\((x,y)\)的修改只会影响第\(x-1\)列维护的答案,等于是线段树中的单点修改,复杂度\(O(rlog(c))\)
code
#include<bits/stdc++.h>
using namespace std;
const int N=5010;
int a[N][N],ans[N][N],R,C,m,to[N][32];
inline int read(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*f;
}
struct SGT{
struct tree{
int t[N];
}T[N<<2];
#define lc (p<<1)
#define rc (p<<1|1)
#define mid (l+r>>1)
inline void pushup(int p){
for(int i=1;i<=R;++i)
T[p].t[i]=T[rc].t[T[lc].t[i]];
}
inline void build(int p,int l,int r){
if(l==r){
for(int i=1;i<=R;++i) T[p].t[i]=ans[l][i];
return ;
}
build(lc,l,mid);
build(rc,mid+1,r);
pushup(p);
}
inline void update(int p,int l,int r,int x){
if(l==r){
for(int i=1;i<=R;++i) T[p].t[i]=ans[l][i];
return ;
}
if(x<=mid) update(lc,l,mid,x);
else update(rc,mid+1,r,x);
pushup(p);
}
#undef lc
#undef rc
#undef mid
}T;
inline void work(int &x,int &y){
int yy=y==C?1:y+1,x1=x>1?x-1:R,x2=x,x3=x==R?1:x+1;
int ans=a[x1][yy],pos=x1;
if(a[x2][yy]>ans) ans=a[x2][yy],pos=x2;
if(a[x3][yy]>ans) ans=a[x3][yy],pos=x3;
x=pos;y=yy;
}
int main(){
R=read();C=read();
for(int i=1;i<=R;++i)
for(int j=1;j<=C;++j)
a[i][j]=read();
for(int i=1;i<=R;++i){
for(int j=1;j<=C;++j){
int x=i,y=j;work(x,y);
ans[j][i]=x;
}
}
T.build(1,1,C);
for(int i=1;i<=R;++i) to[i][0]=T.T[1].t[i];
for(int j=1;j<=30;++j)
for(int i=1;i<=R;++i)
to[i][j]=to[to[i][j-1]][j-1];
m=read();
int nx=1,ny=1;
while(m--){
char s[10];scanf("%s",s);
if(s[0]=='c'){
int x=read(),y=read(),e=read();
a[x][y]=e;y=y>1?y-1:C;
for(int i=1;i<=R;++i){
int t1=i,t2=y;
work(t1,t2);
ans[y][i]=t1;
}
T.update(1,1,C,y);
for(int i=1;i<=R;++i) to[i][0]=T.T[1].t[i];
for(int j=1;j<=30;++j)
for(int i=1;i<=R;++i)
to[i][j]=to[to[i][j-1]][j-1];
}
if(s[0]=='m'){
int k=read();
while(k&&ny!=1) work(nx,ny),k--;
int circle=k/C;k=k%C;
for(int i=30;i>=0;--i) if(circle&(1<<i)) circle^=(1<<i),nx=to[nx][i];
while(k--) work(nx,ny);
printf("%d %d\n",nx,ny);
}
}
return 0;
}
「CERC2017」Donut Drone的更多相关文章
- 「译」JUnit 5 系列:条件测试
原文地址:http://blog.codefx.org/libraries/junit-5-conditions/ 原文日期:08, May, 2016 译文首发:Linesh 的博客:「译」JUni ...
- 「译」JUnit 5 系列:扩展模型(Extension Model)
原文地址:http://blog.codefx.org/design/architecture/junit-5-extension-model/ 原文日期:11, Apr, 2016 译文首发:Lin ...
- JavaScript OOP 之「创建对象」
工厂模式 工厂模式是软件工程领域一种广为人知的设计模式,这种模式抽象了创建具体对象的过程.工厂模式虽然解决了创建多个相似对象的问题,但却没有解决对象识别的问题. function createPers ...
- 「C++」理解智能指针
维基百科上面对于「智能指针」是这样描述的: 智能指针(英语:Smart pointer)是一种抽象的数据类型.在程序设计中,它通常是经由类型模板(class template)来实做,借由模板(tem ...
- 「JavaScript」四种跨域方式详解
超详细并且带 Demo 的 JavaScript 跨域指南来了! 本文基于你了解 JavaScript 的同源策略,并且了解使用跨域跨域的理由. 1. JSONP 首先要介绍的跨域方法必然是 JSON ...
- 「2014-5-31」Z-Stack - Modification of Zigbee Device Object for better network access management
写一份赏心悦目的工程文档,是很困难的事情.若想写得完善,不仅得用对工具(use the right tools),注重文笔,还得投入大把时间,真心是一件难度颇高的事情.但,若是真写好了,也是善莫大焉: ...
- 「2014-3-18」multi-pattern string match using aho-corasick
我是擅(倾)长(向)把一篇文章写成杂文的.毕竟,写博客记录生活点滴,比不得发 paper,要求字斟句酌八股结构到位:风格偏杂文一点,也是没人拒稿的.这么说来,arxiv 就好比是 paper 世界的博 ...
- 「2014-3-17」C pointer again …
记录一个比较基础的东东-- C 语言的指针,一直让人又爱又恨,爱它的人觉得它既灵活又强大,恨它的人觉得它太过于灵活太过于强大以至于容易将人绕晕.最早接触 C 语言,还是在刚进入大学的时候,算起来有好些 ...
- 「2014-3-13」Javascript Engine, Java VM, Python interpreter, PyPy – a glance
提要: url anchor (ajax) => javascript engine (1~4 articles) => java VM vs. python interpreter =& ...
随机推荐
- 没事学些KVM(三)虚拟机基础管理
创建完成虚拟机后,需要对虚拟机进行基础管理学习 virsh list #查看虚拟机列表 改命令只能查看正在运行或挂起的虚拟机 如果需要查看所有的虚拟机需要添加--all 参数 virsh start ...
- deployer服务器端的配置
虽然说是无侵入的部署工具,但是还是需要我们来做一些微小的配置的,因为出于安全考虑,我们一般不会开发 root 用户的 SSH 登录,而是使用其它用户,比如 Ubuntu 默认的 ubuntu 用户.我 ...
- 脑桥Brain-Pons
date: 2014-02-01 15:30:11 updated: 2014-02-01 15:30:11 [一] "2025.7.3.Brain-Pons?Expeiment?Under ...
- CentOS7 Nginx-1.10.3编译安装
cat > nginx.sh <<EOF #停止apache,避免抢占端口号 systemctl stop httpd #创建nginx运行账户,非登录用户,不创建家目录 usera ...
- 计算机CPU是怎么认识代码的?
先说一下半导体,啥叫半导体?就是介于导体和绝缘体中间的一种东西,比如二极管. 电流可以从A端流向C端,但反过来则不行.你可以把它理解成一种防止电流逆流的东西. 当C端10V,A端0V,二极管可以视 ...
- 当浏览器窗口大小发生变化时,重新绘制JsPlumb中的线条、端点
1 window.addEventListener('resize', () => { 2 this.plumbInstance.repaintEverything() 3 }) 参考文章:ht ...
- 20200726_java爬虫_使用HttpClient模拟浏览器发送请求
浏览器获取数据: 打开浏览器 ==> 输入网址 ==> 回车查询 ==> 返回结果 ==> 浏览器显示结果数据 HttpClient获取数据: 创建HttpClient ==& ...
- HTTPDNS开源 Android SDK,赋能更多开发者参与共建
为赋能更多开发者参与共建,阿里云HTTPDNS开源 Android SDK,iOS SDK也在做开源准备,不久也将开放给开发者.HTTPDNS是阿里云移动研发平台面向多端应用(移动端APP,PC客户端 ...
- Java并发原理层面:ReentrantLock中lock()、unlock()全解析
一.前言 Java线程同步两种方式,synchronized关键字和Lock锁机制,其中,AQS队列就是Lock锁实现公平加锁的底层支持. 二.AQS源码对于lock.lock()的实现 2.1 AQ ...
- 十个Pycharm快捷键——提升效率
一些比较实用的Pycharm的快捷键,提升编写开发效率. 1.解除语法限制 默认情况下,Pycharm会对代码进行检查,包括但不仅限于代码是否有语法错误,是否符合PEP8规范. 如命名检查,如下图 变 ...