不错的筛素数题。题意:任何数都可以分解成素数的乘积,现在你的任务是告诉过我一个数的最大素因子在素数表里面的位置,比如2的位置是1,3的位置是2,5的位置是3,等等。当然了,特殊情况1的最大素因子在素数表的位置是0。
我的解题思路:看样例明显是多组测试数据,数据范围又达到了100W,有一种这题要离线的直觉。果然,我第一次做先筛出素数表,然后针对每一个数来根据素数表算最大素因子,TLE了。既然要离线,那就先要直接算出这100W个数的最大素因子在素数表中的位置并存储起来,在筛法求素数表时我们是通过素数来筛除合数的,根据这个思想我们也可以在算合数的时候顺便求出最大素因子。由于我是使用快速线性筛法求素数表,每一个合数只被它的最小素因子筛除一次,那么这个合数的最大素因子要么是这个合数的最小素因子,要么是这个合数除以它最小素因子的商的最大素因子。这算是一种动态规划的思想来解题。
我的解题代码:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
const int N = 1000001;
bool isprime[N];
int primes[N], pn;
int maxfactor[N]; //存储最大素因子在素数表中的位置
int x;
void InitRead();
void FastSieve();
int main()
{
InitRead();
while (~scanf("%d", &x))
{
printf("%d\n", maxfactor[x]);
}
return 0;
}
void InitRead()
{
memset(maxfactor, -1, sizeof(maxfactor));
maxfactor[0] = maxfactor[1] = 0;
memset(isprime, true, sizeof(isprime));
isprime[0] = isprime[1] = false;
pn = 0;
FastSieve();
return;
}
void FastSieve()
{
for (int i=2; i<N; ++i)
{
if (isprime[i])
{
primes[pn++] = i;
maxfactor[i] = pn;
}
for (int j=0; j<pn; ++j)
{
if (i * primes[j] >= N) break;
isprime[i * primes[j]] = false;
maxfactor[i * primes[j]] = max(j + 1, maxfactor[i]);
if (i % primes[j] == 0) break;
}
}
return;
}