poj3048 Max Factor

题目大意：给N个数（1 <= N <= 5,000），每个数的范围是1..20,000，输出他们中素因子最大的那个数。
解决方法：素数打表，找出20000之前的所有素数，存入prime数组，对于每个输入的数a，从prime数组最后一个数往前判断，如果a % prime[当前序号]== 0，那么将该素数就是输入数据a的最大素因子，找出所有输入数据的最大素因子，找出最大的那个对应的数据a即可，当a == 1的时候要特别处理。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;

const int maxn = 5005;
const int maxp = 20005;

int a[maxn];
int maxPfactor[maxn];

int prime[maxp];
int vis[maxp];

void sieve(int n) {
    int m = (int)sqrt(n + 0.5);
    memset(vis, 0, sizeof(vis));

    for(int i = 2; i <= m; i ++) {
        if(!vis[i]) {
            for(int j = i * i; j <= n; j += i) {
                vis[j] = 1;
            }
        }
    }
}

int gen_prime(int n) {
    sieve(n);
    int k = 0;
    for(int i = 2; i <= n; i ++) {
        if(!vis[i])
            prime[k ++] = i;
    }
    return k;
}

int main()
{
    int n;
    int k = gen_prime(maxp);
    scanf("%d", &n);
    for(int i = 0; i < n; i ++) {
        scanf("%d", &a[i]);
        for(int j = k - 1; j >= 0; j --) {
            if(a[i] == 1) {
                maxPfactor[i] = 1;
            }
            if(prime[j] <= a[i]) {
                if(a[i] % prime[j] == 0) {
                    maxPfactor[i] = prime[j];
                    break;
                }
            }
        }
    }
    int id, max = 0;
    for(int i = 0; i < n; i ++) {
        if(maxPfactor[i] > max) {
            max = maxPfactor[i];
            id = i;
        }
    }
    printf("%d\n", a[id]);
    return 0;
}