【HDU-2841-Visible Trees】 容斥原理

HDU-2841-Visible Trees
这道题题意就是给出第一象限的n*m个点，求出站在原点可以看见多少个点
将题意稍微转化一下就变成了，求$(a,b)（1<=a<=n,1<=b<=m,gcd(a,b)=1)$的点对数，由于此题n,m范围均为1e5,所以我们只要枚举a，计算1-m范围内有多少个与a互质的数就可以了，可以发现这是个经典问题，所以就解决了。

#include<stdio.h>
#include<iostream>
#include<algorithm>
using namespace std;
vector<int> v;
int cal(int n,int x)//计算1-n中有多少个与x互质的数
{
    int ans=0;
    v.clear();
    for(int i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            v.push_back(i);
            while(x%i==0) x/=i;
        }
    }
    if(x>1) v.push_back(x);
    int sz=v.size();
    for(int i=0;i<(1<<sz);i++)
    {
        int num=0;
        int tmp=1;
        for(int j=0;j<sz;j++)
        {
            if(i&(1<<j))
            {
                num++;
                tmp*=v[j];
            }
        }
        tmp=n/tmp;
        if(num&1) ans-=tmp;
        else ans+=tmp;
    }
    return ans;
}
int main()
{
    int n,m;
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        long long ans=0;
        for(int i=1;i<=n;i++)
        {
            ans+=cal(m,i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}