BZOJ4870:[SHOI2017]组合数问题——题解
2018-02-21 21:21
323 查看
http://www.lydsy.com/JudgeOnline/problem.php?id=4870
https://www.luogu.org/problemnew/show/P3746
看网上一群人说“傻逼题”,我感觉我傻逼了。
首先我们把式子转换一下变成求有nk件物品,我取的物品数%k==r的方案数有多少。
显然f[i][j]=f[i-1][j]+f[i-1][j-1]。
但就没人教一下f[i][j]=f[i-1][j]+f[i-1][j-1]如何矩乘吗……
那我就引洛谷的题解了:
可以加速的原理,其实就是杨辉三角是一个一维递推,并且可以将递推描述为:复制矩阵到一个新矩阵,然后矩阵右移一格,加到新矩阵中。
#include<cmath> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> using namespace std; typedef long long ll; ll n,p,K,r; struct node{ ll g[51][51]; node(){ memset(g,0,sizeof(g)); } friend node operator *(const node &x,const node &y){ node z; for(int i=0;i<K;i++) for(int j=0;j<K;j++) for(int k=0;k<K;k++) z.g[i][k]=(z.g[i][k]+x.g[i][j]*y.g[j][k]%p)%p; return z; } }f,t,res; int main(){ cin>>n>>p>>K>>r; t.g[0][0]=1; for(int i=0;i<K;i++){ f.g[(i-1+K)%K][i]++; f.g[i][i]++; res.g[i][i]=1; } n*=K; while(n){ if(n&1)res=res*f; f=f*f;n>>=1; } printf("%lld\n",(t*res).g[0][r]); return 0; }
+++++++++++++++++++++++++++++++++++++++++++
+本文作者:luyouqi233。 +
+欢迎访问我的博客:http://www.cnblogs.com/luyouqi233/+
+++++++++++++++++++++++++++++++++++++++++++
相关文章推荐
- 【jzoj5215】【BZOJ4870】【Shoi2017】【GDOI2018模拟7.9】【组合数问题】【矩阵快速幂】
- [DP 倍增] BZOJ 4870 [Shoi2017]组合数问题
- 【BZOJ4870】[Shoi2017]组合数问题 动态规划(矩阵乘法)
- BZOJ4870: [Shoi2017]组合数问题
- BZOJ4870: [Shoi2017]组合数问题
- [bzoj4870] [Shoi2017]组合数问题
- 【BZOJ4870】【SHOI2017】组合数问题
- [BZOJ4870][SHOI2017]组合数问题 DP+矩阵快速幂
- [BZOJ4870][Shoi2017]组合数问题 dp+矩阵乘
- bzoj 4870: [Shoi2017]组合数问题
- bzoj4870 [Shoi2017]组合数问题(dp+矩阵倍增)
- bzoj 4870: [Shoi2017]组合数问题 动态规划
- bzoj 4870: [Shoi2017]组合数问题
- BZOJ 4870 [Shoi2017] 组合数问题
- BZOJ 4870 [Shoi2017]组合数问题 ——动态规划 矩阵乘法
- [bzoj4870][Shoi2017]组合数问题
- BZOJ_4870_[Shoi2017]组合数问题_矩阵乘法
- bzoj4870 [Shoi2017]组合数问题
- BZOJ 4870: [Shoi2017]组合数问题 (递推+矩阵快速幂)
- 【bzoj4870】[Shoi2017]组合数问题 dp+快速幂/矩阵乘法