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

题目传送门

来来来，推式子啦：

∑i=1n∑j=1nϕ(gcd(i,j))=∑i=1n∑j=1n∑d=1n[gcd(i,j)=d]×ϕ(d)=∑d=1n(ϕ(d)×∑i=1⌊nd⌋∑j=1⌊nd⌋[gcd(i,j)=1])

然后我们来证明一个结论： ∑ni=1∑nj=1[gcd(i,j)=1]=2×∑ni=1ϕ(i)−1

证明：首先我们考虑 ∑ni=1∑ij=1[gcd(i,j)=1]=∑ni=1ϕ(i) ，这里显然 i≥j ，显然在 i≤j 的时候也一样，但是 (1,1) 这一对数被算了两次，那么在最后减一即可。

于是可以继续化简式子：

=∑d=1n(ϕ(d)×(2×∑i=1⌊nd⌋ϕ(i)−1))=2×∑d=1n(ϕ(d)×∑i=1⌊nd⌋ϕ(i))−∑i=1nϕ(i)

我们令 sum(n)=∑ni=1ϕ(i) ，那么式子就变成了：

2×∑d=1n(ϕ(d)×sum(⌊nd⌋))−sum(n)

然后就可以数论分块在 O(Tn√) 的时间内解决这题啦。

附上AC代码：

#include 
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;
}