Hankson 的趣味题

题目：
https://ac.nowcoder.com/acm/problem/50558

已知正整数$a_0,a_1,b_0,b_1$，设某未知正整数$x$满足：

1. $gcd(x,a_0)=a_1$
2. $lcm(x,b_0)=b_1$

求有多少个$x$满足条件。

思路：
$a_0={\prod }_i{p}_i^{{\alpha }_i}$，$a_1={\prod }_i{p}_i^{{\beta }_i}$，$b_0={\prod }_i{p}_i^{{\gamma }_i}$，$b_1={\prod }_i{p}_i^{{\delta }_i}$($p_i|a_0$或者$p_i|a_1$或者$p_i|b_0$或者$p_i|b_1$)。所以$x={\prod }_i{p}_i^{b_i}$(因为$x|b_1$)。
对于$p_i$分类讨论：

• 当${\alpha }_i={\beta }_i$，${\gamma }_i={\delta }_i$时，$b_i\ge {\alpha }_i$且$b_i\le {\gamma }_i$
• 当${\alpha }_i={\beta }_i$，${\gamma }_i<{\delta }_i$时，$b_i\ge {\alpha }_i$且$b_i={\delta }_i$
• 当${\alpha }_i>{\beta }_i$，${\gamma }_i={\delta }_i$时，$b_i={\beta }_i$且$b_i\le {\gamma }_i$
• 当${\alpha }_i>{\beta }_i$，${\gamma }_i<{\delta }_i$时，$b_i={\beta }_i$且$b_i={\delta }_i$

#include<bits/stdc++.h>
#define ll long long
using namespace std;
struct node
{
    int x,b[30];
}a[4];
int T,si;
set<int>w;
void inite(int x)
{
    int tmp=x;
    for(int i=2; (ll)i*i<=(ll)x; i++)
    {
        if(tmp%i==0)
        {
            w.insert(i);
            while(tmp%i==0)
            {
                tmp/=i;
            }
        }
    }
    if(tmp!=1)
        w.insert(tmp);
}
void cal(int x)
{
    int cnt=0,tmp=a[x].x;
    for(set<int>::iterator it=w.begin();it!=w.end();it++)
    {
        int num=0;
        while(tmp%(*it)==0)
        {
            num++;
            tmp/=(*it);
        }
        a[x].b[cnt++]=num;
    }
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        int ans=1;
        w.clear();
        for(int i=0; i<4; i++)
            scanf("%d",&a[i].x);
        inite(a[0].x),inite(a[1].x),inite(a[2].x),inite(a[3].x);
        si=w.size();
        for(int i=0;i<4;i++)
            cal(i);
        for(int i=0;i<si;i++)
        {
            if(a[1].b[i]==a[0].b[i]&&a[2].b[i]==a[3].b[i])
                ans*=(a[2].b[i]-a[1].b[i]+1<0?0:a[2].b[i]-a[1].b[i]+1);
            else if(a[1].b[i]==a[0].b[i]&&a[2].b[i]<a[3].b[i])
                ans*=(a[3].b[i]>=a[1].b[i]?1:0);
            else if(a[1].b[i]<a[0].b[i]&&a[2].b[i]==a[3].b[i])
                ans*=(a[1].b[i]<=a[2].b[i]?1:0);
            else
                ans*=(a[1].b[i]==a[3].b[i]?1:0);
        }
        printf("%d\n",ans);
    }
    return 0;
}