首页 > 代码库 > SPOJ 7001

SPOJ 7001

这次求的是三维的。如上篇一样,使用那个特殊的莫比乌斯反演来做。

需要注意的是,求得的只是从(1,1,1)到(n,n,n)的互质的个数,还要注意墙壁三面如(0,1,1)~(0,n,n)等。最后+3。因为是(0,0,1);

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>#define N 1000005using namespace std;typedef long long LL;int mobi[N];bool vis[N];int n;void initial(){    int i,j;    for(i=1;i<N;i++) mobi[i]=1,vis[i]=false;     for(i=2;i<N;i++) {        if(vis[i]) continue;        for(j=i;j<N;j+=i){            vis[j]=true;            if((j/i)%i==0){                mobi[j]=0; continue;            }            mobi[j]=-mobi[j];        }    }}int main(){	initial();	int T;	scanf("%d",&T);	while(T--){		scanf("%d",&n);		LL ans=0;		for(int i=1;i<=n;i++){			ans+=((LL)mobi[i]*(LL)(n/i)*(LL)(n/i)*(LL)(n/i)+3*(LL)mobi[i]*(LL)(n/i)*(LL)(n/i));		}		ans+=3;		printf("%lld\n",ans);	}	return 0;}

  

SPOJ 7001