「FJOI2018」领导集团问题
题意：给你一颗(n)个点的带点权有根树，选择一个点集(S)，使得点集中所有祖先的点权$le (子孙的点权，最大化)|S|$（出题人语死早...）
-
一个显然的(dp)
(dp_{i,j})代表子树(i)中选择的点集中最小的点权为(j)时的最大值。
-
一个朴素的转移需要维护一个后缀最大值
设为(f_{i,j}=maxlimits_{kge i} dp_{i,k})
-
一些显然的事实
- 后缀最大值中仅有最多(siz_v)（指子树大小）个权值不同的位置
- 后缀最大值单调不增
于是我们考虑如何合并两颗子树有限的后缀最大值
在不考虑关键点时，我们只需要对应位置两两加和就可以了，即
[dp_{i,j}=sum_vf_{v,j}
]
下面我们仅考虑关键点
如图所示，蓝色的关键点是子树(u)的关键点，黄色的是子树(v)的关键点，注意到关键点是左边位置的值与它相等。
考虑启发式合并，但是我们注意到关键点是互相影响的，就是黄影响蓝，蓝影响黄。
其实这样可以做，但是比较麻烦。
注意到一个事实，我们只需要求(f_{rt,0})，所以我们考虑维护后缀最大值(f)的差分数组。
这样我们合并两个关键点集合，直接对关键点进行相加就可以了，考虑直接用线段树合并维护。
然后我们还需要计算根(i)自己的贡献，设(i)的权值为(w_i)
考虑到根直接的贡献是使(dp_{i,w_i}=f_{i,w_i}+1)
然后从这个(w_i)往左边走，找到第一个(f_{i,p})比(f_{i,w_i})大(1)的位置，然后更新这一段的贡献，就是一个区间加。具体操作就是在(p)打个(-1)，在(w_i)打一个(1)
形象一点就是把这一段铺平了,就是那(dp)再回去贡献(f)
找第一个位置就是在差分意义下就是找一个非(0)位置，可以在线段树上进行二分操作找到。
Code:
#include <cstdio>
#include <cctype>
#include <vector>
#include <algorithm>
using std::min;
using std::max;
const int SIZE=1<<21;
char ibuf[SIZE],*iS,*iT;
//#define gc() (iS==iT?(iT=(iS=ibuf)+fread(ibuf,1,SIZE,stdin),iS==iT?EOF:*iS++):*iS++)
#define gc() getchar()
template <class T>
void read(T &x)
{
int f=0;x=0;char c=gc();
while(!isdigit(c)) f|=c=='-',c=gc();
while(isdigit(c)) x=x*10+c-'0',c=gc();
if(f) x=-x;
}
const int N=4e5+10;
int n,m,w[N],yuy[N];
int head[N],to[N],Next[N],cnt;
void add(int u,int v)
{
to[++cnt]=v,Next[cnt]=head[u],head[u]=cnt;
}
int root[N],ch[N*40][2],sum[N*40],tot;
#define ls ch[now][0]
#define rs ch[now][1]
void upt(int &now,int l,int r,int p,int d)
{
if(!now) now=++tot;
if(l==r)
{
sum[now]+=d;
return;
}
int mid=l+r>>1;
if(p<=mid) upt(ls,l,mid,p,d);
else upt(rs,mid+1,r,p,d);
sum[now]=sum[ls]+sum[rs];
}
int qry(int now,int l,int r,int p)
{
if(!now||!p) return 0;
if(l==r) return sum[now]?l:0;
int mid=l+r>>1;
if(p<=mid) return qry(ls,l,mid,p);
else
{
int ret=qry(rs,mid+1,r,p);
if(ret) return ret;
else return qry(ls,l,mid,p);
}
}
int Merge(int x,int y,int l,int r)
{
if(!x||!y) return x^y;
if(l==r)
{
sum[x]+=sum[y];
return x;
}
int mid=l+r>>1;
ch[x][0]=Merge(ch[x][0],ch[y][0],l,mid);
ch[x][1]=Merge(ch[x][1],ch[y][1],mid+1,r);
sum[x]=sum[ch[x][0]]+sum[ch[x][1]];
return x;
}
void dfs(int now)
{
for(int v,i=head[now];i;i=Next[i])
{
dfs(v=to[i]);
root[now]=Merge(root[now],root[v],0,m);
}
upt(root[now],1,m,w[now],1);
int pos=qry(root[now],1,m,w[now]-1);
if(pos) upt(root[now],1,m,pos,-1);
}
int main()
{
read(n);
for(int i=1;i<=n;i++) read(w[i]),yuy[i]=w[i];
std::sort(yuy+1,yuy+1+n);
m=std::unique(yuy+1,yuy+1+n)-yuy-1;
for(int i=1;i<=n;i++) w[i]=std::lower_bound(yuy+1,yuy+1+m,w[i])-yuy;
for(int p,i=2;i<=n;i++) read(p),add(p,i);
dfs(1);
printf("%d
",sum[root[1]]);
return 0;
}
2019.5.21