Codeforces 1499 D - The Number of Pairs （数学 + 整数条件推关系）

题意：
给定你$c,d,x$三个数然问你存在多少对二元组$(a,b)$使得等式$c\ast lcm(a,b)-d\ast gcd(a,b)=x$恒成立，$(a,b)$和$(b,a)$算作两个不同的二元组。
纯抄袭大佬思路

思路一：
先把给定的式子转化一下，$c\ast lcm(a,b)=d\ast gcd(a,b)+x$，两边同除$gcd(a,b)$，得$\frac{lcm(a,b)}{gcd(a,b)}=\frac{\frac{x}{gcd(a,b)}+d}{c}$，因为左边肯定是整数，所以这就决定了$gcd(a,b)$必须是$x$的因子。

$\frac{lcm(a,b)}{gcd(a,b)}$就是$a$和$b$中不相同的素因子这里不相等是指两边把相同幂次的素因子约去之后剩余的部分，可知这一部分的因子一定是一部分属于$a$，一部分属于$b$，二者是不相交的，所以本质就相当于我从$n$种素因子中挑出两堆分给$a$和$b$，但是还有$1$和这个数本身也是可以组合的，所以就是$\left(\genfrac{}{}{0ex}{}{n}{0}\right)$ + $\left(\genfrac{}{}{0ex}{}{n}{1}\right)$ + $....$ + $\left(\genfrac{}{}{0ex}{}{n}{n}\right)$ = $2^n$，其中$n$代表素因子个数。

特别地，挑0，那么另一边就是本身那个数，挑n，那么另一边就是0。

思路二：
另一种思路，$c\ast \frac{a\ast b}{gcd(a,b)}=x+d\ast gcd(a,b)$，两侧同除$gcd(a,b)$得$\frac{a\ast b}{gcd(a,b)\ast gcd(a,b)}=\frac{\frac{x}{gcd(a,b)}+d}{c}$,设$A=\frac{a}{gcd(a,b},B=\frac{b}{gcd(a,b)}$，可知$A$和$B$互质，故其实就变为找到两个互质的数使它们相乘等于右侧，而右侧等式又可知$gcd(a,b)$定为$x$的因子，所以我们枚举$x$的因子，然后找右侧这个数的素因子即可，再分配成两个部分相乘，同理也是$2^n$的数量。

其实$gcd(a,b)|x$这个条件也可以有裴蜀定理的推论得到，若$ax+by=z$，对于整数$x$和$y$来说有解，那么$gcd(a,b)|z$，而题目里的$gcd(lcm(a,b),gcd(a,b))=gcd(a,b)$。

需要注意的是，要预处理出每个数的素因子个数，每个数现算的话，会TLE，还有当对于数学式存在分式的同时，给出整数这个条件，可能就是突破口。

#include <bits/stdc++.h>

using namespace std;

#define pb emplace_back
#define MP make_pair
#define pii pair<int,int>
#define pll pair<ll,ll>
#define lson rt<<1
#define rson rt<<1|1
#define CLOSE std::ios::sync_with_stdio(false)
#define sz(x) (int)(x).size()
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-6;
const int N = 2e7;
bool is_prime[N + 10];
int Fac[N + 10],prime[N],cnt;

void Prime() {
	for(int i = 2;i <= N;i ++) is_prime[i] = true;
	for(int i = 2;i <= N;i ++) {
		if(is_prime[i]) {
			prime[++cnt] = i;
			Fac[i] = 1;
		}
		for(int j = 1;j <= cnt && prime[j] * i <= N;j ++) {
			is_prime[prime[j]*i] = false;
			if(i % prime[j] == 0) {
				Fac[prime[j]*i] = Fac[i];//这里其实i的素因子是包含prime[j]的，所以prime[j]是不能放进去的
				break;
			}
			Fac[prime[j]*i] = Fac[i] + 1;
		}
	}
}

int d,c,x;
ll cal(int x) {
	int tmp = d + x;
	if(tmp % c == 0) {
		tmp = tmp / c;
		return 1LL << Fac[tmp];
	} 
	return 0;
}

int main() {
	Prime();
	int T;scanf("%d",&T);
	while(T--) {
		scanf("%d%d%d",&c,&d,&x);
		ll ans = 0;
		for(int i = 1;i * i <= x;i ++) {
			if(x % i == 0) {
				ans += cal(x/i);
				if(x / i != i)
					ans += cal(i);
			}
		}
		printf("%lld\n",ans);
	}
	return 0;
}