hdu--5108--数论

数据因为很大 达到 max=1e9

首先 素数筛选出 sqrt(max)的范围内有多少素数

然后 对于每个n  求出它的所有因子 sort一遍后 从小到大开始遍历过去

特判下 对于n大于 sqrt(max)的情况下 是否为素数的情况 因为N>1e6的情况不超过100组 所以不会特判很多次的

#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;

const int size = 100000;
bool vis[size+10];
vector<int>ve;

void init( )
{
    memset( vis , true , sizeof(vis) );
    vis[0] = vis[1] = false;
    for( int i = 2 ; i<=size ; i++ )
    {
        if( vis[i] )
        {
            int j = i * 2;
            while( j<=size )
            {
                vis[j] = false;
                j += i;
            }
        }
    }
}

void solve( int n )
{
    ve.clear();
    for( int i = 1 ; i*i<=n ; i++ )
    {
        if( n%i==0 )
        {
            ve.push_back( i );
            if( i*i!=n )
                ve.push_back( n/i );
        }
    }
    sort( ve.begin() , ve.end() );
}

bool judge( int x )
{
    if( x%2==0 )
        return false;
    for( int i = 3 ; i*i<=x ; i+=2 )
    {
        if( x%i==0 )
            return false;
    }
    return true;
}

int main()
{
    int n , ans , cnt;
    bool flag;
    init( );
    while( ~scanf("%d",&n) )
    {
        solve( n );
        cnt = ve.size();
        ans = 0;
        flag = false;
        for( int i = 0 ; i<cnt ; i++ )
        {
            if( n/ve[i]>size )
            {
                flag = judge( n/ve[i] );
                if( flag )
                {
                    ans = ve[i];
                    break;
                }
            }
            else if( vis[ n/ve[i] ] )
            {
                ans = ve[i];
                break;
            }
        }
        printf( "%d
",ans );
    }
    return 0;
}