【BZOJ】【2588】COT(Count On a Tree)
可持久化线段树
maya……树么……转化成序列……所以就写了个树链剖分……然后每个点保存的是从它到根的可持久化线段树。
然后就像序列一样查询……注意是多个左端点和多个右端点,处理方法类似BZOJ 1901
然后rausen(Orz!!!)粗来跟我说:你直接减去lca和fa[lca]不就好啦~搞树剖还多一个log……
我恍然大悟!然后两个都交了一下,事实证明:我链剖写的还行,LCA写的太丑……速度反而是多一个log的链剖快QAQ(另:因为边少我就偷懒没写边表,直接vector水过)
链剖:
/**************************************************************
Problem: 2588
User: ProgrammingApe
Language: C++
Result: Accepted
Time:3768 ms
Memory:47652 kb
****************************************************************/ //BZOJ 2588
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
inline int getint(){
int v=,sign=; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}
while(isdigit(ch)) {v=v*+ch-''; ch=getchar();}
return v*sign;
}
const int N=1e5+,INF=~0u>>;
/*******************template********************/
struct tree{
int cnt,l,r;
}t[N*];
int root[N],cnt,num;
int lc,rc,ln[N],rn[N];
#define mid (l+r>>1)
void update(int &o,int l,int r,int pos){
t[++cnt]=t[o]; o=cnt; t[o].cnt++;
if (l==r) return;
if (pos<=mid) update(t[o].l,l,mid,pos);
else update(t[o].r,mid+,r,pos);
}
int query_t(int rank){
int l=,r=num;
int tl=,tr=;
while(l!=r){
tl=tr=;
F(i,,lc) tl+=t[t[ln[i]].l].cnt;
F(i,,rc) tr+=t[t[rn[i]].l].cnt;
if (tr-tl>=rank){
F(i,,lc) ln[i]=t[ln[i]].l;
F(i,,rc) rn[i]=t[rn[i]].l;
r=mid;
}else{
F(i,,lc) ln[i]=t[ln[i]].r;
F(i,,rc) rn[i]=t[rn[i]].r;
l=mid+; rank-=tr-tl;
}
}
return l;
}
/*****************可持久化线段树 ***************/
vector<int>G[N];
int top[N],fa[N],son[N],dep[N],tot,size[N],a[N],b[N],n,m,lastans;
bool vis[N];
void dfs(int x,int f,int d){
vis[x]=;
fa[x]=f; dep[x]=d; size[x]=; son[x]=;
int maxsize=;
rep(i,G[x].size()){
int to=G[x][i];
if (vis[to]) continue;
dfs(to,x,d+);
size[x]+=size[to];
if (size[to]>maxsize) maxsize=size[to],son[x]=to;
}
}
void connect(int x,int f){
vis[x]=;
root[x]=root[fa[x]];
update(root[x],,num,a[x]);
// root[x]=++tot;
top[x]=f;
if (son[x]) connect(son[x],f);
rep(i,G[x].size()){
int to=G[x][i];
if (!vis[to]) connect(to,to);
}
}
void query(int x,int y,int k){
lc=rc=;
while(top[x]!=top[y]){
if (dep[top[x]]<dep[top[y]]) swap(x,y);
ln[++lc]=root[fa[top[x]]]; rn[++rc]=root[x];
x=fa[top[x]];
}
if (dep[x]>dep[y]) swap(x,y);
ln[++lc]=root[fa[x]]; rn[++rc]=root[y];
printf("%d",lastans=b[query_t(k)]);
}
/**********************链剖*********************/
int main(){
// freopen("input.txt","r",stdin);
n=getint(); m=getint();
F(i,,n) b[i]=a[i]=getint();
sort(b+,b+n+);
num=unique(b+,b+n+)-b-;
F(i,,n) a[i]=lower_bound(b+,b+num+,a[i])-b; int x,y,k;
F(i,,n){
x=getint(); y=getint();
G[x].pb(y); G[y].pb(x);
}
dfs(,,);
memset(vis,,sizeof vis);
connect(,); F(i,,m){
x=lastans^getint(); y=getint(); k=getint();
query(x,y,k);
if (i!=m) puts("");
}
return ;
}
倍增LCA:
/**************************************************************
Problem: 2588
User: Tunix
Language: C++
Result: Accepted
Time:4520 ms
Memory:53024 kb
****************************************************************/ //BZOJ 2588
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
inline int getint(){
int v=,sign=; char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') sign=-; ch=getchar();}
while(isdigit(ch)) {v=v*+ch-''; ch=getchar();}
return v*sign;
}
const int N=1e5+,INF=~0u>>;
/*******************template********************/
struct tree{
int cnt,l,r;
}t[N*];
int root[N],cnt,num;
int lc,rc,ln[N],rn[N];
#define mid (l+r>>1)
void update(int &o,int l,int r,int pos){
t[++cnt]=t[o]; o=cnt; t[o].cnt++;
if (l==r) return;
if (pos<=mid) update(t[o].l,l,mid,pos);
else update(t[o].r,mid+,r,pos);
}
int query_t(int rank){
int l=,r=num;
int tl=,tr=;
while(l!=r){
tl=tr=;
F(i,,lc) tl+=t[t[ln[i]].l].cnt;
F(i,,rc) tr+=t[t[rn[i]].l].cnt;
if (tr-tl>=rank){
F(i,,lc) ln[i]=t[ln[i]].l;
F(i,,rc) rn[i]=t[rn[i]].l;
r=mid;
}else{
F(i,,lc) ln[i]=t[ln[i]].r;
F(i,,rc) rn[i]=t[rn[i]].r;
l=mid+; rank-=tr-tl;
}
}
return l;
}
/*****************可持久化线段树 ***************/
vector<int>G[N];
int fa[N][],dep[N],a[N],b[N],n,m,lastans;
void dfs(int x){
F(i,,)
if (dep[x]>=(<<i)) fa[x][i]=fa[fa[x][i-]][i-];
else break;
root[x]=root[fa[x][]];
update(root[x],,num,a[x]);
rep(i,G[x].size()){
int to=G[x][i];
if (to==fa[x][]) continue;
fa[to][]=x; dep[to]=dep[x]+;
dfs(to);
}
}
int LCA(int x,int y){
if (dep[x]<dep[y]) swap(x,y);
int t=dep[x]-dep[y];
for(int i=;(<<i)<=t;i++)
if(t&(<<i)) x=fa[x][i];
D(i,,)
if(fa[x][i]!=fa[y][i])
x=fa[x][i],y=fa[y][i];
if (x==y) return x;
return fa[x][];
}
void query(int x,int y,int k){
lc=rc=;
int lca=LCA(x,y);
ln[++lc]=root[lca]; ln[++lc]=root[fa[lca][]];
rn[++rc]=root[x]; rn[++rc]=root[y];
printf("%d",lastans=b[query_t(k)]);
}
/**********************LCA*********************/
int main(){
n=getint(); m=getint();
F(i,,n) b[i]=a[i]=getint();
sort(b+,b+n+);
num=unique(b+,b+n+)-b-;
F(i,,n) a[i]=lower_bound(b+,b+num+,a[i])-b; int x,y,k;
F(i,,n){
x=getint(); y=getint();
G[x].pb(y); G[y].pb(x);
}
dfs(); F(i,,m){
x=lastans^getint(); y=getint(); k=getint();
query(x,y,k);
if (i!=m) puts("");
}
return ;
}
2588: Spoj 10628. Count on a tree
Time Limit: 12 Sec Memory Limit: 128 MB
Submit: 2581 Solved: 586
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
105 2 9 3 8 5 7 7
1 2
1 3
1 4
3 5
3 6
3 7
4 8
2 5 1
0 5 2
10 5 3
11 5 4
110 8 2
Sample Output
8
9
105
7
HINT
Source
【BZOJ】【2588】COT(Count On a Tree)的更多相关文章
- 【BZOJ1415】【NOI2005】聪聪和可可(动态规划,数学期望)
[BZOJ1415][NOI2005]聪聪和可可(动态规划,数学期望) 题面 BZOJ 题解 先预处理出当可可在某个点,聪聪在某个点时 聪聪会往哪里走 然后记忆化搜索一下就好了 #include< ...
- 【BZOJ2684】【CEOI2004】锯木厂选址(斜率优化,动态规划)
[BZOJ2684][CEOI2004]锯木厂选址(斜率优化,动态规划) 题面 万恶的BZOJ因为权限题的原因而做不了... 我要良心的提供题面 Description 从山顶上到山底下沿着一条直线种 ...
- 【BZOJ4911】[SDOI2017]切树游戏(动态dp,FWT)
[BZOJ4911][SDOI2017]切树游戏(动态dp,FWT) 题面 BZOJ 洛谷 LOJ 题解 首先考虑如何暴力\(dp\),设\(f[i][S]\)表示当前以\(i\)节点为根节点,联通子 ...
- 【BZOJ2324】[ZJOI2011]营救皮卡丘(网络流,费用流)
[BZOJ2324][ZJOI2011]营救皮卡丘(网络流,费用流) 题面 BZOJ 洛谷 题解 如果考虑每个人走的路径,就会很麻烦. 转过来考虑每个人破坏的点集,这样子每个人可以得到一个上升的序列. ...
- 【BZOJ4061】[Cerc2012]Farm and factory(最短路,构造)
[BZOJ4061][Cerc2012]Farm and factory(最短路,构造) 题面 BZOJ 然而权限题QwQ. 题解 先求出所有点到达\(1,2\)的最短路,不妨记为\(d_{u,1}, ...
- 【BZOJ5303】[HAOI2018]反色游戏(Tarjan,线性基)
[BZOJ5303][HAOI2018]反色游戏(Tarjan,线性基) 题面 BZOJ 洛谷 题解 把所有点全部看成一个\(01\)串,那么每次选择一条边意味着在这个\(01\)串的基础上异或上一个 ...
- 【BZOJ1185】[HNOI2007]最小矩形覆盖(凸包,旋转卡壳)
[BZOJ1185][HNOI2007]最小矩形覆盖(凸包,旋转卡壳) 题面 BZOJ 洛谷 题解 最小的矩形一定存在一条边在凸包上,那么枚举这条边,我们还差三个点,即距离当前边的最远点,以及做这条边 ...
- 【BZOJ1047】[HAOI2007]理想的正方形(单调队列,动态规划)
[BZOJ1047][HAOI2007]理想的正方形(单调队列,动态规划) 题面 BZOJ 洛谷 题解 直接一个单调队列维护一下没给点和它前面的\(n\)个位置的最大值,再用一次单调队列维护连续\(n ...
- 【BZOJ3437】小P的牧场(动态规划,斜率优化)
[BZOJ3437]小P的牧场(动态规划,斜率优化) 题面 BZOJ 题解 考虑暴力\(dp\),设\(f[i]\)表示强制在\(i\)处建立控制站的并控制\([1..i]\)的最小代价. 很显然,枚 ...
随机推荐
- SQL Server 数据类型简介
在 SELECT 的查询过程和查询结果中,每个列.变量.表达式和参数都具有一个相关的数据类型.数据类型用于指定某个对象可保存的数据的类型. SQL Server系统的数据类型主要有:数值类型.日期和时 ...
- Ubuntu上OpenStack DashBoard主题修改的方案
最近终于把OpenStack的整个初步环境安装部署并配置成功,后续工作就是安装一些可选项了,如Cinder.Swift部件了.我们的整个项目采用的颜色方案主要为清淡色系,蓝色为主,可OpenStack ...
- (转)Linux概念架构的理解
英文原文:Conceptual Architecture of the Linux Kernel 摘要 Linux kernel成功的两个原因:(1)架构设计支持大量的志愿开发者加入到开发过程中:(2 ...
- 移动Web轮播图IOS卡顿的问题
晚饭前,被测试吐槽说,banner轮播手动左右滑的时候会卡顿.我一看不科学啊,大水果手机怎么会卡顿.我一看测试手中拿的是iPod,我觉得大概是这小玩意性能不强悍,后来又拿来5S,依然会卡顿,有趣的是, ...
- ThinkPHP框架概述
框架概述 1.什么是框架 框架,即framework.其实就是某种应用的半成品,就是一组组件,供你选用完成你自己的系统. 2.框架的应用场景 ① 代码重用 ② 合理分工 ③ 解决团队协作开发问题 3. ...
- Difference between ref and out parameters
Original link: http://www.dotnet-tricks.com/Tutorial/csharp/K0Hb060414-Difference-between-ref-and-ou ...
- 《UNIX编程艺术》读书笔记
最近这段时间比较忙,利用业余时间看完了这本书.虽然书中讲到的很多例子都是上古文物,我没有用过,不过原理都是相通的,对我的启发很大.比如无所不在的KISS原则,实践中慢慢体会到的SPOT原则,无不产生共 ...
- ubuntu 常用参数设置
在Linux下,对于参数的设置,一般来说,都遵循这个规律.每一个功能程序,一定对于一个对于名字的配置文件. 涉及到多用户的功能配置,一定有一个全局的配置文件,对所有用户都生效,而每个用 ...
- 实现textarea自适应的方法
1.用div来模拟实现textarea自适应 <!doctype html> <html lang="en"> <head> <meta ...
- devpress控件属性说明表
XtraEditors 库中所有控件的公共功能 全部都可以绑定数据: 全部都可以独立使用或用于由 Developer Express 提供的容器控件(XtraGrid.XtraVerticalGrid ...