POJ1150_The Last Non-zero Digit_阶乘性质

题意

求组合数C(m, n)的最后一位非零数。

思路

求 n! 中含质因子 x 的个数
举例：求 10! 中质因子 3 的个数:
10! = 10 * 9 * 8 * 7 * 6 * 5 * 4 * 3 * 2 * 1
不能被整除 3 的数不可能贡献出质因子 3，将其删去得到 9 * 6 * 3
从每个数中提取一个因子 3，实际共提取出 10 / 3 = 3 个，然后序列变成了 3 * 2 * 1 = 3!，即变成了子问题：求 3! 中质因子 3 的个数。
归纳上述规律，得到递归公式：
n! 中质因子 x 的个数 fx(n) = n / x + fx(n / x) 【n > 0】 fx(0) = 0.

原题中求的是最后一个非0位，需要除去所有的因子10。但 10 是合数，不能按照上面的公式处理。但是可以把 10 分解成 2 和 5。除去了成对的 2 和 5， 即可除去因子 10。

接下来要求的，是去掉 2 和 5 后剩下的序列乘积最后一位是几。去掉 2 和 5 之后，末尾只有 1 3 7 9，实际只需求解出3 7 9的个数即可。
再次以 10! 为例：
10! = 10 * 9 * 8 * 7 * 6 * 5 * 4 * 3 * 2 * 1
可以将数列按奇偶分开: 1 3 5 7 9 和 2 4 6 8 10
其中偶数列在每个数除去一个因子2后变成了 1 2 3 4 5
由此又得到一个递归式子，n! 剩下的序列中末尾为 x (x = 1 3 7 9) 的数的个数：
f(n) = f(n / 2) + g(n) 其中 g(x) 为奇数列中末尾为 x 的数的个数【x > 0】, f(0) = 0.

然后就是求解 g(n)
观察奇数列 1 3 5 7 9 11 13 15 17 19 21…
也可以分成两部分 1 3 7 9 11 13 17 19 21… 和 5 15 25 …
对于前一部分，在原数列中每 10 个数 1 3 7 9 各出现一次。如果 n % 10 >= x，则需要额外加上
而后一部分除去 5 后就变成了 1 2 3 4 5 …
归纳递归式子
g(n) = x / 10 + (n % 10 >= x) + g(n / 5) 【n > 0】, g(0) = 0.

最后就是将所有的 3 7 9 和 多余的 2 (5) 乘起来模 10 就好了。这一过程可以利用循环减少计算。因为这4个数的幂的末尾是长度为4的循环。特别的，多余的是5，那答案一定是5。

链接

http://poj.org/problem?id=1150

代码

#include<cstdio>
#include<iostream>

using namespace std;

typedef long long LL;

//末尾为3、7、9的数的幂的末尾的循环
const int Three[4] = {1, 3, 9, 7};
const int Seven[4] = {1, 7, 9, 3};
const int Nine[4] = {1, 9, 1, 9};
//2的幂的末尾的循环（2的0次方1不在循环节中）
const int Two[4] = {6, 2, 4, 8};

int n, m;

//n的阶乘中的因子x数
int fun0(int n, int x){
    if(n == 0) return 0;
    return fun0(n / x, x) + n / x;
}

//1~n的奇数列中末尾为x的数的个数
int fun2(int n, int x){
    if(n == 0) return 0;
    return n / 10 + (n % 10 >= x) + fun2(n / 5, x);
}

//数列1~n中末尾为x的数的个数
int fun1(int n, int x){
    if(n == 0) return 0;
    return fun1(n / 2, x) + fun2(n, x);
}

int main(){
    while(scanf("%d %d", &n, &m) == 2){
        m = n - m;

        int two = fun0(n, 2) - fun0(m, 2);
        int five = fun0(n, 5) - fun0(m, 5);
        int d = two - five;

        int three = fun1(n, 3) - fun1(m, 3);
        int seven = fun1(n, 7) - fun1(m, 7);
        int nine = fun1(n, 9) - fun1(m, 9);

        int res = 1;
        if(d < 0){
            res = 5;
        }else{
            res = res * Three[three % 4] * Seven[seven % 4] * Nine[nine % 4] % 10;
        }
        if(d > 0) res = (res * Two[d % 4]) % 10;

        cout << res << endl;
    }

    return 0;
}