TYVJ 1107 Hankson的趣味题 解题报告

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
int a0, a1, b0, b1;
int prim[20000], num[20000], end;
int ans[200000], len;

int gcd(int a, int b)
{
int t;
while(b != 0){
t = a % b;
a = b;
b = t;
}
return a;
}

void dfs(int now, int pro)
{
int i;
if(now == end){
ans[len++] = pro;
return;
}
dfs(now + 1, pro);
for(i = 1; i <= num[now]; i++){
pro *= prim[now];
dfs(now + 1, pro);
}
}

void work(int n)
{
int i = 2;
int lim = sqrt(n);
while(n != 1 && i <= lim){
if(n % i == 0){
if(a0 % i == 0 && a1 % i != 0){
while(n % i == 0){
n /= i;
}
i++;
continue;
}
prim[end] = i;
num[end] = 0;
while(n % i == 0){
num[end]++;
n /= i;
}
end++;
}
i++;
}
if(n != 1){
prim[end] = n;
num[end++] = 1;
}
dfs(0, 1);
}

int main(int argc, char **argv)
{
int i;
int n, t;
scanf("%d", &n);
while(n--){
scanf("%d%d%d%d", &a0, &a1, &b0, &b1);
t = len = end = 0;
work(b1);
for(i = 0; i < len; i++){
if(gcd(ans[i], a0) == a1 && b0 / gcd(b0, ans[i]) * ans[i] == b1){
t++;
}
}
printf("%d\n", t);
}
return 0;
}