hdu2588GCD

来源:互联网 发布:淘宝刷评给礼品微信号 编辑:程序博客网 时间:2024/06/06 03:31
//给出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 ;
}





0 0