codeforces 446E E. Divisors(数论+暴力)

题目链接:

codeforces 446E


题目大意:

给出一个x,k,每次操作都会将x分解因数,得到新的序列,然后每次再分解序列中的每一个数,按照每一个数分解因数从小到大排,整体顺序不做调整。


题目分析:

因为x是10^12以内的数,我们可以 nxlog(n) 105 个答案,就是对于每个因数递归的分解到1,然后输出答案,达到要求的答案数停止,或者不可以再分解了停止,第二种的复杂度是 O(nlog2n),O(105logn)


AC代码:

#include 
#include 
#include 
#include 
#include 

using namespace std;

typedef long long LL;

LL x,k;
vector ans,temp;

void init ( )
{
    ans.clear();
    temp.clear();
    for ( LL i = 1 ; i*i <= x ; i++ )
    {
        if ( x%i ) continue;
        temp.push_back ( i );
        LL y = x/i;
        if ( y == i ) continue;
        temp.push_back ( y );
    }
    sort ( temp.begin() , temp.end() );
}

void dfs ( LL x , LL n )
{
    if ( n == 0 || x == 1 )
    {
        ans.push_back( x );
        return;
    }
    for ( int i = 0 ; i < temp.size() ; i++ )
    {
        LL tt = temp[i];
        if ( tt > x ) break;
        //cout << tt << endl;
        if ( x%tt ) continue;
        dfs ( tt , n-1LL );
        if ( ans.size() >= 100000 ) return;
    }
}

int main ( )
{
    while ( ~scanf ( "%lld%lld" , &x , &k ) )
    {   
        init();
        //cout << temp.size() << endl;
        dfs ( x , k );  
        //cout << ans.size() << endl;
        int lim = 1e5;
        lim = min ( (int)ans.size() , lim );
        for ( int i = 0 ; i < lim ; i++ )
            printf ( "%lld " , ans[i] );
        puts ("");
    }
}

你可能感兴趣的:(codeforces的数论专题)