矩阵快速幂 题目整理 HDU5015、CF392 C、HDU4686、CF385 E

最新推荐文章于 2021-07-16 10:04:10 发布
最新推荐文章于 2021-07-16 10:04:10 发布
阅读量186 收藏
点赞数
分类专栏： 数论 文章标签： 矩阵快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qcwlmqy/article/details/97822252
版权

矩阵快速幂

例行公事先介绍下矩阵乘法~~（相信大家都会）~~

矩阵

  • 在数学中，矩阵（Matrix）是一个按照长方阵列排列的复数或实数集合。

矩阵乘法

  • 矩阵相乘只有在第一个矩阵的列数（column）和第二个矩阵的行数（row）相同时才有意义
  • 一个m×n的矩阵就是m×n个数排成m行n列的一个数阵
  • 由于它把许多数据紧凑的集中到了一起，所以有时候可以简便地表示一些复杂的模型，如电力系统网络模型

设A为 m ∗ p m*p m∗p的矩阵，B为 p ∗ n p*n p∗n的矩阵，那么称 m ∗ n m*n m∗n的矩阵C为矩阵A与B的乘积，记作 C = A B C = AB C=AB
其中矩阵C中的第 i 行第 j 列元素可以表示为：
C [ i ] [ j ] = ∑ k = 1 p a [ i ] [ k ] ∗ b [ k ] [ j ] = a [ i ] [ 1 ] ∗ b [ 1 ] [ j ] + a [ i ] [ 2 ] ∗ b [ 2 ] [ j ] + ⋯ + a [ i ] [ p ] ∗ b [ p ] [ j ] C[i][j]=\sum_{k=1}^{p}{a[i][k]*b[k][j]}=a[i][1]*b[1][j]+a[i][2]*b[2][j]+\cdots+a[i][p]*b[p][j] C[i][j]=k=1∑p​a[i][k]∗b[k][j]=a[i][1]∗b[1][j]+a[i][2]∗b[2][j]+⋯+a[i][p]∗b[p][j]

  • 对于矩阵 C = A B C = AB C=AB ，矩阵A的第m行“乘上”矩阵B的第n列

精简板子

struct Matrix {
	LL n, m;
	LL v[maxn][maxn];
	Matrix(int n = 0, int m = 0) :n(n), m(m) {	//构造函数
		memset(v, 0, sizeof(v));
	}
	void clear() {	//重置函数
		memset(v, 0, sizeof(v));
	}
	void show() {	//debug函数
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= m; j++)
				cout << v[i][j] << ' ';
			cout << '\n';
		}
	}
	friend Matrix operator *(const Matrix& a, const Matrix& b) {//乘法封装
		Matrix res(a.n, b.m);
		for (int i = 1; i <= a.n; i++)
			for (int j = 1; j <= b.m; j++)
				for (int k = 1; k <= a.m; k++)
					res.v[i][j] = (a.v[i][k] * b.v[k][j] % mod + res.v[i][j]) % mod;
		return res;
	}
};

矩阵快速幂

  • 所谓矩阵快速幂，即为使用快速幂的计算方式去计算相同矩阵连乘（即矩阵的幂）的结果
  • 前面已经讲过，可以将矩阵构造成结构体，重载运算符后就可以直接当做正常变量来使用了
void quickpow(Matrix& res, Matrix m, LL p) {
	while (p) {
		if (p & 1)
			res = res * m;
		m = m * m;
		p >>= 1;
	}
}

重点

以下为矩阵快速幂，常见的题目整理和公式

Fib数

( F i b ( 1 ) F i b ( 2 ) ) ∗ ( 0 1 1 1 ) \left(\begin{array}{ccc}Fib(1)&amp;Fib(2)\\\end{array}\right)*\left(\begin{array}{ccc}0&amp;1\\1&amp; 1\\\end{array}\right) (Fib(1)​Fib(2)​)∗(01​11​) = = = ( F i b ( 2 ) F i b ( 1 ) + F i b ( 2 ) ) \left(\begin{array}{ccc} Fib(2)&amp;Fib(1)+Fib(2)\end{array}\right) (Fib(2)​Fib(1)+Fib(2)​)

P3390 【模板】矩阵快速幂模板题

常用推导

在这里插入图片描述

HDU 5015 233 Matrix

题意

