题目

最近比较懒，题目描述都直接截图了。

题目大意

给你一棵树，还有树上的几条路径，一条路径上的点到路径上其它任意点的代价为111。然后是一堆询问，问从一个点到另一个点的最小代价。

思路

一开始做这题时，就自然地往链剖或倍增方面想。
链剖想不出，于是就想倍增。
对于每个点，预处理出它花111的代价能到达的最远祖先。
然后倍增一下~
询问的时候，就将两个点的LCALCALCA求出，然后两个点往上跳，直到再跳就超过LCALCALCA为止。
所以问题就转化成了求两个点是否被一条路径直接连通，如果是，就加111，否则加222。
想不出来……
最后只能打暴力：将一个点作为根，然后粗暴地处理……

正解

事实证明我前面所思考的是正确的。
对于这个问题，其实可以转化成：是否有一条路径的两个端点分别在这两个点的子树内……
好像是很显然的，可我为什么想不到……
于是我们就可以求出它们的dfndfndfn序，这就成了一个平面上的问题。
每一条路径可以看做一个点，每一个询问可以看作一个矩形，求这个矩形当中是否有点出现。
然后就是一个简单的扫描线和树状数组了。

代码

有一点需要补充一下：
有的同学反映，这题递归会爆栈，所以要打人工栈。
反正我是没有爆栈……
后来听说xzb大爷想出一个不用dfs求dfn序的方法。
首先，我们可以通过拓扑排序之类的自底向上递推求出它们的siz等信息。（这题比较良心，由于每个点的父亲编号小于它，所以直接从后往前扫就好了。）
然后自顶向下（可以用bfs，这题直接从前往后扫就好了），
对于每个点分两种情况：

1. 它是父亲的第一个儿子。那么它的dfndfndfn即为父亲加一。
2. 否则，就是父亲上一个儿子的dfndfndfn加它的sizsizsiz。

然后就可以简单地求出dfndfndfn了。
正常的bfs中，上一个点就是父亲的上一个儿子。
不过这题不用bfs，直接记录就好了。
在此膜拜xzb

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 200000
#define M N
#define Q N
int n;
int last[N+1];
int siz[N+1],dfn[N+1];
int fa[N+1][19],dep[N+1];
inline void get_dfn(){//这是一种不需要dfs求dfn的方法，解释见上dfn[1]=1;for (int i=2;i<=n;++i){if (last[fa[i][0]])dfn[i]=dfn[last[fa[i][0]]]+siz[last[fa[i][0]]];elsedfn[i]=dfn[fa[i][0]]+1;last[fa[i][0]]=i;}
}
inline int LCA(int u,int v){if (dep[u]<dep[v])swap(u,v);for (int k=dep[u]-dep[v],i=0;k;k>>=1,++i)if (k&1)u=fa[u][i];if (u==v)return u;for (int i=17;i>=0;--i)if (fa[u][i]!=fa[v][i]){u=fa[u][i];v=fa[v][i];}return fa[u][0];
}
int m;
int up[N+1][19],dep2[N+1];
struct Oper{int x,l,r;int ty;int num;
} o[M*2+Q*2+1];
int cnt;
bool cmp(const Oper &a,const Oper &b){return a.x<b.x || a.x==b.x && a.ty<b.ty;
}
int t[N+1];
#define lowbit(x) ((x)&(-x))
inline void change(int x){do{t[x]++;x+=lowbit(x);}while (x<=n);
}
inline int query(int x){int res=0;while (x){res+=t[x];x-=lowbit(x);}return res;
}
int ans[Q+1],cha[Q+1];
int main(){freopen("car.in","r",stdin);freopen("car.out","w",stdout);scanf("%d",&n);for (int i=2;i<=n;++i){scanf("%d",&fa[i][0]);e[++ne]={i,last[fa[i][0]]};last[fa[i][0]]=e+ne;}for (int i=1;i<=n;++i)dep[i]=dep[fa[i][0]]+1;//因为题目说每个点父亲的编号一定小于它，所以直接递推就好，下同。for (int i=n;i>=1;--i)siz[fa[i][0]]+=++siz[i];get_dfn(1);for (int i=1;i<=17;++i)for (int j=1;j<=n;++j)fa[j][i]=fa[fa[j][i-1]][i-1];//倍增处理for (int i=1;i<=n;++i)up[i][0]=i;//up[i][j]表示花费2^j的代价最远能到达的祖先scanf("%d",&m);for (int i=1;i<=m;++i){int u,v;scanf("%d%d",&u,&v);if (dfn[u]>dfn[v])swap(u,v);//处理花费1代价能到达的最远祖先int lca=LCA(u,v);if (dep[lca]<dep[up[u][0]])up[u][0]=lca;if (dep[lca]<dep[up[v][0]])up[v][0]=lca;//加入操作列表中，相当于点（因为它有对称性，所以正着反着都加进去）o[++cnt]={dfn[u],dfn[v],dfn[v],0,0};o[++cnt]={dfn[v],dfn[u],dfn[u],0,0};}for (int i=n;i>=1;--i)if (dep[up[i][0]]<dep[up[fa[i][0]][0]])up[fa[i][0]][0]=up[i][0];//从底向上递推出每个点花费1代价到达的最远祖先for (int i=1;i<=18;++i)for (int j=1;j<=n;++j)up[j][i]=up[up[j][i-1]][i-1];int q;scanf("%d",&q);for (int i=1;i<=q;++i){int u,v;scanf("%d%d",&u,&v);if (u==v)continue;if (up[u][18]!=up[v][18]){//判断两个点是否可以互相到达（因为N=200000，2^18>N，所以这是它能到达的最远祖先，因为这是一棵树，所以如果可达到的最远祖先不一样，那么他们就不能互相到达）ans[i]=-1;continue;}int lca=LCA(u,v);if (lca==u || lca==v){//在同一条链上的情况if (lca==u)swap(u,v);for (int j=17;j>=0;--j)if (dep[up[u][j]]>dep[v])u=up[u][j],ans[i]+=1<<j;//跳到LCA下的最远点ans[i]++;continue;}for (int j=17;j>=0;--j)if (dep[up[u][j]]>dep[lca])u=up[u][j],ans[i]+=1<<j;//跳到LCA下的最远点，下同for (int j=17;j>=0;--j)if (dep[up[v][j]]>dep[lca])v=up[v][j],ans[i]+=1<<j;//加入操作o[++cnt]={dfn[u]-1,dfn[v],dfn[v]+siz[v]-1,1,i};o[++cnt]={dfn[u]+siz[u],dfn[v],dfn[v]+siz[v]-1,-1,i};}sort(o+1,o+cnt+1,cmp);for (int i=1;i<=cnt;++i)//扫描线求出二维数点问题if (o[i].ty==0)change(o[i].l);else if (o[i].ty==1)cha[o[i].num]=query(o[i].r)-query(o[i].l-1);elseans[o[i].num]+=((query(o[i].r)-query(o[i].l-1))-cha[o[i].num]?1:2);for (int i=1;i<=q;++i)printf("%d\n",ans[i]);return 0;
}

总结

以后，看见有关树的题目，如果普通的倍增、链剖都做不出来，就得要往dfndfndfn序方面想……