HDOJ-2136-Largest prime factor 解题报告

       不错的筛素数题。题意：任何数都可以分解成素数的乘积，现在你的任务是告诉过我一个数的最大素因子在素数表里面的位置，比如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;
}