HDOJ 1239 Calling Extraterrestrial Intelligence Again

Calling Extraterrestrial Intelligence Again

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4871    Accepted Submission(s): 2569

Problem Description

A message from humans to extraterrestrial intelligence was sent through the Arecibo radio telescope in Puerto Rico on the afternoon of Saturday November 16, 1974. The message consisted of 1679 bits and was meant to be translated to a rectangular picture with 23 * 73 pixels. Since both 23 and 73 are prime numbers, 23 * 73 is the unique possible size of the translated rectangular picture each edge of which is longer than 1 pixel. Of course, there was no guarantee that the receivers would try to translate the message to a rectangular picture. Even if they would, they might put the pixels into the rectangle incorrectly. The senders of the Arecibo message were optimistic.
We are planning a similar project. Your task in the project is to find the most suitable width and height of the translated rectangular picture. The term "most suitable" is defined as follows. An integer m greater than 4 is given. A positive fraction a / b less than or equal to 1 is also given. The area of the picture should not be greater than m. Both of the width and the height of the translated picture should be prime numbers. The ratio of the width to the height should not be less than a / b nor greater than 1. You should maximize the area of the picture under these constraints.

In other words, you will receive an integer m and a fraction a / b. It holds that m > 4 and 0 < a / b < 1. You should find the pair of prime numbers p, q such that pq <= m and a / b <= p / q <= 1, and furthermore, the product pq takes the maximum value among such pairs of two prime numbers. You should report p and q as the "most suitable" width and height of the translated picture.

 

Input

The input is a sequence of at most 2000 triplets of positive integers, delimited by a space character in between. Each line contains a single triplet. The sequence is followed by a triplet of zeros, 0 0 0, which indicated the end of the input and should not be treated as data to be processed.

The integers of each input triplet are the integer m, the numerator a, and the denominator b described above, in this order. You may assume 4 < m <= 100000 and 1 <= a <= b <= 1000.

 

Output

The output is a sequence of pairs of positive integers. The i-th output pair corresponds to the i-th input triplet. The integers of each output pair are the width p and the height q described above, in this order.

Each output line contains a single pair. A space character is put between the integers as a delimiter. No other characters should appear in the output.

 

Sample Input

5 1 299999 999 9991680 5 161970 1 12002 4 110 0 0

 

Sample Output

2 2313 31323 7343 4337 53

 

题目大意：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

 

思想流程：a.从1—100000中搜出质数       b.两层循环，试遍所有的组合(p,q可能相等）       c.每种组合去判断是否符合条件，如是，将p*q与当前最大值比较，判断，保存

 

 

剪枝与优化：考虑大于10000的某个质数，不妨设为Q，另一个质数为P，则：  如果P<10，P/Q<0.001；如果P>10，P*Q>100000
而考虑到a,b的取值范围(1<=a<=b<=1000)；可知min(a/b)=0.001  ；同时，要求： p*q<=m<=100000，所以无论如何质数都不能超过10000


AC代码：

 

#include<stdio.h>#include<string.h>#include<math.h>#define max 100000+10//max的值依然要大于m的最大值 int primes[max];void is_primes()//打表，值为1的下标是质数，值为0的不为质数 {int i,j;memset(primes,1,sizeof(primes));primes[0]=primes[1]=0;for(i=2;i<sqrt(10000*1.0);i++)//此处可优化和缩小质数的范围，小于10000 {if(!primes[i])   continue;for(j=i*i;j<10000;j=j+i)   primes[j]=0;}}int main(){int m,a,b,p,q,i,j,sum;is_primes();while(scanf("%d%d%d",&m,&a,&b)&&a!=0||b!=0||m!=0){sum=0;p=q=0;for(i=1;i<=sqrt(m*1.0)+1;i++){for(j=i;j<=m;j++){if(primes[i]&&primes[j]&&(i*j<=m)&&(b*i>=a*j)&&(i*j>sum)){sum=i*j;//i*j>sum用于更新sum的值 p=i;q=j;}} }printf("%d %d\n",p,q);}return 0;}