[JZOJ100019]A--dfn序+扫描线
[JZOJ100019]A--dfn序+扫描线
题目链接
太懒了自行搜索
分析
这道题查了一个下午的错,真的心态崩了
不过这道题确实妙啊
类比于喝喝喝,我们发现任何一条覆盖了非法路径的路径一定不合法,假如非法路径为\((x,kx)\),设\(u\)为\(dfn\)序中较小的那个点,\(v\)为较大的那个;假如覆盖了路径\((u,v)\)的路径为\((a,b)(dfn_a<dfn_b)\);
设数组\(ed[x]\)为\(x\)的子树中\(dfs\)最大的那个点的\(dfs\)序,显然这可以和\(dfs\)序在一遍\(dfs\)中求出
我们还需要知道:\(x\) 的子树中\(dfs\)序是连续的一段
考虑两种情况:
\(u\)是\(lca(u,v)\),分析发现这也有两种情况
先设\(g\)为路径\((u,v)\)上距\(u\)最近的那个点Case#1
\(dfn[v]<=dfn[a]<=ed[v]\) , \(dfn[b]>ed[g]\)
Case#2
\(dfn[v]<=dfn[b]<=ed[v]\) , \(dfn[g]>dfn[a]\)
\(u\)不是\(lca(u,v)\)
这个只有一种情况,\(a,b\)分别在\(u,v\)子树中,即
Case#3
\(dfn[u]<=dfn[a]<=ed[u]<=dfn[v]<=dfn[b]<=ed[v]\)
我们把\((a,b)\)看成一个有序数对的话,发现其实每一种\(case\)都围成了一个矩形,我们只需要求出矩形面积的并之后即可得出所有不合法的路径,用总路径数\(n*(n-1)/2\)减去不合法即所求
求矩形面积的并使用线段树的扫描线法,但是发现一个点我们也要算进贡献,并非严格意义的二维图形;
因为值域很小\((1e5)\),学会了一种船新操作,我们不排序,将纵坐标相同的丢入一个容器,直接从小到大钦定纵坐标.不断取出对应容器内的线段加入,然后计算剩余的点数即可(纵坐标已钦定)
一道图论题转化成了数据结构,也是精妙啊
代码
/*
code by RyeCatcher
*/
inline char gc(){
static char buf[SIZE],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,SIZE,stdin),p1==p2)?EOF:*p1++;
}
template <class T>inline void read(T &x){
x=0;int ne=0;char c;
while((c=gc())>'9'||c<'0')ne=c=='-';x=c-48;
while((c=gc())>='0'&&c<='9')x=(x<<3)+(x<<1)+c-48;x=ne?-x:x;return ;
}
const int maxn=500005;
const int inf=0x7fffffff;
int n;
struct Edge{
int ne,to;
}edge[maxn<<1];
int h[maxn],num_edge=1;
int fa[maxn][17];
inline void add_edge(int f,int to){
edge[++num_edge].ne=h[f];
edge[num_edge].to=to;
h[f]=num_edge;
}
int dfn[maxn],ed[maxn],dep[maxn],tot=0;
void dfs(int now){
int v;dfn[now]=++tot;
for(ri i=h[now];i;i=edge[i].ne){
v=edge[i].to;
if(v==fa[now][0]||dfn[v])continue;
dep[v]=dep[now]+1,fa[v][0]=now;
for(ri i=1;i<=16;i++)fa[v][i]=fa[fa[v][i-1]][i-1];
dfs(v);
}
ed[now]=tot;
return ;
}
int get_g(int x,int y){//x到y路径上距y最近
for(ri i=16;i>=0;i--){
//printf("%d %d\n",fa[x][i],x);
if(dep[fa[x][i]]>dep[y]){
x=fa[x][i];
//printf("%d\n",x);
}
}
return x;
}
int L,R,dta;
struct Seg{
int l,r,h,d;
Seg(){l=r=h=d=0;}
Seg(int _l,int _r,int _h,int _d){l=_l,r=_r,h=_h,d=_d;}
}seg[maxn<<4];
int poi=0;
ll sum[maxn<<2];
ll tag[maxn<<2];
vector <int> dd[maxn];
inline void modify(int now,int l,int r){
if(tag[now]>0)sum[now]=(r-l+1);
else sum[now]=sum[now<<1]+sum[now<<1|1];
}
void update(int now,int l,int r){
if(L<=l&&r<=R){
tag[now]+=dta;
modify(now,l,r);
return ;
}
int mid=(l+r)>>1;
if(L<=mid)update(now<<1,l,mid);
if(mid<R)update(now<<1|1,mid+1,r);
modify(now,l,r);
return ;
}
ll ans=0;
int main(){
int x,y,ex,ey,g;
FO(a);
//freopen("a5.in","r",stdin);
read(n);
for(ri i=1;i<n;i++){
read(x),read(y);
add_edge(x,y);
add_edge(y,x);
}
dep[1]=1,fa[1][0]=0;
dfs(1);
int p,q;
for(ri i=1;i<=n;i++){
for(ri j=i+i;j<=n;j+=i){
p=i,q=j;
if(dfn[p]<dfn[q])swap(p,q);
x=dfn[p],y=dfn[q];
ex=ed[p],ey=ed[q];
if(x>=y&&x<=ey){
g=get_g(p,q);
dd[x].push_back(++poi);
seg[poi]=Seg(1,dfn[g]-1,x,1);
dd[ex+1].push_back(++poi);
seg[poi]=Seg(1,dfn[g]-1,ex,-1);
if(ed[g]==n)continue;
dd[ed[g]+1].push_back(++poi);
seg[poi]=Seg(x,ex,ed[g]+1,1);
dd[n+1].push_back(++poi);
seg[poi]=Seg(x,ex,n,-1);
}
else {
dd[x].push_back(++poi);
seg[poi]=Seg(y,ey,x,1);
dd[ex+1].push_back(++poi);
seg[poi]=Seg(y,ey,ex,-1);
}
}
}
for(ri i=1;i<=n;i++){
for(ri j=0;j<dd[i].size();j++){
x=dd[i][j];
L=seg[x].l,R=seg[x].r,dta=seg[x].d;
update(1,1,n);
}
ans+=sum[1];
}
printf("%lld\n",1ll*n*(n-1)/2-ans);
return 0;
}
[JZOJ100019]A--dfn序+扫描线的更多相关文章
- 洛谷P4689 [Ynoi2016]这是我自己的发明(莫队,树的dfn序,map,容斥原理)
洛谷题目传送门 具体思路看别的题解吧.这里只提两个可能对常数和代码长度有优化的处理方法. I 把一个询问拆成\(9\)个甚至\(16\)个莫队询问实在是有点珂怕. 发现询问的一边要么是一个区间,要么是 ...
- 线段树分治总结(线段树分治,线段树,并查集,树的dfn序,二分图染色)
闲话 stO猫锟学长,满脑子神仙DS 网上有不少Dalao把线段树分治也归入CDQ分治? 还是听听YCB巨佬的介绍: 狭义:只计算左边对右边的贡献. 广义:只计算外部对内部的贡献. 看来可以理解为广义 ...
- 洛谷P3703 [SDOI2017]树点涂色(LCT,dfn序,线段树,倍增LCA)
洛谷题目传送门 闲话 这是所有LCT题目中的一个异类. 之所以认为是LCT题目,是因为本题思路的瓶颈就在于如何去维护同颜色的点的集合. 只不过做着做着,感觉后来的思路(dfn序,线段树,LCA)似乎要 ...
- Kattis - boxes (dfn序)
Boxes There are N boxes, indexed by a number from 1 to N . Each box may (or not may not) be put into ...
- 树的遍历顺序 - dfs序|欧拉序|dfn序(备忘)
(仅作备忘) dfs序是dfs过程中对于某节点进入这个节点的子树和离开子树的顺序 满足每个节点都会在dfs序上出现恰好两次 任意子树的dfs序都是连续的 欧拉序是dfs过程中经过节点的顺序 每个节点至 ...
- bzoj 3779 重组病毒 好题 LCT+dfn序+线段树分类讨论
题目大意 1.将x到当前根路径上的所有点染成一种新的颜色: 2.将x到当前根路径上的所有点染成一种新的颜色,并且把这个点设为新的根: 3.查询以x为根的子树中所有点权值的平均值. 分析 原题codec ...
- bzoj 2819 Nim dfn序+树状数组维护区间异或值
题目大意 著名游戏设计师vfleaking,最近迷上了Nim.普通的Nim游戏为:两个人进行游戏,N堆石子,每回合可以取其中某一堆的任意多个,可以取完,但不可以不取.谁不能取谁输.这个游戏是有必胜策略 ...
- 51nod 1576 Tree and permutation(树的重心+dfn序)
乍一看我不会. 先不考虑加点. 先考虑没有那个除\(2\). 考虑每一条边的贡献,假设某一条边把这棵树分成大小为x,y的两个部分. 那么这条边最多可以被使用\(min(x,y)*2\)次(因为有进有出 ...
- BZOJ2819 Nim 【dfn序 + lca + 博弈论】
题目 著名游戏设计师vfleaking,最近迷上了Nim.普通的Nim游戏为:两个人进行游戏,N堆石子,每回合可以取其中某一堆的任意多个,可以取完,但不可以不取.谁不能取谁输.这个游戏是有必胜策略的. ...
随机推荐
- C++在线编程网站
1.推荐 http://www.dooccn.com/cpp/ 2.https://wandbox.org/ 3.https://www.tutorialspoint.com/compile_cpp_ ...
- QSplitter的使用案例
#include <QApplication> #include <QSplitter> #include <QTextEdit> #include <QTe ...
- JS微信网页使用图片预览(放大缩小)
前言 需求在微信网页中客户点击图片可进行预览放大缩小功能,网上找了各种js方式实现, 唯一的麻烦就是不兼容或者和项目框架不兼容 次函数只只用于部分客户端,否则会出现 WeixinJSBridge is ...
- (十八)Centos之firewall 防火墙命令
如果你的系统上没有安装使用命令安装 #yum install firewalld //安装firewalld 防火墙 开启服务 # systemctl start firewalld.service ...
- (十四)Centos之安装vsftp服务
一.为什么要安装vsftp服务 我们需要向centos操作系统的服务器上上传文件或者下载文件,这时候,ftp有必要安装下,我们选择主流的vsftp 二.安装 第一步:安装vsftp yum insta ...
- pycharm重命名文件
先右键要重命名的文件,然后按照下图操作:
- 最新 好未来java校招面经 (含整理过的面试题大全)
从6月到10月,经过4个月努力和坚持,自己有幸拿到了网易雷火.京东.去哪儿. 好未来等10家互联网公司的校招Offer,因为某些自身原因最终选择了 好未来.6.7月主要是做系统复习.项目复盘.Leet ...
- 多年经验【Parallels Desktop14.0.1 永久激活 】版 推荐苹果mac 虚拟机pmg序列号
parallels desktop 14 mac 激活码 parallels 13免费密钥 parallels desktop 14 激活码 很多用 MAC 的朋友发现平时离不开 W ...
- [转帖]Linux - NetworkManager网络管理工具
Linux - NetworkManager网络管理工具 linux运维菜 发布时间:18-10-3020:32优质原创作者 简介 网络管理器(NetworManager)是检测网络.自动连接网络的程 ...
- [转帖]Oracle报错ORA-26563--当重命名表时碰到物化视图
Oracle报错ORA-26563--当重命名表时碰到物化视图 https://www.toutiao.com/i6739137279115133447/ 原创 波波说运维 2019-09-26 00 ...