noip2009 luogu1072 Hankson的趣味题

变量x。
已知
 gcd(x,a0)=a1  lcm(x,b0)=b1  gcd(x,b0)=xb0lcm(x,b0)=xb0b1
化简得  gcd(xa1,a0a1)=1  gcd(b1x,b1b0)=1 ,显然x为a1倍数且为b1因数。


以上是我看的别人的。事实上,不需推理,十分容易发现x为a1倍数且为b1因数……

#include <iostream>
#include <cstdio>
using namespace std;
int gcd(int a, int b){
    return !b?a:gcd(b, a%b);
}
int main(){
    int n, a0, a1, b0, b1, x;
    cin>>n;
    while(n--){
        scanf("%d %d %d %d", &a0, &a1, &b0, &b1);
        int cnt=0;
        for(int i=1; i*i<=b1; i++)
            if(b1%i==0){
                if(i%a1==0 && gcd(i/a1, a0/a1)==1 && gcd(b1/b0, b1/i)==1)
                    ++cnt;
                x = b1/i;
                if(x==i)    continue;
                if(x%a1==0 && gcd(x/a1, a0/a1)==1 && gcd(b1/b0, b1/x)==1)
                    ++cnt;
            }
        printf("%d\n", cnt);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值