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

最新推荐文章于 2020-12-05 22:21:18 发布

原创最新推荐文章于 2020-12-05 22:21:18 发布 · 458 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#莫比乌斯反演 #数论 #GCD

**********数论********** 同时被 2 个专栏收录

18 篇文章

订阅专栏

莫比乌斯反演

9 篇文章

订阅专栏

本文介绍了一个数学问题的解决方法，即计算N*N*N正方体中有多少整点无法从原点(0,0,0)直接看到。文章通过讨论不同情况下的整点可见性条件，并利用莫比乌斯函数简化计算过程。

此题意思很简单：在 $N*N*N$ 的正方体中，有多少整点不能被 $(0,0,0)$ 看到。如果一个点 $(x,y,z)$ 不能被 $(0,0,0)$ 看到，那么一定有一个节点 $(x',y',z')$ 在 $(0,0,0)$ 到 $(x,y,z)$ 的连线上，那么 $gcd(x,y,z)！=1$ ，反之，如果
$gcd(x,y,z)=1$ ，那么这个点一定能被看到。

考虑退化情况： $x,y,z$ 中有一个为 $0$ 那么只需要满足 $gcd(x,y)=1$ ，如果有两个为0，那么只有可能是 $(0,0,1),(0,1,0),(1,0,0)$ 。

然后情况就变成了一个莫比乌斯水题。

附上代码：

#include <bits/stdc++.h>
#define LL long long
#define FOR(i,x,y)  for(int i = x;i < y;++ i)
#define IFOR(i,x,y) for(int i = x;i > y;-- i)

using namespace std;

const int maxn = 1000010;

int prime[maxn],mu[maxn];
bool check[maxn];

void Mobius(){
    memset(check,false,sizeof(check));
    prime[0] = 0;
    mu[1] = 1;
    FOR(i,2,maxn){
        if(!check[i]){
            prime[++prime[0]] = i;
            mu[i] = -1;
        }
        for(int j = 1;j <= prime[0];++ j){
            if(i*prime[j] >= maxn)  break;
            check[i*prime[j]] = true;
            if(i % prime[j]){
                mu[i*prime[j]] = -mu[i];
            }
            else{
                mu[i*prime[j]] = 0;
                break;
            }
        }
    }
}

int n;

void work(){
    LL res = 0;
    for(LL d = 1;d <= n;++ d){
        res += mu[d]*((n/d)+3LL)*(n/d)*(n/d);
    }
    printf("%lld\n",res+3LL);
}

int main(){
    //freopen("test.in","r",stdin);
    Mobius();
    int T;  scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        work();
    }
    return 0;
}