HDU ACM 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): 3934    Accepted Submission(s): 2065

【题目链接】http://acm.hdu.edu.cn/showproblem.php?pid=1239

【解题思路】题目感觉是水过去，应该数据也是很水的一回事。题目长得也水，但最终水的还是自己，手贱的敲了那么长的时间，首先是筛出素数，找出最接近面积的素数作为一个待定的长然后一直在a/b 到 1 之间得到 待定的高，及时更新最大的面积，更新完待定的高后再继续更新待定的长，两个循环 PS：代码有点乱，因为我也是写给我自己的，目的是回看的时候还能找到点痕迹：这点痕迹就是，#ifndef 后面只能带ONLINE_JUDGE 写成其他比如说Joseph是不可以读取的；第二点是小心处理数据的范围

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#define SIZE 60000

using namespace std;

int prime[SIZE], digit[SIZE];

int init_prime()
{
    memset(digit, 0, sizeof(digit));
    int sum = 0;
    for(int i=2; i<=SIZE; ++i)
    {
        if(!digit[i])
        {
            prime[sum++] = i;
    //        printf("%d ", i);
            for(int j=i+i; j<=SIZE; j+=i) digit[j] = 1;
        }
    }
    return sum;
}


int main()
{
    #ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
    #endif ONLIN_JUDGE
    init_prime();
//    long long n;
    int m, a, b, p, q, area;
    while(scanf("%d%d%d", &m, &a, &b) != EOF && a+b+m)
    {
        area = 0;
        int temp = m/2;
        double rate = a*1.0/b;
        int k;
        for(k=0; prime[k] <= temp; ++k); 
    //    printf("%d\n", prime[k-1]);
        for(--k; k>=0; --k)
        {
    //        n = (long long)prime[k]*prime[k];
    //        if(area > n) { break;}
            for(int t=k; t>=0 && a*prime[k] <= b*prime[t] ; --t)
            {
                if(area < prime[t]*prime[k] && prime[t]*prime[k] <= m)
                {
                    p = prime[t];
                    q = prime[k];
                    area = p*q;
                }
            }    
        }
        printf("%d %d\n", p, q);
    }
    return 0;
}