【题解】[CJOI2019Chebnear]
【题解】[CJOI2019Chebnear]
题目描述
给定平面上有\(n\)个仇人,\((x,y) ,x,y \in R^+\) ,\(n\)个人都是仇人关系,而有仇人关系的一对人的切比雪夫距离若\(\le k\)则会发生“吵架”。可是,对于有些仇人\(i\)和\(j\),先记为三元组\(P(t,i,w)\),他们的的仇人关系并不深,假设给全体\(n\)个人每人\(\ge w\)的钱他们就变为朋友了。而朋友不会吵架,而且朋友的朋友也是你的朋友。
现在问你,你最少给多少钱\(w\)可是让所有人不吵架?
数据范围
\(n\le100000,m\le100000\)
说明
切比雪夫距离:\(d=max(|x_i-x_j|,|y_i-y_j|)\)
\(Solution\)
显然二分答案\(w\)。接下来的任务就是当给每个人\(w\)钱的时候,是否有人吵架了。
然后\(O(m)\)枚举那些朋友可以连接,直接用并查集维护朋友关系。
考虑暴力。
- \(O(n^2)\)枚举每个点对,看他们是否吵架。这样显然超时。
我们仔细思考一下,发现吵架只和以下两个因素有关。
- \(y\)安全但\(x\)距离过近。
- \(x\)安全但\(y\)距离过近。
- 朋友
朋友关系已经用并查集维护了,可以\(O(\alpha(n))=O(1)\)查询。接下来就是\(x,y\)的问题了。
考虑把所有点按\(x\)排序。这样我们就可以快速确定一段区间内的人是否有可能\(x\)过近,然后对于这些点\(O(n_?^2)\)枚举是否会吵架就好了。
但是怎么缩小\(n_?\)呢?
对于平面最近点对问题,可以证明\(n_? \le 6\)。这里我们也可以类比,每次和最近点对一样,直接从中间分开,分治。我们的\(n_?\)一定比最近点对的\(n_?\)还要小,因为我们还有\(k\)的限制呀,帮助我们剪枝。
接下来上代码了,今天的代码是不是码风很棒?
#include<bits/stdc++.h>
#define RP(t,a,b) for(register int (t)=(a),edd_=(b);t<=edd_;++t)
#define DRP(t,a,b) for(register int (t)=(a),edd_=(b);t>=edd_;--t)
#define ERP(t,a) for(int t=head[a];t;t=e[t].nx)
#define pushup(x) seg[(x)]=seg[(x)<<1]+seg[(x)<<1|1]
#define midd register int mid=(l+r)>>1
#define TMP template<class ccf>
#define rgt L,R,mid,r,pos<<1|1
#define lef L,R,l,mid,pos<<1
#define all 1,n,1
using namespace std;typedef long long ll;
TMP inline ccf Max(ccf a,ccf b){return a<b?b:a;}
TMP inline ccf Min(ccf a,ccf b){return a<b?a:b;}
TMP inline ccf Abs(ccf a){return a<0?-a:a;}
TMP inline ccf qr(ccf k){
char c=getchar();ccf x=0;int q=1;
while(c<48||c>57)q=c==45?-1:q,c=getchar();
while(c>=48&&c<=57)x=x*10+c-48,c=getchar();
return q==-1?-x:x;}
//-----------------template&IO------------------
int n,m;double k;
const int maxn=1e5+5;
struct node{
double x,y;int id;
inline bool operator <(node a){return x<a.x;}
inline double operator -(node a){return Max(Abs(x-a.x),Abs(y-a.y));}
}data[maxn];
struct E{
int fr,to;double w;
inline bool operator <(E a){return w<a.w;}
}e[maxn<<1];
int rr[maxn];
inline int q(int x){register int t=x,temp,i=x;
while(rr[t]!=t) t=rr[t];
while(rr[i]!=i) temp=rr[i],rr[i]=t,i=temp;
return t;
}
inline void j(int x,int y){rr[q(x)]=q(y);}
inline bool in(int x,int y){return q(x)==q(y);}
struct DFN{
double pos;int id;
inline bool operator <(DFN x){return pos<x.pos;}
}dfn[maxn];
bool divd(int lb,int rb){register int mid=(lb+rb)>>1;
if(lb>=rb) return 1;
if(!divd(lb,mid)) return 0;
if(!divd(mid+1,rb)) return 0;
register int L=mid,R=mid;
while(data[mid].x-data[L-1].x<=k&&L>lb) --L;
while(data[R+1].x-data[mid].x<=k&&R<rb) ++R;
RP(t,L,R)RP(i,t+1,R)
if(!in(data[t].id,data[i].id)&&data[t]-data[i]<=k)
return 0;
return 1;
}
inline bool chek(double x){
RP(t,1,n) rr[t]=t;
RP(t,1,m) if(x>=e[t].w) j(e[t].fr,e[t].to);
return divd(1,n);
}
int main(){
#ifndef ONLINE_JUDGE
freopen("chebnear.in","r",stdin);
freopen("chebnear.out","w",stdout);
#endif
cin>>n>>m>>k;
RP(t,1,n)
cin>>data[t].x>>data[t].y,data[t].id=t;
RP(t,1,m)
cin>>e[t].fr>>e[t].to>>e[t].w;
sort(data+1,data+n+1);
sort(e+1,e+m+1);
int l=1,r=m,ans=m;
do{
midd;
if(chek(e[mid].w))
r=mid-1,ans=Min(ans,mid);
else
l=mid+1;
}while(l<=r);
if(chek(0)) ans=0;
printf("%.3lf\n",e[ans].w);
return 0;
}
【题解】[CJOI2019Chebnear]的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
随机推荐
- Data-structures-and-algorithms-interview-questions-and-their-solutions
https://techiedelight.quora.com/500-Data-structures-and-algorithms-interview-questions-and-their-sol ...
- 手机端用来上传用户头像的代码canvas
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- Microsoft JET Database Engine(0x80004005)未指定错误的解决方法
今天在给一台新的电脑安装IIS,安装成功,建立虚目录后,运行一个已经在别的机器上的正确的asp文件,就是不成功,提示:Microsoft JET Database Engine (0x80004005 ...
- MFC中 SDI/MDI框架各部分指针获取方式
VC MFC SDI/MDI框架各部分指针获取方式 整理总结一下,希望能帮助到别人. 获得CWinApp 获得CMainFrame 获得CChildFrame 获得CDocument 获得CV ...
- JS: document.getElementBy(), setInerval()
ylbtech-JavaScript-DOM document.getElementBy(),setInerval() 1.A,document.getElementBy()返回顶部 document ...
- AngularJS的稍复杂form验证
代码下载:https://files.cnblogs.com/files/xiandedanteng/angularjsSoccerFormCheck.rar 代码: <!DOCTYPE HTM ...
- C 作用域规则
C 作用域规则 任何一种编程中,作用域是程序中定义的变量所存在的区域,超过该区域变量就不能被访问.C 语言中有三个地方可以声明变量: 在函数或块内部的局部变量 在所有函数外部的全局变量 在形式参数的函 ...
- MQTT---HiveMQ源代码具体解释(十四)Persistence-LocalPersistence
源博客地址:http://blog.csdn.net/pipinet123 MQTT交流群:221405150 简单介绍 HiveMQ的Persistence提供配置包含File和Memory,以解决 ...
- 阿里云官方教程 Linux 系统挂载数据盘
适用系统:Linux(Redhat , CentOS,Debian,Ubuntu) * Linux的云服务器数据盘未做分区和格式化,可以根据以下步骤进行分区以及格式化操作. 下面的操作将会把数据盘划 ...
- 一致性哈希算法(c#版)
最近在研究"一致性HASH算法"(Consistent Hashing),用于解决memcached集群中当服务器出现增减变动时对散列值的影响.后来 在JAVAEYE上的一篇文章中 ...