BZOJ 2190 仪仗队 （莫比乌斯反演）

思路：

从 $(0,0)$ 开始编号，构造一个n+1阶的矩阵，通过观察我们发现，从$(0,0)$ 点看去，点$(4,6)$ 一定被点$(2,4)$ 挡住了。我们发现，只有当$gcd(i，j)==1$ 的点，才能被看到，所以，问题转化成了求$(1，n)$ 区间内，$gcd(i,j)==1$ 的无序点对的个数。然后我们需要加上$(0,1)$和$(1,0)$两个点。 此题解决。

#include<stdio.h>
#include <iostream>
#include<string.h>
#include<math.h>
#include<algorithm>
#define eps 1e-8
typedef long long int lli;
using namespace std;

const int maxn = 4e4+10;
bool isprime[maxn];
int prime[maxn];
char miu[maxn];
void moblus(){
    int cnt = 0;miu[1] = 1;
    for(lli i = 2;i < maxn;i++){
        if(!isprime[i]){
            prime[cnt++] = i,miu[i] = -1;//phi[i] = i-1;
        }
        for(lli j = 0;j < cnt && i*prime[j] < maxn;j++){
            lli x = prime[j];
            isprime[i*x] = 1;
            if(i%x){
                miu[i*x] = -miu[i];
                //phi[i*x] = phi[i] * phi[x];
            }
            else{
                miu[i*x] = 0;
                //phi[i*x] = phi[i] * x;
                break;
            }
        }
    }
}
int ff[maxn],sum[maxn];
int main(){
    moblus();
    for(int i = 1;i <= maxn;i++){
        sum[i] = sum[i-1] + miu[i];
    }
    lli a,b;
    scanf("%lld",&a);
    a--;
    lli ans = 0,l;
    for(int j=1;j<=a;j=l+1){//分块加速
        l= a/(a/j);
        ans += (lli)(sum[l]-sum[j-1])*(a/j)*(a/j);
    }
    printf("%lld\n",ans+2);
    return 0;
}