noip2019集训测试赛(五)A.lcm

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接: https://blog.csdn.net/tylon2006/article/details/99476171

Description

t 组询问, 每组询问给定 n , 求 i = 1 n l c m ( n , i ) \sum_{i=1}^nlcm(n,i)


Input

第一行一个整数 t 。
接下来 t 行, 每行一个整数 n 。


Output

输出共 m 行, 每行一个整数表示答案。


Solution

欧拉函数筛(线性筛求 ϕ \phi ):https://blog.csdn.net/Tc_To_Top/article/details/48025849
我也没法讲得再清楚了,请看学长博客:https://blog.csdn.net/ez_yww/article/details/77922114

只要注意s[1]要取1,因为 1 ×   ϕ ( 1 ) = 0 1\times\ \phi(1) = 0 但实际上1有贡献,为1

%%%yww


Code

#include<bits/stdc++.h>
using namespace std;
const int maxn=1000010;
long long phi[1000010];
long long sum[1000010];
long long s[1000010];
int p[1000010],cnt;
bool vis[1000010];
void work(){
	for(int i=2;i<=maxn;i++){
		if(!vis[i]){
			p[cnt++]=i;
			phi[i]=i-1;//线性欧拉筛求phi
		}
		for(int j=0;j<cnt&&p[j]*i<=maxn;j++){
			vis[p[j]*i]=1;
			if(i%p[j]==0){
				phi[p[j]*i]=p[j]*phi[i];
				break;
			}
			phi[p[j]*i]=phi[p[j]]*phi[i];
		}
	}
	s[1]=1;
	for(int i=2;i<=maxn;i++)
	s[i]=i*1ll*phi[i]/2;
	for(int i=1;i<=maxn;i++)
	for(int j=i;j<=maxn;j+=i)
	sum[j]+=s[i];
}
int main(){
	work();
	int t,n;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		printf("%lld\n",sum[n]*n);
	}
}

猜你喜欢

转载自blog.csdn.net/tylon2006/article/details/99476171