AcWing 1295. X的因子链

题目描述

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

代码

思路：求最大长度即求：各个质因子出现的次数之和。
求满足最大长度的序列个数即求：所有质因子出现的次数的总和的阶乘 / 各个质因子出现次数的阶乘的积。

//思路：求最大长度即求：各个质因子出现的次数之和.
//求满足最大长度的序列个数即求：所有质因子出现的次数的总和的阶乘 / 各个质因子出现次数的阶乘的积
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=(1<<20)+10;
int a[N],b[N],prime[N],minn[N];
//欧式筛素数法，时间复杂度O(n)
void get_prime(int x)
{
    int num=0;
    for(int i=2;i<=x;i++)
    {
        if(!b[i])
        {
            minn[i]=i;//minn用来存储每个数的最小质因子
            prime[num++]=i;//prime用来存储所有数的质因子
        }
        for(int j=0;prime[j]*i<=x;j++)
            {
                b[prime[j]*i]=1;//b用来标记已经遍历的数
                minn[prime[j]*i]=prime[j];
                if(i%prime[j]==0)
                    break;
            }
    }
}
int main()
{
    get_prime(N);
    int n;
    while(cin>>n)
    {
        int count1=0,sum[N],k=0;
        LL res=1;//注意是1，用于计算阶乘
        while(n>1)//该循环是为了得到每个质因子以及出现的次数
        {
            int p=minn[n];//p为得到n的最小质因子
            sum[k]=0;
            while(n%p==0)//该循环是为了得到n分解后存在几个最小质因子p
            {
                count1++;//序列的最大长度即为该数可以分解成质因子的个数
                sum[k]++;//最小质因子的数目
                n/=p;
            }
            k++;//k表示质因子的数目，若n未被除尽，则继续处理下一个质因子。
        }
        for(int i=1;i<=count1;i++)//该循环从1开始！！
            res*=i;//计算count1的阶乘，即得到了所有质因子的排列数
        for(int i=0;i<k;i++)//除去重复的排列
            for(int j=1;j<=sum[i];j++)//该循环从1开始！！
                res/=j;
        cout<<count1<<" "<<res<<endl;
    }
}