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;
}