ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Codeforces 1401 D

2021-03-03 16:57:06  阅读:238  来源: 互联网

标签:int ll Codeforces leq 1401 lld top 边权


传送门

题目大意

给定一棵 n n n个节点, n − 1 n-1 n−1 条边的树。你可以在每一条树上的边标上边权,使得:
每个边权都为 正整数;
这 n − 1 n-1 n−1个边权的乘积等于 k k k;
边权为 1 1 1 的边的数量最少。
定义 f ( u , v ) f(u,v) f(u,v) 表示节点 u u u到节点 v v v的简单路径经过的边权总和。你的任务是让 ∑ i = 1 n − 1 ∑ j = i + 1 n f ( i , j ) \sum_{i=1}^{n-1}\sum_{j=i+1}^{n}{f(i,j)} ∑i=1n−1​∑j=i+1n​f(i,j)最大。
最终答案可能很大,对 1 e 9 + 7 1e9+7 1e9+7取模即可。
k k k有可能很大,输入数据中包含了 m m m个质数 p i p_i pi​ ,那么 k k k为这些质数的乘积。
输入格式
第一行,一个整数 t ( 1 ≤ t ≤ 100 ) t(1\leq t\leq 100) t(1≤t≤100),表示多组测试数据个数。对于每一个测试数据:
第一行,一个整数 n ( 2 ≤ n ≤ 1 0 6 ) n(2\leq n\leq 10^6) n(2≤n≤106),表示树上节点数;
第 2 2 2至 n n n行,每行两个整数 u i u_i ui​和 v i ( 1 ≤ u i , v i ≤ n , u i ≠ v i ) ( 1 ≤ u i , v i ≤ n , u i ≠ v i ) v_i(1\leq u_i,v_i \leq n,u_i\neq v_i)(1≤u_i,v_i ≤n,u_i\neq v_i) vi​(1≤ui​,vi​≤n,ui​​=vi​)(1≤ui​,vi​≤n,ui​​=vi​),描述了一条无向边;
第 n + 1 n+1 n+1行,一个整数 m ( 1 ≤ m ≤ 6 × 1 0 4 ) m(1\leq m\leq 6\times 10^4) m(1≤m≤6×104),表示 k k k分解成质因子的个数;
第 n + 2 n+2 n+2行, m m m个质数 p i ( 2 ≤ p i < 6 × 1 0 4 , 有 k = ∏ i = 1 m p i p_i (2\leq p_i< 6\times 10^4,有 k=\prod\limits_{i=1}^m p_i pi​(2≤pi​<6×104,有k=i=1∏m​pi​
数据保证所有的 n n n总和不超过 1 0 5 10^5 105,所有的 m m m总和不超过 6 × 1 0 4 6\times 10^4 6×104。
数据给出的边保证能够形成一棵树。

思路

考虑每条边对答案有多少次贡献。一条路径被经过,仅当起点是它左边的一个点,终点是它右边的一个点。贡献=左边点个数 × \times ×右边点个数 × \times ×边权。然后把最大的边权分给出现次数最多的就可以了。
当 m < n − 1 m<n-1 m<n−1时,要填上 1 1 1直到 m = n − 1 m=n-1 m=n−1。
当 m > n − 1 m>n-1 m>n−1时,要把多余的最大的合并成一个新的数,值为数的积,直到 m = n − 1 m=n-1 m=n−1。

代码

ll t,n,m,ans;
ll p[maxn],s[maxn];
vector<ll> e[maxn]; 
priority_queue<ll> q;

void init(){
	memset(s,0,sizeof s);
	memset(p,0,sizeof p);
	ans=0;
	for(int i=1;i<=n;i++){
		e[i].clear();
	}
	while(!q.empty()) q.pop();
}

void dfs(ll x,ll f){
	s[x]=1;
	for(int i=0;i<e[x].size();i++){
		ll y=e[x][i];
		if(y==f) continue;
		dfs(y,x);
		s[x]+=s[y];
	}
}

void dfs2(ll x,ll f){
	for(int i=0;i<e[x].size();i++){
		ll y=e[x][i];
		if(y==f) continue;
		q.push(s[y]*(n-s[y]));
		dfs2(y,x);
	}
}

int cmp(ll a,ll b){
	return a>b;
}

int main(){
	int t;
	cin>>t;
	while(t--){
		init();
		scanf("%lld",&n);ll x,y;
		for(int i=1;i<=n-1;i++){
			scanf("%lld%lld",&x,&y);
			e[x].push_back(y);
			e[y].push_back(x);
		}
		scanf("%lld",&m);
		for(int i=1;i<=m;i++){
			scanf("%lld",&p[i]);//质数 
		}
		dfs(1,0);
		dfs2(1,0);
		sort(p+1,p+1+m,cmp);
		if(m<=(n-1)){
			ll top=1;
			while(!q.empty()){
				if(top<=m) ans+=q.top()*p[top];
				else ans+=q.top();
				ans%=mod;
				q.pop();
				top++;
			}
		}
		else{
			ll top=m-n+2;
			ll sum=1;
			for(int i=1;i<=(m-(n-1))+1;i++)
			{
				(sum*=p[i])%=mod;
			}
			p[top]=sum;
			
			while(q.size())
			{
				(ans+=q.top()*p[top])%=mod;
				q.pop();
				top++;
			}
		}
		printf("%lld\n",ans);
	}
}

标签:int,ll,Codeforces,leq,1401,lld,top,边权
来源: https://blog.csdn.net/weixin_45633221/article/details/112630337

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有