P1891 疯狂 LCM

最新推荐文章于 2024-07-08 21:41:38 发布

lzc__123

最新推荐文章于 2024-07-08 21:41:38 发布

阅读量213

点赞数 1

文章标签： c++ 算法开发语言

本文链接：https://blog.csdn.net/liangsheng_meng/article/details/127602161

版权

思路：
请添加图片描述

代码

#include<bits/stdc++.h>
#define ll long long 
using namespace std;
const int N=1e6+7;
int primes[N],cnt;
ll euler[N];
bool st[N];
ll ans[N];
void get_eulers(int n)
{
    euler[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            euler[i] = i - 1;
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            int t = primes[j] * i;
            st[t] = true;
            if (i % primes[j] == 0)
            {
                euler[t] = euler[i] * primes[j];
                break;
            }
            euler[t] = euler[i] * (primes[j] - 1);
        }
    }
    for(int i=1;i<=N;i++)
    {
    	for(int j=1;j*i<=N;j++)
    	{
    		ans[i*j]+=1ll*euler[j]*j/2;
		}
	}
	for(int i=1;i<=N;i++)
	{
		ans[i]++;
		ans[i]*=i;
	}
}
int main()
{
	int t;scanf("%d",&t);
	get_eulers(N);
	while(t--)
	{
		int n;scanf("%d",&n);
		printf("%lld\n",ans[n]);
	}
}