NOIP2011 day2 第一题 计算系数

时间:2023-03-09 15:46:55
NOIP2011 day2 第一题 计算系数

计算系数

NOIP2011 day2 第一题 描述

给定一个多项式(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。

NOIP2011 day2 第一题 计算系数

这里放了一发纯纯的暴力的代码。。。

没有用快速幂,什么什么组合知识。。。

只用到了——->Bling 杨辉三角

(相信大家小学奥数就接触到这个东西了)

(这里面有很多奇奇怪怪的性质我现在都不知道。)

NOIP2011 day2 第一题 计算系数

(摘自百度百科) 里面好多东西看不懂。。。。。

反正就用这个神奇的东西 一层一层递推就好啦,即AC

//By SiriusRen
#include <cstdio>
using namespace std;
int a,b,k,n,m,s[2005],t[2005],ansa=1,ansb=1;
int main(){
scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
a%=10007;b%=10007;
s[1]=s[2]=t[1]=1;
for(int i=2;i<=k;i++){
if(i%2==0)
for(int j=2;j<=i+1;j++)
t[j]=(s[j]+s[j-1])%10007;
else
for(int j=2;j<=i+1;j++)
s[j]=(t[j]+t[j-1])%10007;
}
for(int i=1;i<=n;i++)ansa=(ansa*a)%10007;
for(int i=1;i<=m;i++)ansb=(ansb*b)%10007;
if(k&1)printf("%d\n",(((s[n+1]*ansa)%10007)*ansb)%10007);
else printf("%d\n",(((t[n+1]*ansa)%10007)*ansb)%10007);
}