P5002 专心OI - 找祖先

传送门(洛谷)

这道题就是一个树上乘法原理

当有一个点是询问点p时:

s u m [ p ] + = 2 ∗ s i z e [ p ] − 1 , ( p , p ) 是 一 种 情 况 sum[p]+=2*size[p]-1,(p,p)是一种情况 sum[p]+=2size[p]1,(p,p)

当两个点都是询问点p的子树点的时候

设x有k棵子树
根据乘法原理,p有size[p]个点,则
s u m [ p ] + = ∑ i = 1 k ∑ j = 1 k s i z e [ s o n [ i ] ] ∗ s i z e [ s o n [ j ] ] sum[p]+=\sum_{i=1}^k\sum_{j=1}^ksize[son[i]]*size[son[j]] sum[p]+=i=1kj=1ksize[son[i]]size[son[j]]
= ∑ i = 1 k s i z e [ s o n [ i ] ] ∗ ( s i z e [ p ] − 1 ) =\sum_{i=1}^ksize[son[i]]*(size[p]-1) =i=1ksize[son[i]](size[p]1)
= ( s i z e [ p ] − 1 ) ∗ ( s i z e [ p ] − 1 ) =(size[p]-1)*(size[p]-1) =(size[p]1)(size[p]1)
但是如果是两个相同的点,则我们算重了一部分(子树大小不能乘以他自己的大小,也就是i=j时),应该减去:
∑ i = 1 k s i z e [ s o n [ i ] ] 2 \sum_{i=1}^ksize[son[i]]^2 i=1ksize[son[i]]2
所以总答案为:
a n s = 2 ∗ s i z e [ p ] − 1 + ( s i z e [ p ] − 1 ) ∗ ( s i z e [ p ] − 1 ) − ∑ i = 1 k s i z e [ s o n [ i ] ] 2 ans=2*size[p]-1+(size[p]-1)*(size[p]-1)-\sum_{i=1}^ksize[son[i]]^2 ans=2size[p]1+(size[p]1)(size[p]1)i=1ksize[son[i]]2
= ( s i z e [ p ] ) 2 − ∑ i = 1 k s i z e [ s o n [ i ] ] 2 =(size[p])^2-\sum_{i=1}^ksize[son[i]]^2 =(size[p])2i=1ksize[son[i]]2

#include<bits/stdc++.h>
#define mod 1000000007
#define rep(i,a,b) for(register int (i)=(a);(i)<=(b);(i)++) 
#define don(i,a,b) for(register int (i)=(a);(i)>=(b);(i)--)
using namespace std;
const int maxn=1e5+10;
const int maxm=1e3+10;
int r,n,m,cnt=0;
int a[maxn],ans[maxn],vis[maxn],sum[maxn],head[maxn],siz[maxn];

template <class t> inline void read(t &x)
{
    x=0;int f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=10*x+ch-'0';ch=getchar();}
    x*=f;
}

template <class t> inline void write(t x) {
	if(x<0) {putchar('-');x=~x+1;}
	if(x>9) write(x/10);
	putchar(x%10+48);
}

struct node{
	int v,nex;
}e[maxn];

void add(int u,int v)
{
	e[++cnt].v=v;
	e[cnt].nex=head[u];
	head[u]=cnt;
}

void readdata()
{
	read(n);read(r),read(m);
	rep(i,1,n-1) {
		int x,y;
		read(x),read(y);
		add(x,y);
		add(y,x);
	}
}

inline void dfs(int x,int f) {
    siz[x]=1;
    for(int i=head[x];i;i=e[i].nex){
        int xx=e[i].v;
        if(xx==f) continue;
        dfs(xx,x);
        siz[x]+=siz[xx];
        sum[x]+=siz[xx]*siz[xx];
    }
    ans[x]=(siz[x]*siz[x]-sum[x])%mod;
}

void work()
{
	memset(siz,false,sizeof(siz));
	dfs(r,0);
	rep(i,1,m) {
		int p;
		read(p);
		ans[p]%=mod;
		write(ans[p]);
		putchar('\n');
	}
}

int main()
{
	//freopen("input.txt","r",stdin);
	readdata();
	work();
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
该资源内项目源码是个人的课程设计、毕业设计,代码都测试ok,都是运行成功后才上传资源,答辩评审平均分达到96分,放心下载使用! ## 项目备注 1、该资源内项目代码都经过测试运行成功,功能ok的情况下才上传的,请放心下载使用! 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习,也适合小白学习进阶,当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行,也可在此代码基础上进行修改,以实现其他功能,也可用于毕设、课设、作业等。 下载后请首先打开README.md文件(如有),仅供学习参考, 切勿用于商业用途。 该资源内项目源码是个人的课程设计,代码都测试ok,都是运行成功后才上传资源,答辩评审平均分达到96分,放心下载使用! ## 项目备注 1、该资源内项目代码都经过测试运行成功,功能ok的情况下才上传的,请放心下载使用! 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习,也适合小白学习进阶,当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行,也可在此代码基础上进行修改,以实现其他功能,也可用于毕设、课设、作业等。 下载后请首先打开README.md文件(如有),仅供学习参考, 切勿用于商业用途。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值