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

题目描述

Hanks 博士是 BT（Bio-Tech，生物技术) 领域的知名专家，他的儿子名叫 Hankson。现在，刚刚放学回家的 Hankson 正在思考一个有趣的问题。

今天在课堂上，老师讲解了如何求两个正整数 �1c1​ 和 �2c2​ 的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数 �0,�1,�0,�1a0​,a1​,b0​,b1​，设某未知正整数 �x 满足：

1. �x 和 �0a0​ 的最大公约数是 �1a1​；

2. �x 和 �0b0​ 的最小公倍数是 �1b1​。

Hankson 的“逆问题”就是求出满足条件的正整数 �x。但稍加思索之后，他发现这样的 �x 并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的 �x 的个数。请你帮助他编程求解这个问题。

输入格式

第一行为一个正整数 �n，表示有 �n 组输入数据。接下来的�n 行每行一组输入数据，为四个正整数 �0,�1,�0,�1a0​,a1​,b0​,b1​，每两个整数之间用一个空格隔开。输入数据保证 �0a0​ 能被 �1a1​ 整除，�1b1​ 能被 �0b0​ 整除。

输出格式

共 �n 行。每组输入数据的输出结果占一行，为一个整数。

对于每组数据：若不存在这样的 �x，请输出 00，若存在这样的 �x，请输出满足条件的 �x 的个数；

输入输出样例

输入 #1复制

2 
41 1 96 288 
95 1 37 1776 

输出 #1复制

6 
2

说明/提示

【样例解释】

第一组输入数据，�x 可以是 9,18,36,72,144,2889,18,36,72,144,288，共有 66 个。

第二组输入数据，�x 可以是 48,177648,1776，共有 22 个。

【数据范围】

• 对于 50%50% 的数据，保证有 1≤�0,�1,�0,�1≤100001≤a0​,a1​,b0​,b1​≤10000 且 �≤100n≤100。
• 对于 100%100% 的数据，保证有 1≤�0,�1,�0,�1≤2×1091≤a0​,a1​,b0​,b1​≤2×109 且 �≤2000n≤2000。

NOIP 2009 提高组 第二题

• 首先来分析一下这个题目

证明:

• 把上面的结论推广一下，得到结论�P

对于两个正整数�,�a,b，设���(�,�)=�gcd(a,b)=k，则存在���(�/�,�/�)=1gcd(a/k,b/k)=1

• 应用结论�P

• 整理一下式子

用心体会这两个式子，你会发现�x是�1a1​的整数倍而且是�1b1​的因子

~~好像这个由gcd和lcm也可以得到？~~嗯，就这样

于是得到一种解题思路

�1b​1​枚举�1b1​的因子(也就是�x)，如果这个数是�1a1​的整数倍并且满足那两个式子，则���++ans++

• code

#include<cstdio>
using namespace std;
int gcd(int a,int b) {
    return b==0?a:gcd(b,a%b);
}
int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        int a0,a1,b0,b1;
        scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
        int p=a0/a1,q=b1/b0,ans=0;
        for(int x=1;x*x<=b1;x++) 
            if(b1%x==0){
                if(x%a1==0&&gcd(x/a1,p)==1&&gcd(q,b1/x)==1) ans++;
                int y=b1/x;//得到另一个因子
                if(x==y) continue; 
                if(y%a1==0&&gcd(y/a1,p)==1&&gcd(q,b1/y)==1) ans++;
            }
        printf("%d\n",ans);
    }
    return 0;
}