题意：

给你一个n*m的矩形，在1到m行，和1到n列上都有一棵树，问你站在（0，0）位置能看到多少棵树

题解：

用（x，y）表示某棵树的位置，那么只要x与y互质，那么这棵树就能被看到。不互质的话说明前面已经有树挡住了这棵树

i是[1,m]中的任意一个数

我们可以for循环求在区间[1,n]内有多少数与i互质

求法就是容斥原理，具体见这里：HDU - 4135 容斥原理

代码：

1 /*
2 题意：
3 给你一个n*m的矩形，在1到m行，和1到n列上都有一棵树，问你站在（0，0）位置能看到多少棵树
4
5 题解：
6 用（x，y）表示某棵树的位置，那么只要x与y互质，那么这棵树就能被看到。不互质的话说明前面已经有树挡住了这棵树
7 */
8 #include
9 #include
10 #include
11 #include
12 #include
13 #include
14 using namespace std;
15 typedef long long ll;
16 ll v[10000],index,n;
17 void oula(ll n) //获取n的所有质因数
18 {
19 for(ll i=2; i<=sqrt(n); ++i) 20 { 21 if(n%i==0) 22 { 23 v[index++]=i; 24 n/=i; 25 while(n%i==0) 26 n/=i; 27 } 28 } 29 if(n>1)
30 v[index++]=n;
31 }
32 ll get_result(ll m)//返回1——m这个范围内与n有公因数的数的个数
33 {
34 ll que[10000],i,j,k,t=0,sum=0;
35 que[t++]=-1;
36 for(i=0; i<index; i++)
37 {
38 k=t;
39 for(j=0; j<k; j++)
40 que[t++]=que[j]*v[i]*(-1);
41 }
42 for(i=1; i<t; i++)
43 sum=sum+m/que[i];
44 return sum;
45 }
46 int main()
47 {
48 ll t;
49 scanf("%lld",&t);
50 while(t--)
51 {
52 ll n,m,ans=0;
53 scanf("%lld%lld",&n,&m);
54 ans=0;
55 for(ll i=1;i<=n;++i)
56 {
57 index=0;
58 oula(i);
59 ans+=(m-get_result(m)); //注意这里可不是这样写：ans+=get_result(m);
60 }
61 printf("%lld\n",ans);
62 }
63 return 0;
64 }