【写在前面】
明天是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次,后面再次分析就发现有很多可以减少的地方