#MNS145. gcdmst（AC代码）

最新推荐文章于 2024-10-08 12:54:42 发布

sjy_1234

最新推荐文章于 2024-10-08 12:54:42 发布

阅读量311

点赞数 5

文章标签：算法

本文链接：https://blog.csdn.net/sjy_1234/article/details/141322036

版权

#MNS145. gcdmst

WARNING：本题暂无题解，请自行理解代码

https://marsoj.cn/p/MNS145https://marsoj.cn/p/MNS145

【题目描述】

按照如下规则给定一个完全图：

共 $n - 1$ 个点，点的编号为 $2 \dots n$ ，两点 $i, k$ 之间的连边边权如下：

若 $i, k$ 不互素为： $\gcd(i, k)$ 。

若 $i, k$ 互素为： $\operatorname{lcm}(i, k)$ 。

求这个完全图的最小生成树。

【输入格式】

从gcdmst.in中读入数据。

第 11 行给定一个正整数 $T$ ，代表询问个数。

第 22 行给定 $T$ 个正整数，第 $i$ 个代表第 $i$ 组询问的 $n$ 。

【输出格式】

输出到gcdmst.out中。

$n$ 行，每行一个正整数，第 $i$ 行代表第 $i$ 组询问的答案。

【输入样例 1】

1
4

【输出样例 1】

【样例 1 解释】

$n = 4$ 时，最优的边集为 $\{(2, 3), (2, 4)\}$ ，连边的代价为 $\text{lcm}(2, 3) + \gcd(2, 4) = 6 + 2 = 8$

【数据范围及约定】

对于 30% 的数据： $1 \leq T \leq 300, 2 \leq n \leq 1000$ 。
对于 100% 的数据： $1 \leq T \leq 10^5, 2 \leq n \leq 10^5$ 。

【代码】

#include<bits/stdc++.h>
using namespace std;
const int N=1e5;
long long mp[N+5],sump[N+5];
vector<int>prime;
int main()
{
	//freopen("gcdmst.in","r",stdin);
	//freopen("gcdmst.out","w",stdout);
	for(int i=2;i<=N;i++){
		if(!mp[i])prime.push_back(i);
		for(int p:prime){
			if(p>N/i)break;
			mp[p*i]=p;
			if(i%p==0)break;
		}
	}
	for(int p:prime)sump[p]=2*p;
	for(int i=2;i<=N;i++)mp[i]+=mp[i-1],sump[i]+=sump[i-1];
	int T,n;
	long long ans;
	cin>>T;
	while(T--){
		cin>>n;
		if(n==2)ans=0;
		else if(n==3)ans=6;
		else ans=mp[n]+sump[n]-sump[n/2]/2-2;
		cout<<ans<<endl;
	}
	return 0;
}