bzoj4804: 欧拉心算

最新推荐文章于 2019-08-31 19:34:00 发布

*ACoder*

最新推荐文章于 2019-08-31 19:34:00 发布

阅读量1k

点赞数

分类专栏： # 欧拉函数

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/69945454

版权

欧拉函数专栏收录该内容

14 篇文章 0 订阅

订阅专栏

这篇博客介绍了BZOJ4804题目，作者提供了简洁的题解，指出并不需要复杂的莫比乌斯反演，直接给出了解题代码。

摘要由CSDN通过智能技术生成

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=4804

题解

　　网上的题解都好麻烦啊还用莫比乌斯反演，哪用得着啊。
　　

\sum i = 1 n \sum j = 1 n ϕ (g c d (i, j)) = \sum d = 1 n ϕ (i) \sum i = 1 ⌊ N d ⌋ \sum j = 1 ⌊ N d ⌋ [g c d (i, j) = 1]

$\sum_{i=1}^n\sum_{j=1}^n\phi(gcd(i,j)) 　　\\=\sum_{d=1}^n\phi(i)\sum_{i=1}^{\lfloor\frac{N}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{N}{d}\rfloor}[gcd(i,j)=1]$
　　令

s(x)=∑xi=1ϕ(i) $s(x)=\sum_{i=1}^x\phi(i)$ ，原式
　　

= \sum d = 1 n ϕ (i) (s (⌊ n i ⌋) \times 2 - 1)

$=\sum_{d=1}^n\phi(i)(s(\lfloor\frac{n}{i}\rfloor)\times 2-1)$
　　时间复杂度

O(N+TN−−√) $O(N+T\sqrt{N})$

代码

//欧拉函数
#include <cstdio>
#include <algorithm>
#define maxn 10000001
#define ll long long
using namespace std;
bool mark[maxn];
int prime[maxn/10];
ll phi[maxn], s[maxn];
void init()
{
    ll i, j;
    phi[1]=1;
    for(i=2;i<maxn;i++)
    {
        if(!mark[i])prime[++prime[0]]=i, phi[i]=i-1;
        for(j=1;j<=prime[0] and i*prime[j]<maxn;j++)
        {
            mark[i*prime[j]]=1;
            if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}
            phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
    for(i=1;i<maxn;i++)s[i]=s[i-1]+phi[i];
}
void work(ll N)
{
    ll i, last, ans=0;
    for(i=1;i<=N;i=last+1)
    {
        last=N/(N/i);
        ans+=(s[last]-s[i-1])*((s[N/i]<<1)-1);
    }
    printf("%lld\n",ans);
}
int main()
{
    ll T, N;
    init();
    for(scanf("%lld",&T);T;T--)
    {
        scanf("%lld",&N);
        work(N);
    }
    return 0;
}