标签:style blog http io os ar for sp strong
题目链接: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; }
标签:style blog http io os ar for sp strong
原文地址:http://blog.csdn.net/bigbigship/article/details/40582553