bzoj3051[WC2013]平面图(树上倍增+平面图转对偶图+扫描线)
简要题意:二维平面上n个点,点之间有一些连线,连线不在点之外的地方相交,将平面分为若干个区域。给出一些询问点对,问从这个点所在的区域走到另一个点所在的区域的最小代价。
题解:这道题首先可以把平面图转对偶图,这一点比较容易发现。然后对于从左指向右的线段,运用扫描线的思想,扫到左端点加入平衡树,扫到右端点从平衡树中删除。因为两线互不相交,所以相对位置不变。然后建立平面直角坐标系,y轴可以随意左右平移。对于每一条线段,我们使右端点正好在y轴上,然后选择两线段左端点x坐标比较大的作为比较直线,计算这条直线与两线段的交点的高低。然后就是一堆细节。真实的码农题。我写了3个namespace不然受不了
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+,inf=1e9;
const double eps=1e-;
int n,m,q,num,bel[N<<];
struct edge{int v,nxt,w,id;}e[N<<];
struct Edge{int u,v,w;double o;}d[N<<];
struct point{
int x,y,id;
void init(int i)
{
double xx,yy;scanf("%lf%lf",&xx,&yy);
x=(int)(xx*+eps),y=(int)(yy*+eps),id=i;
}
}p[N],Q[N<<];
struct node{int x,y,t;};
bool operator<(node a,node b){return a.x!=b.x?a.x<b.x:a.t<b.t;}
namespace MST{
int cnt,ecnt,f[N],hd[N],fa[N][],d[N][],dep[N];
Edge E[N<<];
edge e2[N<<];
void adde(int x,int y,int z)
{
e2[++ecnt]=(edge){y,hd[x],z},hd[x]=ecnt;
e2[++ecnt]=(edge){x,hd[y],z},hd[y]=ecnt;
}
int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
bool cmp(Edge a,Edge b){return a.w<b.w;}
void add(int u,int v,int w){E[++cnt]=(Edge){u,v,w};}
void dfs(int x)
{
for(int i=;i<=;i++)
{
fa[x][i]=fa[fa[x][i-]][i-];
d[x][i]=max(d[x][i-],d[fa[x][i-]][i-]);
}
for(int i=hd[x];i;i=e2[i].nxt)
if(e2[i].v!=fa[x][])
dep[e2[i].v]=dep[x]+,fa[e2[i].v][]=x,d[e2[i].v][]=e2[i].w,dfs(e2[i].v);
}
int lca(int x,int y)
{
if(dep[x]<dep[y])swap(x,y);
int ret=,t=dep[x]-dep[y];
for(int i=;(<<i)<=t;i++)if(t&(<<i))ret=max(ret,d[x][i]),x=fa[x][i];
if(x==y)return ret;
for(int i=;~i;i--)
if(fa[x][i]!=fa[y][i])ret=max(ret,max(d[x][i],d[y][i])),x=fa[x][i],y=fa[y][i];
return max(ret,max(d[x][],d[y][]));
}
void work()
{
sort(E+,E+cnt+,cmp);
for(int i=;i<=num;i++)f[i]=i;
for(int i=;i<=cnt;i++)
{
int x=find(E[i].u),y=find(E[i].v);
if(x!=y)adde(E[i].u,E[i].v,E[i].w),f[x]=y;
}
dfs();
for(int i=;i<=q;i++)
{
if(bel[i]==||bel[i+q]==){puts("-1");continue;}
int ans=lca(bel[i],bel[i+q]);
if(ans==inf)puts("-1");else printf("%d\n",ans);
}
}
}
namespace ScanLine{
struct cmp{
bool operator()(int i,int j)
{
if(d[i].u==d[j].u)return d[i].o>d[j].o;
int x=max(p[d[i].u].x,p[d[j].u].x);
double x1=1.0*(p[d[i].v].y-p[d[i].u].y)*(x-p[d[i].v].x)/(p[d[i].v].x-p[d[i].u].x)+p[d[i].v].y;
double x2=1.0*(p[d[j].v].y-p[d[j].u].y)*(x-p[d[j].v].x)/(p[d[j].v].x-p[d[j].u].x)+p[d[j].v].y;
return x1>x2;
}
};
set<int,cmp>S;
node a[N<<];
int cnt=;
void work()
{
for(int i=;i<=(m<<)+;i++)
if(p[d[i].u].x<p[d[i].v].x)a[++cnt]=(node){p[d[i].u].x,i,},a[++cnt]=(node){p[d[i].v].x,i,};
for(int i=;i<=q;i++)a[++cnt]=(node){Q[i].x,i,},a[++cnt]=(node){Q[i+q].x,i+q,};
sort(a+,a+cnt+);
for(int i=;i<=cnt;i++)
if(!a[i].t)S.erase(a[i].y);
else if(a[i].t==)S.insert(a[i].y);
else{
p[n+]=(point){a[i].x,Q[a[i].y].y,};
p[n+]=(point){a[i].x-,Q[a[i].y].y,};
d[(m<<)+]=(Edge){n+,n+,,atan2(p[n+].y-p[n+].y,p[n+].x-p[n+].x)};
S.insert((m<<)+);
set<int,cmp>::iterator it=S.find((m<<)+);
if(it!=S.begin())bel[a[i].y]=e[*--it].id;else bel[a[i].y]=;
S.erase((m<<)+);
}
}
}
namespace Graph{
int cnt=,hd[N],nxt[N<<];
vector<pair<double,int> >G[N];
point s=(point){inf,inf,};
void adde(int x,int y,int z)
{
e[++cnt]=(edge){y,hd[x],z,-},hd[x]=cnt;
e[++cnt]=(edge){x,hd[y],z,-},hd[y]=cnt;
}
void dosort(int x)
{
for(int i=hd[x];i;i=e[i].nxt)
G[x].push_back(make_pair(atan2(p[e[i].v].y-p[x].y,p[e[i].v].x-p[x].x),i));
sort(G[x].begin(),G[x].end());
int sz=G[x].size();
for(int i=;i<sz;i++)nxt[G[x][i].second^]=G[x][(i+)%sz].second;
}
void find(int x){for(int i=x;e[i].id<;i=nxt[i])e[i].id=num;}
void work()
{
cnt=;
for(int i=;i<=n;i++)
{
p[i].init(i);
if(s.x>p[i].x||s.x==p[i].x&&s.y>p[i].y)s=p[i];
}
for(int i=,x,y,z;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
d[i<<]=(Edge){x,y,z,atan2(p[y].y-p[x].y,p[y].x-p[x].x)};
d[(i<<)+]=(Edge){y,x,z,atan2(p[x].y-p[y].y,p[x].x-p[y].x)};
adde(x,y,z);
}
scanf("%d",&q);
for(int i=;i<=q;i++)Q[i].init(i),Q[i+q].init(i+q);
for(int i=;i<=n;i++)dosort(i);
num=;
find(G[s.id][].second);
for(int i=;i<=cnt;i++)if(e[i].id<)++num,find(i);
for (int i=;i<=cnt;i+=)
if(e[i].id!=&&e[i^].id!=)MST::add(e[i].id,e[i^].id,e[i].w);
else MST::add(e[i].id,e[i^].id,inf);
}
}
int main()
{
scanf("%d%d",&n,&m);
Graph::work();
ScanLine::work();
MST::work();
}
bzoj3051[WC2013]平面图(树上倍增+平面图转对偶图+扫描线)的更多相关文章
- Codevs 2370 小机房的树 LCA 树上倍增
题目描述 Description 小机房有棵焕狗种的树,树上有N个节点,节点标号为0到N-1,有两只虫子名叫飘狗和大吉狗,分居在两个不同的节点上.有一天,他们想爬到一个节点上去搞基,但是作为两只虫子, ...
- NOIP2013 货车运输 (最大生成树+树上倍增LCA)
死磕一道题,中间发现倍增还是掌握的不熟 ,而且深刻理解:SB错误毁一生,憋了近2个小时才调对,不过还好一遍AC省了更多的事,不然我一定会疯掉的... 3287 货车运输 2013年NOIP全国联赛提高 ...
- HDU 4822 Tri-war(LCA树上倍增)(2013 Asia Regional Changchun)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4822 Problem Description Three countries, Red, Yellow ...
- [NOIP2013/Codevs3287]货车运输-最小[大]生成树-树上倍增
Problem 树上倍增 题目大意 给出一个图,给出若干个点对u,v,求u,v的一条路径,该路径上最小的边权值最大. Solution 看到这个题第一反应是图论.. 然而,任意路径最小的边权值最大,如 ...
- 树上倍增求LCA及例题
先瞎扯几句 树上倍增的经典应用是求两个节点的LCA 当然它的作用不仅限于求LCA,还可以维护节点的很多信息 求LCA的方法除了倍增之外,还有树链剖分.离线tarjan ,这两种日后再讲(众人:其实是你 ...
- [树上倍增+二分答案][NOIP2012]运输计划
题目背景 公元 2044 年,人类进入了宇宙纪元. 题目描述 公元 2044 年,人类进入了宇宙纪元 L 国有 nn 个星球,还有 n-1n−1 条双向航道,每条航道建立在两个星球之间,这 n-1n− ...
- 两种lca的求法:树上倍增,tarjan
第一种:树上倍增 f[x,k]表示x的2^k辈祖先,即x向根结点走2^k步达到的结点. 初始条件:f[x][0]=fa[x] 递推式:f[x][k]=f[ f[x][k-1] ][k-1] 一次bfs ...
- 【NOIP2013/Codevs3287】货车运输-最小生成树(大)-树上倍增
https://www.luogu.org/problemnew/show/P1967 由题可知,我们走的路的边应尽可能大,所以通过kruscal建最大生成树的图,再树上倍增,注意可能有多棵树; #i ...
- LCA树上倍增
LCA就是最近公共祖先,比如 节点10和11的LCA就是8,9和3的LCA就是3. 我们这里讲一下用树上倍增来求LCA. 大家都可以写出暴力解法,两个节点依次一步一步往上爬,直到爬到了相同的一个节点. ...
随机推荐
- linux文件行首行尾添加或替换
sed -i 's/\(^.*\)/http:\/\/www.blutmagie.de\/img\/flags\//g' cc.txt sed -i 's/\($\)/.gif/g' cc.txt
- Unity2018 Open C# Project Error
错误日志 升级到Unity2018之后,使用 Open C# Project 打开VS工程,出现报错,无法启动VS. 错误日志如下: ArgumentException: Value does not ...
- 自动化测试之路3-selenium3+python3环境搭建
1.首先安装火狐浏览器 有单独文章分享怎么安装 2.搭建python环境 安装python,安装的时候把path选好,就不用自己在配置,安装方法有单独文档分享 安装好以后cmd打开输入python查 ...
- 4.机器学习——统计学习三要素与最大似然估计、最大后验概率估计及L1、L2正则化
1.前言 之前我一直对于“最大似然估计”犯迷糊,今天在看了陶轻松.忆臻.nebulaf91等人的博客以及李航老师的<统计学习方法>后,豁然开朗,于是在此记下一些心得体会. “最大似然估计” ...
- 调用远程主机上的 RMI 服务时抛出 java.rmi.ConnectException: Connection refused to host: 127.0.0.1 异常原因及解决方案
最近使用 jmx 遇到一个问题,client/server 同在一台机器上,jmx client能够成功连接 server,如果把 server 移植到另一台机器上192.168.134.128,抛出 ...
- 推荐一套Angular2的UI模板
Core UI Core UI是一款基于Bootstrap4的UI模板,有html.angular2,react和vue版.我是在使用angular2版本中发现其项目结构不符合angular风格指南推 ...
- css_属性
老师的博客:https://www.cnblogs.com/liwenzhou/p/7999532.htm css的属性 整体属性的:作用于全局 width:表示宽 height:表示长 color: ...
- Python简单的多线程demo:装逼写法
用面向对象来写多线程: import threading class MyThread(threading.Thread): def __init__(self, n): super(MyThread ...
- 【English Email】CIP payouts now in Workday
simplification简化的[ˌsɪmplɪfɪˈkeɪʃn] quota配额[ˈkwoʊtə] regional区域的[ˈriːdʒənl] mechanics技工[məˈkænɪks] ...
- Django-CRM项目学习(一)-admin组件
开始今日份整理 1.admin组件使用 1.1 创建django项目以及开启APP01 略 1.2 创建类 使用django自带的sqlite3的小型文件型的数据库 注:使用sqlite3类型的数据库 ...