前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >LightOJ - 1236 Pairs Forming LCM 算数基本定理

LightOJ - 1236 Pairs Forming LCM 算数基本定理

作者头像
用户2965768
发布2019-07-31 14:55:05
4380
发布2019-07-31 14:55:05
举报
文章被收录于专栏:wymwym

a=p1 ^ a1 * p2 ^ a2 *..........*pn ^ an

b=p1 ^ b1 * p2 ^ b2 *..........*pn ^ bn

gcd(a,b)=p1 ^ min(a1,b1) * p2 ^ min(a2,b2) *..........*pn ^ min(an,bn)

lcm(a,b)=p1 ^ max(a1,b1) * p2 ^ max(a2,b2) *..........*pn ^ max(an,bn)

n = p1 ^ e1* p2 ^ e2 *.....* pn ^en

lcm(a,b) == n , 则 max(a1,b1) ==e1 max(a2,b2) ==e2 ... max(an,bn) == en;

则 当 a 的 因数pi 指数为ei 时, b 的因数 pi 指数 有 [0,ei ] 共 ei + 1种选择,

a , b交换也是一样, 对于 两种都选 ei 重复 ,则 对每个指数 ei, 有 2*(ei +1)- 1 种选择

题中要求 i < = j 的对数,除了 (n,n)这一对,另外都重复了 ans = (ans+1)/2

代码语言:javascript
复制
//LightOJ - 1236 Pairs Forming LCM 
//算数基本定理 
#include <bits/stdc++.h>
#define ll long long
#define LL long long
using namespace std;
const int maxn = 10000010;
bool is_prime[maxn];
int prime[maxn/10];
void get_prime(){
	memset(is_prime,0,sizeof(is_prime));
	for(int i=2;i<maxn;i++){
		if(!is_prime[i])prime[++prime[0]] = i;
		for(int j=1;j<=prime[0]&&prime[j]*i<maxn;j++){
			is_prime[prime[j]*i] = 1;
			if(i%prime[j]==0)break;
		}
	}
}

ll getfact(ll x){
	ll ans = 1,tmp = x;
	ll j = 1;
	while(j<prime[0]&&prime[j]*prime[j]<=tmp){
			if(tmp%prime[j]==0){
				ll cnt = 0;
				while(tmp%prime[j]==0)tmp/=prime[j],cnt++;
				ans*=(1+2*cnt);
			}
			j++;
	}
	if(tmp>1)ans*=3;	
	return (ans+1)/2;
}

int main()
{
	get_prime();
	int t,cnt=1;
	ll a;
	scanf("%d",&t);
	while(t--){
		scanf("%lld",&a);
		printf("Case %d: %lld\n",cnt++,getfact(a));
	}
	return 0;
 } 
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2019年07月13日,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档