poj 1411 Calling Extraterrestrial Intelligence Again

链接：http://poj.org/problem?id=1411

题意：给出三个数：m,a,b(m > 4 and 0 < a/b <= 1) ,求出两个数p、q满足p*q<=m且a/b <= p/q <= 1,输出所有满足条件的pq组合中p*q最大那组。

如果直接打素数表太过直接地暴力的话肯定是要超时的。此题保持10000以内的素数就够了，然后求每两个素数之积，保持两个素数中小的那个，打好标记。表打好后，直接从输入的m开始往下找直到找到第一组能满足条件的。

代码：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#define min(a,b) (a)<(b)?(a):(b)
#define MAXN 100010
int prim[MAXN],up;
int flag[MAXN];
int aa[MAXN];
void getprim()
{
	int i,j;
	up=0;
	memset(flag,-1,sizeof(flag));
	for(i=2;i<=10000;i++)
	{
		if(flag[i]==-1)
		{
			prim[up++]=i;
			for(j=i*2;j<MAXN;j+=i)
				flag[j]=0;
		}
	}
	memset(aa,0,sizeof(aa));
	for(i=0;i<up;i++)
		for(j=0;j<up;j++)
			if((prim[i]*prim[j])>=MAXN)
				continue;
			else
				aa[prim[i]*prim[j]]=min(prim[i],prim[j]);
}
int main()
{
	getprim();
	int i,j;
	int m,a,b;
	while(scanf("%d%d%d",&m,&a,&b)&&b)
	{
		while(m)
		{
			if(aa[m])
			{
				if(aa[m]*1.0/(m/aa[m])>=(a*1.0/b))
				{
					printf("%d %d\n",aa[m],m/aa[m]);
					break;
				}
				else
					m--;
			}
			else
				m--;
		}
	}
	return 0;
}