博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
NOIP2011 计算系数
阅读量:5300 次
发布时间:2019-06-14

本文共 1589 字,大约阅读时间需要 5 分钟。

描述

给定一个多项式(ax+by)^k,请求出多项式展开后x^n*y^m项的系数。

输入格式

共一行,包含5 个整数,分别为 a ,b ,k ,n ,m,每两个整数之间用一个空格隔开。 

输出格式

输出共1 行,包含一个整数,表示所求的系数,这个系数可能很大,输出对10007 取模后的结果。 

测试样例1

输入

1 1 3 1 2

输出

3

备注

对于30% 的数据,有 0 ≤k ≤10 ; 
对于50% 的数据,有 a = 1,b = 1; 
对于100%的数据,有 0 ≤k ≤1,000,0≤n, m ≤k ,且n + m = k ,0 ≤a ,b ≤1,000,000。 
#include
#include
#include
#include
#include
using namespace std;const int mod = 10007;int a,b,k,n,m,fac[1005];short c[1005][1005];void getc(){ for(int i = 1;i <= 1000;i++){ c[i][0] = c[i][i] = 1; c[i][1] = c[i][i-1] = i; } for(int i = 2;i <= 1000;i++){ for(int j = 2;j < i;j++){ c[i][j] = (c[i-1][j-1] + c[i-1][j]) % mod; } }}void getfac(){ fac[0] = 1; for(int i = 1;i <= 1000;i++){ fac[i] = (fac[i-1] * i) % mod; }}int q_mul(int a,int b){ int ans = 0; while(b){ if(b&1){ ans = (ans + a) % mod; b--; } b >>= 1; a = (a + a) % mod; } return ans;}int q_pow(int a,int b){ int ans= 1; while(b){ if(b&1){ ans = q_mul(ans,a); } b >>= 1; a = q_mul(a,a); } return ans;}int C(int n,int m){ if(m > n) return 0; else return q_mul(fac[n],q_pow(q_mul(fac[m],fac[n-m]),mod-2));}int lucas(int n,int m){ if(m == 0) return 1; else return q_mul(C(n%mod,m%mod),lucas(n/mod,m/mod));}int main(){ getfac(); cin>>a>>b>>k>>n>>m; cout<

 

转载于:https://www.cnblogs.com/hyfer/p/5846021.html

你可能感兴趣的文章
pandas 修改指定列中所有内容
查看>>
「 Luogu P2285 」打鼹鼠
查看>>
lua语言入门之Sublime Text设置lua的Build System
查看>>
vue.js基础
查看>>
电脑的自带图标的显示
查看>>
[转载] redis 的两种持久化方式及原理
查看>>
C++ 删除字符串的两种实现方式
查看>>
ORA-01502: 索引'P_ABCD.PK_WEB_BASE'或这类索引的分区处于不可用状态
查看>>
Java抽象类和接口的比较
查看>>
php连接postgresql数据库
查看>>
开发进度一
查看>>
MyBaits学习
查看>>
管道,数据共享,进程池
查看>>
CSS
查看>>
[LeetCode] 55. Jump Game_ Medium tag: Dynamic Programming
查看>>
[Cypress] Stub a Post Request for Successful Form Submission with Cypress
查看>>
程序集的混淆及签名
查看>>
thinkphp框架 中 ajax 的应用
查看>>
JAVA排序(一) Comparable接口
查看>>
判断9X9数组是否是数独的java代码
查看>>