2090. 「ZJOI2016」旅行者 分治,最短路
2090. 「ZJOI2016」旅行者
链接
思路
\((l,mid)(mid+1,r)\).考虑跨过mid的贡献。
假设选的中间那条线的点为gzy,贡献为\(dis(x,gzy)+dis(gzy,y)\)
那就计算n遍最短路,一次分治为\(n^2mlog{nm}\)
设S=n*m.矩阵的长度是不定的,每次取最长的边进行分治是最好的,n最坏为\(\sqrt{n}\)。
\(f(n)=2*f(\frac{n}{2})+S\sqrt{S}logS。所以总的复杂度就是\)\(S\sqrt{S}logS\)
都在同侧的也需要跨一跨
代码
#include <bits/stdc++.h>
#define FOR(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
const int _=1e5+7,INF=0x3f3f3f3f;
int read() {
int x=0,f=1;char s=getchar();
for(;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
for(;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
return x*f;
}
int n,m,q,ans[_],vis[_];
struct node {
int x,y,X,Y,u,v,id;
}Q[_],tmp[_];
bool cmp(node a,node b) {return a.id<b.id;}
struct edge {int v,nxt,q;}e[_];
int head[_],tot;
void add(int u,int v,int q) {
e[++tot].v=v;
e[tot].q=q;
e[tot].nxt=head[u];
head[u]=tot;
}
int id(int x,int y) {return (x-1)*m+y;}
struct T_T {
int u,val;
T_T(int a=0,int b=0) {u=a,val=b;}
bool operator < (const T_T &b) const {return val>b.val;}
};
int dis[_];
void dij(int S) {
dis[S]=0;
priority_queue<T_T> q;
q.push(T_T(S,0));
while(!q.empty()) {
T_T u=q.top();q.pop();
if(dis[u.u]!=u.val) continue;
for(int i=head[u.u];i;i=e[i].nxt) {
int v=e[i].v;
if(vis[v]&&dis[v]>u.val+e[i].q) {
dis[v]=u.val+e[i].q;
q.push(T_T(v,dis[v]));
}
}
}
}
void solve(int x,int y,int X,int Y,int l,int r) {
if(l>r) return;
if(x==X&&y==Y) {
for(int i=l;i<=r;++i) ans[Q[i].id]=0;
return;
}
if(Y-y>X-x) {
int mid=(Y+y)>>1;
FOR(i,x,X) {
FOR(j,x,X) FOR(k,y,Y) dis[id(j,k)]=INF;
dij(id(i,mid));
FOR(j,l,r)
ans[Q[j].id]=min(dis[Q[j].u]+dis[Q[j].v],ans[Q[j].id]);
}
FOR(i,x,X) vis[id(i,mid)]=0;
int p=l,q=r;
FOR(i,l,r) {
if(max(Q[i].y,Q[i].Y)<=mid) tmp[p++]=Q[i];
if(min(Q[i].y,Q[i].Y)>mid) tmp[q--]=Q[i];
}
FOR(i,l,r) Q[i]=tmp[i];
solve(x,y,X,mid,l,p-1);
solve(x,mid+1,X,Y,q+1,r);
} else {
int mid=(X+x)>>1;
FOR(i,y,Y) {
FOR(j,x,X) FOR(k,y,Y) dis[id(j,k)]=INF;
dij(id(mid,i));
FOR(j,l,r)
ans[Q[j].id]=min(dis[Q[j].u]+dis[Q[j].v],ans[Q[j].id]);
}
FOR(i,y,Y) vis[id(mid,i)]=0;
int p=l,q=r;
FOR(i,l,r) {
if(max(Q[i].x,Q[i].X)<=mid) tmp[p++]=Q[i];
if(min(Q[i].x,Q[i].X)>mid) tmp[q--]=Q[i];
}
FOR(i,l,r) Q[i]=tmp[i];
solve(x,y,mid,Y,l,p-1);
solve(mid+1,y,X,Y,q+1,r);
}
}
int main() {
n=read(),m=read();
FOR(i,1,n) FOR(j,1,m-1) {
int val=read();
add(id(i,j),id(i,j+1),val);
add(id(i,j+1),id(i,j),val);
}
FOR(i,1,n-1) FOR(j,1,m) {
int val=read();
add(id(i,j),id(i+1,j),val);
add(id(i+1,j),id(i,j),val);
}
FOR(i,1,n) FOR(j,1,m) vis[id(i,j)]=1;
q=read();
FOR(i,1,q) {
Q[i].x=read(),Q[i].y=read(),
Q[i].u=id(Q[i].x,Q[i].y);
Q[i].X=read(),Q[i].Y=read(),
Q[i].v=id(Q[i].X,Q[i].Y);
Q[i].id=i,ans[i]=INF;
}
solve(1,1,n,m,1,q);
FOR(i,1,q) printf("%d\n",ans[i]);
return 0;
}
2090. 「ZJOI2016」旅行者 分治,最短路的更多相关文章
- @loj - 2090@ 「ZJOI2016」旅行者
目录 @description@ @solution@ @accepted code@ @details@ @description@ 小 Y 来到了一个新的城市旅行.她发现了这个城市的布局是网格状的 ...
- 【LOJ】#2090. 「ZJOI2016」旅行者
题解 每次按较长边把矩形分成两半,找一个中间轴,轴上的每个点跑一边最短路更新所有的答案 然后把矩形分成两半,递归下去 代码 #include <bits/stdc++.h> #define ...
- 「ZJOI2016」旅行者 解题报告
「ZJOI2016」旅行者 对网格图进行分治. 每次从中间选一列,然后枚举每个这一列的格子作为起点跑最短路,进入子矩形时把询问划分一下,有点类似整体二分 至于复杂度么,我不会阿 Code: #incl ...
- loj2090 「ZJOI2016」旅行者
分治+最短路,很套路的 #include <algorithm> #include <iostream> #include <cstring> #include & ...
- 「ZJOI2016」解题报告
「ZJOI2016」解题报告 我大浙的省选题真是超级神仙--这套已经算是比较可做的了. 「ZJOI2016」旅行者 神仙分治题. 对于一个矩形,每次我们从最长边切开,最短边不会超过 \(\sqrt{n ...
- 「ZJOI2016」大森林 解题报告
「ZJOI2016」大森林 神仙题... 很显然线段树搞不了 考虑离线操作 我们只搞一颗树,从位置1一直往后移动,然后维护它的形态试试 显然操作0,1都可以拆成差分的形式,就是加入和删除 因为保证了操 ...
- 「bzoj1003」「ZJOI2006」物流运输 最短路+区间dp
「bzoj1003」「ZJOI2006」物流运输---------------------------------------------------------------------------- ...
- [BZOJ4456] [Zjoi2016]旅行者 分治+最短路
4456: [Zjoi2016]旅行者 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 777 Solved: 439[Submit][Status] ...
- 【BZOJ4456】[Zjoi2016]旅行者 分治+最短路
[BZOJ4456][Zjoi2016]旅行者 Description 小Y来到了一个新的城市旅行.她发现了这个城市的布局是网格状的,也就是有n条从东到西的道路和m条从南到北的道路,这些道路两两相交形 ...
随机推荐
- Oracle Hint用法整理笔记
目录 1./+ result_cache / 2./+ connect_by_filtering / 3./+ no_unnset / 4./+ index(表别名 索引名) / 5./+ INDEX ...
- python做一个简易图片下载工具
代码有点乱,先这样 # -*- coding:utf-8 -*- #__author__ :kusy #__content__:文件说明 #__date__:2018/11/01 11:01 impo ...
- 利用SQL计算两个地理坐标(经纬度)之间的地表距离
两个地理坐标(经纬度)地表距离计算公式: 公式解释如下: Long1,Lat1表示A点经纬度,Long2,Lat2表示B点经纬度: a=Lat1–Lat2 为两点纬度之差,b=Long1-Long2为 ...
- Ubuntu19 安装 pylearn2
环境: /etc/issue # Ubuntu 19.10 \n \l python -V # Python 其中,python环境是我自己建立的虚拟 venv 方便测试( 相当于你的python 实 ...
- js 固定div 不随着滚动条滚动
css .fixed { position: fixed; top:; } javascript function my$(id) { return document.getElementById(i ...
- leetcode求峰值,js实现
原题: 最开始是照着提示的思路进行,中规中矩,用时64ms /** * @param {number[]} nums * @return {number} */var findPeakElement ...
- delegate里的Invoke和BeginInvoke
Invoke和BeginInvoke都是调用委托实体的方法,前者是同步调用,即它运行在主线程上,当Invode处理时间长时,会出现阻塞的情况,而BeginInvod是异步操作,它会从新开启一个线程,所 ...
- pom.xml管理jar包——安全性框架配置文件
<properties> <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding> &l ...
- web的应用模式
在开发web应用中,有两种模式: 1.前后端分离. 2.前后端不分离. 一.前后端不分离 在前后端不分离的应用模式中,前端页面看到的效果是有后端控制的,由后端渲染页面或重定向,也就是后端需要控制前端的 ...
- Extjs 树菜单的自动展开数据的请求
今天在做extjs开发的时候,在树菜单上遇到了一个坑,也许是我刚接触extjs 不熟的缘故 问题描述:后台设置的树自动展开,但是在前端总是只显示一条数据,但是数据确实都请求到了. 经过几个小时不屑的努 ...