给出矩阵的第0行 ( 233 , 2333 , 23333 , ⋯ &ThinSpace; ) (233,2333,23333,\cdots) (233,2333,23333,⋯)和第0列 a 1 , a 2 , ⋯ &ThinSpace; , a n a1,a2,\cdots,an a1,a2,⋯,an ( n &lt; = 10 ， m &lt; = 1 0 9 ) (n&lt;=10，m&lt;=10^9) (n<=10，m<=109)，给出式子： A [ i ] [ j ] = A [ i − 1 ] [ j ] + A [ i ] [ j − 1 ] A[i][j] = A[i-1][j] + A[i][j-1] A[i][j]=A[i−1][j]+A[i][j−1]，要求 A [ n ] [ m ] A[n][m] A[n][m]

思路

较为复杂的一个推导

  • 可以假设未给出的第0行第0列值为23
  • 对于每一列，可以将所有元素一起构造成初始矩阵。
    在这里插入图片描述

CF392C Yet Another Number Sequence

题意

给定 n , k , F i b ( 1 ) = 1 , F i b ( 2 ) = 2 n,k,Fib(1)=1,Fib(2)=2 n,k,Fib(1)=1,Fib(2)=2,令 A k ( i ) = F i b ( i ) ∗ i k A_{k}(i)=Fib(i)*i^k Ak​(i)=Fib(i)∗ik
求 ∑ i = 1 n A k ( i ) \sum_{i=1}^{n}A_{k}(i) ∑i=1n​Ak​(i) ( k ≤ 40 ) (k\leq40) (k≤40)

思路

  • 本题的难点在于 i k i^k ik的推导
  • 以下表格展示使用二项式定理推导 i 3 i^3 i3的过程
1 0 1^0 10 1 1 1^1 11 1 2 1^2 12
C [ 0 ] [ 0 ] ∗ 1 0 = 2 0 C[0][0]*1^0=2^0 C[0][0]∗10=20 C [ 1 ] [ 0 ] ∗ 1 0 + C [ 1 ] [ 1 ] ∗ 1 1 = 2 1 C[1][0]*1^0+C[1][1]*1^1=2^1 C[1][0]∗10+C[1][1]∗11=21 C [ 2 ] [ 0 ] ∗ 1 0 + C [ 2 ] [ 1 ] ∗ 1 1 + C [ 2 ] [ 2 ] ∗ 1 2 = 2 2 C[2][0]*1^0+C[2][1]*1^1+C[2][2]*1^2=2^2 C[2][0]∗10+C[2][1]∗11+C[2][2]∗12=22
C [ 0 ] [ 0 ] ∗ 2 0 = 3 0 C[0][0]*2^0=3^0 C[0][0]∗20=30 C [ 1 ] [ 0 ] ∗ 2 0 + C [ 1 ] [ 1 ] ∗ 2 1 = 3 1 C[1][0]*2^0+C[1][1]*2^1=3^1 C[1][0]∗20+C[1][1]∗21=31 C [ 2 ] [ 0 ] ∗ 2 0 + C [ 2 ] [ 1 ] ∗ 2 1 + C [ 2 ] [ 2 ] ∗ 2 2 = 3 2 C[2][0]*2^0+C[2][1]*2^1+C[2][2]*2^2=3^2 C[2][0]∗20+C[2][1]∗21+C[2][2]∗22=32

HDU4686 Arc of Dream

题意

a i = a i − 1 ∗ a x + a y a_{i}=a_{i-1}*ax+ay ai​=ai−1​∗ax+ay
b i = b i − 1 ∗ b x + b y b_{i}=b_{i-1}*bx+by bi​=bi−1​∗bx+by
求 ∑ i = 0 i &lt; n a i ∗ b i \sum_{i=0}^{i &lt; n} a_{i}*b_{i} ∑i=0i<n​ai​∗bi​

思路

本题难度在于 a i ∗ b i a_{i}*b_{i} ai​∗bi​的推导
显然由 a i a_{i} ai​和 b i b_{i} bi​使用矩阵快速幂推导比较困难
但是由 a i ∗ b i , a i , b i a_{i}*b_{i},a_{i},b_{i} ai​∗bi​,ai​,bi​推导就可以了
( a i ∗ a x + a y ) ∗ ( b i ∗ b x + b y ) = a x ∗ b x ∗ a i ∗ b i + b x ∗ a y ∗ b i + a x ∗ b y ∗ a i + a y ∗ b y (a_{i}*ax+ay)*(b_{i}*bx+by)=ax*bx*a_{i}*b_{i}+bx*ay*b_{i}+ax*by*a_{i}+ay*by (ai​∗ax+ay)∗(bi​∗bx+by)=ax∗bx∗ai​∗bi​+bx∗ay∗bi​+ax∗by∗ai​+ay∗by
以 ( a i ∗ b i , a i , a y , b i , b y , s u m ) (a_{i}*b_{i},a_{i},ay,b_{i},by,sum) (ai​∗bi​,ai​,ay,bi​,by,sum)为初始矩阵推导即可

