需要將問題分解一下。

我們需要求的是這個立方體從(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; }

總結

以上是生活随笔為你收集整理的SPOJ-VLATTICE Visible Lattice Points-莫比乌斯反演的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。