C语言-蓝桥杯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<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const ll mod=999101;
ll C(ll n,ll m)//组合数 
{
    if(n<n-m)
    n=n-m;
    ll ans=1;
    for(int i=m+1;i<=n;i++)    ans*=i;
    for(int i=1;i<=n-m;i++)    ans/=i;
    return ans;
}
int main()
{
    ll n,m,k,ans=0;
    scanf("%lld%lld%lld",&n,&m,&k);
//    if(n>2*pow(10,10))
//    {
//        printf("0\n");
//        return 0;
//    }
    for(int i=0;i<=n;i++)
    {
        ll temp=(ll)pow(i,k);
        ans+=(C(n,i)%mod*C(n,m)%mod*temp%mod)%mod;
        ans%=mod;
    }
    printf("%lld",ans);
    return 0;
}
相关推荐
qq_29655327几秒前
【LeetCode】最大子数组乘积:三种解法从暴力到最优
数据结构·算法·leetcode·职场和发展·动态规划·柔性数组
不知名的老吴几秒前
关于C++中的placement new
数据结构·c++·算法
平行侠1 分钟前
023Pollard-ρ 因子分解算法
数据结构·算法
谭欣辰6 分钟前
C++倍增算法详解
数据结构·c++·算法
MATLAB代码顾问6 分钟前
差分进化算法(DE)原理与Python实现
开发语言·python·算法
MicroTech20257 分钟前
微算法科技(NASDAQ :MLGO)基于后量子密码学的动态BFT共识机制:QDBFT架构
科技·算法·密码学
Brilliantwxx15 分钟前
【C++】认识 list(初步认识+模拟实现)
开发语言·数据结构·c++·笔记·算法·list
南宫萧幕18 分钟前
锂电池二阶 RC 模型仿真实战:从理论解析到 Simulink 闭环搭建全流程
开发语言·人工智能·算法·机器学习
故事和你9119 分钟前
洛谷-数据结构2-1-二叉堆与树状数组2
开发语言·javascript·数据结构·算法·ecmascript·动态规划·图论
智者知已应修善业22 分钟前
【51单片机流水灯中断嵌套,低优先级中断完成后如何返回主程序】2023-10-15
c++·经验分享·笔记·算法·51单片机