Hdu 1239 Calling Extraterrestrial Intelligence Again

简单模拟题。只要理解了题意，基本上都写得出的。

思路1：首先知道求的是素数，那么通过筛选法打表。由于p/q<=1 --> p <= sqrt(m); a/b <= p/q --> a*q <= b*p;p*q的值要最大，那么就枚举所有的值那么就可以求出正确的数字啦。（优化过的代码，顺便复习下最基础的东西，其实直接枚举也可以过的~）

CODE：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;

const int maxn = 1200;      //sqrt(100000);
int prime[maxn];
int vis[maxn];
int tot;

void init()    //素数打表 
{
    int i, j;
    tot = 0;
    memset(vis, 0, sizeof(vis));
    for(i = 2 ; i < maxn ; i++) if(!vis[i])
    {
        prime[tot++] = i;
        for(j = i*i; j < maxn; j+=i) vis[j] = 1;
    }
    return ;
}

int bsearch(int v)               //2分查找离p最近的值 
{
    int l = 0, h = tot-1;
    while(l < h)
    {
        int m = (l+h)/2;
        if(prime[m] == v) return prime[m];
        else if(prime[m] > v) h = m-1;
        else l = m+1;
    }
    return prime[h-1];
}


int main()
{
    int m, a, b;
    int p, q, index;
    init();
    while(~scanf("%d%d%d", &m, &a, &b), m, a, b)
    {
        int i, j;
        index = 0;
        int index = bsearch((int)sqrt(m));
        int max = -1;
        for(i = index ; i >= 0 ; i--) //p
        {
            for(j = tot-1; j >= i; j--) //q
            {
                if(b*prime[i] >= a*prime[j] && prime[i]*prime[j] <= m)
                {
                    if(prime[i]*prime[j] > max)
                    {
                        p = prime[i];
                        q = prime[j];
                        max = prime[i]*prime[j];
                    }
                }
            }
        }
        printf("%d %d\n", p, q);
    }
    return 0;

} 

 

 

思路2： 后来想了想，其实不需要知道p <= sqrt(m)。只需知道p <= q就行了，然后一个一个的枚举。也过了，哈哈~

 CODE：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
using namespace std;

const int maxn = 1200;      //sqrt(100000);
int prime[maxn];
int vis[maxn];
int tot;

void init()    //素数打表 
{
    int i, j;
    tot = 0;
    memset(vis, 0, sizeof(vis));
    for(i = 2 ; i < maxn ; i++) if(!vis[i])
    {
        prime[tot++] = i;
        for(j = i*i; j < maxn; j+=i) vis[j] = 1;
    }
    return ;
}


int main()
{
    int m, a, b;
    int p, q, index;
    init();
    while(~scanf("%d%d%d", &m, &a, &b), m, a, b)
    {
        int i, j;
        int max = -1;
        for(i = tot-1 ; i >= 0 ; i--) //q
        {
            for(j = i; j >= 0; j--) //p
            {
                if(a*prime[i] <= b*prime[j] && prime[i]*prime[j] <= m)
                {
                    if(prime[i]*prime[j] > max)
                    {
                        q = prime[i];
                        p = prime[j];
                        max = prime[i]*prime[j];
                    }
                }
            }
        }
        printf("%d %d\n", p, q);
    }
    return 0;

}