题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2841
题意:
一个人在(0,0)点,然后前面有一个m*n的格子 ,每个格子的节点上有一棵树,问这个人站在原地能看到多少棵树
如果两棵树在一条直线上那么只能看到最前面的一棵树。
分析:
如果一个数的坐标为(a,b),那么坐标为(a*k,b*k)的都不能看见,如果a,b有公因子c那么我们肯定只能看到(a/c,b/c);
因此我们得出结论,能看到的树的横纵坐标一定互质。那么我们就可以把问题转化为,从[1,n]有多少个数与m里的数互质。
那么只用把[1,n]里的数都素因子分解了,然后容斥一下就搞定了。
代码如下:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
using namespace std;
const int maxn = 100010;
typedef long long LL;
int Prim[maxn],num[maxn][20];
void init()
{
int i,j;
memset(Prim,0,sizeof(Prim));
for(i=1; i<=100000; i++) num[i][0]=0;
for(i=2; i<=100000; i++)
if(Prim[i]==0)
{
num[i][1]=i;
num[i][0]++;
for(j=i*2; j<=100000; j+=i)
{
num[j][++num[j][0]]=i;
Prim[j]=1;
}
}
}
LL dfs(int id,int b,int now)//求不大于b的数中,与now不互质的数的个数;
{
LL ans=0;
for(int i=id;i<=num[now][0];i++)
ans+=b/num[now][i]-dfs(i+1,b/num[now][i],now);
return ans;
}
int main()
{
int m,n,t;
cin>>t;
while(t--){
cin>>m>>n;
init();
long long sum=0;
for(int i=2;i<=m;i++)
sum+=n-dfs(1,n,i);
printf("%I64d\n",sum+n);
}
return 0;
}