题目链接


Solution

分块+\(Dijkstra\).

难点在于建边,很明显 \(O(n^2)\) 建边会挂一堆 .

那么考虑一下, \(n^2\) 建边多余的是哪些东西 \(???\)

很显然是冗杂的边,即两个点在之前已经可以互达了,但是在这一次仍然又连接一遍.

所以我们对于 \(n\) 个点都开 \(sqrt(n)\) 个辅助点.代表第 \(i\) 个点可以走出 \(j\) .

辅助点之间也需要与相邻的连上一条边权为 \(1\) 的边.

然后对于 \(m\) 个点分类讨论.

  • 如果 \(p_i<sqrt(n)\)

    那么在这 \(sqrt(n)\) 里面对应连边.

  • 如果 \(p_i>sqrt(n)\)

    那么很显然他所连向的边一般不会有冗杂(大概率).所以直接暴力连边即可.

然后最后面跑一遍 \(Diskstra\) 即可.

注意距离要开 \(long~long\).


Code

#include<bits/stdc++.h>
#define il inline
#define ll long long
#define rg register
#define num(x,y) x*n+y
using namespace std;
const int maxn=30008;
const int inf=0x3f3f3f3f; struct sj{int to,next; ll w;}a[maxn*500];
int n,m,head[maxn*105],size,s,t,tmp;
int b[maxn],p[maxn];ll dis[maxn*105];
il void add(int x,int y,ll w)
{
a[++size].to=y;
a[size].next=head[x];
head[x]=size;
a[size].w=w;
} il int read()
{
char ch=getchar();int f=1,w=0;
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch<='9'&&ch>='0'){w=w*10+ch-'0';ch=getchar();}
return f*w;
} struct node{int u; ll d;
bool operator <(const node& kkk)const
{return d>kkk.d;}
}; il void Dijkstra()
{
priority_queue<node>q;
memset(dis,127,sizeof(dis));
q.push((node){s,0});
dis[s]=0;
while(!q.empty())
{
node x=q.top();q.pop();
int u=x.u;
for(int i=head[u];i;i=a[i].next)
{
int tt=a[i].to;
if(dis[tt]>dis[u]+a[i].w)
{
dis[tt]=dis[u]+a[i].w;
q.push((node){tt,dis[tt]});
}
}
}
return;
} int main()
{
n=read();m=read();
for(int i=1;i<=m;i++) b[i]=read()+1,p[i]=read();
s=b[1];t=b[2];
tmp=min((int)sqrt(n),100);
for(int i=1;i<=tmp;i++)
for(int j=1;j<=n;j++)
add(num(i,j),j,0);
for(int i=1;i<=tmp;i++)
for(int j=1;j<=n-i;j++)
add(num(i,j),num(i,j+i),1),
add(num(i,j+i),num(i,j),1);
for(int i=1;i<=m;i++)
{
if (p[i]<=tmp) add(b[i],num(p[i],b[i]),0);
else
{
for(ll j=1;b[i]+j*p[i]<=n;j++) add(b[i],b[i]+j*p[i],j);
for(ll j=1;b[i]-j*p[i]>=1;j++) add(b[i],b[i]-j*p[i],j);
}
} Dijkstra();
cout<<(dis[t]>192608173?-1:dis[t])<<endl;
return 0;
}

