POJ 3421 X-factor Chains（构造）

这条链依次乘一个因子。因为n<2^20，sqrt(n)分解因子，相同的因子相对顺序取一个。

组合公式计算一下就好。

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<queue>
#include<vector>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<cmath>
//#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
#define PB push_back

const int maxn = <<;

ll fac[];
void calFac(int n = )
{
    fac[] = ;
    for(int i = ; i <= n; i++){
        fac[i] = fac[i-]*i;
    }
}

int primeDecompose(int n,ll &div)
{
    int re = ; div = ;

    for (int i = ; i*i <= n; i++){
        int cnt = ;
        while (n % i == ){
            n /= i;
            cnt++;
            re++;
        }
        if(cnt) div *= fac[cnt];
    }
    if(n>) re++;
    return re;
}

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    calFac();
    int X;
    while(~scanf("%d",&X)){
        ll div;
        int len = primeDecompose(X,div);
        ll shc = fac[len]/div;
        printf("%d %I64d\n",len,shc);
    }
    return ;
}

我尝试用预先筛sqrt(n)以内的素数来加速判断，结果却慢了26ms..

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<queue>
#include<vector>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<cmath>
//#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
#define PB push_back

const int maxn = <<;
const int sqrn = <<;//(int)sqrt(maxn+0.5); + 5
const int PrmSz = (<<)+; // (sqrn>>2) + 5

bool isNot[sqrn];
int Prm[PrmSz], tot;

void seive(int n = sqrn) //sqrt(sqrt(maxn))
{
    int m = sqrt(n+0.5);
    isNot[] = isNot[] = true;
    for(int i = ; i <= m; i++){
        if(!isNot[i])
        for(int j = i*i; j <= n; j += i){
            isNot[j] = true;
        }
    }
    for(int i = ; i <= n; i++){
        if(!isNot[i]) //Prm.PB(i);
            Prm[tot++] = i;
    }
}

ll fac[];
void calFac(int n = )
{
    fac[] = ;
    for(int i = ; i <= n; i++){
        fac[i] = fac[i-]*i;
    }
}

inline int primeDecompose(int n,ll &div)
{
    int re = ; div = ;
    for (int i = ; i < tot && Prm[i] <= n; i++){
        int p = Prm[i], cnt = ;
        while (n % p == ){
            n /= p;
            cnt++;
            re++;
        }
        if(cnt) div *= fac[cnt];
    }
    if(n>) re++;
    return re;
}

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    seive();
    calFac();
    int X;
    while(~scanf("%d",&X)){
        ll div;
        int len = primeDecompose(X,div);
        ll shc = fac[len]/div;
        printf("%d %I64d\n",len,shc);
    }
    return ;
}