CF385 E Bear in the Field

题意

有一片 n × n n×n n×n的草莓地，每个位置的初始草莓量为横坐标和纵坐标的和，然后每过一秒增长一个草莓
然后给出熊的初始位置 ( s x , s y ) (sx,sy) (sx,sy)，以及移动的速度 ( d x , d y ) (dx,dy) (dx,dy)，每一秒发生的事：

  • 速度增加k（k为该位置的草莓数）
  • 熊的位置发生移动
  • 每个位置上草莓数+1
    求t秒后熊的位置

思路

为了方便求余n，将坐标由 [ 1 , n ] [1,n] [1,n]变成 [ 0 , n − 1 ] [0,n-1] [0,n−1]
在这里插入图片描述

233 Matrix 代码

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace  std;
#pragma warning (disable:4996)
typedef long long LL;
const int maxn = 13;
const LL mod = 10000007;
struct Matrix
{
	int n, m;
	int v[maxn][maxn];
	Matrix(int n = 0, int m = 0):n(n),m(m) {
		memset(v, 0, sizeof(v));
	}
	void clear() {
		memset(v, 0, sizeof(v));
	}
	void show() {
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= m; j++)
				cout << v[i][j] << ' ';
			cout << '\n';
		}
	}
	friend Matrix operator *(const Matrix& a, const Matrix& b) {
		Matrix res(a.n, b.m);
		for (int i = 1; i <= a.n; i++)
			for (int j = 1; j <= b.m; j++)
				for (int k = 1; k <= a.m; k++)
					res.v[i][j] = (res.v[i][j] + LL(a.v[i][k]) * b.v[k][j] % mod) % mod;
		return res;
	}
};
void quickpow(Matrix &res, Matrix m, LL p) {
	while (p) {
		if (p & 1)
			res = res * m;
		m = m * m;
		p >>= 1;
	}
}
int main() {
	int n, m;
	while (~scanf("%d%d", &n, &m)) {
		Matrix a(1, n + 2), b(n + 2, n + 2);
		a.v[1][1] = 23;
		for (int i = 2; i <= n + 1; i++)
			scanf("%d", &a.v[1][i]);
		a.v[1][n + 2] = 3;
		for (int i = 1; i <= n + 2; i++)
			b.v[1][i] = 10, b.v[n + 2][i] = 1;
		b.v[1][n + 2] = 0;
		for (int i = 2; i <= n + 1; i++)
			for (int j = i; j <= n + 1; j++)
				b.v[i][j] = 1;
		quickpow(a, b, m);
		printf("%d\n", a.v[1][n + 1]);
	}
}

