银河战舰 [启发式合并+dp]
题面
思路
我们首先考虑传统的链上LIS做法:保存每个长度的LIS末端的最小值,二分查找
那么这道题其实就只是搬到树上来做了而已
我们考虑一个节点,假设它的儿子已经处理完毕了
那么我们选择LIS最长的儿子进行继承,然后每次和其它儿子合并
合并的方法就是逐位取min值
合并完了以后再把当前节点二分一下加进去
记录的LIS,因为本题中的路径可以往根走也可以往叶子走,所以我们要记录从根开始的LIS和从叶子开始的LIS(其实就是从根开始的LDS),维护方法是一样的
然后对于答案的统计,有三种情况:
以当前节点为一端,向下递增
以当前节点为一端,向下递减
以当前子树中某一个点为一端,经过当前点,进入当前子树另一个点
前两种直接算就好了,第三种可以枚举LIS的长度,在LDS数组上面二分
合并的时候用启发式合并,这样总复杂度是$O(n\log^2 n)$
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cassert>
#include<vector>
using namespace std;
inline int read(){
int re=0,flag=1;char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') flag=-1;
ch=getchar();
}
while(isdigit(ch)) re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
return re*flag;
}
int n,a[200010],dep[200010],cnte,ans=0;
int first[200010];
vector<int>f[200010],g[200010];
struct edge{
int to,next;
}e[400010];
inline void add(int u,int v){
e[++cnte]=(edge){v,first[u]};first[u]=cnte;
e[++cnte]=(edge){u,first[v]};first[v]=cnte;
}
bool anti(int l,int r){return l>r;}
void fg(int u,int v){
int i,pos,su=f[u].size(),sv=g[v].size();
if(su<sv){
for(i=0;i<su;i++){
pos=lower_bound(g[v].begin(),g[v].end(),f[u][i],anti)-g[v].begin();
ans=max(ans,i+1+pos);
}
}
else{
for(i=0;i<sv;i++){
pos=lower_bound(f[u].begin(),f[u].end(),g[v][i])-f[u].begin();
ans=max(ans,i+1+pos);
}
}
}
void only(int u,int v){
int l1,l2;
l1=lower_bound(f[u].begin(),f[u].end(),a[u])-f[u].begin();
l2=lower_bound(g[v].begin(),g[v].end(),a[u],anti)-g[v].begin();
ans=max(ans,l1+l2+1);
l1=lower_bound(f[v].begin(),f[v].end(),a[u])-f[v].begin();
l2=lower_bound(g[u].begin(),g[u].end(),a[u],anti)-g[u].begin();
ans=max(ans,l1+l2+1);
}
void merge(int u,int v){
int i,su,sv;
su=f[u].size();sv=f[v].size();
if(su>=sv) for(i=0;i<sv;i++) f[u][i]=min(f[u][i],f[v][i]);
else{
for(i=0;i<su;i++) f[v][i]=min(f[u][i],f[v][i]);
swap(f[u],f[v]);
}
su=g[u].size();sv=g[v].size();
if(su>=sv) for(i=0;i<sv;i++) g[u][i]=max(g[u][i],g[v][i]);
else{
for(i=0;i<su;i++) g[v][i]=max(g[u][i],g[v][i]);
swap(g[u],g[v]);
}
}
void join(int u){
int len,pos;
len=f[u].size();
pos=upper_bound(f[u].begin(),f[u].end(),a[u])-f[u].begin();
if(!len||!pos||a[u]>f[u][pos-1]){
if(pos==len) f[u].push_back(a[u]);
else f[u][pos]=a[u];
}
len=g[u].size();
pos=upper_bound(g[u].begin(),g[u].end(),a[u],anti)-g[u].begin();
if(!len||!pos||a[u]<g[u][pos-1]){
if(pos==len) g[u].push_back(a[u]);
else g[u][pos]=a[u];
}
}
void dfs(int u,int f){
int i,v;
for(i=first[u];~i;i=e[i].next){
v=e[i].to;if(v==f) continue;
dfs(v,u);
only(u,v);
fg(u,v);fg(v,u);
merge(u,v);
}
join(u);
}
int main(){
memset(first,-1,sizeof(first));
n=read();int i,t1,t2;
for(i=1;i<=n;i++) a[i]=read();
for(i=1;i<n;i++){
t1=read();t2=read();
add(t1,t2);
}
dfs(1,0);
printf("%d\n",ans);
}
银河战舰 [启发式合并+dp]的更多相关文章
- Codeforces 1455G - Forbidden Value(map 启发式合并+DP)
Codeforces 题面传送门 & 洛谷题面传送门 首先这个 if 与 end 配对的结构显然形成一个树形结构,考虑把这棵树建出来,于是这个程序的结构就变为,对树进行一遍 DFS,到达某个节 ...
- ARC 086 E - Smuggling Marbles(dp + 启发式合并)
题意 Sunke 有一棵 \(N + 1\) 个点的树,其中 \(0\) 为根,每个点上有 \(0\) 或 \(1\) 个石子, Sunke 会不停的进行如下操作直至整棵树没有石子 : 把 \(0\) ...
- [2016北京集训试题7]thr-[树形dp+树链剖分+启发式合并]
Description Solution 神仙操作orz. 首先看数据范围,显然不可能是O(n2)的.(即绝对不是枚举那么简单的),我们考虑dp. 定义f(x,k)为以x为根的子树中与x距离为k的节点 ...
- AtCoder AGC007E Shik and Travel (二分、DP、启发式合并)
题目链接 https://atcoder.jp/contests/agc007/tasks/agc007_e 题解 首先有个很朴素的想法是,二分答案\(mid\)后使用可行性DP, 设\(dp[u][ ...
- P5979 [PA2014]Druzyny dp 分治 线段树 分类讨论 启发式合并
LINK:Druzyny 这题研究了一下午 终于搞懂了. \(n^2\)的dp很容易得到. 考虑优化.又有大于的限制又有小于的限制这个非常难处理. 不过可以得到在限制人数上界的情况下能转移到的最远端点 ...
- [多校 NOIP 联合模拟 20201130 T4] ZZH 的旅行(斜率优化dp,启发式合并,平衡树)
题面 题目背景 因为出题人天天被 ZZH(Zou ZHen) 吊打,所以这场比赛的题目中出现了 ZZH . 简要题面 数据范围 题解 (笔者写两个log的平衡树和启发式合并卡过的,不足为奇) 首先,很 ...
- 2018.10.14 loj#516. DP 一般看规律(启发式合并)
传送门 注意到一种颜色改了之后就不能改回去了. 因此可以启发式合并. 每次把小的合并给大的. 这样每个数最多被合并logloglog次. 如果维护一棵比较下标的平衡树的话,对于答案有贡献的就是每个数与 ...
- BZOJ4919 [Lydsy1706月赛]大根堆 【dp + 启发式合并】
题目链接 BZOJ4919 题解 链上的\(LIS\)维护一个数组\(f[i]\)表示长度为\(i\)的\(LIS\)最小的结尾大小 我们可以用\(multiset\)来维护这个数组,子树互不影响,启 ...
- loj516 DP一般看规律(set启发式合并)
题目: https://loj.ac/problem/516 分析: 每次将一个颜色更改为另一个颜色相当于将两个集合合并 然后对于答案的更新,一个点插入到一个集合中,那么可能更新答案的就是其前驱节点或 ...
随机推荐
- CSS3一些特殊属性
(一)-webkit-tap-highlight-color 这个属性只用于iOS (iPhone和iPad).当你点击一个链接或者通过Javascript定义的可点击元素的时候,它就 ...
- 关于Pycharm基本操作笔记
创建 project(工程,译音:破拽科特) 1.Create New project(创建一个新的工程,译音:科瑞特 纽 破摘科特) 2.pure python(纯派森,译音:皮忧儿 派森) 3.l ...
- thinkphp5,单图,多图,上传
/** * 上传单图 */ function upload($path, $filename) { $file = request()->file($filename); $info = $fi ...
- ES6笔记01-声明变量
ES6只有六种声明变量的方法:var命令和function命令,let和const命令,import命令和class命令.所以,ES6一共有6种声明变量的方法. const声明一个只读的常量.一旦声明 ...
- python——“/”运算符和“//”运算符的区别
首先先看单斜杆的用法:举几个例子 >>> print(5/3),type(5/3)1.6666666666666667(None, <class 'float'>) &g ...
- [Bzoj4408]神秘数(主席树)
Description 一个可重复数字集合S的神秘数定义为最小的不能被S的子集的和表示的正整数. 例如S={1,1,1,4,13}, 1 = 1 2 = 1+1 3 = 1+1+1 4 = 4 5 = ...
- 用for循环计算(1-3+5-7...99)的结果(两种方法)
1) sum=0 count=1 for i in range(1,100,2): if count % 2==0: sum = sum - i else: sum = sum + i count = ...
- Kotlin 二分法算法游戏--猜价格
本人最新想学习算法,算法是提高程序性能的关键! 程序就是数据结构和算法! 写了一个二分法的游戏,供大家参考: 当然,语言基于kotlin import java.util.* /** * Create ...
- 《Cracking the Coding Interview》——第13章:C和C++——题目5
2014-04-25 19:59 题目:C的关键字volatile有什么用? 解法:搞硬件设计的人好像更关注这个关键字.volatile本身是易变的意思,应该和persistent有反义词关系吧.说一 ...
- 《Cracking the Coding Interview》——第12章:测试——题目6
2014-04-25 00:53 题目:你要如何测试一个分布式银行系统的ATM机? 解法:ATM是Automatic Teller Machine,取钱的.我想了半天,没找到什么很清晰的思路,也许是因 ...