暴力搜索=80分,对于暴力搜索这种方法来讲,已经不错了!
注意:下面是80分的代码。
正宗解法也是暴力搜索,不过搜索对象不一样 ;)
/*
* noip-2009-b2 hankson的趣味题
* mike-w
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<assert.h>
#define PSIZE 10000
long prm[PSIZE];
long ma0[PSIZE],ma1[PSIZE],mb0[PSIZE],mb1[PSIZE];
long a0,a1,b0,b1;
inline int min(int e1,int e2)
{
return e1<e2?e1:e2;
}
inline int max(int e1,int e2)
{
return e1>e2?e1:e2;
}
long gprm(long maxv)
{
long id,i,j;
prm[0]=2;
prm[1]=3;
for(i=5,id=2;i<=maxv;i++)
{
for(j=0;j<id;j++)
if(i%prm[j]==0)
j=id;
if(j==id)
prm[id++]=i;
if(id==PSIZE)
return id;
}
return id;
}
int fmul(long *s,long e)
{
int i;
for(i=0;prm[i]<=e;i++)
while(e && e%prm[i]==0)
s[i]++,e/=prm[i];
return i;
}
long solve(void)
{
int lim=0;
memset(ma0,0,sizeof(ma0));
memset(ma1,0,sizeof(ma1));
memset(mb0,0,sizeof(mb0));
memset(mb1,0,sizeof(mb1));
lim=max(lim,fmul(ma0,a0));
lim=max(lim,fmul(ma1,a1));
lim=max(lim,fmul(mb0,b0));
lim=max(lim,fmul(mb1,b1));
int i,j,t,up;
long result=1;
for(i=0;i<=lim;i++)
{
t=0;
up=max(max(ma0[i],ma1[i]),max(mb1[i],mb0[i]));
for(j=0;j<=up;j++)
if(min(j,ma0[i])==ma1[i] && max(j,mb0[i])==mb1[i])
t++;
if(t==0)
return 0;
else
result*=t;
}
return result;
}
int main(void)
{
int n,i;
scanf("%d",&n);
gprm(2000000001L);
for(i=0;i<n;i++)
{
scanf("%ld%ld%ld%ld",&a0,&a1,&b0,&b1);
printf("%ld\n",solve());
}
return 0;
}