poj1141

#include<stdio.h>
#include<math.h>
#include<algorithm>
using namespace std;

double r[10002]={0.0};


//bool comp(double x,double y){if(a/b<=x/y&&x/y<=1) return true;else return false;}
int main(){
double m,n,a,b;
 while(1){
  scanf("%lf%lf%lf",&m,&a,&b);
  if(m==0&&a==0&&b==0) break;
bool prim[10002]={0};//标记i是否为素数!!!
  int i,j,k=0;
  for(i=3;i<=100;i+=2)//这里的100是为了平摊j而设的!!!
    if(prim[i]==0)   
   for(j=i;j*i<=10000;j+=2)
    prim[i*j]=1;//奇数*奇数一定是合数;
  for(i=4;i<=10000;i+=2)
   prim[i]=1;
  k=0;
  for(i=2;i<=10000;i++)
   if(prim[i]==0)
    r[k++]=i+0.0;
   
 

 /* if(m>10000) n=10000;else n=m;

  for(i=3;i<n;i+=2){
   for(j=2;j<sqrt(i+0.0)+1;j++){
    if(i%j==0)
     break;
    r[k++]=i;}}//此处的素数筛选是一种普通方法应根据题目条件优化,不然就会超时*/
 
  k--; while(r[k]*2>=m){k--;}//求出最大的素数!!!


 double ans=0.0;double p=0.0,q=0.0;

 /* for(i=k-1;i>=0;i--)
   for(j=i;j<=k;j++){
   double ans1=r[i]*r[j];
    if(r[i]*r[j]<=m&&comp(r[i],r[j])&&ans1>ans){
     p=r[i];q=r[j];ans=ans1;continue;}
     }

   if(p>q) swap(p,q);此处搜索时剪枝不够导致超时*/

 for(i=k;i>=0;i--)
  for(j=i;j<=k&&r[j]*r[i]<=m&&r[i]/r[j]>=0.001;j++)//已隐含p[i]/p[j]<=1;
   if(r[i]/r[j]>=a/b)
    if(r[j]*r[i]>ans)
   {p=r[i];q=r[j];ans=r[i]*r[j];}


 printf("%.0lf %.0lf/n",p,q); 
 }
 return 0;}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值