LWDB
题意:
给一棵 $n$ 个节点的树,维护两种操作:
1.将距离 $x$ $distance \leq d$ 的点染成 $c$
2.询问 $x$ 的颜色。
解法:
首先将染色可以转换为每个时间对应一个颜色,问题转化为区间取 $max$
动态树分治,即可。
考虑朴素点分治中的每一个重心,将其管辖的所有点按照到重心的距离从大到小排序,然后用线段树维护。
修改时,从 $x$ 向上走,将每一个重心 $t$ 中 $dist<x,t> + dist<t,y>\leq d$ 的点修改。
查找时,从 $x$ 向上走,查询每一个重心中 $x$ 所对应的值,求 $max$ 即可。
$O(nlog^2n)$
#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm> #define p E[i].x
#define N 100010
#define M 7000010
#define l(x) ch[x][0]
#define r(x) ch[x][1]
#define MP(x, y) make_pair(x, y)
#define LL long long using namespace std; typedef pair<int,int> PII; struct edge
{
int x, to, v;
}E[N<<]; int totE, totn, n, m, hev, lef_now;
int g[N], h[N], de[N], siz[N], root[N], fa[N][], col[N], d[N];
int ch[M][], setv[M];
vector<PII> nod[N], id[N];
vector<int> hevs[N];
int dis[N][];
bool v[N]; void addedge(int x, int y, int v)
{
E[++totE] = (edge){y, g[x], v}; g[x] = totE;
} int build(int l, int r)
{
int x = ++totn;
setv[x] = ;
if(l == r) return x;
int mid = (l+r)>>;
l(x) = build(l, mid);
r(x) = build(mid+, r);
return x;
} void push(int x)
{
if(!setv[x]) return;
setv[l(x)] = setv[x];
setv[r(x)] = setv[x];
setv[x] = ;
} int ask(int x, int l, int r, int qx)
{
if(l<r) push(x);
if(l == r) return setv[x];
int mid = (l+r)>>;
if(qx <= mid) return ask(l(x), l, mid, qx);
else return ask(r(x), mid+, r, qx);
} void change(int x, int l, int r, int ql, int qr, int qv)
{
if(l<r) push(x);
if(ql<=l && r<=qr)
{
setv[x] = qv;
return;
}
int mid=(l+r)>>;
if(ql <= mid) change(l(x),l,mid,ql,qr,qv);
if(mid < qr) change(r(x),mid+,r,ql,qr,qv);
} int dist(int x, int y)
{
if(de[x] < de[y]) swap(x, y);
int ans = ;
for(int i = ;~i;i--)
if(de[fa[x][i]] >= de[y]) ans += dis[x][i], x = fa[x][i];
if(x == y) return ans;
for(int i = ;~i;i--)
if(fa[x][i] != fa[y][i])
{
ans += dis[x][i]; x = fa[x][i];
ans += dis[y][i]; y = fa[y][i];
}
return ans+dis[x][]+dis[y][];
} void dfs1(int x, int tmp)
{
siz[x] = ;
h[x] = ;
for(int i = g[x];i;i = E[i].to)
if(!v[p] && p != tmp)
{
dfs1(p, x);
h[x] = max(h[x], siz[p]);
siz[x] += siz[p];
}
h[x] = max(h[x], lef_now - siz[x]);
if(!hev || h[x] < h[hev]) hev = x;
} void dfs(int x)
{
for(int i = g[x];i;i = E[i].to)
if(p != fa[x][])
{
fa[p][] = x;
dis[p][] = E[i].v;
de[p] = de[x] + ;
dfs(p);
}
} void dfs2(int x, int tmp, int now)
{
nod[now].push_back(MP(d[x], x));
for(int i = g[x];i;i = E[i].to)
if(!v[p] && p != tmp)
{
d[p] = d[x] + E[i].v;
dfs2(p, x, now);
}
} void solve(int x)
{
v[x] = ;
d[x] = ;
dfs2(x, x, x);
sort(nod[x].begin(), nod[x].end());
for(int i = ;i < (int)nod[x].size();i++)
{
id[x].push_back( MP(nod[x][i].second, i) );
hevs[nod[x][i].second].push_back(x);
}
sort(id[x].begin(), id[x].end());
root[x] = build(, nod[x].size()-);
for(int i = g[x];i;i = E[i].to)
if(!v[p])
{
hev = ;
lef_now = siz[p];
dfs1(p, x);
solve(hev);
}
} int find(int x,int len)
{
int l=,r=nod[x].size()-;
while(r-l>)
{
int mid = (l+r)>>;
if(nod[x][mid].first <= len) l=mid;
else r=mid;
}
for(int i=r;i>=l;i--)
if(nod[x][i].first <= len) return i;
return -;
} void com_change(int x, int len, int v)
{
for(int i = ;i < (int)hevs[x].size();i++)
{
int tmp = hevs[x][i];
int t = find(tmp, len - dist(x, tmp));
if(t == -) continue;
change(root[tmp], , id[tmp].size()-, , t, v);
}
} int com_ask(int x)
{
int ans = ;
for(int i = ;i < (int)hevs[x].size();i++)
{
int tmp = hevs[x][i];
int t = (*lower_bound(id[tmp].begin(), id[tmp].end(), MP(x, -))).second;
ans = max(ans, ask(root[tmp], , id[tmp].size()-, t));
}
return col[ans];
} int main()
{
scanf("%d", &n);
for(int i = , x, y, v;i < n;i++)
{
scanf("%d %d %d", &x, &y, &v);
addedge(x, y, v);
addedge(y, x, v);
}
de[]=;
dfs();
for(int j = ;j <= ;j++)
{
for(int i = ;i <= n;i++)
{
fa[i][j] = fa[fa[i][j-]][j-];
dis[i][j] = dis[fa[i][j-]][j-] + dis[i][j-];
}
}
totn = ;
hev = ;
lef_now = n;
dfs1(, );
solve(hev);
scanf("%d", &m);
int x,cmd,len,tot=;
while(m--)
{
scanf("%d%d",&cmd,&x);
if(cmd==)
{
scanf("%d%d",&len,&col[++tot]);
com_change(x,len,tot);
}
else printf("%d\n",com_ask(x));
}
return ;
}
LWDB的更多相关文章
- 【Codeforces】Gym100633 D. LWDB
题解 点分治,然后每个点上挂着一个距离不超过\(a_{i}\)的颜色改成\(c\) 用一个单调栈维护距离单调递减,每次查询在每个包括这个点的分治中心的单调栈上二分,找到修改最靠前的颜色作为这个点的颜色 ...
- XV Open Cup named after E.V. Pankratiev. GP of Tatarstan
A. Survival Route 留坑. B. Dispersed parentheses $f[i][j][k]$表示长度为$i$,未匹配的左括号数为$j$,最多的未匹配左括号数为$k$的方案数. ...
- 转weblogic 10.3新建域
一.安装前准备 1.解决linux中文乱码问题 修改/etc/sysconfig/i18n文件 #LANG="en_US.UTF-8"#SUPPORTED="en_US. ...
- SQL Server优化技巧之SQL Server中的"MapReduce"
日常的OLTP环境中,有时会涉及到一些统计方面的SQL语句,这些语句可能消耗巨大,进而影响整体运行环境,这里我为大家介绍如何利用SQL Server中的”类MapReduce”方式,在特定的统计情形中 ...
- 写给iOS程序员的命令行使用秘籍
http://www.jianshu.com/p/44d3b8f713f2 Mac OS是Unix系统的分支,有着强大的命令行功能.很多事情在命令行下处理会事半功倍,所以我就iOS程序员可能会用到的功 ...
- VxWorks 引导程序
前言:vxworks 的一些文件,如 usrconfig.c 在 config,comp目录中均有出现,因编译方式而选择某一个文件,命令行方式采用 config 目录文件,tornado 图形界面配置 ...
随机推荐
- python(9)- python基础知识刷题
1. 执行 Python 脚本的两种方式 交互方式:命令行 Windows操作系统下,快捷键cmd,输入“python”启动交互式python解释器. 文件方式:python文件 2. 简述位.字 ...
- Linux安装Java/Maven
所需文件:jdk 下载 安装Java INSTALL_PATH=/opt/soft TAR_FILE=/mnt/d/resources/soft/jdk-8u152-linux-x64.tar.gz ...
- 使用Caffe完成图像目标检测 和 caffe 全卷积网络
一.[用Python学习Caffe]2. 使用Caffe完成图像目标检测 标签: pythoncaffe深度学习目标检测ssd 2017-06-22 22:08 207人阅读 评论(0) 收藏 举报 ...
- JMeter中使用Put请求方式请求接口
前言 现在有如下接口,是以PUT的方式请求的: 请求URL:IP+Port+/api/v1/apps/{appId} 请求参数: 参数名 必选 类型 nameCn 是 string nameEn 是 ...
- 使用 lstat 函数获取文件信息
前言 在之前的文章中,描述过如何用 fcntl 函数改变文件的状态标记.但,文件还有很多信息,如文件类型,权限设置,设备编号,访问时间等等.如果要获取这些信息,则使用函数 lstat 可以轻松达到这个 ...
- 最新wap手机agent
名称 agent 铃声格式 和弦数 数据量 删除 LGE-CU8080 LGE-CU8080/1.0 UP.Browser/4.1.26l UP.Link/5.1.2.9 pmd2.0 40 ...
- EasyDarwin Streaming Server对Task的调用方法
我们在EasyDarwin流媒体服务器的二次开发过程中,经常会需要定义自己的Task类,例如在EasyDarwin中,RTSPSessioin.HTTPSession.RTCPTask等,都是Task ...
- java设计模式之综述
一.什么是设计模式 设计模式是一套被反复使用的.多数人知晓的.经过分类编目的.代码设计经验的总结.使用设计模式是为了重用代码.让代码更容易被他人理解.保证代码可靠性. 毫无疑问,设计模式于己于他人于系 ...
- [SCOI2009] 最长距离
题目描述 windy有一块矩形土地,被分为 NM 块 11 的小格子. 有的格子含有障碍物. 如果从格子A可以走到格子B,那么两个格子的距离就为两个格子中心的欧几里德距离. 如果从格子A不可以走到格子 ...
- 【C++基础学习】类型声明
1.初始化 在C++中,初始化与赋值操作是完全不同的两个操作.初始化不是赋值,初始化的含义是创建变量时赋予其一个初始值,而赋值的含义是把对象的当前值擦除,而以一个新值来代替. 初始化的方式有: 1 i ...