标签：10000 因数分解 int 3421 factorial 因子 POJ include LL
原题
题目分析
依题意可以知道,X每次只分解出一个质因子即可得到最长的链.像20=22*5=(2*5,5,1)就是20最长因子链.因此最长因子链的长度等于质因子的幂的和.最长因子链的组合方法由组合数学可知有(质因子的幂的和的全排列)/(每个质因子的幂的全排列)种.由题目数据可知用组合数学求解的过程中并不会爆longlong,因此不需要额外处理.
代码
1 #include <cstdio>
2 #include <cstdlib>
3 #include <iostream>
4 #include <algorithm>
5 #include <utility>
6 #include <ctime>
7 #include <cmath>
8 #include <cstring>
9 #include <string>
10 #include <stack>
11 #include <queue>
12 #include <vector>
13 #include <set>
14 #include <map>
15
16 using namespace std;
17 typedef long long LL;
18 const int INF_INT=0x3f3f3f3f;
19 const LL INF_LL=0x3f3f3f3f3f3f3f3f;
20
21 LL factorial[21];
22 bool is_prime[10000];
23 int prime[10000];
24 int cnt;
25 int num[10000];
26 int k;
27
28 void pre()
29 {
30 factorial[0]=factorial[1]=1;
31 for(int i=2;i<=20;i++) factorial[i]=factorial[i-1]*i;
32 for(int i=2;i<=9000;i++) is_prime[i]=true;
33 for(int i=2;i<=9000;i++)
34 {
35 if(is_prime[i])
36 {
37 prime[cnt++]=i;
38 for(int j=(i<<1);j<=9000;j+=i) is_prime[j]=false;
39 }
40 }
41 }
42
43 void prime_factor(int x)
44 {
45 k=0,num[k]=0;
46 LL sum=0;
47 for(int i=0;i<cnt;i++)
48 {
49 while(!(x%prime[i])) x/=prime[i],num[k]++;
50 sum+=num[k];
51 k++,num[k]=0;
52 }
53 if(x!=1) sum++;
54 printf("%lld ",sum);
55 sum=factorial[sum];
56 for(int i=0;i<k;i++) sum/=factorial[num[i]];
57 printf("%lld\n",sum);
58 }
59
60 int main()
61 {
62 // freopen("black.in","r",stdin);
63 // freopen("black.out","w",stdout);
64 pre();
65 int x;
66 while(~scanf("%d",&x)) prime_factor(x);
67 return 0;
68 }
标签：10000,因数分解,int,3421,factorial,因子,POJ,include,LL 来源： https://www.cnblogs.com/VBEL/p/11441686.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。