POJ 1411
1 #include<iostream> 2 #include<stdio.h> 3 #include<math.h> 4 #define MAXN 50000 5 using namespace std; 6 7 void give_list(); 8 double prime[5200] = {0}; 9 int num = 5200; 10 int main() 11 { 12 //freopen("acm.acm","r",stdin); 13 int i; 14 int j; 15 double m; 16 double a; 17 double b; 18 give_list(); 19 i = 0; 20 j = 0; 21 int x; 22 int y; 23 int k; 24 int tem; 25 int max = 0; 26 while(cin>>m>>a>>b) 27 { 28 if(!m) 29 break; 30 max = 0; 31 while(prime[i]*prime[i] <= m) ++ i; 32 -- i; 33 k = i; 34 for(; i >= 0; -- i) 35 { 36 for(j = k; j < num; ++ j) 37 { 38 tem = prime[i]*prime[j]; 39 if(tem > m||prime[i]/prime[j] < a/b) 40 break; 41 if(tem > max) 42 { 43 max = tem; 44 x = prime[i]; 45 y = prime[j]; 46 } 47 } 48 } 49 cout<<x<<" "<<y<<endl; 50 } 51 } 52 void give_list() 53 { 54 int i; 55 int j; 56 bool prim1[MAXN]; 57 memset(prim1,true,sizeof(bool)*MAXN); 58 prim1[1] = false; 59 for(i = 2; i < sqrt(long double(MAXN))+1; i ++) 60 { 61 if(prim1[i]) 62 { 63 j = 2*i; 64 while(j < MAXN) 65 { 66 prim1[j] = false; 67 j += i; 68 } 69 } 70 } 71 j = 0; 72 for(i = 2; i < MAXN; i ++) 73 { 74 if(prim1[i]) 75 { 76 prime[j] = i; 77 ++ j; 78 } 79 } 80 }
关注我的公众号，当然，如果你对Java, Scala, Python等技术经验，以及编程日记，感兴趣的话。
技术网站地址： vmfor.com
