NOIP 2011 提高组 复赛 day2 factor 计算系数

NOIP 2011 提高组 复赛 day2 factor 计算系数

1.只看标题，就可以猜到，这是一个数学题。

2.准备好稿纸，一看，学过的二项式定理用上了，系数C(k,n)*a^n*b^m。

3.本题，一不小心，数据就会越界，取余必不可少，怎么取余才不会影响结果，看来要好好学习初等数论。

4.按自己理解，乘之前，先将系数全部取余，再进行乘运算，但很遗憾，提交0分，马上按照最后结果才取余，提交20分，这样才稍微平复了心情。

5.翻看书籍《信息学奥赛之数学一本通》，取余发现没有除运算，马上想到，“除”处理过程中有问题，若无整除，保留原数，若整除，值改变为商。信心满满，提交估计50分以上。提交还是20分。

6.想尽办法，没辙。取余遇到除运算，没辙。上网搜索。

7.看了http://blog.sina.com.cn/s/blog_606a23dd010128lo.html这篇文章，发现还是二项式，只是采用杨辉三角的表达形式，递推（或是动态规划DP），有了思路，自个编。跟本人是两种截然不同的思路，虽然都是二项式。

这道题虽然出现在提高组，却并不一定只能用高中知识解决。
其实，这道题可以用递推解决。
设f[i][j]为(ax+by)^i的x^j*y^(n-j)的系数。显然可以得到公式：

   f[i][j]=(f[i-1][j-1]*a+f[i-1][j]*b)%10007。
时间复杂度O（N^2）

8.编写过程中，发现memset(d,0,sizeof(d));有效，memset(d,-1,sizeof(d));有效，memset(d,1,sizeof(d));无效。

9.经过一番编码，提交AC。

10.此题收获，取余一定要学，学习初等数论。

附上AC代码，编译环境Dev-C++4.9.9.2

//2011 factor2 计算系数
//递推公式：d[n][m]=d[n-1][m-1]+d[n-1][m];动态规划 
 
#include <stdio.h>
#include <string.h>
#define mod 10007
int d[1000+10][1000+10];
int main(){
    int a,b,k,n,m;
    int i,j;
    int ans=1;
    scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
    for(i=0;i<=k;i++){//初始化 
        d[i][0]=1;
        d[i][i]=1;
    }
    //i从0开始，j从0开始 
    for(i=2;i<=k;i++)//行
        for(j=1;j<i;j++){//列 
            d[i][j]=(d[i-1][j-1]+d[i-1][j])%mod;
        }
    
    for(i=1;i<=n;i++){
        a%=mod;
        ans*=a;
        ans%=mod;
    }    
    
    for(i=1;i<=m;i++){
        b%=mod;
        ans*=b;
        ans%=mod;
    }
        
    ans*=d[k][n];
    ans%=mod;
    
    printf("%d\n",ans);
    return 0;
} 

附上20分代码，编译环境Dev-C++4.9.9.2

//2011 factor
#include <stdio.h>
#define mod 10007
int main(){
    int a,b,k,n,m;
    long long ans;
    int i,j;
    int result;
    scanf("%d%d%d%d",&a,&b,&k,&n,&m);
    ans=1;
    for(i=1;i<=n;i++){
        //a%=mod;
        ans*=a;
        //ans%=mod;
    }
    for(i=1;i<=m;i++){
        //b%=mod;
        ans*=b;
        //ans%=mod;
    }
    for(i=0;i<n;i++){
        ans*=(k-i);
        //ans%=mod;
    }
    for(i=n;i>0;i--){
        ans/=i;//此处写成了ans/=n，查了会 
    }
    ans%=mod;
    result=ans;
    printf("%d\n",result);
    return 0;
}