hdu2841Visible Trees 容斥原理

//给定一个n*m的方格，农场主从(0 , 0 )开始看 ， 只能看到
//一条直线上的第一个格子，问农场主能看到几个格子
//对于任意的坐标(x,y) ,与其在同一条直线上的坐标有(k*x , k*y)
//故而可以枚举所有的x，找其在(1,m)有多少个互质的数
#include<cstring>
#include<iostream>
#include<cstdio>
#include<vector>
using namespace std ;
typedef __int64 ll ;
const int maxn = 100010 ;
int isp[maxn] ;
vector<int> vec[maxn] ;
void set()
{
    memset(isp , 0 , sizeof(isp)) ;
    for(int i = 2;i < maxn;i+=2)
    vec[i].push_back(2) ;
    for(int i = 3;i < maxn;i+=2)
    {
        if(isp[i])continue ;
        for(int j = i;j < maxn;j+=i)
        {
            isp[j] = 1;
            vec[j].push_back(i) ;
        }
    }
}
int dfs(int pos , int x , int num)
{
    int ans = 0 ;
    for(int i = pos ;i < vec[x].size() ;i++)
    ans += num/vec[x][i] - dfs(i+1 , x , num/vec[x][i]) ;
    return ans ;
}
int main()
{
    int n , m ;
    set() ;
    int T;scanf("%d" ,&T) ;
    while(T--)
    {
        scanf("%d%d" ,  &n, &m) ;
        ll ans = m;   
        for(int i = 2;i <= n;i++)
            ans +=( m - dfs(0 , i , m)) ;
        printf("%I64d\n" , ans) ;
    }
    return  0 ;
}