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;
}