loj2090. 「ZJOI2016」旅行者
时间:2019-08-24
本文章向大家介绍loj2090. 「ZJOI2016」旅行者,主要包括loj2090. 「ZJOI2016」旅行者使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
loj2090. 「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;
}
原文地址:https://www.cnblogs.com/dsrdsr/p/11405996.html
- Logback 简明使用手册
- Python编写渗透工具学习笔记一 | 0x06 Zip包破解程序
- Java 并发排序
- Python编写渗透工具学习笔记一 | 0x07 Python实现键盘记录器
- Python编写渗透工具学习笔记一 | 0x08字典生成程序
- 漫步VR——Unity语音聊天室开发
- Python编写渗透工具学习笔记一 | 0x01 目录扫描程序
- JDK容器学习之HashMap (一) : 底层存储结构分析
- Python编写渗透工具学习笔记一 | 0x02实现一个反弹shell
- 动手实现MVC: 1. Java 扫描并加载包路径下class文件
- 动手实现MVC: 2. bean加载, IoC依赖注入
- Python编写渗透工具学习笔记一 | 0x04 nmap实现端口扫描(准确性更高)
- spring-boot & ffmpeg 搭建一个音频转码服务
- java 实现二维码生成工具类
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 推荐 3 款超好用的 Docker 图形化管理工具
- python标准库之glob介绍
- Python 为什么只需一条语句“a,b=b,a”,就能直接交换两个变量?
- 使用List中的remove方法遇到的坑,不信你没有踩过!
- python opencv 图像尺寸变换
- OpenCv保存图像
- 机器学习|KNN
- docker 查看容器日志
- consul配置ACL
- CentOS7.5更改python版本后及yum不能用的解决办法,非编译!
- nginx 配置websocket
- PostgreSQL10分区表性能研究报告
- linux文件目录管理基本命令总结
- Linux中的硬链接与软链接?
- linux中实现灵活的特殊权限管理ACL,以及ACL权限判断顺序!