Visible Trees HDU - 2841 （容斥 莫比乌斯函数）

https://cn.vjudge.net/problem/HDU-2841

题意 给你一个n*m的格子 问在格子外（0,0）点能看到多少个点。

思路 如果改点能被表示成 k*（x,y），则这个点会被（x,y）遮挡，所有gcd(x,y) == 1.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 101000;
bool is[maxn];
ll n,ans,sum[maxn],tot = 0;
ll pri[maxn],miu[maxn];
void init() //首先把莫比乌斯函数筛出来
{
    miu[1]=1;
    for(int i=2; i<=maxn; i++)
    {
        if(!is[i])
        {
            pri[++tot]=i;
            miu[i]=-1;
        }
        for(int j=1; j<=tot; j++)
        {
            int k=pri[j]*i;
            if(k>maxn)
                break;
            is[k]=1;
            if(i%pri[j]==0)
            {
                miu[k]=0;
                break;
            }
            else
                miu[k]=-miu[i];
        }
    }
}
int main()
{
    init();
    ll t,m;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        ans = 0;
        ll c = min(n,m);
        for(int i = 1; i<=c; i++)
        {
            ans += (ll)miu[i]*(n/i)*(m/i);
        }
        printf("%lld\n",ans);
    }
}