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 ;

}