zoj 3435 Ideal Puzzle Bobble (莫比乌斯反演)

题目链接：哆啦A梦传送门
题意：这题跟 bzoj2190 类似。
在三维坐标(a,b,c)下，原点点(1,1,1)能看到多少个点？
这题也同理，我们把点(1,1,1)看成是点(0,0,0)就可以了。

那么我们就可这几种情况：
原点在平面 (x,y) 能看多少个点？$1\le x\le (a-1),1\le y\le (b-1)$
原点在平面 (x,z) 能看多少个点？$1\le x\le (a-1),1\le z\le (c-1)$
原点在平面 (y,z) 能看多少个点？$1\le y\le (b-1),1\le z\le (c-1)$

原点在平面 (x,y,z) 能看多少个点？$1\le x\le (a-1)1\le y\le (b-1),1\le z\le (c-1)$

原点分别在$x,y,z$ 轴上能看多少个点？显然总共能看3个点。

那么我们这里也就跟bzoj2190做法类似。

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N=1e6+10;

int prime[N],tot,mu[N];
LL sum[N];

bool vis[N];

void init()
{
    mu[1]=sum[1]=1;

    for(int i=2;i<N;i++)
    {
        if(!vis[i]){
            prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=tot&&i*prime[j]<N;j++){
            int x=i*prime[j];

            vis[x]=1;
            if(i%prime[j]==0){
                mu[x]=0;
                break;
            }
            mu[x]=-mu[i];

        }
        sum[i]=sum[i-1]+mu[i];
    }

}

LL solve(int a,int b)
{
    LL ret=0;

    if(a>b) swap(a,b);
    int r;
    for(int l=1;l<=a;l=r+1)
    {
        r=min(a/(a/l),b/(b/l));


        ret=ret+1LL*(sum[r]-sum[l-1])*(a/l)*(b/l);
    }
    return ret;
}

LL solve1(int a,int b,int c)
{
    LL ret=0;

    if(a>b) swap(a,b);
    if(a>c) swap(a,c);

    int r;
    for(int l=1;l<=a;l=r+1)
    {
        r=min(a/(a/l),min(b/(b/l),c/(c/l)));

        ret=ret+1LL*(sum[r]-sum[l-1])*(a/l)*(b/l)*(c/l);

    }
    return ret;

}

int main()
{

    init();
    int a,b,c;
    while(~scanf("%d%d%d",&a,&b,&c)){


    a--;b--;c--;
    LL ret=solve(a,b)+solve(a,c)+solve(b,c)+solve1(a,b,c);
    printf("%lld\n",ret+3);

    }
    return 0;
}