hdu 2588 GCD 欧拉函数
来源:互联网 发布:vb,action属性值 编辑:程序博客网 时间:2024/06/05 01:52
#include <cstdio>#include <cstring>#include <cmath>#include <queue>#include <iostream>#include <algorithm>using namespace std;const int maxn=1000;int e[maxn];int euler_phi(int n){ int m=(int)sqrt(n+0.5); int ans=n,i; for(i=2;i<=m;i++) { if(n%i==0) { ans=ans/i*(i-1); while(n%i==0)n/=i; } } if(n>1)ans=ans/n*(n-1); return ans;}int main(){ int T,n,m; cin>>T; while(T--) { cin>>n>>m; int i,j,k,t=0,num,ans=0; num=(int)sqrt(n+0.5); for(i=1;i<num;i++) { if(n%i==0) { if(i>=m)e[t++]=n/i; if(n/i>=m)e[t++]=i; } } if(num*num==n&&num>=m)e[t++]=num; for(i=0;i<t;i++) ans+=euler_phi(e[i]); cout<<ans<<endl; } return 0;}/* 题意:求有多少x(1<=x<=n),使得gcd(x,n)>=m; 先求n的所有大于等于m的因子,ei 答案ans=∑phi[n/ei];phi[i]为欧拉函数,为不大于i且与i互质的正整数个数 why? 对于一个与ei互质且小于等于n/ei的正整数p来说,p*ei<=n,gcd(p*ei,n)=ei;则phi[n/ei]就是1~n中的与n最大公约数是ei的个数。而n与1~n的最大公约数必定是n的因子。所以符合gcd(x,n)>=m的x为n所有大于等于m因子的倍数,用phi即可避免重复。*/