POJ 3233 Matrix Power Series 矩阵快速幂+二分…
2018-06-17 23:39:30来源:未知 阅读 ()
矩阵快速幂,请参照模板 http://www.cnblogs.com/pach/p/5978475.html
直接sum=A+A2+A3...+Ak这样累加肯定会超时,但是
sum=A+A2+...+Ak/2+A(k/2)*(A+A2+...+Ak/2) k为偶数时;
sum=A+A2+...+A(k-1)/2+A((k-1)/2)*(A+A2+...+A(k-1)/2)+Ak k为奇数时。
然后递归二分求和
PS:刚开始mat定义的是__int64,于是贡献了n次TLE。。。
#include <iostream> #include <cstring> #include <cstdio> using namespace std; int n,m; const int N=55; struct Mat { int mat[N][N]; }; Mat Multiply(Mat a, Mat b) { Mat c; memset(c.mat, 0, sizeof(c.mat)); for(int k = 0; k < n; ++k) for(int i = 0; i < n; ++i) if(a.mat[i][k]) for(int j = 0; j < n; ++j) if(b.mat[k][j]) c.mat[i][j] = (c.mat[i][j] +a.mat[i][k] * b.mat[k][j])%m; return c; } Mat QuickPower(Mat a, int k) { Mat c; memset(c.mat,0,sizeof(c.mat)); for(int i = 0; i < n; ++i) c.mat[i][i]=1; for(; k; k >>= 1) { if(k&1) c = Multiply(c,a); a = Multiply(a,a); } return c; } Mat Add(Mat a,Mat b) { for(int i=0; i<n; i++) for(int j=0; j<n; j++) a.mat[i][j]=(a.mat[i][j]+b.mat[i][j])%m; return a; } Mat Solve(Mat a,int k) { if(k==1) return a; Mat e,ret; memset(e.mat,0,sizeof(e.mat)); for(int i=0; i<n; i++) e.mat[i][i]=1; ret=Multiply(Add(e,QuickPower(a,k>>1)),Solve(a,k>>1)); if(k%2) return Add(ret,QuickPower(a,k)); return ret; } int main() { //freopen("in.txt","r",stdin); int k; scanf("%d%d%d",&n,&k,&m); Mat a; for(int i=0; i<n; i++) for(int j=0; j<n; j++) scanf("%d",&a.mat[i][j]); Mat ans=Solve(a,k); for(int i=0; i<n; i++) { for(int j=0; j<n-1; j++) printf("%d ",ans.mat[i][j]); printf("%d\n",ans.mat[i][n-1]); } return 0; }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
下一篇:const in C/C++
- POJ-3278 2020-04-01
- Asteroids!_poj2225 2020-02-09
- poj-1753题题解思路 2020-01-26
- POJ1852 2019-11-11
- POJ2431 优先队列+贪心 - biaobiao88 2019-11-03
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash