蓝桥杯2013年第四届真题-公式求值

题目描述

输入n,  m,  k，输出下面公式的值。

其中C_n^m是组合数，表示在n个人的集合中选出m个人组成一个集合的方案数。组合数的计算公式如下：

输入格式

输入的第一行包含一个整数n；第二行包含一个整数m，第三行包含一个整数k。

数据规模和约定
对于100%的数据，n在十进制下不超过1000位，即1≤n< 10^1000，1≤k≤1000，同时0≤m≤n，k≤n。
提示
999101是一个质数；
当n位数比较多时，绝大多数情况下答案都是0，但评测的时候会选取一些答案不是0的数据；

输出格式

计算上面公式的值，由于答案非常大，请输出这个值除以999101的余数。 

样例输入

3 
1 
3 

样例输出

162

---

参考代码

#include <iostream>
#include <cmath>
using namespace std;
int zuhe(int n, int m){
    int s1 = 1, s2 = 1, s3, i, f = 0;
    for (i = n; ; i--){
        s1 *= i;
        f++;
        if (f == m)
            break;
    }
    for (i = 1; i <= m; i++)
        s2 *= i;
    s3 = s1 / s2;
    return s3;
}
int main(){
    int n, m, k, i, s1=0, s2, s3, t, d;
    cin >> n >> m >> k;
    s2 = zuhe(n,m);
    for (i = 1; i <= n; i++){
        t=pow(i, k);
        d = zuhe(n, i);
        s1+= t * d;
        t = 0; d = 0;
    }
    s3 = s1 * s2;
    cout << s3 << endl;
    system("pause");
    return 0;
}