HDU 1239 思路值得学习

刘春英课件说的很清楚。

#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;
}