BZOJ 2005 NOI2010 能量采集 莫比乌斯反演

BZOJ又被卡了

---

这道题水过，作为第一道莫比乌斯反演的习题，我还是理解了好长时间分块的想法的，可能是我太菜了吧…

---

首先一个推导我们需要求的是

$$\sum_{i=0}^n\sum_{j=0}^m({gcd(i,j)*2-1})$$
进行推导
$$=2*\sum_{i=0}^n\sum_{j=0}^m(gcd(i,j))-n*m$$
$$根据n=\sum_{d|n}φ(d),将上式中的n置换成gcd(i,j),则$$
$$=2*\sum_{i=0}^n\sum_{j=0}^m\sum_{d|gcd(i,j)}φ(d)-n*m$$
$$=2*\sum_{i=0}^n\sum_{j=0}^m\sum_{d=1}^n[d|i][d|j]φ(d)-n*m$$
$$=2*\sum_{d=1}^n\sum_{i=0}^n[d|i]\sum_{j=0}^m[d|j]φ(d)-n*m$$
$$=2*\sum_{d=1}^n\lfloor n/d \rfloor \lfloor m/d\rfloorφ(d)-n*m$$
一组询问可以O（n）水过去，多组询问就分块搞
注意其中有一步，是直接把d转换到1到n了，其实这里有个取较小值操作，我们可以把n当作小的那个，然后来搞

---

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int MAXN=100000+50;
long long n,m,is_not_prime[MAXN],phi[MAXN],prime[MAXN];
inline void getphi(){
    phi[1]=1;
    for(register int i=2;i<=n;i++){
        if(!is_not_prime[i]){
            prime[++prime[0]]=i;
            phi[i]=i-1;
        }
        for(register int j=1;j<=prime[0]&&i*prime[j]<=n;j++){
            is_not_prime[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            } 
        }
        phi[i]+=phi[i-1];
    }
}
int main(){
    scanf("%lld%lld",&n,&m);
    if(n<m) swap(n,m);
    getphi();
    long long end,final=0;
    for(register long long start=1;start<=m;start=end+1){
        end=min(n/(n/start),m/(m/start));
        final+=(phi[end]-phi[start-1])*(n/start)*(m/start);
    }
    printf("%lld\n",final*2-m*n);
    return 0;
}