时间限制:1s 内存限制:256MB
题目描述:
为什么 1 小时有 60 分钟,而不是 100 分钟呢?这是历史上的习惯导致。
但也并非纯粹的偶然:60 是个优秀的数字,它的因子比较多。
事实上,它是 1 至 6 的每个数字的倍数。即 1,2,3,4,5,6 都是可以除尽 60。
我们希望寻找到能除尽 1 至 nn 的的每个数字的最小整数。
不要小看这个数字,它可能十分大,比如 nn = 100, 则该数为:
69720375229712477164533808935312303556800
输入描述:
输入一个数字 N (N<100)。
输出描述:
输出出 1 ~ n 的最小公倍数。
示例
输入:
6
输出:
60
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
using namespace std;
void trans(vector<int> &a, int b) //低精度b转换成高精度a
{
a.clear();
if(b == 0)
{
a.push_back(0);
return;
}
while(b)
{
a.push_back(b % 10);
b /= 10;
}
return;
}
vector<int> mul(vector<int> A, int b) //高精度A乘低精度b
{
vector<int> C;
int t = 0;
for(int i = 0; i < A.size() || t; i++)
{
if(i < A.size()) t += A[i] * b;
C.push_back(t % 10);
t /= 10;
}
return C;
}
vector<int> div(vector<int> A, int b, int &r) //高精度A除以低精度b,余数r
{
vector<int> C;
r = 0;
for(int i = A.size() - 1; i >= 0; i--)
{
r = r * 10 + A[i];
C.push_back(r / b);
r %= b;
}
reverse(C.begin(), C.end());
while(C.size() > 1 && C.back() == 0) C.pop_back(); //去除前导0
return C;
}
int gcd(vector<int> a, int b) //求高精度与低精度的最大公约数
{
int r;
div(a, b, r);
while(r)
{
int t = b % r;
b = r;
r = t;
}
return b;
}
int main()
{
int n; cin >> n;
vector<int> A, B;
A.push_back(1);
for(int i = 1; i <= n; i++)
{
trans(B, i);
int r;
A = div(mul(A, i), gcd(A, i), r); //求高精度A和低精度i的最小公倍数
}
for(int i = A.size() - 1; i >= 0; i--) cout << A[i];
return 0;
}
注意事项:
最终结果的值可能位数非常多,故采用高精度数模拟运算。
如有问题,欢迎提出。