0902-数论之组合数-POJ 3421 X-factor Chains

【写在前面】

明天是2018年9月3号，是博主开启高中生涯的第一天。进入高中的自己请一定要加倍努力，平衡好文化课和竞赛之间的关系，尽自己的努力不留遗憾就好了，你知道的竞赛是高风险的一项任务，但同时也是高回报的。为什么不去试一试呢？基础再薄弱只要肯练习，就总会补上去的，亲爱的自己，你常常会因为自己不如别人那般总能想到正解而沮丧，懊恼，甚至怀疑人生，可真的没有必要。你一步一步来吧，不管别人怎么说，怎么怀疑，只要你还愿意执着地搞竞赛，那就放手一搏！明天是全新的一天，是充满机会的一天！每天的自己都要元气满满，加油加油加油！！！！！！！！好了，开始吧，啊啊啊啊啊啊

传送门

大致题意

构造一个序列 a1,a2,a3,……,an，满足 ai < a(i+1) , ai | a(i+1) , an = x，其中 | 表示a(i+1)是a(i)的整数倍， x 给定

求这个序列最长的长度，以及在这个长度下有多少种方案数

 

分析

我们将  x  分解质因数，每一个数 ai 都是在 a( i-1 )的基础上随便乘以 x 的一个质因数，但注意只能乘1个，因为这样才能保证最后的长度是最优的情况，那么这个序列最长的长度就是 x 的质因数个数

现在来求方案数：先假设 x 一共有 num 个质因数，则可以得到 num! 种方案数，但里面肯定是有重复的，因为质因数并不是互不相同，我们再除去重复的部分。而重复的部分就是每种质因数的指数（也就是相同质因数的个数）的阶乘之乘，类似组合数的思想进行分析

写成公式：

 

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#define ll long long
#define ull unsigned long long
using namespace std;
int p[2000],x,num=0;
ull fac[25];
void fact(int x){
	for(int i=2;i*i<=x;++i){
		while(x%i==0){
			p[i]++;num++;
			x/=i;
		}
	}
	if(x!=1) num++;
}
void init(int x){
	fac[1]=1;
	for(int i=2;i<=x;++i) fac[i]=fac[i-1]*1ll*i;
}
int main(){
	init(22);//为什么是22呢？
//因为在这道题的范围下一个数的质因数个数不可能大于20，2^20就是上限了嘛
	while(scanf("%d",&x)!=EOF){
		if(x==1){printf("1 1\n");continue;} //1的话，单独搞一下
		memset(p,0,sizeof(p));num=0;
		fact(x);
		printf("%d ",num);
		ll ans2=fac[num];
		for(int i=2;i<2000;++i){//这个地方为什么又只取到2000呢？
//其实只要稍微大于1024即可，因为一个数的质因数肯定是小于它的sqrt的，而2^20的sqrt就是2^10=1024
			if(p[i])	ans2/=fac[p[i]];
		}
		printf("%lld\n",ans2);
	}
}

一开始总是超时，就是因为傻傻的循环x次，后面再次分析就发现有很多可以减少的地方