Brave Balloonists

题目来自:http://acm.timus.ru/problem.aspx?space=1&num=1049

 

题目的大意是: 

   输入10个整数:a1,a2,...., a10.

   假设他们的乘积:k = a1*a2*....*a10.

   求k的公约数的个数N的个位数上是什么数字。

   

  解这个题目的关键在于,怎样求出一个数的公约数的个数.(如6的公约数的个数为4(1,2,3,6)).

  很多人一看,很简单啊,直接循环1~k,总可以求出来。

  这话虽没错,但只会使蛮力,不懂取巧,除了让自己的形象看起来更像一头牛,通常也会遇到效率上的瓶颈。

  这里输入的10个数,有可能会很大,因此直接暴力搜索约数的个数是不行的,第4个测试数据会被卡下来。

  

   事实上这类的题目用上一点点数论的知识,就可以得出很巧的答案。

   通常,在对一个整数进行分解的时候,我们常常把这个数写成素数的形式,这里我们也试试这样来分析一下。

 

  假设: 

      k = (p1^a1)*(p2^a2)*.....*(pi^ai)  (p1,...,pi是素数, ai > 0)

  

  显然,k的所有约数,肯定都是这些素数们的乘积,比如:

  6 = 2*3

  因此,6的约数为:2,3,2*3.

  刚好是: (1+1)*(1+1) - 1 = 3.

 

  问题其实就等价于:在a1个P1,a2个p2,....,ai个Pi,共(a1+a2+...+ai)个数中,通过相乘,可以找出多少个互不相等的数。

  这其实就排列组合的问题了。

  对于i个集合,每个集合各有aj个相同的元素( 1 =< j <= i),现在要从各个集合中抽出一定数量的元素,组成另一个集合,有多少种取法?

  每个集合有aj个元素,我们可以共有(aj+1)种取法(分别取0,1,...,aj个元素出来)

  因此对i个结合,共有(a1+1)*(a2+1)*....*(ai+1) - 1种取法(减1是排除掉各个集合里都只抽出0个元素的情况)

 

  现在回到我们的题目中来。

 对于, k = (p1^a1)*(p2^a2)*.....*(pi^ai)。k的约数的个数为:(a1+1)*(a2+2)*....*(ai+1) - 1

  但上面并没有包括1进去(1必然是k的约数)

 因此最后的结果是:(a1+1)*(a2+2)*....*(ai+1) - 1 + 1.

 

 由此可见,我们只要对各个输入数求一下他们的素数的表示形式,就很容易可以计算出他们的约数的个数了。

 而素数相对来说,是相当稀疏的,因此效率上必然比暴力搜索要好很多。

 

 

#include<iostream>
#include<math.h>
#include<map>
using namespace std;

#define MAXV (10002)  //整数的最大上限。

static map<int,int> si;   //用保存素数的个数,如si[3] = 5,表示3^5

//下面两个数组用于筛法求素数表。
static bool isPrime[MAXV];
static int  prime[MAXV];

//求素数。
static void MakePrime()
{
    memset(isPrime,1,sizeof(isPrime));
    memset(prime,0,sizeof(prime));

    for(int i=2; i < MAXV; ++i)
    {
        if(isPrime[i])
        {
            prime[++prime[0]] = i;
            for(int j = 2;j * i < MAXV;++j)
            {
                isPrime[i*j] = false;
            }
        }
    }
}


//获取整数num的素数表示形式
//如 6 = 2*3   12 = (2^2)*(3)
void GetDivisor(int num)
{
     
    for(int i = 1;i <= prime[0] && prime[i] <= num;++i)
    {
        int t = prime[i];
        int s = t;

        while(num % s == 0)
        {
                si[t]++;
                s *= t;
        }
    }
}



int main()//
{
    int i = 0;
    int k ;
     
    MakePrime();
    si.clear();

    while(i++ < 10)
    {
        cin >> k;
        GetDivisor(k);
    }

    __int64 num = 1;

    map<int,int>::iterator iter = si.begin();
    
    for(;iter != si.end();++iter)
        num *=1+ iter->second;


    cout << num%10;

    return 0;
}

 

 

 

 

 

     

 

 

你可能感兴趣的:(OO)