noip 2011 提高组 计算系数
来源:互联网 发布:python sha1加密 编辑:程序博客网 时间:2024/06/10 12:46
题目描述:
给定一个多项式(ax + by)^k,请求出多项式展开后(x^n)*(y^m)项的系数。(保证n+m==k)
输入:
共一行,包含 5 个整数,分别为a,b,k,n,m,每两个整数之间用一个空格隔开。
输出:
输出共 1 行,包含一个整数,表示所求的系数,这个系数可能很大,输出对10007 取模后的结果。
输入输出样例:
输入:1 1 3 1 2输出:3
题解:
我们来看几组简单数据:(a*x+b*y)^2=(a*x)^2+2*a*b*x*y+(b*y)^2(a*x+b*y)^3=(a*x)^3+3*(a^2)*b*(x^2)*y+3*a*(b^2)*x*(y^2)+(b*y)^3(a*x+b*y)^4=(a*x)^4+4*(a^3)*b*(x^3)*y+6*(a^2)*(b^2)*(x^2)*(y^2)+4*a*(b^3)*x*(y^3)+(b*y)^3 易得(x^n)*(y^m)的系数(答案)为k*(a^n)*(b^m),而k值有以下规律: 1 2 1 1 3 3 1 1 4 6 4 1 ...... 由此可见,这是一个杨辉三角(dalao都用组合数), 因为(a^n)*(b^m)是可以算的(直接循环,边乘边取余或者快速幂计算) 因此,在构造的杨辉三角中找到对应系数 k,再计算出(a^n)*(b^m),即可得到答案。
代码:
#include<stdio.h>long long h[1010][1010],k,a,b,n,m;int Tpow(int x,int p){ //快速幂计算(a^n)*(b^m) long long tmp=x,ans=1; while(p>0){ if(p&1)ans=(ans*tmp)%10007; tmp=(tmp*tmp)%10007; p>>=1; } return ans;}int main(){ scanf("%lld%lld%lld%lld%lld",&a,&b,&k,&n,&m); h[1][1]=1; for(int i=2;i<=k+1;i++) //暴力构造杨辉三角,别学我 for(int j=1;j<=i;j++) h[i][j]=(h[i-1][j-1+h[i-1][j])%10007; long long s=(h[k+1][n+1]%10007*(Tpow(a,n)%10007)*(Tpow(b,m)%10007))%10007; printf("%lld",s);}
阅读全文
1 0
- noip 2011 提高组 计算系数
- NOIP 2011 提高组 复赛 day2 factor 计算系数
- NOIP 2011 计算系数
- Noip 2011 计算系数
- Codevs 1137 计算系数 2011年NOIP全国联赛提高组
- NOIP 2011 Day2 T1 计算系数
- noip2011提高组——计算系数
- rqnoj[NOIP2011提高组]计算系数
- 【递推】【NOIP2011提高组】计算系数
- 计算系数(NOIP2011提高组Day2T1)
- [NOIP 2011]计算系数 乘法逆元+扩展欧几里得
- codevs 1137||NOIP 2011 计算系数 二项式定理 解题报告
- 洛谷 1313||NOIP 2011 计算系数 二项式定理 解题报告
- NOIP提高组 旷野大计算
- luogu1313【2011提高】计算系数(求组合数取模)
- NOIP2011复赛提高组day2(A:计算系数 B:聪明的质监员 C:观光公交)
- 2011年NOIP提高组 铺地毯
- NOIP提高组2011 观光公车
- ACM集训日记-8月15日
- git log 查找文件变更
- FFMPEG之编译生成
- 使用TestNG-xslt+ant来美化selenium测试报告
- caffe的两个测试案例
- noip 2011 提高组 计算系数
- 你告诉我,读书没有用
- Python 的 mysql 操作
- 通过原型的方式实现遍历
- Proxy 代理器(1)
- HDU 3549 Flow Problem (Ford-Fulkerson&Dinic)
- 数塔问题
- Erlang搭建cowboy服务器
- go sockt 编程(初识)