蓝桥杯决赛试题：求1到n的最小公倍数

为什么1小时有60分钟，而不是100分钟呢？这是历史上的习惯导致。
但也并非纯粹的偶然：60是个优秀的数字，它的因子比较多。
事实上，它是1至6的每个数字的倍数。即1,2,3,4,5,6都是可以除尽60。
1  2  3  4  5  6

我们希望寻找到能除尽1至n的的每个数字的最小整数。

不要小看这个数字，它可能十分大，比如n=100, 则该数为：
6972 0375 2297 1247 7164 5338 0893 5312 3035 5680 0

请编写程序，实现对用户输入的 n （n<100）求出1~n的最小公倍数。

例如：
用户输入：
6
程序输出：
60

用户输入：
10
程序输出：
2520

 

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define MAX 100000
#define N 101
void init(double *r)
{
    int i;
    for (i=0; i<N; i++)
        r[i]=0;
}
void cpy(double *r,double *t)
{
    int i;
    for (i=0; i<N; i++)
        t[i]=r[i];
}
void add(double *r,double *t)
{
    int i;
    double carry=0;
    for (i=N-1; i>=0; i--)
    {
        r[i] = r[i]+t[i]+carry;
        if (r[i] >= MAX)
        {
            r[i] -= MAX;
            carry = 1;
        }
        else
        {
            carry = 0;
        }
    }
}
int mod(double *r,double b)
{
    int i;
    double tmp,remain=0,xx;
    for (i=0; i<N; i++)
    {
        tmp = r[i]+remain;
        xx = modf(tmp/b,&r[i]);
        remain = fmod(tmp,b)*MAX;
    }
    return xx?1:0;
}

void main()
{
    double n,r[N],a[N],b[N],i;
    int j=0;
    init(r);
    scanf("%lf",&n);
    r[N-1]=n;
    for (i=n-1;i>=n/2;i--)
    {
        cpy(r,a);cpy(r,b);
        while(mod(a,i)!=0)
        {
            add(r,b);
            cpy(r,a);
        }
    }
    while (r[j] == 0)         j++;
    printf("%.0lf",r[j++]);
    while (j<N-1)        printf("%05.0lf",r[j++]);
    printf("\n");
}