题解

每次按较长边把矩形分成两半,找一个中间轴,轴上的每个点跑一边最短路更新所有的答案

然后把矩形分成两半,递归下去

代码

#include <bits/stdc++.h>
#define enter putchar('\n')
#define space putchar(' ')
#define pii pair<int,int>
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define eps 1e-8
//#define ivorysi
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 - '0' + c;
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) out(x / 10);
putchar('0' + x % 10);
}
int N,M,Q;
int id(int x,int y) {return (x - 1) * M + y;}
pii point(int s) {return mp((s - 1) / M + 1,(s - 1) % M + 1);}
struct qry_node {
int s,t,id;
}qry[100005],tmp[100005];
struct node {
int to,next,val;
}E[200005];
int head[20005],sumE;
int R[20005],C[20005],ans[100005],dis[20005];
bool vis[20005];
priority_queue<pii > que;
void add(int u,int v,int c) {
E[++sumE].to = v;
E[sumE].next = head[u];
E[sumE].val = c;
head[u] = sumE;
}
bool in_Matrix(int x,int X1,int Y1,int X2,int Y2) {
pii p = point(x);
return p.fi >= X1 && p.fi <= X2 && p.se >= Y1 && p.se <= Y2;
}
void Dijkstra(int st,int X1,int Y1,int X2,int Y2) {
for(int i = X1 ; i <= X2 ; ++i) {
for(int j = Y1 ; j <= Y2 ; ++j) {
dis[id(i,j)] = 0x7fffffff;
vis[id(i,j)] = 0;
}
}
dis[st] = 0;
que.push(mp(-dis[st],st));
while(!que.empty()) {
pii now = que.top();que.pop();
if(vis[now.se]) continue;
vis[now.se] = 1;
int u = now.se;
for(int i = head[u] ; i ; i = E[i].next) {
int v = E[i].to;
if(!vis[v] && in_Matrix(v,X1,Y1,X2,Y2) && dis[v] > dis[u] + E[i].val) {
dis[v] = dis[u] + E[i].val;
que.push(mp(-dis[v],v));
}
}
}
}
void update(int &x,int y) {
x = min(x,y);
}
void Solve(int l,int r,int X1,int Y1,int X2,int Y2) {
if(Y1 > Y2 || X1 > X2) return;
if(Y1 == Y2 && X1 == X2) return;
if(l > r) return;
if(X2 - X1 <= Y2 - Y1) {
int mid = (Y2 + Y1) >> 1;
for(int i = X1 ; i <= X2 ; ++i) {
Dijkstra(id(i,mid),X1,Y1,X2,Y2);
for(int k = l ; k <= r ; ++k) {
update(ans[qry[k].id],dis[qry[k].s] + dis[qry[k].t]);
}
}
int p = l - 1,m = l - 1;
for(int k = l ; k <= r ; ++k) {
if(in_Matrix(qry[k].s,X1,Y1,X2,mid - 1) && in_Matrix(qry[k].t,X1,Y1,X2,mid - 1)) {
tmp[++p] = qry[k];
}
}
m = p;
for(int k = l ; k <= r ; ++k) {
if(in_Matrix(qry[k].s,X1,mid + 1,X2,Y2) && in_Matrix(qry[k].t,X1,mid + 1,X2,Y2)) {
tmp[++p] = qry[k];
}
}
for(int i = l ; i <= p ; ++i) qry[i] = tmp[i];
Solve(l,m,X1,Y1,X2,mid - 1);
Solve(m + 1,p,X1,mid + 1,X2,Y2);
}
else {
int mid = (X1 + X2) >> 1;
for(int j = Y1 ; j <= Y2 ; ++j) {
Dijkstra(id(mid,j),X1,Y1,X2,Y2);
for(int k = l ; k <= r ; ++k) {
update(ans[qry[k].id],dis[qry[k].s] + dis[qry[k].t]);
}
}
int p = l - 1,m = l - 1;
for(int k = l ; k <= r ; ++k) {
if(in_Matrix(qry[k].s,X1,Y1,mid - 1,Y2) && in_Matrix(qry[k].t,X1,Y1,mid - 1,Y2)) {
tmp[++p] = qry[k];
}
}
m = p;
for(int k = l ; k <= r ; ++k) {
if(in_Matrix(qry[k].s,mid + 1,Y1,X2,Y2) && in_Matrix(qry[k].t,mid + 1,Y1,X2,Y2)) {
tmp[++p] = qry[k];
}
}
for(int i = l ; i <= p ; ++i) qry[i] = tmp[i];
Solve(l,m,X1,Y1,mid - 1,Y2);
Solve(m + 1,p,mid + 1,Y1,X2,Y2);
}
}
void Init() {
read(N);read(M);
for(int i = 1 ; i <= N ; ++i) {
for(int j = 1 ; j <= M - 1; ++j) {
read(R[id(i,j)]);
add(id(i,j),id(i,j + 1),R[id(i,j)]);
add(id(i,j + 1),id(i,j),R[id(i,j)]);
}
}
for(int i = 1 ; i <= N - 1 ; ++i) {
for(int j = 1 ; j <= M ; ++j) {
read(C[id(i,j)]);
add(id(i,j),id(i + 1,j),C[id(i,j)]);
add(id(i + 1,j),id(i,j),C[id(i,j)]);
}
}
read(Q);
int x,y;
for(int i = 1 ; i <= Q ; ++i) {
read(x);read(y);qry[i].s = id(x,y);
read(x);read(y);qry[i].t = id(x,y);
qry[i].id = i;
ans[i] = 0x7fffffff;
if(qry[i].s == qry[i].t) ans[i] = 0;
} }
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve(1,Q,1,1,N,M);
for(int i = 1 ; i <= Q ; ++i) {out(ans[i]);enter;}
}

【LOJ】#2090. 「ZJOI2016」旅行者的更多相关文章

  1. @loj - 2090@ 「ZJOI2016」旅行者

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 小 Y 来到了一个新的城市旅行.她发现了这个城市的布局是网格状的 ...

  2. 2090. 「ZJOI2016」旅行者 分治,最短路

    2090. 「ZJOI2016」旅行者 链接 loj 思路 \((l,mid)(mid+1,r)\).考虑跨过mid的贡献. 假设选的中间那条线的点为gzy,贡献为\(dis(x,gzy)+dis(g ...

  3. 「ZJOI2016」旅行者 解题报告

    「ZJOI2016」旅行者 对网格图进行分治. 每次从中间选一列,然后枚举每个这一列的格子作为起点跑最短路,进入子矩形时把询问划分一下,有点类似整体二分 至于复杂度么,我不会阿 Code: #incl ...

  4. loj2090 「ZJOI2016」旅行者

    分治+最短路,很套路的 #include <algorithm> #include <iostream> #include <cstring> #include & ...

  5. @loj - 2093@ 「ZJOI2016」线段树

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 小 Yuuka 遇到了一个题目:有一个序列 a1,a2,..., ...

  6. @loj - 2092@ 「ZJOI2016」大森林

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 小 Y 家里有一个大森林,里面有 n 棵树,编号从 1 到 n. ...

  7. @loj - 2091@ 「ZJOI2016」小星星

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 小 Y 是一个心灵手巧的女孩子,她喜欢手工制作一些小饰品.她有 ...

  8. 「ZJOI2016」解题报告

    「ZJOI2016」解题报告 我大浙的省选题真是超级神仙--这套已经算是比较可做的了. 「ZJOI2016」旅行者 神仙分治题. 对于一个矩形,每次我们从最长边切开,最短边不会超过 \(\sqrt{n ...

  9. Loj #2192. 「SHOI2014」概率充电器

    Loj #2192. 「SHOI2014」概率充电器 题目描述 著名的电子产品品牌 SHOI 刚刚发布了引领世界潮流的下一代电子产品--概率充电器: 「采用全新纳米级加工技术,实现元件与导线能否通电完 ...

随机推荐

  1. MT【94】由参数前系数凑配系数题1

    评:根据$b,c$前系数凑配系数,也是比较常见的思路.

  2. uva11235 FrequentValues (ST表)

    既然他是非降的,那我们可以把这个序列每一位转化成到这位位置连续相同的个数,比如001111233444变成121234112123,然后一个区间内的最大值就是众数的位数.但有个问题,就是这个区间的左端 ...

  3. jQuery源码之 empty与html('')的区别

    empty: function() { var elem, i = 0; for ( ; (elem = this[i]) != null; i++ ) { // Remove element nod ...

  4. 解决VMware虚拟机网络时长中断的问题

    1. 操作环境 VMware VMware® Workstation 14 Pro Windows Win7旗舰版 2. 操作过程   VMware虚拟机在使用一段时间后,经常会出现时常断网的情况,而 ...

  5. 怎么使用 JavaScript 将网站后台的数据变化实时更新到前端

    实时这个工作现在大体有两种方法一.前端不断地向后台轮询请求数据查询的接口(不管你是用AJAX还是什么)然后将返回的数据重绘在页面上,这以前端页面为主动的方式.二.如果浏览器支持Websocket 那么 ...

  6. JavaScript中作用域回顾(避免使用全局变量)(瀑布流的实现)(scroll事件)以及Django自定义模板函数回顾

    页面显示照片样式为瀑布流: 上面的div个数可以按照自己安排进行划分.img的分布可以使用模板标签以及自定义模板函数进行排布: 自定义模板函数实现可以看,最后几列:python---django中模板 ...

  7. js实现表单提交submit(),onsubmit

    通常表单的提交有两种方式,一是直接通过html的form提交,代码如下: <form action="" method="" id="forms ...

  8. 拼接html a标签字符串,onClick传递两个字符串类型参数写法

    在参数传递过程中字符串类型的参数要有引号,我一开始拼接的完成后,没有想到要加引号,后来想到了这一问题,可是怎么拼都不对,于是就搜了很多拼接的例子,发现并没有几个能借鉴的,最后终于在一个人的博客中看到, ...

  9. 流媒体技术学习笔记之(八)海康、大华IpCamera RTSP地址和格式

    海康: rtsp://[username]:[password]@[ip]:[port]/[codec]/[channel]/[subtype]/av_stream 说明: username: 用户名 ...

  10. charles https抓包 (安卓安装证书)

    的Android APP使用的都是http请求,之后改成了https,就出现了以下情况,无法正常读取抓取的内容 下面阐述一下,正确的安装步骤,为出现类似情况的朋友提供一个参考: 1.第一步: 最后点击 ...