问题 L: 超超的中等意思
Posted dwvictor
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了问题 L: 超超的中等意思相关的知识,希望对你有一定的参考价值。
问题 L: 超超的中等意思
时间限制: 1 Sec 内存限制: 128 MB提交: 366 解决: 27
[提交] [状态] [命题人:jsu_admin]
题目描述
已知p,q,k和一个难搞得多项式(px+qy)^k。想知道在给定a和b的值下计算多项式展开后x^a*y^b得系数s。
输入
多组输入,每组数据一行输入p,q,k,a,b。其中保证k = a+b,所有输入都为非负数且不大于50000
输出
输出系数s对2^61-1取模后的值
样例输入 Copy
1 1 2 1 1
样例输出 Copy
2
式子我们很容易就可以得知我们要求的是那个系数, 答案也就是求 组合数*x 的系数的 a 次方*y 的系数的 b 次方,
根据这个
我们在计算次方的时候可以用快速幂节约时间,组合数也可以由 C(m+1,n+1)
=C(m,n)+C(m+1,n)递推打表
1 #include<iostream> 2 #include<cmath> 3 //using namespace std; 4 //#include <iostream> 5 #include<string.h> 6 #include <stdio.h> 7 #define ll long long 8 typedef __int128 lll; 9 const lll MOD = 2305843009213693951; 10 using namespace std; 11 const int N = 10000; 12 lll nn[N],mm[N]; 13 14 15 lll dp[50020][50020]; 16 void D() 17 { 18 lll n, k; 19 dp[1][1] = 1; 20 for(int i = 1; i < 50010; i++) 21 { 22 dp[i][0] = 1; 23 } 24 for(int i = 2; i < 50010; i++) 25 { 26 for(int j = 1; j <= i; j++) 27 dp[i][j] = dp[i - 1][j] + dp[i - 1][j - 1]; 28 } 29 //int n, k; 30 //while(scanf("%d %d", &n, &k) == 2) 31 // { 32 //printf("%d ", dp[n][k]); 33 // } 34 return ; 35 // return 0; 36 } 37 lll C(lll n,lll m) 38 { 39 if(m==0||n==m) return 1; 40 lll sb=min(m,n-m); 41 lll f=1,f1; 42 for(lll i=1;i<=sb;i++) 43 { 44 f1=f*(n-i+1)/(i); 45 f=f1; 46 } 47 return f1; 48 } 49 lll pow64(lll x, lll y) { 50 51 if(!y) return 1; 52 if(x >= MOD) x %= MOD; 53 lll ans = 1; 54 while(y) { 55 56 if(y & 1) ans = ans * x % MOD; 57 x = x * x % MOD; 58 y >>= 1; 59 60 } 61 return ans; 62 63 } 64 int main() 65 { 66 int p,q,k,a,b; 67 // D(); 68 while(scanf("%d%d%d%d%d",&p,&q,&k,&a,&b)!=EOF){ 69 70 // long long m,n; 71 // while(cin>>m>>n)///C(m,n) 72 // { 73 // cout<<C[m][n]<<endl; 74 // } 75 //int c = (int)C(m,n); 76 //int l = (int)pow(p,b); 77 //int f = (int)pow(q,n); 78 //prlong longf("%lld",(long long)pow(p,b)); 79 printf("%lld ",ll(pow64(p, a) * pow64(q, b) % MOD * C(k, a) % MOD)); 80 // return 0; 81 } 82 }
以上是关于问题 L: 超超的中等意思的主要内容,如果未能解决你的问题,请参考以下文章
Gym 101490KSafe Racing (dp转换, 超超超级详细,包你看懂)