NOIP-2009-B2 HANKSON 的趣味题

暴力搜索=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;
}