牛客 Celestial Resort 质因数分解求最小公倍数 除法取模

题目描述

Madeline来到了Celeste山的天空度假山庄（Celestial Resort）。工作人员Oshiro先生希望Madeline住在这里，而这里却因为早就没有客人来光顾而脏乱不堪。

虽然Madeline不想留在这里，她的目标是登上Celeste山，但她不忍心看到Oshiro先生失落的样子，于是想要帮助Oshiro先生来整理天空度假山庄。

Madeline发现这里有n种书，每一种书都有无限本，第i种书的厚度是ai ,她想把这 n 种书在摆成一排，每一种书叠 b i 本，使得每一种书叠起来的高度相同。(b i ≥1) 即所有 1<=i,j<=n, ai * bi=aj *bj，她想知道b1+b2+b3+…+bn
的最小值是多少。

比如当前有3种书，这三本书的厚度分别为2、3、4，那么这三本书各自需要6、4、3本书可以使得每一种书叠起来的高度相同，那么总共需要13本。可以证明这种情况下需要的书的总数是最少的。

这个结果可能会很大，求结果对10^9+7取模的结果.

输入描述:

只有一组输入数据。
第一行输入一个数n，代表书的种类数。( 1<=n<=10^4 )
第二行输入{n}n个整数ai，代表第i本书的高度。( 1<=ai<=10^9 )

输出描述:

b1+b2+b3+…+bn的最小值对10^9+7取模的结果。

示例1

输入
3
1000000 999999 999998

输出
996989508

链接：https://ac.nowcoder.com/acm/contest/6591/C
来源：牛客网

关键是求n个数的最小公倍数，利用质因数分解求最小公倍数
利用STL中的map存质因数的个数，因为用数组存的话数据太大
除法取模 ( a / b )%mod = a * ( b ^ ( mod - 2 ) ) % mod
乘方用快速幂
刚开始用scanf输入，超时了，转换成快读，就能过了

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+10,M=1e9;
const int inf=0x3f3f3f3f;
const int mod=1e9+7;
inline int rd()
{
	register int w=1,s=0;
	register char ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-') w=-1;
		ch=getchar();
	} 
	while(ch>='0'&&ch<='9')
	{
		s=s*10+ch-'0';
		ch=getchar();
	}
	return s*w;
}
int w[N];
map<int,int> p;
ll ans=1;
ll poww(int a,int b)
{
	ll base=a%mod,an=1;
	while(b)
	{
		if(b&1) an=(an%mod*(base%mod))%mod;
		base=(base%mod*(base)%mod)%mod;
		b>>=1;
	}
	return an;
}
void get_s(int x,int y)
{
	if(p[x]>=y) return ;
	int z=y-p[x];
	ans=(ans%mod*poww(x,z)%mod)%mod;
	p[x]=y;
}
void ff(int x)
{
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0)
		{
			int t=0;
			while(x%i==0)
			{
				t++;
				x=x/i;
			}
			get_s(i,t);	
		}	
	}
	if(x!=1) get_s(x,1);
}
int main()
{
	int n;
	ll cnt=0;
	n=rd();
	for(int i=1;i<=n;i++)
	{
		w[i]=rd();
		ff(w[i]);
	}
	for(int i=1;i<=n;i++)
	{
		cnt=(cnt+(ans*(poww(w[i],mod-2)%mod))%mod)%mod;
	}
	printf("%lld",cnt);
	return 0;
}