hdu 2841 Visible Trees（容斥原理）

#include <bits/stdc++.h>
using namespace std;

vector<int> p;
//计算1-n有多少数字与num互质
int solve(int num, int n)
{
p.clear();
for(int i = 2; i <= num/i; ++i)
{
if(num%i == 0)
{
p.push_back(i);
while(num%i == 0) num /= i;
}
}
if(num > 1)
p.push_back(num);
int len = p.size();
int ans = 0;
for(int i = 1; i < (1<<len); ++i)
{
int num = 0;
int mult = 1;
for(int j = 0; j < len; ++j)
{
if(i&(1<<j))
{
++num;
mult *= p[j];
}
}
if(num&1) ans += n/mult;
else ans -= n/mult;
}
return n-ans;
}

int main()
{
ios::sync_with_stdio(false);
int T,m,n;
long long res;
cin >> T;
while(T--)
{
cin >> m >> n;
if(m > n) swap(m,n);//让m作为小的
res = n;
for(int i = 2; i <= m; ++i)
res += solve(i,n);
cout << res << endl;
}
return 0;
}