筛素数并且需要优化技巧的题,当然暴力貌似也能过,不推荐暴力。题意:很简单,给你三个正整数m,a,b(其中4 < m <= 100000 且 1 <= a <= b <= 1000)。现在让你找到两个素数p和q,使得pq <= m 且 a/b <= p/q <= 1,并且要让pq的乘积尽可能大。输出p和q。
我的解题思路:如果是暴力的话筛出m最大值以内的素数然后开始从小到大枚举p的值就是了,很明显q必定大于等于p。我初步优化的话通过a/b <= p/q这个不等式可以推出q <= pb/a,这样又限制q的值,大于pb/a的值就不用枚举了。暴力在POJ上教的600多ms,初步优化后200多ms,提交榜上还有0ms,不行了,实在想不出还要怎么优化了。
我的解题代码:
#include <stdio.h>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long Long;
const int N = 100001;
const double eps = 1e-8;
bool isprime[N];
int primes[N], pn;
int m, a, b;
void InitRead();
void DataProcess();
void FastSieve(int maxn);
int main()
{
InitRead();
while (~scanf("%d %d %d", &m, &a, &b))
{
if (a == b && b == m && m == 0) break;
DataProcess();
}
return 0;
}
void InitRead()
{
memset(isprime, true, sizeof(isprime));
isprime[0] = isprime[1] = false;
pn = 0;
FastSieve(N-1);
return;
}
void DataProcess()
{
Long ansp, ansq, anspq = -1;
double temp = (double)a / b;
for (int i=0; i<pn; ++i) //枚举p
{
for (int j=i; j<pn; ++j) //枚举q
{
if ((Long)primes[i] * primes[j] > m) break;
if ((Long)a * primes[j] > (Long)b * primes[i]) break; //初步优化添加的代码
if ((Long)primes[i] * primes[j] > anspq && (double)primes[i] / primes[j] >= temp - eps)
{
ansp = primes[i];
ansq = primes[j];
anspq = (Long)primes[i] * primes[j];
}
}
}
printf("%lld %lld\n", ansp, ansq);
return;
}
void FastSieve(int maxn)
{
for (int i=2; i<=maxn; ++i)
{
if (isprime[i]) primes[pn++] = i;
for (int j=0; j<pn; ++j)
{
if (i * primes[j] > maxn) break;
isprime[i * primes[j]] = false;
if (i % primes[j] == 0) break;
}
}
return;
}