Yet Another Number Sequence 代码

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace  std;
#pragma warning (disable:4996)
typedef long long LL;
const int maxn = 85;
const LL mod = 1e9 + 7;
LL n; int k;
struct Matrix
{
	LL n, m;
	LL v[maxn][maxn];
	Matrix(int n = 0, int m = 0) :n(n), m(m) {
		memset(v, 0, sizeof(v));
	}
	void clear() {
		memset(v, 0, sizeof(v));
	}
	void show() {
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= m; j++)
				cout << v[i][j] << ' ';
			cout << '\n';
		}
	}
	friend Matrix operator *(const Matrix& a, const Matrix& b) {
		Matrix res(a.n, b.m);
		for (int i = 1; i <= a.n; i++)
			for (int j = 1; j <= b.m; j++)
				for (int k = 1; k <= a.m; k++)
					res.v[i][j] = (res.v[i][j] + a.v[i][k] * b.v[k][j] % mod) % mod;
		return res;
	}
};
int comb[41][41];
void Pre_Work() {
	for (int i = 0; i <= 40; i++)
		comb[i][i] = comb[i][0] = 1;
	for (int i = 2; i <= 40; i++) {
		for (int j = 1; j < i; j++)
			comb[i][j] = (LL(comb[i - 1][j]) + comb[i - 1][j - 1]) % mod;
	}
}
Matrix a;
void Matrix_a() {
	a.n = 1; a.m = 2 * k + 3;
	for (int i = 1; i <= k + 1; i++)
		a.v[1][i] = (LL(1) << (i - 1)) % mod;
	for (int i = k + 2; i <= 2 * k + 2; i++)
		a.v[1][i] = (LL(2) << (i - k - 2)) % mod;
	a.v[1][2 * k + 3] = (LL(1) << k + 1) + 1;
	//a.show();
}
Matrix b;
void Matrix_b() {
	b.n = 2 * k + 3;b.m = 2 * k + 3;
	for (int i = 1; i <= k + 1; i++)
		for (int j = 1; j <= i; j++)
			b.v[j + k + 1][i] = comb[i - 1][j - 1];
	for (int i = k + 2; i <= 2 * k + 2; i++) {
		for (int j = 1; j <= i - k - 1; j++) {
			b.v[j][i] = comb[i - k - 2][j - 1];
			b.v[j + k + 1][i] = comb[i - k - 2][j - 1];
		}
	}
	for (int i = 1; i <= k + 1; i++)
		b.v[i][2 * k + 3] = b.v[i + k + 1][2 * k + 3] = comb[k][i - 1];
	b.v[2 * k + 3][2 * k + 3] = 1;
	//b.show();
}
void quickpow(Matrix& res, Matrix m, LL p) {
	while (p) {
		if (p & 1)
			res = res * m;
		m = m * m;
		p >>= 1;
	}
}
int main() {
	scanf("%I64d%d", &n, &k);
	if (n == 1) {
		printf("1\n");
		return 0;
	}
	Pre_Work();
	Matrix_a();
	Matrix_b();
	quickpow(a, b, n - 2);
	printf("%I64d\n", a.v[1][2 * k + 3]);
}

C - Arc of Dream 代码

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#pragma warning (disable:4996)
typedef long long LL;
const LL mod = 1e9 + 7;
const int maxn = 7;
struct Matrix
{
	LL n, m;
	LL v[maxn][maxn];
	Matrix(int n = 0, int m = 0) :n(n), m(m) {
		memset(v, 0, sizeof(v));
	}
	void clear() {
		memset(v, 0, sizeof(v));
	}
	void show() {
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= m; j++)
				cout << v[i][j] << ' ';
			cout << '\n';
		}
	}
	friend Matrix operator *(const Matrix& a, const Matrix& b) {
		Matrix res(a.n, b.m);
		for (int i = 1; i <= a.n; i++)
			for (int j = 1; j <= b.m; j++)
				for (int k = 1; k <= a.m; k++)
					res.v[i][j] = (a.v[i][k] * b.v[k][j] % mod + res.v[i][j])%mod;
		return res;
	}
};
LL a0, ax, ay, b0, bx, by;
Matrix a;
void Matrix_a() {
	a.clear();
	a.n = 1; a.m = 6;
	a.v[1][1] = a0 * b0; a.v[1][1] %= mod;
	a.v[1][2] = a0;
	a.v[1][3] = ay;
	a.v[1][4] = b0;
	a.v[1][5] = by;
	a.v[1][6] = 0;
	//a.show();
}
Matrix b;
void Matrix_b() {
	b.clear();
	b.n = 6; b.m = 6;
	b.v[1][1] = ax * bx; b.v[1][1] %= mod;
	b.v[2][1] = ax * by; b.v[2][1] %= mod;
	b.v[3][1] = by;
	b.v[4][1] = bx * ay; b.v[4][1] %= mod;
	b.v[2][2] = ax;
	b.v[3][2] = 1;
	b.v[3][3] = 1;
	b.v[4][4] = bx;
	b.v[5][4] = 1;
	b.v[5][5] = 1;
	b.v[1][6] = b.v[6][6] = 1;
	//b.show();
}
void quickpow(Matrix& res, Matrix m, LL p) {
	while (p) {
		if (p & 1)
			res = res * m;
		m = m * m;
		p >>= 1;
	}
}
int main() {
	LL n;
	while (~scanf("%lld", &n)) {
		scanf("%lld%lld%lld", &a0, &ax, &ay);
		a0 %= mod; ax %= mod; ay %= mod;
		scanf("%lld%lld%lld", &b0, &bx, &by);
		b0 %= mod; bx %= mod; by %= mod;
		Matrix_a();
		Matrix_b();
		quickpow(a, b, n);
		printf("%lld\n", a.v[1][6]);
	}
}

