poj 3233 矩阵乘方
来源:互联网 发布:美工自学 编辑:程序博客网 时间:2024/05/19 00:08
转载于:http://xuyemin520.is-programmer.com/posts/26387.html
// A^n.cpp : 定义控制台应用程序的入口点。////2分法求矩阵冥次//Poj 3233#include "stdafx.h"#include <iostream>#include <algorithm>using namespace std;//如果矩阵过大或者冥过多都会造成栈溢出const int MAXN = 31;int Mod;typedef struct Matrix{int s[MAXN][MAXN];}M;int n;//a+bM Add(M a, M b){M c;for(int i=1; i<=n; i++){for(int j=1; j<=n; j++){c.s[i][j] = (a.s[i][j] + b.s[i][j]) % Mod;}}return c;}//a*bM Multiply(M a, M b){M c;memset(c.s,0,sizeof(c.s));for(int i=1; i<=n; i++){for(int j=1; j<=n; j++){for(int k=1; k<=n; k++){c.s[i][j] = (c.s[i][j] + (a.s[i][k] * b.s[k][j] % Mod))%Mod ;}}}return c;}//a^tM Pow(M a, int t){M b;if(t == 0){memset(b.s,0,sizeof(b.s));for(int i=1; i<=n; i++){b.s[i][i] = 1;}return b;}else{M k = Pow(a,t/2);if(t&1) //t为奇数{return Multiply(Multiply(k,k),a);}else{return Multiply(k,k);}}}//A + A^2 + A^3 + ... + A^tM Sum(M a, int t){M b;if(t == 1){return a;}else{M k = Sum(a,t/2);if(t&1){M p = Pow(a,t/2+1);return Add (Add(k,p) , Multiply(k,p) );}else{return Add(k , Multiply(Pow(a,t/2),k));}}}int main(){M a;int t;cin>>n>>t>>Mod;for(int i=1; i<=n; i++){for(int j=1; j<=n; j++){cin>>a.s[i][j];a.s[i][j] %= Mod;}}a = Sum(a,t);for(int i=1; i<=n; i++){for(int j=1; j<=n; j++){cout<<a.s[i][j];if(j != n)cout<<" ";}cout<<endl;}return 0;}
- poj 3233 矩阵乘方
- 矩阵乘方和(POJ 3233)
- POJ-3233+3070(矩阵加、乘、乘方)()
- 矩阵乘方
- 快速幂+矩阵乘方
- 算法训练 矩阵乘方
- 算法训练 矩阵乘方
- 算法训练 矩阵乘方
- [POJ]1001 乘方计算
- poj 2406:字符串乘方
- 蓝桥杯 算法提高 矩阵乘方
- poj 1001 高精度乘方 java
- poj 1001 Exponentiation 高精度乘方
- HDU 4350--矩阵n乘方的利用!!
- 标量的数组和矩阵乘方
- HDOJ——5171(矩阵乘方)
- POJ 2109-Power of Cryptography(double乘方)
- POJ 3761 Bubble Sort (乘方取模)
- mp4文件格式解析
- 常用内存数据库介绍
- 四种领域模型
- mp4文件格式系列
- 一些人生感悟
- poj 3233 矩阵乘方
- 一个被‘\0'分割成多段的字符串无法全部显示
- Android平台刷机包 修改制作方法
- memcached之java客户端:spymemcached与spring整合
- 写了个进程监控程序,有简易MFC界面
- 标准C语言中的预定义宏(部分)
- POJ 1503 Integer Inquiry 高精度加法
- USACO 1.2 Dual Palindromes
- 从U盘安装Ubuntu的最简易方式