hdu2588GCD

//给出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 ;

}

时间: 2024-10-14 14:16:56

hdu2588GCD的相关文章