2018.10.26 poj3421X-factor Chains（数论+排列组合）

最新推荐文章于 2019-11-30 22:11:03 发布

SC.ldxcaicai

最新推荐文章于 2019-11-30 22:11:03 发布

阅读量196

点赞数 1

分类专栏： # 组合数学 # 数论文章标签：数学

本文链接：https://blog.csdn.net/dreaming__ldx/article/details/83420317

版权

数学专题同时被 3 个专栏收录

72 篇文章 1 订阅

订阅专栏

组合数学

35 篇文章 0 订阅

订阅专栏

数论

28 篇文章 0 订阅

订阅专栏

这是一道排列组合入门题目，要求求解X=p1a1p2a2..pkak中答案1（即各质因数指数之和）和答案2（各质因数阶乘的乘积除以所有质因数指数阶乘的乘积）。解题策略是先进行质因数分解，然后统计指数并计算结果。

摘要由CSDN通过智能技术生成

传送门
排列组合入门题。

令 $X=p_1^{a_1}p_2^{a_2}..p_k^{a_k}$
那么答案1就等于 $\sum_{i=1}^ka_i$
答案2等于 $\frac {(\sum_{i=1}^ka_i)!} {\prod_{i=1}^ka[i]!}$
边分解质因数边统计就行了。
代码：

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=1e6+5;
int tot=0,n,pri[N];
bool vis[N];
long long ans1,ans2,fac[25];
inline void init(){
	fac[0]=1;
	for(int i=1;i<=20;++i)fac[i]=fac[i-1]*i;
	for(int i=2;i<=1e6;++i){
		if(!vis[i])pri[++tot]=i;
		for(int j=1;j<=tot&&i*pri[j]<=1e6;++j){
			int k=i*pri[j];
			vis[k]=1;
			if(!(i%pri[j]))break;
		}
	}
}
int main(){
	init();
	while(~scanf("%d",&n)){
		ans1=0,ans2=1;
		for(int i=1;pri[i]*pri[i]<=n;++i){
			if(!(n^1))break;
			int a=pri[i],b=0;
			if(n!=n/a*a)continue;
			while(n==n/a*a)n/=a,++b;
			ans1+=b,ans2*=fac[b];
		}
		if(n^1)++ans1;
		cout<<ans1<<' '<<fac[ans1]/ans2<<'\n';
	}
	return 0;
}