【杜教筛】51Nod1239[欧拉函数之和]题解

最新推荐文章于 2020-08-30 12:04:03 发布

ZigZagK

最新推荐文章于 2020-08-30 12:04:03 发布

阅读量430

点赞数

分类专栏： 51Nod题解杜教筛

本文链接：https://blog.csdn.net/zzkksunboy/article/details/78902549

版权

51Nod题解同时被 2 个专栏收录

17 篇文章 0 订阅

订阅专栏

杜教筛

3 篇文章 0 订阅

订阅专栏

题目概述

求 $\sum_{i=1}^{n}\varphi(i)$ 。

解题报告

因为 $n=\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(i,n)=j]=\sum_{d|n}\sum_{i=1}^{n\over d} [gcd(i,{n\over d})=1]=\sum_{d|i}\varphi({n\over d})$ ，所以 $Id=\varphi*1$ 。

那么这道题只要像莫比乌斯函数之和一样用 $1$ 卷 $\varphi$ 就可以了，得到：

S (n) = \sum i = 1 n I d (i) - \sum i = 2 n S (⌊ n i ⌋) = n ( n + 1 ) 2 - \sum i = 2 n S (⌊ n i ⌋)

$S(n)=\sum_{i=1}^{n}Id(i)-\sum_{i=2}^{n}S(\lfloor{n\over i}\rfloor)={n(n+1)\over 2}-\sum_{i=2}^{n}S(\lfloor{n\over i}\rfloor)$

示例程序

#include<cstdio>
#include<map>
using namespace std;
typedef long long LL;
const int maxn=4700000,MOD=1e9+7,INV2=MOD+1>>1;

int p[maxn+5],phi[maxn+5];bool pri[maxn+5];
LL n;map<LL,int> f;

#define Mod(x) ((x)%MOD)
inline void AMOD(int &x,int tem) {if ((x+=tem)>=MOD) x-=MOD;}
inline void Make()
{
    pri[1]=true;phi[1]=1;
    for (int i=2;i<=maxn;i++)
    {
        if (!pri[i]) p[++p[0]]=i,phi[i]=i-1;
        for (int j=1,t;j<=p[0]&&(t=i*p[j])<=maxn;j++)
            if (i%p[j]) pri[t]=true,phi[t]=phi[i]*phi[p[j]]; else
            {pri[t]=true;phi[t]=phi[i]*p[j];break;}
    }
    for (int i=2;i<=maxn;i++) AMOD(phi[i],phi[i-1]);
}
int Sum(LL n)
{
    if (n<=maxn) return phi[n];if (f.count(n)) return f[n];
    int ans=Mod(Mod(Mod(n+1)*Mod(n))*INV2);
    for (LL l=2,r;l<=n;l=r+1)
        r=n/(n/l),AMOD(ans,MOD-Mod(Mod(r-l+1)*Sum(n/l)));
    return f[n]=ans;
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    return Make(),scanf("%lld",&n),printf("%d\n",Sum(n)),0;
}