这场CF是以unrated选手的身份参加的,所以打的时候节奏比较慢,然而尽管慢,总共还是爆了5发OJ(A题忘判左界大于右界*1,B题看错d的含义*2 ,出题人真是的,样例给那么水 ,C题写了一个显然(是错的)的DP*2,D题反倒*0)
等到开E题时已经只剩30分钟,看到群里Claris说
辣鸡蒟蒻Claris(524936896) 0:34:20
我啥都不维护
辣鸡蒟蒻Claris(524936896) 0:34:30
只需要维护任意一条直径,这个随便做,
辣鸡蒟蒻Claris(524936896) 0:35:17
就是你随意维护一条直径
辣鸡蒟蒻Claris(524936896) 0:35:20
就知道一个端点x
辣鸡蒟蒻Claris(524936896) 0:35:27
就是问离x最远的点有几个
辣鸡蒟蒻Claris(524936896) 0:35:31
点分治随便做
感觉好奇妙,根本不会啊。可是知道离x最远的点有几个,并不等同于知道有几个点是直径,可能是我理解错了。
官方题解看的一知半解,大致是维护直径中点(然而也可能是“中边”,可以通过加虚点来解决)。
然后去扒fatego的代码,好像是用线段树维护dfn序,区间(l,r)表示有多少个点离中心距离等于半径,似乎很短,不到100行。
#include<cstdio>
const int N=300005;
int n,a[N],i,c,d,j,x;
struct segtree{
struct node{
int l,r,m,s;
bool v;
}t[N*18];
void build(int i,int l,int r){
t[i].l=l;
t[i].r=r;
t[i].m=(l+r)>>1;
if(l^r)build(i<<1,l,t[i].m),build(i<<1|1,t[i].m+1,r);
}
inline void pushdown(int i){
//printf("%d\n",i);
if(t[i].v)t[i<<1].v=t[i<<1|1].v=1,t[i].v=0,t[i<<1].s=t[i<<1|1].s=0;
}
void set(int i,int l,int r){
if(t[i].l==l && t[i].r==r){
t[i].s=0;
t[i].v=1;
}else{
pushdown(i);
if(l>t[i].m)set(i<<1|1,l,r);
else if(r<=t[i].m)set(i<<1,l,r);
else set(i<<1,l,t[i].m),set(i<<1|1,t[i].m+1,r);
t[i].s=t[i<<1].s+t[i<<1|1].s;
}
}
inline void add(int l){
for(register int i=1;t[i].l;i=i<<1|(l>t[i].m))pushdown(i),++t[i].s;
}
}tt;
inline void swap(int&a,int&b){
register int t=a;a=b;b=t;
}
struct tree{
static const int M=N<<1;
struct edge{
int to,next;
}e[N<<2];
int h[M],id[M],right[M],sz[M],ma[M],top[M],dad[M],xb,dep[M],dfn[M],sl;
inline void addedge(int u,int v){
e[++xb]=(edge){v,h[u]};
h[u]=xb;
}
void dfs(int x,int fa){
sz[x]=1;
dep[x]=dep[dad[dfn[id[x]=++xb]=x]=fa]+1;
for(int i=h[x];i;i=e[i].next)
if(e[i].to!=fa){
dfs(e[i].to,x);
sz[x]+=sz[e[i].to];
if(sz[e[i].to]>sz[ma[x]])ma[x]=e[i].to;
}
right[x]=xb;
}
inline void prepare(){
dfs(1,xb=0);
for(register int i=1;i<n<<1;++i)
top[dfn[i]]=dfn[i]==ma[dad[dfn[i]]]?top[dad[dfn[i]]]:dfn[i];
}
inline int dis(int a,int b){
register int u=a,v=b;
while(top[u]^top[v]){
if(dep[top[u]]>dep[top[v]])swap(u,v);
v=dad[sl=top[v]];
}
if(dep[u]>dep[v])swap(u,v);
if(u^v)sl=ma[u];
return dep[a]+dep[b]-(dep[u]<<1);
}
}t;
int main(){
scanf("%d",&n);
for(++n,i=2;i<=n;++i)scanf("%d",a+i),t.addedge(a[i],n+i-1),t.addedge(n+i-1,i);
t.prepare();
tt.build(1,c=1,n*2-1);
tt.add(1);
for(i=2;i<=n;++i){
j=t.dis(c,i);
x=t.sl;
if(j>d){
if(j==t.dep[i]-t.dep[c]){
tt.set(1,t.id[x],t.right[x]);
c=x;
}else{
if(1<t.id[c])tt.set(1,1,t.id[c]-1);
if(t.right[c]+1<n<<1)tt.set(1,t.right[c]+1,(n<<1)-1);
c=t.dad[c];
}
j=d=t.dis(c,i);
}
if(d==j)tt.add(t.id[i]);
printf("%d\n",tt.t[1].s);
}
return 0;
}