矩阵和POJ 3233

矩阵快速幂的应用，时间复杂度n³logn
假设s_k=I+A+…+A^K-1那么就有，s_k=A^k-1+s_k-1由这个表达我们就得到了通项公式。
矩阵套矩阵便可以做出来。
记录一下模板。
wa来源于开的是2*n的矩阵，但是i，j都是1~n相乘。

#include<iostream>
#include<string>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#include<stack>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<cstdio>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
using namespace std;
typedef pair<int, int> pir;
int n, k, m;
struct mat {int  a[70][70];
};
mat mat_mul(mat x, mat y)
{mat res;memset(res.a, 0, sizeof(res.a));up(i, 0, 2*n){up(j, 0, 2*n){up(k, 0, 2*n){res.a[i][j] = (res.a[i][j] + x.a[i][k] * y.a[k][j])%m;}/*cout << res.a[i][j];if (j+1 != 2*n)cout << " ";else cout << endl;*/}}return res;
}
mat mat_pow(mat x, int i)
{mat res;memset(res.a, 0, sizeof(res.a));up(j, 0, 2*n)res.a[j][j] = 1;while (i > 0){if (i & 1)res = mat_mul(res, x);x = mat_mul(x, x);i >>= 1;}return res;
}
void solve(mat model)
{mat c;mat res;memset(c.a, 0, sizeof(c.a));up(i, 0, n){up(j, 0, n)c.a[i][j] = model.a[i][j];c.a[n + i][i] = 1;c.a[n + i][n + i] = 1;}res = mat_pow(c, k+1);int a;up(i, 0, n){up(j, 0, n){a = res.a[n + i][j];if (i == j)a = (a + m - 1) % m;printf("%d%c", a, (j == n - 1) ? '\n' : ' ');}}
}
int main()
{cin >> n >> k >> m;mat model;up(i, 0, n){up(j,0,n)scanf("%d", &model.a[i][j]);}solve(model);return 0;
}

矩阵和POJ 3233相关推荐

I-Matrix Power Series POJ - 3233 矩阵快速幂+分治
I-Matrix Power Series POJ - 3233 矩阵快速幂+分治 Problem Description Given a n × n matrix A and a positive ...
矩阵儿快速幂 - POJ 3233 矩阵力量系列
不要管上面的标题的bug 那是幂的意思,不是力量... POJ 3233 Matrix Power Series 描述 Given a n × n matrix A and a positive in ...
【矩阵乘法】Matrix Power Series（poj 3233）
Matrix Power Series poj 3233 题目大意给你一个矩阵A,让你求S=A+A2+A3+-+AkS = A + A^2 + A^3 + - + A^kS=A+A2+A3+-+Ak ...
poj 3233 矩阵乘法（分块矩阵）
POJ 3233 题解:Sn为所求矩阵, 则这样, 此题就变成了求矩阵幂和矩阵乘法, 分块矩阵乘法和普通矩阵一样的. code: /* adrui's submission Language : C ...
POJ 3233 Matrix Power Series （矩阵分块，递推）
矩阵乘法是可以分块的,而且幂的和也是具有线性的. 不难得到 Si = Si-1+A*Ai-1,Ai = A*Ai-1.然后矩阵快速幂就可以了. /*************************** ...
POJ 3233 Matrix Power Series 矩阵快速幂 + 二分
题意:求矩阵的次方和解题思路:最容易想到方法就是两次二分因为我们可以把一段 A^1 + A^2 + .......A^K 变成 A^1 + ..A^(K/2) +( A^1 + ..A^( ...
POJ 3233 Matrix Power Serie （矩阵快速幂）
Description Given a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + - + Ak. ...
POJ 3233 矩阵快速幂等比数列
//http://www.cnblogs.com/jiangjing/archive/2013/05/28/3103336.html/* 分析:求a^1+..a^n这是矩阵乘法中关于等比矩阵的求法:| ...
poj 3233 Matrix Power Series
最近一直在做矩阵的题目,这算是一道简单题目吧. 建立一个矩阵的结构体,便于赋值和返回值. S = A + A2 + A3 + - + Ak 如果 S=A+A2 + A3 +A4 S= (A + A ...
poj 3233 Matrix Power Series
Matrix Power Series 思路题意比较简单,就是要求S(n)=∑i=1nAiS(n) = \sum _{i = 1} ^{n} A^ {i}S(n)=∑i=1nAi,显然有S(n)= ...

矩阵和POJ 3233

矩阵和POJ 3233相关推荐

最新文章

热门文章