题目描述:
输入正整数 X,求 X的大于 1 的因子组成的满足任意前一项都能整除后一项的严格递增序列的最大长度,以及满足最大长度的序列的个数。
输入格式
输入包含多组数据,每组数据占一行,包含一个正整数表示 X。
输出格式
对于每组数据,输出序列的最大长度以及满足最大长度的序列的个数。
每个结果占一行。
数据范围
1≤X≤2^20
输入样例:
2
3
4
10
100
输出样例:
1 1
1 1
2 1
2 2
4 6
解题思路:
1.将x的所有因子进行排序,切bi+1 / bi == 0;
由算术基本定理得出 X = p1^a1 + p2^a2 + ... +pkak;
所以X的最长因子链为a1 + a2 + ... +ak;
但是由于有重复的质因子,所以要将重复的进行分组排序得到不重复的答案数 = (a1+a2+...+ak)! / a1! + a2! + ... ak!
2.如何筛选出一个数的最小质因子:线性筛法 (筛出所有合数)prims[j] 表示质数,i 表示 1-n任意数。
primes j * i的最小质因数 == primes j 并保存下来。
3.对每个输入的x进行质因数分解得到质因数总个数的阶乘,并 / 每个质因数的阶乘得到答案。
参考代码:
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int N = (1<<20) + 10;
int x;
int primes[N],minp[N];
int cnt;
bool st[N];
void get_primes(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i]) minp[i] = i,primes[cnt++] = i;
for(int j=0;primes[j] * i <=n;j++)
{
int t = primes[j] * i;
st[t] = true;
minp[t] = primes[j];
if(i % primes[j]==0) break;
}
}
}
int main()
{
get_primes(N-1);
int sum[N];
while(cin>>x && x != -1)
{
int k = 0,tot = 0;
while(x>1)
{
int p = minp[x];
sum[k] = 0;
while(x%p==0)
{
x/=p;
sum[k]++;
tot++;
}
k++;
}
LL res = 1;
for(int i=1;i<=tot;i++) res *= i;
for(int i=0;i<k;i++)
for(int j=1;j<=sum[i];j++)
res /= j;
printf("%d %lld\n",tot,res);
}
return 0;
}