SPOJ 7001-Visible Lattice Points （莫比乌斯反演）

题目链接：https://cn.vjudge.net/contest/180599#problem/A

题意：给出n，求出有多少小于n的三元组（自然数），使得gcd(x, y, z)为1，

 即从三维空间(0, 0, 0) 可以看到的点数。

题解：莫比乌斯反演的应用，首先回顾莫比乌斯的两种形式

本题采用了第二种形式，由于n可以取到0，故要考虑0的情况。。。。

#include<map>     
#include<stack>            
#include<queue>            
#include<vector>    
#include<string>  
#include<math.h>            
#include<stdio.h>            
#include<iostream>            
#include<string.h>            
#include<stdlib.h>    
#include<algorithm>   
#include<functional>    
using namespace std;
typedef long long  ll;
#define inf  1000000000       
#define mod 1000000007             
#define maxn  1000005  
#define lowbit(x) (x&-x)            
#define eps 1e-9  
int a[maxn] = { 1,1 }, b[maxn/10], mu[maxn], cnt;
void init()
{
	ll i, j;mu[1] = 1;
	for (i = 2;i < maxn;i++)
	{
		if (a[i] == 0)
			b[++cnt] = i, mu[i] = -1;
		for (j = 1;j <= cnt && b[j] * i < maxn;j++)
		{
			a[b[j] * i] = 1;
			if (i % b[j] == 0)
			{
				mu[b[j] * i] = 0;
				break;
			}
			else
				mu[b[j] * i] = -mu[i];
		}
	}
}
int main(void)
{
	init();
	ll i, j, n, T;
	scanf("%lld", &T);
	while (T--)
	{
		ll sum = 3;//有两个坐标为0的情况
		scanf("%lld", &n);
		for (i = 1;i <= n;i++)//其中一个左边为0的情况
			sum += (ll)mu[i] * (n / i)*(n / i) * 3;
		for (i = 1;i <= n;i++)
			sum += (ll)mu[i] * (n / i)*(n / i)*(n / i);
		printf("%lld\n", sum);
	}
	return 0;
}