藍橋杯決賽試題：求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");
}