HDOJ 1239

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

  1 <= a <= b <= 1000)

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

      p,q均为质数；

    p*q<=m;

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

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

题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=1239

#include <iostream>
#include <string>
#include <cstdio>
#include <cmath>
#include <vector>
#include <algorithm>
#include <sstream>
#include <cstdlib>
using namespace std;
int num=0;
bool isprime[10000];
int prime[10000];
void primemake(){
	memset(isprime,0,sizeof(isprime));
	double d=sqrt(double(10000));
	for(int i=2;i<=d;i++)
		for(int j=2;i*j<10000;j++)
			isprime[i*j]=1;
	for(int i=2;i<10000;i++)
		if(!isprime[i]){
		prime[num++]=i;
	}
}
	
int main()
{
	int m,a,b;
	primemake();
	while(cin>>m>>a>>b)
	{
		if(m==0&&a==0&&b==0)break;
		double s=double(a)/double(b);
		int p,q,mmax=0;
		for(int i=num-1;i>=0;i--)
		{
			if(prime[i]>m)continue; //注意剪枝
			for(int j=i;j>=0;j--)
			{
				if(prime[j]>m||prime[i]*prime[j]>m)continue;//注意箭枝
				if((double(prime[j])/double(prime[i]))>=s){
					if(prime[i]*prime[j]>mmax)
					{
						p=prime[j];
						q=prime[i];
						mmax=prime[i]*prime[j];
						break;
					}
				}
			}
		}
		cout<<p<<" "<<q<<endl;
	}
	return 0;

}