HDU 1239 Calling Extraterrestrial Intelligence Again

是个小搜索题吧！不过也要注意一下剪枝。

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

1、 p,q均为质数；

2、p*q<=m;

3、a/b <= p/q <= 1;

要找到输出所有满足以上条件的p,q中乘积最大的一对p,q。

剪枝：

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

a、如果P<10，P/Q<0.001。

b、如果P>10，P*Q>100000。

而考虑到a,b的取值范围(1<=a<=b<=1000)： 可知min(a/b)=0.001； 同时要求： p*q<=m<=100000。所以无论如何质数都不能超过10000。

AC代码：

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int q[10008],w[2013];
int main()
{
int n,ant=0,sum,i,j,maxx;
int x,y,a,b;
double ans,cot;
memset(q,0,sizeof(q));
for(i=2;i<=10000;i++)
{
if(q[i] == 0)
{
w[ant++]=i;
for(j=2*i;j<=10000;j+=i) q[j]=1;
}
}
while(scanf("%d%d%d",&n,&a,&b) && n && a && b)
{
ans=(double )a/b;
sum=maxx=0;
for( i=0;i<ant;i++)
{
for(j=i;j<ant;j++)
{
cot=(double )w[i]/w[j];
if(cot < ans || n < w[i]*w[j]) break;

sum=w[i]*w[j];
if(maxx < sum)
{
maxx=sum;
x=w[i];y=w[j];
}
}
}
printf("%d %d\n",x,y);
}
return 0;
}

路途中。。。。