【CF】1499D-The Number of Pairs 题解

最新推荐文章于 2024-10-06 14:42:09 发布

M铭风F

最新推荐文章于 2024-10-06 14:42:09 发布

阅读量687

点赞数 16

分类专栏： codeforces题解文章标签： c++

本文链接：https://blog.csdn.net/qq_37272986/article/details/141597173

版权

codeforces题解专栏收录该内容

32 篇文章 0 订阅

订阅专栏

传送门：1499D
标签：数学

题目大意

给你三个正整数（大于零） $c, d$ 和 $x$ 。你需要找出满足等式 $c \cdot l c m (a, b) - d \cdot g c d (a, b) = x$ 的正整数对 $(a, b)$ 的数量。其中 $l c m (a, b)$ 是 $a$ 和 $b$ 的最小公倍数， $g c d (a, b)$ 是 $a$ 和 $b$ 的最大公约数。

输入：第一行包含一个整数 $t$ （ $1≤t≤10^4$ ）—— 测试用例的数量。每个测试用例都包含一行，包含三个整数 $c, d$ 和 $x$ （ $1≤c,d,x≤10^7$ ）。

输出：对于每个测试用例，打印一个整数 —— 满足上述等式的 $(a, b)$ 对的数量。

算法分析

如果将 $a$ 表示为 $A g$ ，将 $b$ 表示为 $B g$ ，则其中 $g = g c d (a, b)$ 且 $g c d (A, B) = 1$ 。根据定义 $lcm(a,b)=\frac{ab}{g}$ ，所以我们可以将 $l c m (a, b)$ 表示为 $\frac{AgB}{g}=ABg$ 。然后可以重写陈述中的方程如下： $\implies g(c·AB-d)=x$ 。由于左侧可被 $g$ 整除，右侧也应如此。所以我们可以遍历 $g$ 作为 $x$ 的约数。如果右侧的 $c·AB=\frac{x}{g}+d$ 不可被 $c$ 整除，则可以跳过这样的 $g$ 。 $AB=\frac{x}{c}+\frac{d}{c}$ （我们称它为 $k$ ）。如果 $k$ 具有素数因子 $p$ ，则 $A$ 或 $B$ 应该可被 $p$ 整除，因为 $g c d (A, B) = 1$ （ $A$ 和 $B$ 没有公共因子）。因此对于当前的 $g$ 值，有 $2^z$ 对 $(A$ 和 $B)$ ，其中x为k的素数因子个数。
我们可以预计算每个数字的最大素因数（可能需要的最大值为 $k$ ）上的最小素因数，在 $O(2·10^7loglog(2·10^7))$ 时间内使用欧拉筛法。那么在 $O(\sqrt{x} log x)$ 时间内可以为每个测试用例解决这个问题，但这还不够快。为了加快这个方法，我们可以预计算每个数字的最大素因数的数量。令 $mind_i$ 为 $i$ 的最小素因数， $val_i$ 为 $i$ 的素因数数量。那么 $val_i=val_i/mind_i$ 加上 1，如果 $mind_i≠mind_i/mind_i$ 。现在要解决问题，我们只需遍历 $x$ 的约数，所以每个测试用例的时间复杂度为 $O(\sqrt{x})$ 。

代码实现

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 1e7+5;
int tt;
int c,d,x;
int pr[2*N],iki[22];
ll ans;

void f(int nw){
	if((x/nw+d) % c) return;
	int k = (x/nw+d)/c;
	ans += iki[pr[k]];
}

int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);	
	iki[0] = 1;
	for(int i = 1;i <= 20;i++) iki[i] = iki[i-1] * 2;
	
	for(int i = 2;i < 2*N;i++){
		if(pr[i]) continue;
		for(int j = i;j < 2*N;j += i) pr[j]++;
	}
	
	cin >> tt;
	while(tt--){
		cin >> c >> d >> x;
		ans = 0;
		for(int i = 1;i*i <= x;i++){
			if(x % i) continue;
			f(i);
			if(i != x/i) f(x/i);
		}
		cout<<ans<<'\n';
	}
}