算数基本定理 + 例题

算数基本定理即质因子分解定理,顾名思义,将一个数的所有质因子,及每个质因子的数量分解出来。(算数基本定理是数论的一个公理)
在代码的实现的时候不能直接暴力,太慢,一般配合欧拉筛,这里欧拉筛的作用不是素数打表,而是求所有小于N的自然数的最小质因子,及minp。及对任意一个x都让他除尽minp[x],到x % minp[x] 不为0为止,之后再让新的x 除尽 新的minp[x],x为1为止。
看道板子题
1295. X的因子链
输入正整数 X,求 X 的大于 1 的因子组成的满足任意前一项都能整除后一项的严格递增序列的最大长度,以及满足最大长度的序列的个数。

输入格式
输入包含多组数据,每组数据占一行,包含一个正整数表示 X。

输出格式
对于每组数据,输出序列的最大长度以及满足最大长度的序列的个数。

每个结果占一行。

数据范围
1≤X≤220
输入样例:
2
3
4
10
100
输出样例:
1 1
1 1
2 1
2 2
4 6

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>

using namespace std;
typedef long long ll;
const int N = (1 << 20) + 15;
int x,p[N],minp[N],cnt,c[30][30],fact[30],sum[30];
bool book[N];
void oula()
{
	for(int i = 2; i < N; i ++)
	{
		if(!book[i])
		{
			p[cnt ++] = i;
			minp[i] = i;
		}
		for(int j = 0; j < cnt && i * p[j] < N; j ++)
		{
			int temp = i * p[j];
			book[temp] = true;
			minp[temp] = p[j];
			if(i % p[j] == 0) break;
		}
	}
}
int main()
{
	oula();
	//
	for(int i = 0; i < 30; i ++)
	{
		c[i][0] = 1;
		c[i][i] = 1;
		c[i][1] = i;
	}
	for(int i = 3; i < 30; i ++)
	{
		for(int j = 2; j < i; j ++)
		{
			c[i][j] = c[i-1][j-1] + c[i-1][j];
		}
	}
	/
	while(scanf("%d",&x) != EOF)
	{
		int tot = 0,ans = 0,k = 0;
		memset(sum,0,sizeof sum);
		memset(fact,0,sizeof fact);
		while(x > 1)
		{
			int temp = minp[x];
			fact[ans ++] = temp;
			while(x % temp == 0)
			{
				x /= temp;
				sum[k] ++;
				tot ++;
			}
			//printf("x = %d temp = %d sum = %d\n",x,temp,sum[tot]);
			k ++;
		}
		ll res = 1;
		int sniper = tot;
		for(int i = 0; i < tot; i ++)
		{
			//printf("c[%d][%d] = %d\n",sniper,sum[i],c[sniper][sum[i]]);
			res *= (ll)c[sniper][sum[i]];
			sniper -= sum[i];
		}
		printf("%d %lld\n",tot,res);
	}
	return 0;
} 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值