码迷,mamicode.com
首页 > 其他好文 > 详细

hdu2588GCD

时间:2015-05-28 21:33:05      阅读:140      评论:0      收藏:0      [点我收藏+]

标签:

//给出n , m 求出小于n且与n的最大公因数大于m的个数
//枚举所有n的大于m的因数i为最大公因数
//那么就只需要求n/i的欧拉函数就行
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 100010 ;
typedef __int64 ll ;
ll f[maxn] ;
ll get_factor(ll n)
{
    ll len = 0 ;
    for(ll i = 1;i*i <= n;i++)
    {
        if(n%i == 0)
        {
            f[++len] = i ;
            if(i != n/i)
            f[++len] = n/i ;
        }
    }
    return len ;
}
ll Euler(ll n)
{
   ll rea = n ;
   for(ll i = 2;i*i <= n;i++)
   {
       if(n%i == 0)
       rea -= (rea/i) ;
       while(n%i == 0)
       n /= i ;
   }
   if(n > 1)rea -= rea/n ;
   return rea ;
}
int main()
{
    int T ;
    scanf("%d" , &T) ;
    while(T--)
    {
        ll n , m ;
        scanf("%I64d%I64d" , &n , &m) ;
        ll len = get_factor(n) ;
        ll ans = 0 ;
        for(ll i = 1;i <= len;i++)
        if(f[i] >= m)
        ans += Euler(n/f[i]) ;
        printf("%I64d\n" , ans);
    }
    return  0 ;
}





hdu2588GCD

标签:

原文地址:http://blog.csdn.net/cq_pf/article/details/46128635

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!