【蓝桥云】最小公倍数（高精度数运算模拟）

时间限制: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;
}

注意事项：

最终结果的值可能位数非常多，故采用高精度数模拟运算。

如有问题，欢迎提出。