当前位置: 代码迷 >> 综合 >> HOJ 2588 GCD(约数,欧拉函数)
  详细解决方案

HOJ 2588 GCD(约数,欧拉函数)

热度:68   发布时间:2023-12-13 19:19:39.0

约数,欧拉函数
题目意思:
给出 n, m 两个数(1 <= m <= n < 1e9), 求有多少个数 x ,满足 1 <= x <= n, 使得 gcd(x, n) >= m

本题要点:
1、先求出 gcd(x, n)(假设等于 g), 那么 gcd(x / g, n / g) = 1, 联想到 欧拉函数
x / g <= n / g, 欧拉函数 phi(n / g) 表示 小于 n / g 的数,有多少个与 之互质的。
2、解答:
先求出 n的所有的约数 y,看看那些约数大于m,然后 用 n 除以 y , 得到 n / y,
然后计算 phi(n / y), 全部累计起来就是答案。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int T, n, m;int phi(int x)	//求x 的欧拉函数
{
    int ans = x, k = x;for(int i = 2; i * i <= k; ++i){
    if(x % i == 0){
    ans = ans / i * (i - 1);while(x % i == 0)x /= i;}}if(x > 1)ans = ans / x * (x - 1);return ans;
}void solve()
{
    int ans = 0;for(int i = 1; i * i <= n; ++i)	//只存放大于m的约数{
    if(n % i == 0){
    if(i >= m)ans += phi(n / i);if(i != n / i && n / i >= m){
    ans += phi(i);}}}printf("%d\n", ans);
}int main()
{
    scanf("%d", &T);while(T--){
    scanf("%d%d", &n, &m);	solve();}return 0;
}/* 3 1 1 10 2 10000 72 *//* 1 6 260 */