题目
题目链接
题解
数学 + 高精度。
如果直接按照计算多个数连续计算最小公倍数,那么显然要经过高精度乘法、高精度除法,两个高精度过于麻烦了。
换个思路,我们将每个数都分解质因数,全部数的最小公倍数必然由分解得到的质因数相乘得到,而且构成最小公倍数的每种质因子的个数等于某个数分解得到该质因子最多的个数,举个例子:n=6
,1,2,3,4,5,6
分解质因子得到{1}, {2}, {3}, {2,2}, {5}, {2,3}
,那么构成这些数的最小公倍数的质因子2
的个数应该为2
,因为4
存在两个质因子2
,提供的2
最多,我们为了保证最小公倍数能被4
整除,所以最小公倍数的质因子2
的数量就不能少于4
的质因子数量,但是质因子3
就只用取1
个,因为3
只有一个质因子3
,6
也只有一个质因子3
。这就是为什么每个质因子都要取每个数数量最多的质因子个数的原因。
代码
#include<bits/stdc++.h>
using namespace std;
int c[1000];
string mul (string s, int b) {
string res = "";
int m = s.size();
reverse (s.begin (), s.end ());
memset (c, 0, sizeof c);
for (int i = 0;i < m;i ++) {
c[i] += (s[i] - '0') * b;
c[i+1] += c[i] / 10;
c[i] %= 10;
}
while (c[m]) c[m+1] += c[m] / 10, c[m] %= 10, m ++;
for (int i = m-1;i >= 0;i --)
res += c[i] + '0';
return res;
}
int n;
int cnt[110], res[110];
string ans = "1";
int main()
{
cin >> n;
for (int i = 2;i <= n;i ++) {
int x = i;
memset (cnt, 0, sizeof cnt); // cnt 统计每个数的每种质因子个数
for (int j = 2;j <= x/j;j ++) {
while (x % j == 0) {
x /= j;
cnt[j] ++;
}
res[j] = max (res[j], cnt[j]); // res 记录全部数对于每种质因子的最大贡献数量
}
if (x > 1) res[x] = max (res[x], 1); // 不要忘记
}
for (int i = 2;i <= n;i ++)
for (int j = 0;j < res[i];j ++)
ans = mul (ans, i); // 高精度 * 低精度
cout << ans << endl;
return 0;
}