Bear in the Field 代码

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#pragma warning (disable:4996)
typedef long long LL;
const int maxn = 7;
LL mod;
struct Matrix
{
	LL n, m;
	LL v[maxn][maxn];
	Matrix(int n = 0, int m = 0) :n(n), m(m) {
		memset(v, 0, sizeof(v));
	}
	void clear() {
		memset(v, 0, sizeof(v));
	}
	void show() {
		for (int i = 1; i <= n; i++) {
			for (int j = 1; j <= m; j++)
				cout << v[i][j] << ' ';
			cout << '\n';
		}
	}
	friend Matrix operator *(const Matrix& a, const Matrix& b) {
		Matrix res(a.n, b.m);
		for (int i = 1; i <= a.n; i++)
			for (int j = 1; j <= b.m; j++)
				for (int k = 1; k <= a.m; k++)
					res.v[i][j] = (a.v[i][k] * b.v[k][j] % mod + res.v[i][j]) % mod;
		return res;
	}
};
LL sx, sy, dx, dy, n;
Matrix a;
void Matrix_a() {
	a.clear();
	a.n = 1; a.m = 6;
	a.v[1][1] = sx;
	a.v[1][2] = sy;
	a.v[1][3] = dx;
	a.v[1][4] = dy;
	a.v[1][5] = 0;
	a.v[1][6] = 1;
	//a.show();
}
Matrix b;
void Matrix_b() {
	b.clear();
	b.n = 6; b.m = 6;
	b.v[1][1] = b.v[6][1] = 2;
	b.v[2][1] = b.v[3][1] = b.v[5][1] = 1;
	b.v[2][2] = b.v[6][2] = 2;
	b.v[1][2] = b.v[4][2] = b.v[5][2] = 1;
	b.v[6][3] = 2;
	b.v[1][3] = b.v[2][3] = b.v[3][3] = b.v[5][3] = 1;
	b.v[6][4] = 2;
	b.v[1][4] = b.v[2][4] = b.v[4][4] = b.v[5][4] = 1;
	b.v[5][5] = b.v[6][5] = 1;
	b.v[6][6] = 1;
	//b.show();
}
void quickpow(Matrix& res, Matrix& m, LL p) {
	while (p) {
		if (p & 1)
			res = res * m;
		m = m * m;
		p >>= 1;
	}
}
int main() {
	scanf("%I64d", &mod);
	scanf("%I64d%I64d%I64d%I64d%I64d", &sx, &sy, &dx, &dy, &n);
	sx--; sy--;
	Matrix_a();
	Matrix_b();
	quickpow(a, b, n);
	a.v[1][1]++; a.v[1][2]++;
	if (a.v[1][1] <= 0)a.v[1][1] += mod;
	if (a.v[1][2] <= 0)a.v[1][2] += mod;
	printf("%I64d %I64d\n", a.v[1][1], a.v[1][2]);
}
qcwlmqy
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
矩阵快速幂专题（一）
shengtao96的专栏
03-25 2619
最近闲来无事，准备集中精力刷一波数论与图论。矩阵快速幂是数论里面的重要组成部分，值得我好好学习一下。因为题目比较多，分析也比较多，所以将此专题分成几个部分。做完这一专题，可能会暂时转向图论部分，然后等我组合数学学得差不多了，再回过头来继续做数论题。 矩阵快速幂算法的核心思想是将问题建模转化为数学模型（有一些简单题目是裸的矩阵模型，但是大部分难题就是难在要构造矩阵，用矩阵方法解决问题），推倒递推式
【专题】—【数学】—【矩阵快速幂】
ACer-TomatoFish
08-16 4651
参考文章： http://blog.csdn.net/shiwei408/article/details/8818386
矩阵快速幂专题【完结】
caobihole
08-25 397
第一题 hdu 1757 A Simple Math Problem 点击打开链接 思路：矩阵快速幂 分析: 1 最简单的矩阵快速幂的题目，直接利用矩阵求解即可 点击打开查看代码 第二题 hdu 1575Tr A 点击打开hdu 1575 思路: 矩阵快速幂 分析: 1 题目给定一个n*n的矩阵要求矩阵的k次幂之后的矩阵的对角线的和 2 矩阵快速幂的裸题 点...
矩阵快速幂题集
wangtao的博客
07-10 632
矩阵快速幂 矩阵必须全部初始化,注意矩阵大小 POJ - 3070 最简单的矩阵快速幂 f[n]=f[n-1]+f[n-2] ac: #include <iostream> #include <cstdio> #include <cstring> #define ll long long using namespace std; /* (f(n-...
矩阵快速幂的一些题目
ypeijasd的专栏
07-16 625
矩阵快速幂求斐波那契数列第n项，递推公式 进而转换为矩阵幂次 矩阵快速幂代码 //矩阵快速幂计算 a^n #include <bits/stdc++.h> using namespace std; int mod=1e9+7; //模数 long long n; struct matrix //定义结构体来表示矩阵 { long long m[3][3]; }; matrix a, an, fib; //fib是最终斐波那契矩阵, a是矩阵, an为计算完a
leetcode和hdoj-Algorithm-Packages:算法包
07-06
用于存放参与和vp的cf比赛，比赛文件夹以比赛序号和div描述，平时练习的题目放在problems下 hdoj 存放hdoj题目代码 hdu-winter-2020 存放2020寒假hdu相关的训练，problems下存放集训的专题，contests下存放寒假比赛...
OJOJO:ACMer 的 Alfred 2 工作流程
07-16
然后，您可以使用这些命令： OJOJOCLEAR清除缓存TC推出 TopCoder Arena OJNAME ProblemID ProblemID 在 OJNAME 中显示问题 ProblemID 例子： CF 123A显示 Codeforces 中的问题 123A HDU 4567在 HDU OJ 中显示问题 ...
万能头文件stdc++.h
04-26
HDU部分支持（G++支持，C++不支持）。 其他国外的oj，还有台湾的oj都支持，CF，Topcoder也都支持。 当然，其实这是一个偷懒的写法，但是会降低编译速度（为何会降低编译速度，我还不能知道，等到之后学编译原理再...
矩阵快速幂（快速幂）模板题目Decoding Genome
名字不能相同
05-22 1004
这类题目的难点在于，数据量大，需要用到高精度，也就是快速幂 https://codeforces.com/group/vFwRVj9WjO/contest/328371/problem/H #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <cmath> #include <algorithm> #define M
矩阵快速幂（模板+例题）
lr7682的博客
05-15 645
矩阵快速幂推导过程：https://blog.csdn.net/u012061345/article/details/52224623 矩阵快速幂求解数列第n项的关键在于计算系数矩阵A。之后就是套模板了。 模板（求解斐波那契数列第n项）： #include <cstdio> #include <cstring> #define LL long long #define mo...
hdu1575 矩阵快速幂 模板题
LATTE
10-13 391
传送门 题意：求给定n维方阵的k次幂后的对角线元素和 思路：模板题 代码如下： #include &lt;bits/stdc++.h&gt; int N;//N个系数，N维矩阵 typedef long long ll; using namespace std; struct matrix { int m[20][20]; }; matrix ans, base, m; matr...
一类常见的矩阵快速幂题型 总结
herongwei 的 BLOG
09-21 1969
Googlecode jam 2008 Round 1A(c.Numbers) 【题意】计算的小数点前三位数，不足三位补0,正整数n的最大值为20亿。 【前提】：满足 的值在【0，1】范围 首先将 展开之后可以发现 的形式，同样的，有 因此， 是个整数， 其中 这正是解题的关键！ 由于 所以的整
hdoj 1575 Tr A (矩阵快速幂)
R
11-13 1855
Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 4549    Accepted Submission(s): 3428 Problem Description A为一个方阵，则Tr A表示A的迹（就是主
快速幂+矩阵快速幂（总结+例题）
热门推荐
red_red_red的博客
05-14 1万+
1.快速幂 以求a的b次方来介绍： 首先把b转换成二进制数 该二进制数第i位的权为 2^i - 1 . 比如 : 11的二进制是1011 11 = 2³×1 + 2²×0 + 2¹×1 + 2º×1 所以假设我们要求a^b，那么其实b是可以拆成二进制的，该二进制数第i位的权为2^(i-1)，例如当b==11时 　　　　　　　　　　a^11=a^(2^0+2^1+2^3)模板： ll p...
矩阵快速幂专题（矩阵快速幂入门、矩阵构造法、数论规律题）
林下的码路
08-04 2321
Link1：http://acm.hdu.edu.cn/showproblem.php?pid=1575 Tr A Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 3628    Accepted Submissio
数据分析 day02(二)
沉觞的博客
06-29 396
2.9 矩阵与矩阵运算 2.10 数组排序
矩阵的初等变换 hdu
最新发布
07-27
初等变换是指对矩阵进行一系列基本操作，包括交换两行（列），将某行（列）乘以一个非零标量，将某行（列）的倍数加到另一行（列）。这些操作可以用来简化矩阵，求解线性方程组，计算矩阵的秩等。 在求解线性方程组...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值