矩阵快速幂

矩阵乘法

const int N = 3;
struct Matrix{
	int x[N][N];
	Matrix(){memset(x,0,sizeof(x));}
	Matrix(const Matrix& t){memcpy(x,t.x,sizeof(x));}
	Matrix operator*(const Matrix& t)const{
	Matrix ret;
	for(int i = 0;i < N;i++)
		for(int j = 0;j < N;j++)
		{
			ret[i][j] = 0;
			for(int k = 0;k < N;k++)
				ret[i][j] += x[i][k] * t.x[k][j];
		}
		return ret;
	}
	
}

Matrix quick_pow (Matrix a,int x)
{
	Matrix ret;
	//定义单位矩阵以 2*2矩阵为例
	ret.x[0][1] = ret.x[1][0] = 0;
	ret.x[0][0] = ret.x[1][1] = 1;
	while(x){
	if(x&1) ret = ret * a;
	a = a*a;
	x >>= 1;
	}
	return ret;
}

利用矩阵的幂进行运算

线性递推

利用矩阵快速幂求解

自然数幂和

矩阵快速幂优化DP

连续格子染色

状态表示 A[i] 表示 i 个格子， 红绿均为偶数的染色方案 B[i] 表示 i 个格子，红绿均为奇数的方案数 C[i] 表示 i 个格子 红绿只有一个为奇数的方案数

衬衫快乐值

状态定义 f[i][j] 表示第 i 天， 穿的衣服为 j 所获得的快乐值得总和最大值