bzoj3631 松鼠的新家[LCA][树上差分]

3631: [JLOI2014]松鼠的新家

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1643 Solved: 776
[Submit][Status][Discuss]
Description

松鼠的新家是一棵树，前几天刚刚装修了新家，新家有n个房间，并且有n-1根树枝连接，每个房间都可以相互到达，且俩个房间之间的路线都是唯一的。天哪，他居然真的住在“树”上。松鼠想邀请小熊维尼前来参观，并且还指定一份参观指南，他希望维尼能够按照他的指南顺序，先去a1，再去a2，……，最后到an，去参观新家。
可是这样会导致维尼重复走很多房间，懒惰的维尼不听地推辞。可是松鼠告诉他，每走到一个房间，他就可以从房间拿一块糖果吃。维尼是个馋家伙，立马就答应了。
现在松鼠希望知道为了保证维尼有糖果吃，他需要在每一个房间各放至少多少个糖果。因为松鼠参观指南上的最后一个房间an是餐厅，餐厅里他准备了丰盛的大餐，所以当维尼在参观的最后到达餐厅时就不需要再拿糖果吃了。
Input

第一行一个整数n，表示房间个数
第二行n个整数，依次描述a1-an
接下来n-1行，每行两个整数x，y，表示标号x和y的两个房间之间有树枝相连。
Output

一共n行，第i行输出标号为i的房间至少需要放多少个糖果，才能让维尼有糖果吃。
Sample Input

5

1 4 5 3 2

1 2

2 4

2 3

4 5

Sample Output

1

2

1

2

1

HINT

2<= n <=300000

题意：对给定的树按给定的序列遍历树上的每一个点，求最后每个点被走了多少次。
分析：数据范围过大，线段树过比较艰难，因为是对一条链进行同一修改，考虑树上差分，即对两端点加1，对两端点lca和lca的父节点减一，对链上的点进行讨论可以知道这个方法可以解决问题。
然后我就把LCA写错了很多次……
嗯RMQ算法要定大一点因为是时间戳不是点啊！
f要先初始化为很大值不然取小值会错啊！
我的RMQ原来怎么会对啊！

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
const int maxn=3e5+5;
int fr[maxn*2],tov[maxn*2],des[maxn*2],f[maxn*4][20],tii[maxn*2],sig[maxn],val[maxn],tot[maxn],fa[maxn];
int n,cnt,sta,fin,a[maxn];
void addedge()
{
    scanf("%d %d",&sta,&fin);
    tov[++cnt]=fr[sta];fr[sta]=cnt;des[cnt]=fin;
    tov[++cnt]=fr[fin];fr[fin]=cnt;des[cnt]=sta;
}
void dfs(int u,int fath)
{
    fa[u]=fath;
    tii[++cnt]=u;sig[u]=cnt;f[cnt][0]=cnt;
    for(int i=fr[u];i;i=tov[i])if(des[i]!=fath){
    dfs(des[i],u);
    tii[++cnt]=u;
    f[cnt][0]=sig[u];
    }
}
int poww(int tmp)
{
    int tmp2=1;
    for(int i=1;i<=tmp;i++)tmp2*=2;
    return tmp2;
}
void bz()
{
    for(int i=1;i<=19;i++)for(int j=1;j<=2*n-2;j++)
        f[j][i]=min(f[j][i-1],f[j+poww(i-1)][i-1]);
}
int lca(int sta,int fin)
{
    int tmp1=sig[fin],tmp2=sig[sta];
    if(tmp1<tmp2)swap(tmp1,tmp2);
    int tmp=floor(log(tmp1-tmp2+1)/log(2));
    return tii[min(f[tmp2][tmp],f[tmp1+1-poww(tmp)][tmp])];
}
void init()
{
    scanf("%d",&n);
    memset(f,127,sizeof(f));
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<n;i++)addedge();
    cnt=0;dfs(a[1],0);
    bz();
    for(int i=2;i<=n;i++){
        int llca=lca(a[i-1],a[i]);
        val[llca]--;
        val[fa[llca]]--;
        val[a[i-1]]++;
        val[a[i]]++;
    }
}
void dfs2(int u,int fath)
{
    tot[u]+=val[u];
    for(int i=fr[u];i;i=tov[i])if(des[i]!=fath){
        dfs2(des[i],u);
        tot[u]+=tot[des[i]];
    }
}
void work()
{
    dfs2(a[1],0);
    for(int i=2;i<=n;i++)tot[a[i]]--;
    for(int i=1;i<=n;i++)
        printf("%d\n",tot[i]);
}
int main()
{
    //freopen("bzoj3631.in","r",stdin);
   // freopen("bzoj3631.out","w",stdout);
    init();
    work();
    return 0;
}