题意
求组合数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;
}