[APIO2015] 雅加达的摩天楼 (分块,最短路)的更多相关文章

  1. 【BZOJ4070】[Apio2015]雅加达的摩天楼 set+最短路

    [BZOJ4070][Apio2015]雅加达的摩天楼 Description 印尼首都雅加达市有 N 座摩天楼,它们排列成一条直线,我们从左到右依次将它们编号为 0 到 N−1.除了这 N 座摩天楼 ...

  2. luogu P3645 [APIO2015]雅加达的摩天楼 分块 根号分治

    LINK:雅加达的摩天楼 容易想到设\(f_{i,j}\)表示第i个\(doge\)在第j层楼的最小步数. 转移显然是bfs.值得一提的是把初始某层的\(doge\)加入队列 然后转移边权全为1不需要 ...

  3. 洛谷P3645 [APIO2015]雅加达的摩天楼(最短路+分块)

    传送门 这最短路的建图怎么和网络流一样玄学…… 一个最朴素的想法是从每一个点向它能到达的所有点连边,边权为跳的次数,然后跑最短路(然而边数是$O(n^2)$除非自创复杂度比spfa和dijkstra还 ...

  4. BZOJ 4070 [Apio2015]雅加达的摩天楼 ——分块 SPFA

    挺有趣的分块的题目. 直接暴力建边SPFA貌似是$O(nm)$的. 然后考虑分块,$\sqrt n$一下用虚拟节点辅助连边, 以上的直接暴力连边即可. 然后卡卡时间,卡卡空间. 终于在UOJ上T掉辣. ...

  5. 【题解】P3645 [APIO2015]雅加达的摩天楼(分层图最短路)

    [题解]P3645 [APIO2015]雅加达的摩天楼(分层图最短路) 感觉分层图是个很灵活的东西 直接连边的话,边数是\(O(n^2)\)的过不去 然而我们有一个优化的办法,可以建一个新图\(G=( ...

  6. BZOJ 4070:[APIO2015]雅加达的摩天楼 最短路

    4070: [Apio2015]雅加达的摩天楼 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 464  Solved: 164[Submit][Sta ...

  7. bzoj 4070 [Apio2015]雅加达的摩天楼 Dijkstra+建图

    [Apio2015]雅加达的摩天楼 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 644  Solved: 238[Submit][Status][D ...

  8. BZOJ4070 [Apio2015]雅加达的摩天楼 【分块 + 最短路】

    题目链接 BZOJ4070 题解 考虑暴力建图,将每个\(B_i\)向其能到的点连边,复杂度\(O(\sum \frac{n}{p_i})\),当\(p\)比较小时不适用 考虑优化建图,每个\(dog ...

  9. 洛谷$P3645\ [APIO2015]$雅加达的摩天楼 最短路

    正解:最短路 解题报告: 传送门$QwQ$ 考虑暴力连边,发现最多有$n^2$条边.于是考虑分块 对于长度$p_i$小于等于$\sqrt(n)$的边,建立子图$d=p_i$.说下关于子图$d$的定义? ...

随机推荐

  1. django 数据库中中文转化为韩语拼音

    1.安装模块 django-uuslug pip install django-uuslug 2.导入模块 from uuslug import slugify 3.使用模块 slugify('天龙八 ...

  2. java基础——快速排序

    今天又把以前学的快速排序拿出来回忆一下 高快省的排序算法 有没有既不浪费空间又可以快一点的排序算法呢?那就是“快速排序”啦!光听这个名字是不是就觉得很高端呢. 假设我们现在对“6 1 2 7 9 3 ...

  3. 01_1_准备ibatis环境

    01_1_准备ibatis环境 1. 搭建环境:导入相关的jar包 mysql-connector-java-5.1.5-bin.jar(mysql)或者ojdbc6.jar(oracle).ibat ...

  4. CF-1143D. The Beatles

    题意:有间隔为k的n个点在数轴上,下标为 \(1,k+1, 2*k+1,\cdots (n-1)*k+1\) 首尾相接.设起点为s,步长为L,而现在只知道s距离最近的点的距离为a,和(s+L)距离最近 ...

  5. 02等待单个线程返回WaitForSingleObject

    windows 多线程之等待线程返回 ​ 多线程编程中,有时我们需要等待某一线程完成了特定的操作之后再继续做其他事情,要实现这个目的,可以使用 Windows API 函数 WaitForSingle ...

  6. 【php】expose_php 作用

    php.ini expose_php = On expose_php = Off

  7. 初学Python01

    1.文本编辑器区别于交互模式的Python,它可以保存Python代码文件,再次打开还是存在.文件保存时要注意是.py的模式. 2.Windows系统下,应使用命令行模式打开.py 首先进入文件所在磁 ...

  8. Python学习笔记:字符串

    字符串 字符串定义:字符串可以使用一对单引号.双引号或三引号来定义,即便是单个字符也会当做字符串来处理(Python中没有字符类型,单个字符也就是只有一个字符的字符串而已). 原始字符串:字符串中反斜 ...

  9. python能干什么?

    python能干什么? 网络爬虫 爬虫,指的是从互联网采集数据的程序脚本 . 爬天爬地爬空气 ,无聊的时候爬一爬吃鸡数据.b站评论,能得出很多有意思的结论.知乎有个很有意思的问题——"利用爬 ...

  10. Aizu - 1378 Secret of Chocolate Poles (DP)

    你有三种盘子,黑薄,白薄,黑厚. 薄的盘子占1,厚的盘子占k. 有一个高度为L的桶,盘子总高度不能超出桶的总高度(可以小于等于).相同颜色的盘子不能挨着放. 问桶内装盘子的方案数. 如 L = 5,k ...