HDU 2588 GCD（欧拉函数）

最新推荐文章于 2021-02-04 00:30:41 发布

ACMiao_

最新推荐文章于 2021-02-04 00:30:41 发布

阅读量405

点赞数

分类专栏：数论文章标签：欧拉函数数论

本文链接：https://blog.csdn.net/Yeluorag/article/details/48262207

版权

数论专栏收录该内容

8 篇文章 0 订阅

订阅专栏

1e9，显然不能用常规做法。

思路：先找出n的约数t，则有gcd（t，n）= t 。令t >= m则满足了题意。

接下来要找出所有小于n / t且与n / t互素的数的个数并令其与t和n / t相乘，设这个互素的数为p，p与n / t互素，且t又为n的约数，所以显然有gcd（t * p，n）= t >= m。所以对于一个t，用欧拉函数找出与n / t互素的个数，然后找出所有满足条件的t求和，即得到答案。

这题除了思路是思路之外，就是模板题了。。

代码：

// Header.
#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <string>
#include <bitset>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <set>
#include <map>
using namespace std;

// Macro
typedef long long LL;
#define TIME cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << "s." << endl;
#define IN freopen("/Users/apple/input.txt", "r", stdin);
#define OUT freopen("/Users/apple/out.txt", "w", stdout);
#define mem(a, n) memset(a, n, sizeof(a))
#define rep(i, n) for(int i = 0; i < (n); i ++)
#define REP(i, t, n) for(int i = (t); i < (n); i ++)
#define FOR(i, t, n) for(int i = (t); i <= (n); i ++)
#define ALL(v) v.begin(), v.end()
#define Min(a, b) a = min(a, b)
#define Max(a, b) a = max(a, b)
#define put(a) printf("%d\n", a)
#define ss(a) scanf("%s", a)
#define si(a) scanf("%d", &a)
#define sii(a, b) scanf("%d%d", &a, &b)
#define siii(a, b, c) scanf("%d%d%d", &a, &b, &c)
#define VI vector<int>
#define pb push_back
const int inf = 0x3f3f3f3f, N = 1e2 + 5, MOD = 1e9 + 7;
// Macro end

int T, cas = 0;
int n, m;

// Imp
int Euler_Phi(int n) {
	int ret = n, m = (int)sqrt(n + 0.5);
	FOR(i, 2, m) if(n % i == 0) {
		ret = ret / i * (i - 1);
		while(n % i == 0) n /= i;
	} if(n > 1) ret = ret / n * (n - 1);
	return ret;
}

void Phi_Table(int n, int * phi) {
	mem(phi, 0);
	phi[1] = 1;
	FOR(i, 2, n) {
		if(!phi[i]) {
			for(int j = i; j <= n; j += i) {
				if(!phi[j]) phi[j] = j;
				phi[j] = phi[j] / i * (i - 1);
			}
		}
	}
}

int solve(int n, int m) {
	int ul = (int)sqrt(n + 0.5), ret = 0;
	FOR(i, 1, ul) if(n % i == 0) {
		if(i >= m && i < ul) ret += Euler_Phi(n / i);
		if(n / i >= m) ret += Euler_Phi(i);
	}
	return ret;
}

int main(){
#ifdef LOCAL
    IN // OUT
#endif
	
    si(T);
    while(T --) {
    	sii(n, m);
    	put(solve(n, m));
    }
    
    return 0;
}