UVALive 7148 LRIP【树分治+线段树】
题意就是要求一棵树上的最长不下降序列,同时不下降序列的最小值与最大值不超过D。
做法是树分治+线段树,假设树根是x,y是其当前需要处理的子树,对于子树y,需要处理出两个数组MN,MX,MN[i]表示以x为第一个数字的不下降子序列中第i个数的最小值,MX[i]表示以x为第一个数字的不上升子序列中第i个数的最大值。如果当前子树有一个以x为首的不下降序列,那么我们就需要在之前处理的子树中找一条以x为首的满足约束条件不上升序列,可以用线段树来查询。同时每做完一颗子树的时候,用MN,MX对线段树进行更新。对于不经过x的情况可以递归下去处理。
代码:
#include<cstdio>
#include<algorithm>
#include<cstring>
#define N 1000010
using namespace std;
int dp,p[N],s[N],pre[N],tt[N],flag[N],father[N],tmproot,n,d,ttt;
int i,v[N],a,b,mn,mx,MN[N],MX[N],U,V,ans;
int l[N],r[N],smn[N],smx[N],vsmn[N],vsmx[N];
void build(int x,int a,int b)
{
int m;
l[x]=a;r[x]=b;
if (x>ttt) ttt=x;
if (b-a>)
{
m=(a+b)>>;
build(*x,a,m);
build(*x+,m,b);
}
}
void clean(int x)
{
if (vsmn[x])
{
smn[x]=;
if (*x<=ttt)
vsmn[*x]=;
if (*x+<=ttt)
vsmn[*x+]=;
vsmn[x]=;
} if (vsmx[x])
{
smx[x]=;
if (*x<=ttt)
vsmx[*x]=;
if (*x+<=ttt)
vsmx[*x+]=;
vsmx[x]=;
}
}
void change(int x,int a,int b,int c,int typ)
{
int m;
clean(x);
if ((a<=l[x])&&(r[x]<=b))
{
if (typ==)
smn[x]=max(smn[x],c);
else
smx[x]=max(smx[x],c);
return;
}
m=(l[x]+r[x])>>;
if (a<m) change(*x,a,b,c,typ);
if (m<b) change(*x+,a,b,c,typ);
clean(*x);clean(*x+);
if (!typ)
smn[x]=max(smn[*x],smn[*x+]);
else
smx[x]=max(smx[*x],smx[*x+]);
}
int query(int x,int a,int b,int typ)
{
int m,ans=;
clean(x);
if ((a<=l[x]&&(r[x]<=b)))
{
if (typ==)
return smn[x];
else
return smx[x];
}
m=(l[x]+r[x])>>;
if (a<m) ans=max(ans,query(*x,a,b,typ));
if (m<b) ans=max(ans,query(*x+,a,b,typ));
return ans;
}
void link(int x,int y)
{
dp++;pre[dp]=p[x];p[x]=dp;tt[dp]=y;
}
int getroot(int x,int fa,int sum)
{
int i,f=;
i=p[x];
father[x]=fa;
s[x]=;
while (i)
{
if ((tt[i]!=fa)&&(!flag[tt[i]]))
{
getroot(tt[i],x,sum);
s[x]=s[x]+s[tt[i]];
if (s[tt[i]]>sum/) f=;
}
i=pre[i];
}
if (sum-s[x]>sum/) f=;
if (!f) tmproot=x;
}
void dfs(int x,int fa,int mn,int mx)
{
int i;
i=p[x];
if (v[x]>=v[fa])
{
if (mn)
{
mn++;
MN[mn]=min(MN[mn],v[x]);
}
if (v[x]>v[fa])
mx=;
}
if (v[x]<=v[fa])
{
if (mx)
{
mx++;
MX[mx]=max(MX[mx],v[x]);
}
if (v[x]<v[fa])
mn=;
}
U=max(U,mn);
V=max(V,mx);
while (i)
{
if ((tt[i]!=fa)&&(!flag[tt[i]]))
dfs(tt[i],x,mn,mx);
i=pre[i];
}
}
void work(int x,int sum)
{
int i,j,k,root;
getroot(x,,sum);
root=tmproot;
i=p[root];
flag[root]=;
while (i)
{
if (flag[tt[i]]==)
{
if (tt[i]==father[root])
work(tt[i],sum-s[root]);
else
work(tt[i],s[tt[i]]);
}
i=pre[i];
}
//------------------------------
i=p[root];
vsmn[]=;
vsmx[]=;
while (i)
{
if (!flag[tt[i]])
{
for (j=;j<=U;j++)
MN[j]=0x37373737;
for (j=;j<=V;j++)
MX[j]=;
U=;V=;
dfs(tt[i],root,,);
for (j=;j<=U;j++)
{
if (MN[j]-d<=v[root])
ans=max(ans,j);
k=MN[j]-d;
if (k<) k=;
if (k<=v[root])
ans=max(ans,j+query(,k-,v[root],)-);
}
for (j=;j<=V;j++)
{
if (MX[j]+d>=v[root])
ans=max(ans,j);
k=MX[j]+d;
if (k>) k=;
if (k>=v[root])
ans=max(ans,j+query(,v[root]-,k,)-);
} for (j=;j<=U;j++)
change(,MN[j]-,MN[j],j,);
for (j=;j<=V;j++)
change(,MX[j]-,MX[j],j,);
}
i=pre[i];
}
flag[root]=;
}
int main()
{
build(,,);
for (i=;i<=;i++)
MN[i]=0x37373737;
int test,ii=;
scanf("%d",&test);
while (test)
{
test--;
ii++;
dp=;memset(p,,sizeof(p));
scanf("%d%d",&n,&d);
for (i=;i<=n;i++)
scanf("%d",&v[i]);
for (i=;i<=n-;i++)
{
scanf("%d%d",&a,&b);
link(a,b);
link(b,a);
}
ans=;
work(,n);
printf("Case #%d: %d\n",ii,ans);
}
}
UVALive 7148 LRIP【树分治+线段树】的更多相关文章
- 【BZOJ4372】烁烁的游戏 动态树分治+线段树
[BZOJ4372]烁烁的游戏 Description 背景:烁烁很喜欢爬树,这吓坏了树上的皮皮鼠.题意:给定一颗n个节点的树,边权均为1,初始树上没有皮皮鼠.烁烁他每次会跳到一个节点u,把周围与他距 ...
- LOJ#6463 AK YOI 树分治+线段树合并
传送门 既然是树上路径统计问题,不难想到要使用树分治,这里以点分治为例 由点分治的性质,每层只需要考虑经过重心的路径 因为需要维护路径长度在一定范围内的最大权值和,所以要用一个数据结构维护一下到根节点 ...
- 【BZOJ3730】震波 动态树分治+线段树
[BZOJ3730]震波 Description 在一片土地上有N个城市,通过N-1条无向边互相连接,形成一棵树的结构,相邻两个城市的距离为1,其中第i个城市的价值为value[i].不幸的是,这片土 ...
- 7.18 NOI模拟赛 因懒无名 线段树分治 线段树维护直径
LINK:因懒无名 20分显然有\(n\cdot q\)的暴力. 还有20分 每次只询问一种颜色的直径不过带修改. 容易想到利用线段树维护直径就可以解决了. 当然也可以进行线段树分治 每种颜色存一下直 ...
- bzoj3730 [震波][动态树分治+线段树+LCA]
震波 Time Limit: 15 Sec Memory Limit: 256 MBSubmit: 1573 Solved: 358[Submit][Status][Discuss] Descri ...
- BZOJ4317Atm的树&BZOJ2051A Problem For Fun&BZOJ2117[2010国家集训队]Crash的旅游计划——二分答案+动态点分治(点分树套线段树/点分树+vector)
题目描述 Atm有一段时间在虐qtree的题目,于是,他满脑子都是tree,tree,tree…… 于是,一天晚上他梦到自己被关在了一个有根树中,每条路径都有边权,一个神秘的声音告诉他,每个点到其他的 ...
- 【loj6145】「2017 山东三轮集训 Day7」Easy 动态点分治+线段树
题目描述 给你一棵 $n$ 个点的树,边有边权.$m$ 次询问,每次给出 $l$ .$r$ .$x$ ,求 $\text{Min}_{i=l}^r\text{dis}(i,x)$ . $n,m\le ...
- 【bzoj4372】烁烁的游戏 动态点分治+线段树
题目描述 给一颗n个节点的树,边权均为1,初始点权均为0,m次操作:Q x:询问x的点权.M x d w:将树上与节点x距离不超过d的节点的点权均加上w. 输入 第一行两个正整数:n,m接下来的n-1 ...
- 【bzoj3730】震波 动态点分治+线段树
题目描述 在一片土地上有N个城市,通过N-1条无向边互相连接,形成一棵树的结构,相邻两个城市的距离为1,其中第i个城市的价值为value[i].不幸的是,这片土地常常发生地震,并且随着时代的发展,城市 ...
随机推荐
- Python之什么是dict
我们已经知道,list 和 tuple 可以用来表示顺序集合,例如,班里同学的名字: ['Adam', 'Lisa', 'Bart'] 或者考试的成绩列表: [95, 85, 59] 但是,要根据名字 ...
- css背景图片定位练习(一)
首先准备一张雪碧图,Like this 背景图片的定位方法有3种,比较常用的两种为 关键字:background-position: top left; (top/bottom/cennter/lef ...
- There has been an error processing your request magento
如果使用magento的过程中,出现以下页面: 说明出现了错误,但是亲,不用紧张,请根据"Error record number:xxxxxxxxx"的数字在网站根目录下的var/ ...
- Boolean强制转换
数据类型 转换为true的值 转换为false的值 Boolean true false String 任何非空字符串 空字符串 Number 任何非0数字值(包括无穷大) 0和NaN Object ...
- LeetCode H-Index II
原题链接在这里:https://leetcode.com/problems/h-index-ii/ 题目: Follow up for H-Index: What if the citations a ...
- RedHat6.6更新Centos6yum源
RedHat6.6更新Centos6yum源 一. 删除自带的RedHat6.6yum源 1- rpm -qa|grep yum|xargs rpm -e --nodeps(不检查依 ...
- 高性能Linux服务器 第10章 基于Linux服务器的性能分析与优化
高性能Linux服务器 第10章 基于Linux服务器的性能分析与优化 作为一名Linux系统管理员,最主要的工作是优化系统配置,使应用在系统上以最优的状态运行.但硬件问题.软件问题.网络环境等 ...
- Android 使用PullToRefresh实现下拉刷新和上拉加载(ExpandableListView)
PullToRefresh是一套实现非常好的下拉刷新库,它支持: 1.ListView 2.ExpandableListView 3.GridView 4.WebView 等多种常用的需要刷新的Vie ...
- iptables调试方法
iptables调试时,使用到raw表.ipt_LOG内核模块.日志记录在kern.log中. 具体的步骤如下: 1.准备ipt_LOG内核模块 modprobe ipt_LOG 2.使用raw表,加 ...
- javascript:void(0)和javascript:;的用法
一.JavaScript:void(0) 我们经常会使用到 javascript:void(0) 这样的代码,那么在 JavaScript 中 javascript:void(0) 代表的是什么意思呢 ...