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

题目大意：

从坐标（0,0,0）处观察到所有在（n,n,n)范围内的点的个数，如果一条直线上出现多个点，除了第一个，后面的都视为被遮挡了

 

这题目稍微推导一下可得知 gcd(x,y,z) = 1的点是可观察到的，若三者的gcd>1，则这个点之前必然出现了一个(x/gcd(x,y,z) , y/gcd(x,y,z) , z/gcd(x,y,z))的点

那么因为 0 是无法计算gcd的 ， 所以先不考虑0

1. x>0 , y>0 , z>0  ans = ∑x<=n∑y<=n∑z<=n[gcd(x,y,z)=1]  -> ∑x<=n∑y<=n∑z<=n∑d|gcd(x,y,z) mu[d]

2. 只存在一个0 ans = 3*(∑x<=n∑y<=n[gcd(x,y)=1]) -> 3*(∑x<=n∑y<=n∑d|gcd(x,y) mu[d])

3.3个数两个为0 ans = 3

那么只要将上面的3个值相加就是我们要求的答案了

 1 #include <cstdio>

 2 #include <cstring>

 3 using namespace std;

 4 const int MAXN = 1000000;

 5 #define ll long long

 6 int mu[MAXN+5] , prime[MAXN/10];

 7 bool check[MAXN+5];

 8 

 9 void getMu(int n)

10 {

11     memset(check , 0 ,  sizeof(check));

12     int tot=0;

13     check[1] = true , mu[1]=1 , mu[2]=-1;

14     for(int i=2 ; i<=n ; i++){

15         if(!check[i]) prime[tot++]=i , mu[i]=-1;

16         check[i] = true;

17         for(int j=0 ; j<tot ; j++){

18             if(i*prime[j]>n) break;

19             check[i*prime[j]]=true;

20             if(i%prime[j] == 0){

21                 mu[i*prime[j]] = 0;

22                 break;

23             }else{

24                 mu[i*prime[j]] = -mu[i];

25             }

26         }

27     }

28 }

29 

30 int main()

31 {

32    // freopen("a.in" , "r" , stdin);

33     getMu(MAXN);

34     int T , n;

35     scanf("%d" , &T);

36     while(T--)

37     {

38         scanf("%d" , &n);

39         ll ans = 0;

40         for(int i=1 ; i<=n  ; i++){

41             ll tmp1 = (ll)(n/i)*(n/i)*(n/i)*mu[i]; //坐标中一个零都不存在的情况

42             ll tmp2 = (ll)3*(n/i)*(n/i)*mu[i];  //坐标中有一个零的情况

43             ll tmp3 = (ll)3*(n/i)*mu[i]; //坐标中有两个零的情况

44             ans += tmp1+tmp2+tmp3;

45         }

46         printf("%lld\n" , ans);

47     }

48     return 0;

49 }