题目链接: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;
}