C++X的因子链（数论、素数筛法、混合排列）

输入正整数 X，求 X 的大于 1的因子组成的满足任意前一项都能整除后一项的严格递增序列的最大长度，以及满足最大长度的序列的个数。
输入格式
输入包含多组数据，每组数据占一行，包含一个正整数表示 X。
输出格式
对于每组数据，输出序列的最大长度以及满足最大长度的序列的个数。每个结果占一行。
数据范围
1≤X≤220
输入样例：
2
3
4
10
100
输出样例：
1 1
1 1
2 1
2 2
4 6

想要因子链最长，且前一项是后一项倍数，先只考虑第一个要求，想要因子链最长，倍数尽量取小就行，倍数取质因子时是最小，再考虑第二个要求，前一项是后一项倍数，本来每项就都是因子，再结合之前的思考，把X的质因数都求出来，组合成若干项且每项除前一项都为一个质因子，就可以得到最长因子链。因为每项除前一项的商是一个质因子，于是这个数列可以等效为X的质因子的排列问题，我们要求出X的所有质因数，然后再求它们不同的排列有几个。

AC代码：

#include<stdio.h>

const int N=(1<<20)+10;

int primes[N],cnt;//存所有质数
int st[N];//当前这个数有没有被筛过
int fact[30],sum[N];//存所有质因数以及其出现次数
int minp[N];//存某个数最小质因数

int get_primes(int n)
{
    for(int i=2;i<=n;++i)
    {
    	//没被筛过，必然是素数，且其最小质因数一定是它自身
        if(!st[i]) primes[cnt++]=i,minp[i]=i;
        for(int j=0;primes[j]*i<=n;++j)
        {
            int t=primes[j]*i;
            st[t]=true;
            minp[t]=primes[j];
            //确保以最小质因数筛掉t
            if(i%primes[j]==0) break;
        }
    }
}

int main()
{
    get_primes(N-1);

    int x;
    while(scanf("%d",&x)!=-1)
    {
        int k=0,tot=0;//素数因子下标、素数因子总数
        while(x>1)
        {
            int p=minp[x];//取最小质因数
            fact[k]=p;//记录新分解质因数
            sum[k]=0;//统计该质因数出现次数
            while(x%p==0)
            {
                ++sum[k];
                x/=p;
                ++tot;
            }
            ++k;
        }
        long long ans=1;//计算答案
        //以每个元素都是独一无二的眼光，进行全排列
        for(int i=1;i<=tot;++i) ans*=i;
        //相同的质因数的排列顺序无所谓，除以它们的全排列
        for(int i=0;i<k;++i)
            for(int j=1;j<=sum[i];++j)
                ans/=j;
        printf("%d %lld\n",tot,ans);
    }
    return 0;
}