poj 3233 Matrix Power Series 矩阵快速幂or二分
来源:互联网 发布:10.1.1越狱抹掉数据 编辑:程序博客网 时间:2024/06/06 15:42
题目:
Description
Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.
Input
The input contains exactly one test case. The first line of input contains three positive integers n (n ≤ 30), k (k ≤ 109) and m (m < 104). Then follow n lines each containing n nonnegative integers below 32,768, giving A’s elements in row-major order.
Output
Output the elements of S modulo m in the same way as A is given.
Sample Input
2 2 40 11 1
Sample Output
1 22 3
分析:
现有题解主要是两种做法 构造矩阵+矩阵快速幂 或者 二分加速
第二种方法不难理解,式子化简一下就可以了
第一种方法挑战中直接给了一个矩阵,并没有细讲怎么想到构造这个矩阵的
http://blog.csdn.net/wangjian8006/article/details/7868864 怎么想到用这个矩阵可以参考这篇博客
代码:
方法一代码参挑战
方法二:二分
/*二分可过 等比数列和的二分加速k为偶数:sum(k) = (1+A^(k/2)) *( A+A^2+……+A^(k/2)) = (1+A^(k/2)) * sum(k/2)k为奇数:sum(k) = (1+A^(k/2)) * sum(k/2) + A^k*/#include<iostream>#include<cstring>#define size 35using namespace std;typedef long long ll;int n,k,mod;struct Mat{ int mat[size][size]; void clear(){ memset(mat,0,sizeof(mat)); } Mat operator *(const Mat &e) const{ Mat tmp; tmp.clear(); for(int k=0;k<size;k++) for(int i=0;i<size;i++){ if(mat[i][k]==0) continue; for(int j=0;j<size;j++){ if(e.mat[k][j]==0) continue; tmp.mat[i][j]+=mat[i][k]*e.mat[k][j]%mod; tmp.mat[i][j]%=mod; } } return tmp; } Mat operator +(const Mat &e) const{ Mat tmp; tmp.clear(); for(int i=0;i<size;i++){ for(int j=0;j<size;j++){ tmp.mat[i][j]=(mat[i][j]%mod+e.mat[i][j]%mod)%mod; } } return tmp; }};Mat m,E;Mat pow(Mat ma,ll num){ Mat ans; ans.clear(); for(int i=0;i<n;i++) ans.mat[i][i]=1; while(num){ if(num&1) ans=ans*ma; num/=2; ma=ma*ma; } return ans;}Mat Sum(int x){ if(x==1) return m; if(x&1) return Sum(x-1)+pow(m,x); else return (pow(m,x/2)+E)*Sum(x/2);}int main() { ios::sync_with_stdio(false); cin.tie(0); cin>>n>>k>>mod; E.clear(); for(int i=0;i<n;i++){ E.mat[i][i]=1; for(int j=0;j<n;j++) cin>>m.mat[i][j]; } Mat ans=Sum(k); for(int i=0;i<n;i++){ for(int j=0;j<n;j++) cout<<ans.mat[i][j]<<" "; cout<<endl; } return 0;}
阅读全文
0 0
- poj 3233 Matrix Power Series 矩阵快速幂or二分
- poj 3233 Matrix Power Series 矩阵 快速幂 两次二分
- POJ 3233 Matrix Power Series(矩阵快速幂+二分求和)
- POJ 3233 Matrix Power Series 矩阵快速幂+二分
- POJ 3233 Matrix Power Series 矩阵快速幂+二分
- POJ 3233 Matrix Power Series (矩阵快速幂 + 二分)
- POJ 3233 Matrix Power Series(矩阵快速幂+二分)
- poj 3233 Matrix Power Series(矩阵二分,快速幂)
- POJ 3233 Matrix Power Series (矩阵快速幂+二分)
- POJ - Matrix Power Series 【矩阵快速幂+二分求和】
- POJ Matrix Power Series (矩阵快速幂+二分)
- poj 1995 Matrix Power Series 二分+矩阵快速幂
- 【矩阵快速幂+二分】Matrix Power Series POJ
- poj 3233 --- Matrix Power Series (二分,矩阵)
- POJ 3233 Matrix Power Series (矩阵 + 二分)
- POJ 3233 Matrix Power Series(矩阵+二分)
- POJ 3233-Matrix Power Series(矩阵快速幂+二分求矩阵和)
- POJ 3233 Matrix Power Series(求矩阵幂的和——分块矩阵快速幂 or 二分递归+矩阵快速幂)
- 3.MariaDB笔记——MariaDB在Windows安装
- Android_AutoRoom 1
- String与StringBuffer的区别
- 第4周项目6 多项式求和
- 4.MariaDB笔记——MariaDB在Windows源码编译
- poj 3233 Matrix Power Series 矩阵快速幂or二分
- 5.MariaDB笔记——MariaDB开发者
- linux c程序中内核态与用户态内存存储问题
- CentOS7搭建FTP服务器
- 实验吧 这个FLAG有点怪 By Assassin
- spark安装教程
- leetcode.array--48. Rotate Image
- Android自定义控件(一) 自定义组合控件
- 6.MariaDB笔记——cmake使用介绍一