P1072 [NOIP2009 提高组] Hankson 的趣味题

最新推荐文章于 2023-03-16 22:14:01 发布

guapi2333

最新推荐文章于 2023-03-16 22:14:01 发布

阅读量275

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/guapi2333/article/details/115706087

版权

数论专栏收录该内容

25 篇文章 2 订阅

订阅专栏

Description

给定 $n(n\leq 2\times 10^3)$ 组正整数 $a,b,c,d(1\leq a,b,c,d\leq 2\times 10^9)$ ，对于每一组 $a, b, c, d$ ，请求出：
$\sum_{i=1}^{d}[gcd(a,i)==b\land lcm(c,i)==d]$

Solution

$g c d$ 与 $l c m$ 在算术基本定理下的表示:

设 $A=\prod_{i=1}^{n}p_i^{a_i},B=\prod_{i=1}^{n}p_i^{b_i}$ ，其中 $p_i$ 为 $A$ 或 $B$ 的质因子，则： $gcd(A,B)=\prod_{i=1}^{n}p_i^{min(a_i,b_i)},lcm(A,B)=\prod_{i=1}^{n}p_i^{max(a_i,b_i)}$

将 $a, b, c, d$ 均改写为算术基本定理下的表示，进而考虑所有可能的 $i$ 的算术基本定理表示形式，之后便可以根据相应的质因子次数推断是否有解或答案了。

Summary

1、 $g c d$ 与 $l c m$ 在算术基本定理下的表示不失为思考相关题的一个好角度；

2、算术基本定理分解目标数并借此转化条件是一个不错的套路。

Code

#include<cstdio>
#include<iostream>
#define ri register int
using namespace std;

int T,a,b,c,d,check,ans;

void divide(int x)
{
	ri ai=0,bi=0,ci=0,di=0;
	while(a%x==0) { ++ai; a/=x; }
	while(b%x==0) { ++bi; b/=x; }
	while(c%x==0) { ++ci; c/=x; }
	while(d%x==0) { ++di; d/=x; }
	if(ai<bi||ci>di) { check=1; return; }
	if(ai==bi&&ci<di)
	{
		if(ai<=di) ans*=1;
		else ans*=1;
	}
	if(ai==bi&&ci==di)
	{
		if(ai<=ci)	ans*=(ci-ai+1);
		else	check=1;
	}
	if(ai>bi&&ci<di)
	{
		if(bi==di)	ans*=1;
		else	check=1;		
	}
	if(ai>bi&&ci==di)
	{
		if(bi<=ci)	ans*=1;
		else	check=1;		
	}
}

int main()
{
	scanf("%d",&T);
	for(ri opt=1;opt<=T;++opt)
	{
		ans=1; check=0;
		scanf("%d%d%d%d",&a,&b,&c,&d);
		for(ri i=2;i*i<=d;++i)
			if(d%i==0)
			{
				divide(i);
				if(check)	break;
			}
		if(d>1)	divide(d);
		if(check) { cout<<"0"<<'\n'; continue; }
		cout<<ans<<'\n';
	}
	return 0;
}