hdu 1239 Calling Extraterrestrial Intelligence Again 枚举

http://acm.hdu.edu.cn/showproblem.php?pid=1239

题意：

a. 给定整数 m,a,b( 4 < m <= 100000 and 1 <= a <= b <= 1000) 

b . 需要找到两个数 ( 不妨设为 p,q) 满足以下条件：

l       p,q 均为质数；

l     p*q<=m;

l     a/b <= p/q <= 1;

c. 输出所有满足以上条件的 p,q 中乘积最大的一对 p,q

分析：

考虑大于10000的某个质数，不妨设为Q，另一个质数为P，则：

l 如果 P<10 ， P/Q<0.001

l 如果 P>10 ， P*Q>100000

而考虑到 a,b 的取值范围 (1<=a<=b<=1000)

   可知 min(a/b)=0.001

   同时，要求： p*q<=m<=100000

所以无论如何质数都不能超过 10000 。（事实上，不会超 过 9091 ）

 

 

 

 

#include<stdio.h>
#include<stdlib.h>
int main()
{
   int m,a,b,p[10000],i,j,max,ans1,ans2,num;
    for(i=0;i<10000;i++)
    p[i]=i;
    for(i=2;i<10000;i++)
    if(p[i])
    {
        for(j=i+1;j<10000;j++)
        if(p[j]%p[i]==0)
        p[j]=0;
    }
    j=0;
    for(i=1;i<10000;i++)
    if(p[i])
    p[j++]=p[i];
    num=j;
    while(scanf("%d%d%d",&m,&a,&b),m+a+b)
    {
  max=0;
  for(i=num-1;i>=0;i--)
  {
   if(p[i]>m)
   continue;
   for(j=i;j>=0;j--)
   {
    if(a*1.0/b>p[j]*1.0/p[i])
    break;
    if(p[i]*p[j]>m)
    continue;
    if(p[i]*p[j]>max)
    {
     max=p[i]*p[j];
     ans1=p[j];
     ans2=p[i];
    }
    break;
   }
  }
  printf("%d %d\n",ans1,ans2);
 }
 return 0;
}