[51Nod1239]欧拉函数之和

最新推荐文章于 2020-01-11 13:02:01 发布

清疚

最新推荐文章于 2020-01-11 13:02:01 发布

阅读量166

点赞数

分类专栏：数论文章标签：杜教筛 51nod

本文链接：https://blog.csdn.net/qq_38064038/article/details/79317814

版权

数论专栏收录该内容

13 篇文章 0 订阅

订阅专栏

题意

给定 $n$ ，求 $\phi(n)=\sum_{i=1}^{n}{\varphi(i)}$

解法

这道题就要用到一种神奇的黑科技：杜教筛了！
首先来推公式：
$\sum_{d|n}{\varphi(d)}=n · · · ①$
由①可知： $\varphi(n)=n-\sum_{d|n,d<n}\varphi(d) · · · ②$
将②带入到所求的式子中去，有：
$ϕ (n) = \sum i = 1 n [i - \sum d | i, d < i φ (d)] = n * ( n + 1 ) 2 - \sum i = 2 n \sum d | i, d < i φ (d)$ $\phi(n)=\sum_{i=1}^n[i-\sum_{d|i,d<i}\varphi(d)]=\dfrac{n*(n+1)}{2}-\sum_{i=2}^n\sum_{d|i,d<i}\varphi(d)$
$= n * ( n + 1 ) 2 - \sum i = 2 n \sum d = 1 ⌊ n i ⌋ φ (d) = n * ( n + 1 ) 2 - \sum i = 2 n ϕ (⌊ n i ⌋)$ $=\dfrac{n*(n+1)}{2}-\sum_{i=2}^n\sum_{d=1}^{\lfloor \dfrac{n}{i} \rfloor}\varphi(d)=\dfrac{n*(n+1)}{2}-\sum_{i=2}^n\phi(\lfloor \dfrac{n}{i} \rfloor )$
然后就只需要不停递归下去，就可以达到比较优的复杂度： $T(n)=O(\sqrt{n})+\sum_{i=1}^{\sqrt{n}}{T(i)+T(\frac{n}{i})}$ ，利用主定理便可得知复杂度为 $O(n^\frac{3}{4})$
如果能够预处理出一部分的 $\phi(n)$ ，那么复杂度可以进一步优化，当预处理出前 $n^{\frac{2}{3}}$ 项，复杂度最优
PS： $有人问：\sum_{i=2}^n\sum_{d|i,d<i}\varphi(d)=\sum_{i=2}^n\sum_{d=1}^{\lfloor \dfrac{n}{i} \rfloor}\varphi(d)$ 是怎么得到的，其实很简单，左边是一个一个枚举 $\varphi(d)$ ，而右边则是算出 $d$ 可以成为多少个数的因子，那么 $d$ 就可以产生那么多的贡献

复杂度

O( $n^{\frac{2}{3}}$ )

代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<map>
#define Rint register int
#define Lint long long int
using namespace std;
const int M=1e9+7;
const int N=10000010;
const int m2=500000004;
bool vis[N];
int pri[N/10],cnt;
Lint phi[N];
map<Lint,int> f;
void Prepare()
{
    phi[1]=1;
    for(int i=2;i<N;i++)
    {
        if( !vis[i] )   pri[++cnt]=i,phi[i]=i-1;
        for(int j=1;j<=cnt;j++)
        {
            int x=pri[j]*i;
            if( x>=N )   break ;
            vis[x]=1;
            if( i%pri[j] )   phi[x]=phi[i]*phi[pri[j]];
            else
            {
                phi[x]=phi[i]*pri[j];
                break ;
            }
        }
    }
    for(int i=1;i<N;i++)   phi[i]=(phi[i]+phi[i-1])%M;
}
int cal(Lint n)
{
    if( n<N )   return phi[n];
    if( f.count(n) )   return f[n];
    Lint x=2,ret=n%M*(n%M+1)%M*m2%M;
    while( x<=n )
    {
        Lint y=n/(n/x);
        ret=(ret-1ll*(y-x+1)%M*cal(n/x)%M+M)%M,x=y+1;
    }
    return f[n]=ret;
}
int main()
{
    Lint n;
    Prepare();
    scanf("%lld",&n);
    printf("%d\n",cal(n));
    return 0;
}