SPOJ-VLATTICE Visible Lattice Points-莫比乌斯反演

需要将问题分解一下。

我们需要求的是这个立方体从(0,0,0)能看到的点的个数。可是在三个含有(0,0,0)的面上我们没有办法和其他的一起进行分析（因为含有坐标是0，而我们的数论工具都是从1开始的），所以我们可以将那三个面分开考虑，剩下的就是一个立方体，这个立方体中能看到的点的坐标就是gcd(x,y,z)=1，然后用莫比乌斯反演进行处理。对于那三个面，首先是三个坐标轴上我们只能看到三个点，剩下的每个面中我们只能看到gcd(x,y)=1的点，也用莫比乌斯反演处理。

AC代码

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<ctime>
#include<climits>
#include<queue>
#include<vector>
#include<set>
#include<map>
using namespace std;

typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=1e6+5;
int prime[MAXN],mobius[MAXN],sum[MAXN];
bool check[MAXN]; int tot;

void pre()
{
    tot=0; mobius[1]=1; sum[1]=1;
    for(int i=2;i<MAXN;i++)
    {
        if(!check[i])
        {
            prime[tot++]=i; mobius[i]=-1;
        }
        for(int j=0;j<tot && prime[j]*i<MAXN;j++)
        {
            check[prime[j]*i]=true;
            if(i%prime[j]) mobius[prime[j]*i]=-mobius[i];
            else
            {
                mobius[prime[j]*i]=0; break;
            }
        }
        sum[i]=sum[i-1]+mobius[i];
    }
}

int main()
{
    pre();
    int T,n;
    scanf("%d",&T);
    while(T--)
    {
        ll ans=3;
        scanf("%d",&n);
        int l,r;
        for(l=1;l<=n;l=r+1)
        {
            r=n/(n/l);
            ans+=(ll)(sum[r]-sum[l-1])*(n/l)*(n/l)*(n/l);
            ans+=(ll)(sum[r]-sum[l-1])*(n/l)*(n/l)*3;
        }
        printf("%lld\n",ans);
    }
    return 0;
}