LightOj1024 - Eid (求n个数的最小公约数+高精度)

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
typedef long long LL;
#define N 10001
using namespace std;
const double eps = 1e-6;

int f[101], p[101], k = 0;

void Prime()
{
    for(int i=2; i<101; i++)
    {
        if(!f[i])p[k++] = i;
        for(int j=i; j<101; j+=i)
            f[j] = 1;
    }
}

int ans[1001], cnt[N];

void Mul(int a[], int num)
{
    for(int i=0; i<1000; i++)
        a[i] = a[i]*num;
    for(int i=0; i<1000; i++)
    {
        a[i+1] += a[i]/10000;
        a[i] = a[i]%10000;
    }
}

void PUTS(int a[])
{
    int i=1000;
    while(i>=0 && a[i]==0) i--;
    printf("%d", a[i--]);
    while(i>=0) printf("%04d", a[i--]);
    printf("\n");
}

int main()
{
    Prime();
    int n, T, t = 1;
    scanf("%d", &T);
    while(T--)
    {
        memset(cnt, 0, sizeof(cnt));
        memset(ans, 0, sizeof(ans));

        int num;
        scanf("%d", &n);
        for(int i=1; i<=n; i++)
        {
            scanf("%d", &num);
            for(int j=0; j<k && p[j]*p[j] <= num; j++)
            {
                int s = 0;
                while(num%p[j] == 0)
                {
                    s ++;
                    num /= p[j];
                }
                cnt[p[j]] = max(cnt[p[j]], s);
            }
            if(num > 1)
                cnt[num] = max(cnt[num], 1);
        }
        ans[0] = 1;
        for(int i=0; i<N; i++)
        {
            if(!cnt[i]) continue;
            int ret = 1;
            for(int j=0; j<cnt[i]; j++)
                ret *= i;
            Mul(ans, ret);
        }
        printf("Case %d: ", t++);
        PUTS(ans);
    }
    return 0;
}