UVa 10892 LCM Cardinality (分解质因数+数学)

本题挺经典的,sqrt(n)时间复杂度分解质因数。

先看一个公式,数字a=p1^a1*p2^a2*p3^a3.........,b=p1^b1*p2^b2*p3^b3(其中pi为质因数,ai,bi为幂次),那么他们的LCM c=p1^max(a1,b1)*p2^max(a2,b2)..........,同理gcd是取每个质因数幂次最小值。


这样,如果lcm:c确定了,且c=p1^c1*p2^c2*p3^c3.........,那么得到他的两个数a,b中,对于每个pi,a,b中一个数的pi的幂次等于c1,另一个小于等于c1,这样对于每个pi,a的幂次为ci时b可以是0~ci-1,b为ci时a同理,再加上a,b同为ci的情况,对于每个pi方案数就是2*c1+1,把这些数乘法原理乘起来,除了(c,c)这组数,每个数都出现了2次,所以最后结果除2加1。


剩下的问题就是如何对10^9大小的数分解质因数。

遍历到根号n,对每个能整除n的数都一直除,这样根号n一下的质因子肯定都处理完了。可以证明,大于根号n的质因子至多有一个,且它的幂次为1。

证明:设p1,p2是n的质因子,且都大于根号n,那么n>=p1*p2,所有p1,p2不可能都大于根号n。


这样遍历到根号n后,只要判断下除剩下的n是不是等于1,大于一就说明还有一个大质数,再把结果乘3。


代码:

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#include <algorithm>
#define LL long long

int N;
int main(){

	while(~scanf("%d",&N)){
		int k=N;
		LL res=1;
		if(!N) break;
		for(int i=2;i*i<=N;i+=2){
			int tmp=1;
			while(N%i==0){
				tmp++;
				N/=i;
			}
			res*=(2*tmp-1);
			if(i==2) i--;
		}
		if(N>1) res*=3;
		res/=2;
		res+=1;
		printf("%d %lld\n",k,res);
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值