HDU 1019 Least Common Multiple

原创 2012年03月22日 12:41:56
/*
--------------------------------------------------------------------
    stratege:素数筛选 + 整数拆分
            一个整数x = p1^a1 + p2 ^ a2 + ... + pn ^ an 
            其中pi为素数。本题关键是找到存在pi^ai 中 ai的
            最大值, 并将pi^ai的结果不断相乘,即为最小公倍数
    Author : Johnsondu
    
    2012-03-21 21:51:34 Accepted1007 0 MS976 KB GNU C++
---------------------------------------------------------------------
*/


#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;


typedef __int64 LL ;
const int maxn = 78500 ;


int prime[7900] ;
int a[maxn], num ;
int mark[maxn] ;


LL res ;


void init ()
{
    int i, j, k ;
    memset (a, 0, sizeof(a)) ;
    num = 1 ;
    prime[0] = 2 ;
    for (i = 4; i < maxn; i += 2)
        a[i] = 1 ;
    for (i = 3; i < maxn; i += 2)
    {
        if (!a[i])
        {
            prime[num++] = i ;
            int tmp = 2 * i ;
            while (tmp < maxn)
            {
                a[tmp] = 1 ;
                tmp += i ;
            }
        }
    }
}


int max (int a, int b)
{
    return a > b ? a : b ;
}


int main()
{
    int i, j, k, tmp ;
    unsigned int tcase, n, cnt ;
    LL res, ans ;


    init () ;


    memset (mark, 0, sizeof(mark)) ;


    cin >> tcase ;
    while (tcase --)
    {
        cin >> n ;
        res = 1 ;
        ans = 1 ;
        for (i = 0; i < n; i ++)
        {
            cin >> cnt ;
            for (j = 0; j < num; j ++)
            {
                if (cnt < prime[j])
                    break ;
                tmp = 0 ;
                if (cnt % prime[j] == 0)
                {
                    tmp ++ ;
                    cnt /= prime[j] ;


                    while (cnt % prime[j] == 0)
                    {
                        tmp ++ ;
                        cnt /= prime[j] ;
                    }
                    mark[prime[j]] = max (tmp, mark[prime[j]]) ;
                }
            }
            if (cnt != 1)
            ans = cnt ;


        }


        for (i = 0; i < num; i ++)
            if (mark[prime[i]])
            {
                while (mark[prime[i]] --)
                    res *= prime[i] ;
            }
        res *= ans ;
        cout << res << endl ;
        memset(mark, 0, sizeof(mark)) ;
    }
    return 0 ;
}



                    
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

Least Common Multiple&&http://acm.hdu.edu.cn/showproblem.php?pid=1019

Problem Description The least common multiple (LCM) of a set of positive integers is the smallest p...

ACM-简单题之Least Common Multiple——hdu1019

ACM 简单题 Least Common Multiple hdu1019 最小公倍数

HDU 1019 Least Common Multiple

Problem Description The least common multiple (LCM) of a set of positive integers is the smallest p...

HDU1019_Least Common Multiple_GCD

Least Common Multiple Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/O...
  • yuege38
  • yuege38
  • 2017年02月12日 13:46
  • 79

HDU 1019 Least Common Multiple(最小公倍数&最大公约数熟悉)

Least Common Multiple Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/O...

HDU1019 The Least Common Multiple

题目: Problem Description The least common multiple (LCM) of a set of positive integers is the small...

HDU 1019.Least Common Multiple【多个数的最小公倍数】【自己的算法】【8月18】【8月19更新】

Least Common Multiple Problem Description The least common multiple (LCM) of a set of positive in...

杭电 1019 Least Common Multiple(最小公倍数)

Least Common Multiple Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java...

【HDOJ】【1019】Least Common Multiple

此题计算多个数值的最小公倍数,利用先求最大公约数间接求解最小公倍数的时候,由于通常会先将两个数相乘,然后再除的形式,导致先做乘法的位置溢出,经过多次提交,发现此问题。另外,此题还存在数值中有0的情况,...

HDOJ——1019:Least Common Multiple(求n个数的最小公倍数)

题目就不说了,就是求n个数的最小公倍数。其实这个问题已经困扰我好久了,今天终于有机会把它解决了。 这道题用到了递归和辗转相除法。 n个数的最小公倍数,可以先求前两个数的最小公倍数(lcm),然后再用前...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 1019 Least Common Multiple
举报原因:
原因补充:

(最多只允许输入30个字)