刘春英课件说的很清楚。
#include<iostream>
#include<cmath>
using namespace std;
int prime[1000]={0,2,3};
void getprime(int n)//试除法求素数
{
int k=2;
for(int i=5;i<=n;i+=2)
{
bool f=true;
for(int j=1;prime[j]*prime[j]<=i;j++)//用前面已经求出来的素数试除,且小于等于sqrt(i)
if(i%prime[j]==0)
{
f=false;
break;
}
if(f)prime[++k]=i;//如果是素数就加入末尾
}
prime[0]=k;//保存素数的长度
}
int main()
{
int m,a,b,p=0,q=0;
getprime(7000);
while(scanf("%d%d%d",&m,&a,&b),m||a||b)
{
double bit=a*1.0/b;//注意类型
int max=0;
for(int i=1;i<=66;i++)
for(int j=i;bit*prime[j]<=prime[i];j++)//结束条件这换了个写法
{
int muti=prime[i]*prime[j];
if(muti>m)break;//乘积不大于m
if(muti>max)
{
max=muti;
p=prime[i];
q=prime[j];
}
}
printf("%d %d\n",p,q);
}
return 0;
}