HDU 1239 Calling Extraterrestrial Intelligence Again
来源:互联网 发布:js数组实例 编辑:程序博客网 时间:2024/06/06 08:41
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.
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 53p q <= m and a / b <= p / q <= 1;就是求满足的这条公式 写了这么长………… 还有就是求P*Q为最大 所以要跑完 不能中间就BREAK上面这个代码是优化过的 HD的数据太水了 我们学校的交上去超时了。。。优化一下才能过#include <stdio.h> #include <string.h> int p[10005],pr[2000],i,j; void pri() { bool pre[10000]; memset(p,1,sizeof p); memset(pre,1,sizeof pre); pre[1]=0; for(i=2;i<=5000;i++) { if(pre[i]==1) { for(j=2;j*i<=10000;j++) pre[j*i] = 0; } } for(i=1,j=0;i<=10000;i++) if( pre[i] == 1) pr[j++] = i;} int main() { int m,a,b,max,p,q,n; pri(); n=j; while(~scanf("%d %d %d",&m,&a,&b),m||a||b) { pq <= m and a / b <= p / q <= 1, max=0; for(i=n-1;i>=0;i--) { for(j=i;j<n;j++) { if(pr[i]*pr[j]<=m) { if(1.0*a/b<=1.0*pr[i]/pr[j]) { if(pr[i]*pr[j]>max) { p=pr[i];q=pr[j]; max=pr[i]*pr[j]; } } } } } printf("%d %d\n",p,q); } return 33; }#include <stdio.h>#include <string.h>int p[10005],sum,pr[2000],i,j;void pri(){ memset(p,1,sizeof p); for(i=2;i<=10000;i++) for(j=2;j*i<=10000;j++) p[i*j]=0; sum=0; for(i=2;i<=10000;i++) { if(p[i]) pr[sum++]=i; }}int main(){ int m,a,b,max,p,q; pri(); while(~scanf("%d %d %d",&m,&a,&b),m||a||b) { max=0; double c=1.0*a/b; for(i=sum-1;i>=0;i--) { for(j=sum-1;j>=0;j--) { if(pr[i]<=pr[j]) { if(pr[i]*pr[j]<=m) { double d=1.0*pr[i]/pr[j]; if(c<=d) { if(pr[i]*pr[j]>max) { p=pr[i];q=pr[j]; max=pr[i]*pr[j];//printf("%d %d %lf %lf\n",pr[i],pr[j],c,d); } } } } } } printf("%d %d\n",p,q); } return 0;}
0 0
- hdu 1239Calling Extraterrestrial Intelligence Again
- Hdu 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- HDU 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- HDU 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- HDU 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again
- 【HDU】 1239 Calling Extraterrestrial Intelligence Again
- HDU Calling Extraterrestrial Intelligence Again
- hdu 1239 Calling Extraterrestrial Intelligence Again (暴力枚举)
- HDU 1239 Calling Extraterrestrial Intelligence Again 素数表
- HDU 1239 Calling Extraterrestrial Intelligence Again(找素数)
- HDU 1239 Calling Extraterrestrial Intelligence Again 【打表+暴力搜索】
- HDU - 1239 - Calling Extraterrestrial Intelligence Again (素数相关~)
- Asp.net mvc4 发布你的程序到IIS7.0
- Android中的线程与进程之间的关系简单解释
- hdu 5073 Galaxy 2014 Asia AnShan Regional Contest
- 头文件中#ifndef后面的如何写及什么作用
- 判断电话号码是否正确
- HDU 1239 Calling Extraterrestrial Intelligence Again
- Grails开发环境的快速搭建
- Jni 简单实战 获取c++返回值
- eclipse luna 安装EMF插件
- MFC的工具条和状态栏
- Android-无法隐藏标题栏
- Spark脚本提交/运行/部署
- HTTP的KeepAlive是开启还是关闭?
- 目前开发手机app的几种方式