计算系数 （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.

杨辉三角形