HDU2588-GCD

题意大概：给定N,M(2<=N<=1000000000, 1<=M<=N), 求1<=X<=N 且gcd(X,N)>=M的个数。

解法：数据量太大，用常规方法做是行不通的。后来看了别人的解题报告说，先找出N的约数x，

          并且gcd（x，N）>= M，结果为所有N/x的欧拉函数之和。

          因为x是Ｎ的约数，所以ｇｃｄ（ｘ，Ｎ）＝ｘ >= M；

　　　设ｙ＝Ｎ／ｘ，ｙ的欧拉函数为小于ｙ且与ｙ互质的数的个数。

　　　设与ｙ互质的的数为ｐ１，ｐ２，ｐ３，…，ｐ４

　　　那么ｇｃｄ（ｘ* pi，N）= x >= M。

          也就是说只要找出所有符合要求的y的欧拉函数之和就是答案

题意：求解小于n的数i且gcd（i，n）大于m的i的个数

分析：

对于所有小于n的最大公约数值一定是n的因子，所以，从这个方面下手找i为n的因子且i>m时求解车i的素数倍数且小于n的个数累加起来就好了。以为这个倍数最大为n/i所以求得n/i的欧拉函数值累加起来就好了。

---

#include<iostream>
#include<cmath> 
using namespace std;
typedef long long ll;
ll phi(ll x) {
	ll i, res = x;
	for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)    {
		if (!(x % i))        {
			res = res / i * (i - 1);
			while (!(x % i))            {
				x /= i;     
			}
		}
	}
	if (x > 1)    {
		res = res / x * (x - 1);
	}
	return res;
}

int main(){
	int t; 
	cin>>t;
	while(t--){
		ll n,m,ans=0;
		cin>>n>>m;
        for(ll i = 1; i * i <= n; i++)
        {   
            if(n%i!=0) continue;
            if(i>=m&&i*i!=n)
            	ans+=phi(n/i);			
			if(n/i>=m)
				ans+=phi(i);
        }
        cout<<ans<<endl;	
	}
	return 0;
}