计算系数 (codevs 1137) 题解

【问题描述】

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

【样例输入】

    1 1 3 1 2

【样例输出】

    3

【解题思路】

     本题为NOIP2011 提高组第一题,看到这玩意我懵了,本人初三狗,完全不会做啊……在老师的指导下终于觉悟,这玩意为二项式定理,其实就是求杨辉三角形的第k行,第n项,在老师推完之后还发现了,其实就是求C(k,m),C为组合数公式,我两种方法都试了,但似乎组合数取模与一半的不一样,于是WA了,顺便求一下哪位神犇看到了请指教一下……这里贴杨辉三角形的代码。求了这个以后,再乘上a^n*b^m即可,注意:a^n与b^m都要边乘边取模,杨辉三角形建立的时候也要取模,最后三者相乘还要取模……

【代码实现】

 1 var f:array[0..1001,0..1001] of longint;
 2     a,b,m,n,k,i,j:longint;
 3 function fac(a,b:int64):int64;
 4 var t:int64;
 5     y:longint;
 6 begin
 7  t:=1;
 8  for y:=1 to b do
 9   t:=(t*a) mod 10007;//幂运算,注意乘一次就要取一次模
10  exit(t);
11 end;
12 begin
13  readln(a,b,k,n,m);
14  for i:=0 to k do
15   f[i,0]:=1;
16  for i:=1 to k do
17   for j:=1 to k do
18    f[i,j]:=(f[i-1,j]+f[i-1,j-1])mod 10007;//建立杨辉三角形,注意取模
19  if k=0 then writeln(1)
20  else
21   writeln((f[k,n]*(fac(a,n)mod 10007)*(fac(b,m)mod 10007))mod 10007);//最后再取模
22 end.
杨辉三角形

 

你可能感兴趣的:(计算系数 (codevs 1137) 题解)