HDU 2588 GCD【欧拉函数的运用】

http://acm.hdu.edu.cn/showproblem.php?pid=2588

题意：输入s个数

输入n m 表示从1到n的数与n的公约数大于m的数的个数

思路：

首先找出n的所有大于m的公约数k，然后求出每个对应的n/k的euler（欧拉函数）即小于n/k的数与n/k互质的个数，那么这些数与n/k互质且小于n/k，那么这些与n/k互质的数 乘以k之后那么就变成了与n公约数为k的数（k>m） 把所有的euler（n/k）相加即是答案 这是参考别人的思路的。

#include<stdio.h>
int arr[40000];
int euler(int n)
{
int ret=n,i;
for(i=2;i*i<=n;i++)
if(n%i==0)
{
ret=ret-ret/i;
while(n%i==0)
n/=i;
}
if(n>1)		ret=ret-ret/n;
return ret;
}

int main()
{
int s,n,m,i;
scanf("%d",&s);
while(s--)
{
int j=0,count=0;
scanf("%d%d",&n,&m);
for(i=1;i*i<n;i++)
{
if(n%i==0)
{
if(i >= m)
arr[j++] = i;
if(n/i >= m)
arr[j++] = n/i;
}
}

if(i*i==n&&i>=m)	arr[j++]=i;

for(i=0;i<j;i++)
count += euler(n/arr[i]);

printf("%d\n",count);
}
return 0;
}