数论专题 K-blocks（矩阵快速幂）

题意

用四种颜色涂长度为n的方块，求其中两种颜色均为偶数的方案数

思路

这题可以用递推dp的方式求出一个数列。

设a,b,c分别为均为偶数，一个为奇数，均为奇数三种情况

两个均为偶数的状态由其中一个为奇数转移来

a[i]=2*a[i-1]*b[i-1]

b[i]=2*a[i-1]+2*b[i-1]+2*c[i-1]

c[i]=b[i-1]+2c[i-1]

于是可以用矩阵快速幂来加快方程的递推；

构造矩阵为

2 1 0

2 2 2

0 1 2

进行快速幂

代码如下

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int mod=10007;
struct node{
	int a[3][3];
};
node cal(node x,node y){
	node b;
	memset(b.a,0,sizeof b.a);
	for(int i=0;i<3;i++){
		for(int j=0;j<3;j++){
			for(int k=0;k<3;k++){
				b.a[i][j]=(b.a[i][j]+x.a[i][k]*y.a[k][j])%mod;
			}
		}
	}
	return b;
}//矩阵相乘函数
int fastpow(int n){//快速幂
	node x;
    x.a[0][0]=2,x.a[0][1]=1,x.a[0][2]=0;
    x.a[1][0]=2,x.a[1][1]=2,x.a[1][2]=2;
    x.a[2][0]=0,x.a[2][1]=1,x.a[2][2]=2;
    node r;
    r.a[0][0]=2,r.a[0][1]=1,r.a[0][2]=0;
    r.a[1][0]=2,r.a[1][1]=2,r.a[1][2]=2;
    r.a[2][0]=0,r.a[2][1]=1,r.a[2][2]=2;
    while(n)
    {
        if(n%2) r=cal(r,x);
        x=cal(x,x);
        n=n/2;
    }
    long long ans=r.a[0][0];
    return ans;
}
int main(){
	int n,t;
	cin>>t;
	while(t--){
	  cin>>n;
	  int ans=fastpow(n-1);
	  cout<<ans<<endl;
    }
	return 0;
}