Hankson的趣味题[NOIP2009]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

因为$\mathrm{lcm}(x,b_0)=b_1$，所以$x$是$b_1$的因数，所以$x$的质因子一定也是时$b_1$的质因子。我们可以对$b_1$的每个质因子p，计算$x$可能包含多少个$p$。

设$a_0,b_0,a_1,b_1,x$分别包含$m_a,m_b,m_c,m_d,m_x$个质因子$p$，其中$m_x$是未知数。

根据最大公约数的定义，在$\gcd (a_0,x)=a_1$中:

1. 若$m_a>m_c$，则$m_x=m_c$
2. 若$m_a<m_c$，则无解
3. 若$m_a=m_c$，则$m_x\ge m_c$即可

同理，根据最小公倍数的定义，在$\mathrm{lcm}(b_0,x)=b_1$中

1. 若$m_b<m_d$，则$m_x=m_d$.
2. 若$m_b>m_d$，则无解
3. 若$m_b=m_d$，则$m_x\le m_d$即可 .

结合两种情况，有以下结论:

1. 若$m_a>m_c,m_b<m_d,m_c=m_d$，$m_x=m_c=m_d$
2. 若$m_a>m_c,m_b=m_d,m_c\le m_d$，$m_x=m_c$
3. 若$m_a=m_c,m_b<m_d,m_c\le m_d$，$m_x=m_d$
4. 若$m_a=m_c,m_b=m_d,m_c\le m_d$，$m_x$可取$m_c$~$m_d$之间的任意值，共有$m_d-m_c+1$种取法。
5. 其他情况，$m_x$ 无解.

我们把$m_x$的取法数记为$cn{t}_p$，也就是$x$包含质因子$p$的方案有$cn{t}_p$种，根据乘法原理，根据题意的$x$数量即为连乘积:
$$\prod _{质数p\mid d}cn{t}_p$$

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define ll long long 
using namespace std;
const int N=45010;
const int inf=45000;
const int M=40000;
bool v[N];int prime[M],m;
inline void g_p()
{
	m=0;memset(v,false,sizeof(v));
	for(int i=2;i<=inf;i++)
	{
		if(!v[i])prime[++m]=i;
		for(int j=1;j<=m&&i*prime[j]<=inf;j++)
		{
			v[i*prime[j]]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int get(int &a,int b)
{
	int ans=0;
	while(a%b==0)a/=b,ans++;
	return ans;
}
int a0,a1,b0,b1;
int calc(int x)
{
	int ans=0;
	int ma=get(a0,x),mb=get(b0,x),mc=get(a1,x),md=get(b1,x);
	if(ma==mc&&mb==md&&mc<=md)ans=md-mc+1;
	else if(ma>mc&&mb==md&&mc<=md)ans=1;
	else if(ma==mc&&mb<md&&mc<=md)ans=1;
	else if(ma>mc&&mb<md&&mc==md)ans=1;
	return ans;
	
}
int main()
{
	g_p();
	int n;scanf("%d",&n);
	while(n--)
	{
		int ans=1;
		scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
		for(int i=1;i<=m&&(ll)prime[i]*prime[i]<=b1;i++)
		{
			if(b1%prime[i]==0)
				ans*=calc(prime[i]);
		}
		if(b1>1)ans*=calc(b1);
		printf("%d\n",ans);
	}
	return 0;
}