【树形DP】COCI2018/2019 Deblo

分析:

非常水的树形DP+二进制分组。

每条路径在其最高处计算。
D p [ i ] [ j ] [ 0 / 1 ] Dp[i][j][0/1] Dp[i][j][0/1]表示在i的子树中,到达i的路径在第j位为0/1的方案数。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#define SF scanf
#define PF printf
#define MAXN 100010
#define y1 warspitewarspitewarspiteprinceofwalesprinceofwalesprinceofwales 
using namespace std;
typedef long long ll;
int n;
int val[MAXN];
vector<int> a[MAXN];
ll ans; 
int dp[MAXN][25][2];
void dfs(int x,int fa=0){
	for(int i=0;i<24;i++)
		dp[x][i][(val[x]>>i)&1]++;
	for(int i=0;i<int(a[x].size());i++){
		int u=a[x][i];
		if(u==fa)
			continue;
		dfs(u,x);
		for(int j=0;j<24;j++){
			int &x0=dp[x][j][0];
			int &x1=dp[x][j][1];
			int &y0=dp[u][j][0];
			int &y1=dp[u][j][1];
			ans=ans+(1ll*x0*y1+1ll*x1*y0)*(1<<j);
			if((val[x]>>j)&1)
				x0+=y1,x1+=y0;
			else
				x0+=y0,x1+=y1;
		}
	}
//	for(int i=0;i<3;i++)
//		PF("[%d %d :%d,%d]\n",x,i,dp[x][i][0],dp[x][i][1]);
}
int main(){
	freopen("hyperspace.in","r",stdin);
	freopen("hyperspace.out","w",stdout);
	SF("%d",&n);
	for(int i=1;i<=n;i++){
		SF("%d",&val[i]);
		ans+=val[i];
	}
	int u,v;
	for(int i=1;i<n;i++){
		SF("%d%d",&u,&v);
		a[u].push_back(v);
		a[v].push_back(u);
	}
	dfs(1);
	PF("%lld",ans);
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值