【BZOJ】4804 欧拉心算莫比乌斯函数+欧拉函数+数论分块

最新推荐文章于 2024-09-01 15:19:02 发布

Chester_King

最新推荐文章于 2024-09-01 15:19:02 发布

阅读量328

点赞数

分类专栏： BZOJ 难题莫比乌斯函数数论分块欧拉函数

本文链接：https://blog.csdn.net/lyfsb/article/details/78806468

版权

BZOJ 同时被 3 个专栏收录

128 篇文章 0 订阅

订阅专栏

难题

80 篇文章 1 订阅

订阅专栏

欧拉函数

4 篇文章 0 订阅

订阅专栏

题目传送门

来来来，推式子啦：

\sum i = 1 n \sum j = 1 n ϕ (g c d (i, j)) = \sum i = 1 n \sum j = 1 n \sum d = 1 n [g c d (i, j) = d] \times ϕ (d) = \sum d = 1 n (ϕ (d) \times \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_{i=1}^n\sum_{j=1}^n\sum_{d=1}^n[gcd(i,j)=d]\times\phi(d)=\sum_{d=1}^n(\phi(d)\times\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{n}{d}\rfloor}[gcd(i,j)=1])$
然后我们来证明一个结论：

∑ni=1∑nj=1[gcd(i,j)=1]=2×∑ni=1ϕ(i)−1 $\sum_{i=1}^n\sum_{j=1}^n[gcd(i,j)=1]=2\times\sum_{i=1}^n\phi(i)-1$

证明：首先我们考虑 $\sum_{i=1}^n\sum_{j=1}^i[gcd(i,j)=1]=\sum_{i=1}^n\phi(i)$ ，这里显然 $i \ge j$ ，显然在 $i \le j$ 的时候也一样，但是 $(1,1)$ 这一对数被算了两次，那么在最后减一即可。

于是可以继续化简式子：

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

$=\sum_{d=1}^n(\phi(d)\times(2\times\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\phi(i)-1))=2\times\sum_{d=1}^n(\phi(d)\times\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\phi(i))-\sum_{i=1}^n\phi(i)$
我们令

sum(n)=∑ni=1ϕ(i) $sum(n)=\sum_{i=1}^n\phi(i)$ ，那么式子就变成了：

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

$2\times\sum_{d=1}^n(\phi(d)\times sum(\lfloor\frac{n}{d}\rfloor))-sum(n)$
然后就可以数论分块在

O(Tn√) $O(T\sqrt n)$ 的时间内解决这题啦。

附上AC代码：

#include <cstdio>
using namespace std;

typedef long long ll;
const int N=1e7+10;
int t,p[N],phi[N],num,n;
ll sum[N];
bool b[N];

inline void prime(int n){
    phi[1]=1;
    for (int i=2; i<=n; ++i){
        if (!b[i]) p[++num]=i,phi[i]=i-1;
        for (int j=1; j<=num&&p[j]*i<=n; ++j){
            b[p[j]*i]=1,phi[p[j]*i]=phi[i]*(p[j]-1);
            if (i%p[j]==0) {phi[p[j]*i]+=phi[i]; break;}
        }
    }
    for (int i=1; i<=n; ++i) sum[i]=sum[i-1]+phi[i];
    return;
}

inline ll calc(int n){
    ll ret=0;
    for (int l=1,r; l<=n; l=r+1) r=n/(n/l),ret+=(sum[r]-sum[l-1])*sum[n/l];
    return 2ll*ret-sum[n];
}

int main(void){
    for (prime(1e7),scanf("%d",&t); t; --t) scanf("%d",&n),printf("%lld\n",calc(n));
    return 0;
}