【搜索剪枝】HOJ 1239 Calling Extraterrestrial Intelligence Again

此题就是求符合条件的数。

首先要打出10000以内的素数表。

 

为什么是10000以内？？

应为根据题目条件，我们只需要

 

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

而且在搜索过程中要注意剪枝

也就是当if (Arr[i]*Arr[j]>m || D Arr[j]/D Arr[i]>1) break;

下面看代码：

#include <stdio.h> #include <string.h> #include <math.h> #define MAX 10000 #define D (double) void GetPrime(int Prime[]) { int i,j,qrt; for (i=4;i<=MAX;i+=2) Prime[i] = 1; qrt = sqrt(MAX); for (i=3;i<=qrt;i++) for(j=i*i;j<=MAX;j+=i) { if (Prime[j] == 1) continue; Prime[j] = 1; } } int Prime[MAX+1]; int main() { int i,j,k; int Arr[1250]; memset(Prime, 0, sizeof(Prime)); memset(Arr, 0, sizeof(Arr)); GetPrime(Prime); for (i=2,k=0;i<=MAX;i++) if (Prime[i] == 0) Arr[k++] = i; int m,a,b; int p,q; int max; while (scanf("%d%d%d", &m, &a, &b)!=EOF,m!=0) { max = 0; for (i=k-1;i>=0;i--) { for (j=0;j<k;j++) { if (Arr[i]*Arr[j]>m || D Arr[j]/D Arr[i]>1) break; if (max < D Arr[j]*D Arr[i]&&D Arr[j]/D Arr[i]<=1&&D Arr[j]/D Arr[i]>=D a/D b&&D Arr[j]*D Arr[i]<=m) { p = Arr[j];q = Arr[i];max = Arr[j]*Arr[i]; } } } printf("%d %d/n", p, q); } return 0; }