洛谷 P3258 [JLOI2014]松鼠的新家 树上差分

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

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

输出格式:

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

输入输出样例
输入样例#1:
5
1 4 5 3 2
1 2
2 4
2 3
4 5
输出样例#1:
1
2
1
2
1
说明
2<= n <=300000

分析:
每次相当于给路径上除了起点的每个点都加 1 1 ,但是中间没有询问,就不用树剖,而且还会T,所以直接上树上差分。差分后的树每个点的子树和就是当前点原来的值。所以,给路径(l,r)都加上 1 1 ,相当于给l r r 1 lca(l,r) l c a ( l , r ) fa[lca(l,r)] f a [ l c a ( l , r ) ] 都减去 1 1 。然后给差分树求子树和即可。因为除了第一步的起点和路径终点都要少放一个,所以除了a1外所有点减一即可。

代码:

#include <iostream>
#include <cstdio>
#include <cmath>

const int maxn=3e5+7;

using namespace std;

int x,y,n,cnt;
int a[maxn],sum[maxn],dep[maxn],ls[maxn];
int f[maxn][20];

struct edge{
    int y,next;
}g[maxn*2];

void add(int x,int y)
{
    g[++cnt]=(edge){y,ls[x]};
    ls[x]=cnt;
}

void dfs(int x,int fa)
{
    f[x][0]=fa;
    for (int i=ls[x];i>0;i=g[i].next)
    {
        int y=g[i].y;
        if (y==fa) continue;
        dep[y]=dep[x]+1;
        dfs(y,x);
        sum[x]+=sum[y];
    }
}

int lca(int x,int y)
{
    if (dep[x]>dep[y]) swap(x,y);
    int d=dep[y]-dep[x],t=1<<19,k=19;
    while (d)
    {
        if (t<=d)
        {
            d-=t;
            y=f[y][k];
        }
        t/=2;
        k--;
    }
    if (x==y) return x;
    k=19;
    while (k>=0)
    {
        if (f[x][k]!=f[y][k])
        {
            x=f[x][k];
            y=f[y][k];
        }
        k--;
    }
    return f[x][0];
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    for (int i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        add(x,y);
        add(y,x);
    }
    dep[1]=1;
    dfs(1,0);
    for (int j=1;j<20;j++)
    {
        for (int i=1;i<=n;i++)
        {
            f[i][j]=f[f[i][j-1]][j-1];
        }
    }
    for (int i=1;i<n;i++)
    {
        int d=lca(a[i],a[i+1]);
        sum[a[i]]+=1;
        sum[a[i+1]]+=1;
        sum[d]-=1;
        if (d!=1) sum[f[d][0]]-=1;
    }
    dfs(1,0);
    for (int i=1;i<=n;i++)
    {
        if (i==a[1]) printf("%d\n",sum[i]);
                else printf("%d\n",sum[i]-1);
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值