bzoj 4805: 欧拉函数求和杜教筛

最新推荐文章于 2023-09-21 23:55:55 发布

_beginend

最新推荐文章于 2023-09-21 23:55:55 发布

阅读量1k

点赞数

分类专栏：杜教筛

本文链接：https://blog.csdn.net/qq_33229466/article/details/70140643

版权

杜教筛专栏收录该内容

11 篇文章 0 订阅

订阅专栏

题意

给出一个数字N，求sigma(phi(i)),1<=i<=N。
正整数N。N<=2*10^9

分析

同bzoj 3994

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=5000005;
const int MOD=1000000007;

int tot,prime[N/10];
bool not_prime[N];
LL phi[N];

void get_prime(int n)
{
    phi[1]=1;
    for (int i=2;i<=n;i++)
    {
        if (!not_prime[i]) prime[++tot]=i,phi[i]=i-1;
        for (int j=1;j<=tot&&i*prime[j]<=n;j++)
        {
            not_prime[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
    for (int i=1;i<=n;i++) phi[i]+=phi[i-1];
}

LL solve(int n)
{
    if (n<=5000000) return phi[n];
    LL ans=(LL)n*(n+1)/2;
    for (int i=2,last;i<=n;i=last+1)
    {
        last=n/(n/i);
        ans-=(LL)(last-i+1)*solve(n/i);
    }
    return ans;
}

int main()
{
    get_prime(5000000);
    int n;
    scanf("%d",&n);
    printf("%lld",solve(n));
    return 0;
}