题目链接

\(Description\)

有一个\(n\times n\)的正方形网格,上面有若干障碍点。\(q\)次询问,每次询问把一个正方形箱子从\((x1,y1)\)推到\((x2,y2)\) ,不能接触障碍点且不超出边界,箱子的边长最大能是多少。

\(n\leq1000,q\leq3\times10^5\)。

\(Solution\)

我们可以直接将非障碍点看做顶点,向周围四个非障碍点连边。这样是要求一棵最大生成树。

从每个障碍点做八连通BFS,可以得到通过每个非障碍点位置的最大直径。

然后按最大直径从大到小依次加入点并标为访问过,每次与周围之前被访问过的点连边(也就是枚举权值小的与大的连边,优先连长边,和Kruskal一样)。

其实就是Kruskal重构树。连边时因为权值在点上,所以不需要新建点来表示边了,直接将权值大的点连到权值小的点作为其子节点即可。

两点间路径询问就是直接求LCA的权值了。

复杂度\(O(n^2+q\log n^2)\)。

注意初始化dis为INF(否则无障碍点时所有dis都为0)。

//126688kb	9332ms
#include <queue>
#include <cstdio>
#include <cctype>
#include <vector>
#include <cstring>
#include <algorithm>
#define mp std::make_pair
#define pr std::pair<int,int>
//#define gc() getchar()
#define MAXIN 1000000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int N=1003,M=1e6+5;
const int dd[]={1,0,-1,0,1};
const int dx[]={0,0,-1,-1,-1,1,1,1};
const int dy[]={-1,1,-1,0,1,-1,0,1}; int n,id[N][N],A[M],fa[M];
std::vector<int> v[N];
char IN[MAXIN],*SS=IN,*TT=IN; namespace HLD
{
int Enum,H[M],nxt[M],to[M],fa[M],sz[M],son[M],dep[M],top[M]; inline void AE(int u,int v)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
}
inline int LCA(int u,int v)
{
while(top[u]!=top[v]) dep[top[u]]>dep[top[v]]?u=fa[top[u]]:v=fa[top[v]];
return dep[u]<dep[v]?u:v;
}
void DFS1(int x)
{
int mx=0; sz[x]=1;
for(int i=H[x],v; i; i=nxt[i])
if((v=to[i])!=fa[x])
{
fa[v]=x, dep[v]=dep[x]+1, DFS1(v), sz[x]+=sz[v];
if(sz[v]>mx) mx=sz[v], son[x]=v;
}
}
void DFS2(int x,int tp)
{
top[x]=tp;
if(son[x])
{
DFS2(son[x],tp);
for(int i=H[x]; i; i=nxt[i])
if(to[i]!=fa[x]&&to[i]!=son[x]) DFS2(to[i],to[i]);
}
}
} inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
int Find(int x)
{
return x==fa[x]?x:fa[x]=Find(fa[x]);
}
void Merge(int x,int y)//x<-y
{
y=Find(y);
if(x!=y) fa[y]=x, HLD::AE(x,y);
}
void BFS()
{
static int dis[N][N];
static bool vis[N][N];
static std::queue<pr> q; memset(dis,0x3f,sizeof dis);
for(int i=1; i<=n; ++i)
for(int j=1; j<=n; ++j)
if(!id[i][j]) q.push(mp(i,j)),vis[i][j]=1,dis[i][j]=0;
while(!q.empty())
{
int x=q.front().first,y=q.front().second;
q.pop();
for(int i=0,xn,yn,d=dis[x][y]; i<8; ++i)
if((xn=x+dx[i])>0&&xn<=n&&(yn=y+dy[i])>0&&yn<=n&&!vis[xn][yn])
vis[xn][yn]=1, dis[xn][yn]=d+1, q.push(mp(xn,yn));
}
for(int i=1; i<=n; ++i)
for(int j=1,ds,t; j<=n; ++j)
if(t=id[i][j])
{
ds=std::min(dis[i][j],std::min(std::min(i,n-i+1),std::min(j,n-j+1)));
ds=(ds<<1)-1;//maxd = 2r-1
A[t]=ds, v[ds].push_back(t);
}
} int main()
{
static pr ref[M];
static bool vis[N][N]; n=read(); int cnt=0;
for(int i=1; i<=n; ++i)
{
register char c=gc(); while(c!='.'&&c!='#') c=gc();
if(c=='.') id[i][1]=++cnt, ref[cnt]=mp(i,1);
for(int j=2; j<=n; ++j) if(gc()=='.') id[i][j]=++cnt, ref[cnt]=mp(i,j);
} BFS();
for(int i=1; i<=cnt; ++i) fa[i]=i;
for(int i=n; i; --i)
{
const std::vector<int> &vec=v[i];
for(int j=0,l=vec.size(); j<l; ++j)
{
int t=vec[j],x=ref[t].first,y=ref[t].second;
vis[x][y]=1;
for(int i=0,xn,yn; i<4; ++i)
if((xn=x+dd[i])>0&&xn<=n&&(yn=y+dd[i+1])>0&&yn<=n&&vis[xn][yn])
Merge(t,id[xn][yn]);
}
}
for(int i=1; i<=cnt; ++i) if(fa[i]==i) Merge(0,i);
HLD::DFS1(0), HLD::DFS2(0,0); for(int q=read(),p1,p2,x1,y1,x2,y2; q--; )
{
p1=id[read()][read()], p2=id[read()][read()];
if(!p1||!p2) puts("0");
else printf("%d\n",A[HLD::LCA(p1,p2)]);
}
return 0;
}

