Arc065_E Manhattan Compass
平面上有$N$个点$(X_i\space, Y_i)$,定义$D(a,b)=|X_a-X_b|+|Y_a-Y_b|$。
如果你当前在$(p,q)$,这个无序二元组(即$(p,q)$和$(q,p)$被认为是相同的),并且存在$r$使得$D(p, q)=D(p, r)$,那么你就可以到达$(p, r)$这个无序二元组。
一开始你在$(a,b)$,这个无序二元组,求你能到达多少不同的无序二元组。
题解
这道题其实没有任何大的思维跨越,大体思想就是先求出$(a,b)$的曼哈顿距离,然后$bfs$看可以连出去多少条边即可。
然而这样朴素的做法是$O(N^2)$的,想办法优化这一过程。
显然枚举出发点不能进行深一步的优化,所以我们考虑优化枚举到达点的方法。
不难发现,从固定曼哈顿距离确定每一个点出发到达的点一定在一个斜$45$度的正方形上,可以对于每一条正反斜$45$度的直线开一个$vector$按$x$升序存下内部的点,这样就可以在$vector$上进行二分。
至于如何将这些点依次加入队列,可以直接开一个并查集维护$vector$中右侧最近的没有被删除的点,由于每个点只会被删除一次,所以总复杂度$O(N)$
考虑到每次在$vector$上找点还需要二分,所以常数较大,但复杂度仍是$O(N\log N)$
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
#include<cstdio>
#include<cmath>
#include<map>
#define LL long long
#define M 100020
using namespace std;
int read(){
int nm=0,fh=1; char cw=getchar();
for(;!isdigit(cw);cw=getchar()) if(cw=='-') fh=-fh;
for(;isdigit(cw);cw=getchar()) nm=nm*10+(cw-'0');
return nm*fh;
}
int n,m,cnt,A,B,q[M],hd,tl,D,X[M],Y[M];
LL ans;
bool vis[M];
vector<int> t[M<<1],f[M<<1];
map<int,int> ID1,ID2;
int fd(int kd,int x){return f[kd][x]==x?x:f[kd][x]=fd(kd,f[kd][x]);}
bool cmp(int x,int y){return X[x]<X[y];}
void gtin(int x){
int tar; X[x]=read(),Y[x]=read();
if(ID1.count(X[x]+Y[x])) tar=ID1[X[x]+Y[x]];
else cnt++,tar=ID1[X[x]+Y[x]]=cnt,f[tar].push_back(0);
t[tar].push_back(x),f[tar].push_back(f[tar].size());
if(ID2.count(X[x]-Y[x])) tar=ID2[X[x]-Y[x]];
else cnt++,tar=ID2[X[x]-Y[x]]=cnt,f[tar].push_back(0);
t[tar].push_back(x),f[tar].push_back(f[tar].size());
}
int getpos(int kd,int minx){
int l=0,r=t[kd].size(),fin=0,md;
while(l<=r){
md=((l+r)>>1);
if(md<t[kd].size()&&X[t[kd][md]]<minx) l=md+1;
else fin=md,r=md-1;
} return fin;
}
void del(int x){
vis[x]=true; q[tl++]=x;
int kd=ID1[X[x]+Y[x]],fin; fin=getpos(kd,X[x]);
f[kd][fin]=fd(kd,fin+1),kd=ID2[X[x]-Y[x]];
fin=getpos(kd,X[x]),f[kd][fin]=fd(kd,fin+1);
}
void calc(int kd,int lx,int rx,bool down){
int t1=getpos(kd,lx),t2=getpos(kd,rx),tot=t[kd].size();
if(down&&t1<tot&&X[t[kd][t1]]<=lx) t1++;
if(down&&t2<tot&&X[t[kd][t2]]<=rx) t2++; ans+=t2-t1;
for(t1=fd(kd,t1);t1<t2&&t1<tot;t1=fd(kd,t1+1)) del(t[kd][t1]); }
void getans(int x){
if(ID2.count(X[x]-Y[x]-D)) calc(ID2[X[x]-Y[x]-D],X[x]-D,X[x],false);
if(ID1.count(X[x]+Y[x]+D)) calc(ID1[X[x]+Y[x]+D],X[x],X[x]+D,false);
if(ID1.count(X[x]+Y[x]-D)) calc(ID1[X[x]+Y[x]-D],X[x]-D,X[x],true);
if(ID2.count(X[x]-Y[x]+D)) calc(ID2[X[x]-Y[x]+D],X[x],X[x]+D,true);
}
void solve(){
del(A),del(B),D=abs(X[A]-X[B])+abs(Y[A]-Y[B]);
while(hd<tl){int now=q[hd++];getans(now);}
}
int main(){
n=read(),A=read(),B=read();
for(int i=1;i<=n;i++) gtin(i);
for(int i=1;i<=cnt;i++) sort(t[i].begin(),t[i].end(),cmp);
solve(),printf("%lld\n",ans/2); return 0;
}
Arc065_E Manhattan Compass的更多相关文章
- [arc065E]Manhattan Compass[曼哈顿距离和切比雪夫距离转换]
Description 传送门 Solution 题目要求的是曼达顿距离,对于每个点(x,y),我们把它变为(x-y,x+y),就可以转换成求切比雪夫距离了. 证明如下:$max(\left | (x ...
- AT2160 へんなコンパス / Manhattan Compass
传送门 乍一看像是一个计算几何,然后想到了BFS,但是苦于无奈\(O(n^2)\)不会优化 然后以下参考zjq_shadow大佬的思路 显然发现曼哈顿距离很麻烦,除了暴力枚举貌似没什么很好的办法 考虑 ...
- AtCoder 瞎做
目录 ARC 058 E - 和風いろはちゃん / Iroha and Haiku 题意 题解 技巧 代码 ARC 059 F - バイナリハック / Unhappy Hacking 题意 题解 技巧 ...
- NOIp2018模拟赛三十八
爆〇啦~ A题C题不会写,B题头铁写正解: 随手过拍很自信,出分一看挂成零. 若要问我为什么?gtmdsubtask! 神tm就一个subtask要么0分要么100,结果我预处理少了一点当场去世 难受 ...
- 【AtCoder】ARC065
ARC065 C - 白昼夢 / Daydream 直接递推就好 #include <bits/stdc++.h> #define fi first #define se second # ...
- Sass之坑Compass编译报错
前段时间在使用Compass时遇到了其为难处理的一个坑,现记录到博客希望能帮助到各位. 一.问题: 利用Koala或者是gulp编译提示如下,截图为koala编译提示错误: 二.解决办法 从问题截图上 ...
- compass reset和layout [Sass和compass学习笔记]
reset 可以重置浏览器的的html的默认样式,因为各个浏览器对有些元素解析差别很大 通过重置样式可以让样式的浏览器兼容 更简单 使用方法简单 @import "compass/reset ...
- 分享15款很实用的 Sass 和 Compass 工具
Sass 是 CSS 的扩展,增加了嵌套规则,变量,混入功能等很多更多.它简化了组织和维护 CSS 代码的成本.Compass 是一个开源的 CSS 框架,使得使用 CSS3 和流行的设计模式比以往任 ...
- SASS+COMPASS 自适应 学习笔记
来源 http://snugug.github.io/RWD-with-Sass-Compass/#/ 1 安装 COMPASS 扩展 安装 方式 gem 'extension', '~>X.Y ...
随机推荐
- window下安装php7的memcache扩展
安装memcache:http://www.runoob.com/memcached/memcached-connection.html1.4.4 c:\memcached\memcached.exe ...
- ios推送服务,php服务端
本文转载至http://my.oschina.net/AStar/blog/176531 生成证书 证书生成参考:https://parse.com/tutorials/ios-push-noti ...
- 在UIWebView中设置cookie
本文转载至 http://blog.csdn.net/chengyakun11/article/details/8863878 项目中,需要在打开3g网页时,通过cookie传递一些信息. 实现代码 ...
- python 深复制与浅复制------copy模块
模块解读: 浅复制: x = copy.copy(y)深复制: x = copy.deepcopy(y)(注:模块特有的异常,copy.Error) 深copy与浅copy的差别主要体现在当有混合对象 ...
- Ubuntu 下安装JDK1.8
好困,不行了,我要睡觉了,先上图吧!
- 让Xcode支持高版本系统设备真机测试
最新支持11.2 (15C107) Xcode只可以支持iPhone手机对应iOS系统以下的真机测试.一般想要支持最新的iPhone手机系统,有两个方法. 第一.就需要更新Xcode,这一个方法有一个 ...
- React-Native 安装改变镜像
3.安装完node后建议设置npm镜像以加速后面的过程,否则后面插件安装巨慢 npm config set registry https://registry.npm.taobao.org --glo ...
- Go 学习笔记
官网: https://golang.org/ 环境: $GOROOT: GOROOT环境变量指定了Go的安装目录. $GOPATH: GOPATH 环境变量指定workspace的目录. 命令行: ...
- php自定义函数: 加密下载地址
function getdownurl($downurl, $extime = "3600", $serverid = 1) { if (empty($downurl)) { re ...
- 我的Android进阶之旅------>解决 Error: ShouldNotReachHere() 问题
在Android项目中创建一个包含main()方法的类,直接右键运行该类时会报如下错误: # # An unexpected error has been detected by Java Runti ...