BZOJ.4793.[CERC2016]Hangar Hurdles(Kruskal重构树 BFS)的更多相关文章

  1. Gym - 101173H Hangar Hurdles (kruskal重构树/最小生成树+LCA)

    题目大意:给出一个n*n的矩阵,有一些点是障碍,给出Q组询问,每组询问求两点间能通过的最大正方形宽度. 首先需要求出以每个点(i,j)为中心的最大正方形宽度mxl[i][j],可以用二维前缀和+二分或 ...

  2. BZOJ 3551: [ONTAK2010]Peaks加强版 [Kruskal重构树 dfs序 主席树]

    3551: [ONTAK2010]Peaks加强版 题意:带权图,多组询问与一个点通过边权\(\le lim\)的边连通的点中点权k大值,强制在线 PoPoQQQ大爷题解传送门 说一下感受: 容易发现 ...

  3. BZOJ 4242: 水壶(Kruskal重构树 + Bfs)

    题意 一块 \(h ∗ w\) 的区域,存在障碍.空地.\(n\) 个建筑,从一个建筑到另一个建筑的花费为:路径上最长的连续空地的长度. \(q\) 次询问:从建筑 \(s_i\) 到 \(t_i\) ...

  4. 【BZOJ】3732: Network【Kruskal重构树】

    3732: Network Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2812  Solved: 1363[Submit][Status][Dis ...

  5. BZOJ 3551: [ONTAK2010]Peaks加强版 Kruskal重构树+dfs序+主席树+倍增

    建出来 $Kruskal$ 重构树. 将询问点向上跳到深度最小,且合法的节点上. 那么,得益于重构树优美的性质,这个最终跳到的点为根的所有子节点都可以与询问点互达. 对于子树中求点权第 $k$ 大的问 ...

  6. BZOJ.3551.[ONTAK2010]Peaks加强版(Kruskal重构树 主席树)

    题目链接 \(Description\) 有n个座山,其高度为hi.有m条带权双向边连接某些山.多次询问,每次询问从v出发 只经过边权<=x的边 所能到达的山中,第K高的是多少. 强制在线. \ ...

  7. [bzoj 3732] Network (Kruskal重构树)

    kruskal重构树 Description 给你N个点的无向图 (1 <= N <= 15,000),记为:1-N. 图中有M条边 (1 <= M <= 30,000) ,第 ...

  8. 【BZOJ 3732】 Network Kruskal重构树+倍增LCA

    Kruskal重构树裸题, Sunshine互测的A题就是Kruskal重构树,我通过互测了解到了这个神奇的东西... 理解起来应该没什么难度吧,但是我的Peaks连WA,,, 省选估计要滚粗了TwT ...

  9. bzoj 3551 kruskal重构树dfs序上的主席树

    强制在线 kruskal重构树,每两点间的最大边权即为其lca的点权. 倍增找,dfs序对应区间搞主席树 #include<cstdio> #include<cstring> ...

随机推荐

  1. shell-检测服务是否运行,并记日志

    目的:每隔*分钟检测服务是否运行:若运行中,则记录执行的进程名称:若不运行,记录当前时间 shell: #!/bin/bash date=`date +%Y%m%d` log=/home/mono_$ ...

  2. Linux 入门记录:十二、Linux 权限机制【转】

    转自:https://www.cnblogs.com/mingc/p/7591287.html 一.权限 权限是操作系统用来限制资源访问的机制,权限一般分为读.写.执行. 系统中每个文件都拥有特定的权 ...

  3. Linux驱动开发必看详解神秘内核(完全转载)

    Linux驱动开发必看详解神秘内核 完全转载-链接:http://blog.chinaunix.net/uid-21356596-id-1827434.html   IT168 技术文档]在开始步入L ...

  4. HTML学习笔记01-HTML简介

    主要是为了做接口测试,试着自己写爬虫,所以学习一下HTML一些基础的东西,方便用来解析网页.学习内容主要来自菜鸟教程的HTML教程,W3school的HTML 超文本标记语言(英语:HyperText ...

  5. android 知识点汇总

    1.activity 它是 android 应用程序的基本功能单元.一个Activity是一个应用程序组件,提供一个屏幕,用户可以用来交互为了完成某项任务,例如拨号.拍照.Activity 本身是没有 ...

  6. [Vue warn]: Do not use built-in or reserved HTML elements as component id: header

    因为header在HTML5里面是个原生的标签,所以在开发的时候会提示错误,解决方法:修改components里面左边的header

  7. [PHP] 链表数据结构(单链表)

    链表:是一个有序的列表,但是它在内存中是分散存储的,使用链表可以解决类似约瑟夫问题,排序问题,搜索问题,广义表 单向链表,双向链表,环形链表 PHP的底层是C,当一个程序运行时,内存分成五个区(堆区, ...

  8. Python 列表推导、迭代器与生成器

    1.列表推导 1 2 3 4 5 6 7 8 9 10 11 numbers = [i for i in range(10) if i % 2 == 0] print(numbers)   seq = ...

  9. python作业(day1)

    1.输入用户名密码 ,认证成功后显示欢迎信息 ,输错三次后锁定 ,锁定后再次登录时直接提示已锁定 2.多级菜单 三级菜单 可依次进入各子菜单 所需新知识点:列表.字典

  10. IDA Pro的patch插件 KeyPatch

    本来这个是没什么可写的,但是安这个插件的时候真是气到爆炸,安装文档写的不明不白,几万行的代码都写了就差那么点时间写个几十字的详细说明吗? 1.下载keypatch.py放到\IDA\plugins里 ...