矩阵快速幂

/*
	假设 A 是 m*p 的矩阵 , B 是 p*n 的矩阵
	记 C = AB (C 是 矩阵 A与B的乘积)
	那么 C 是 m*n 的矩阵 
*/
for (int i = 1;i <= m;++i)//A的行 
{
	for (int j = 1;j <= n;++j)//B的列 
	{
		for (int k = 1;k <= p;++k)//通过公式求C 
		{
			C[i][j] += A[i][k]*B[k][j];
		}
	}
}

求一个矩阵的n次幂

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define LL long long
char str[100005];
struct A
{
    int m[3][3];
}bb,aa;

A multiply(A a,A b)//矩阵相乘
{
    A temp;
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++){
            temp.m[i][j]=0;
            for(int k=0;k<2;k++)
                temp.m[i][j]=(temp.m[i][j]+a.m[i][k]*b.m[k][j])%7;
        }
    return temp;
}
void AA()
{
    aa.m[0][0]=aa.m[1][1]=1;
    aa.m[1][0]=aa.m[0][1]=0;
}
A quick_mod(LL n)//快速幂
{
    while(n)
    {
        if(n&1)
            aa=multiply(aa,bb);
        n>>=1;
        bb=multiply(bb,bb);
    }
    return aa;
}//模板部分

int main()
{
    int num;
    while(~scanf("%s",str))
    {
        num=0;
        int len=strlen(str);
        for(int i=0;i<len;i++)
            num=(num*10+str[i]-'0')%2016;
        scanf("%d%d",&bb.m[0][0],&bb.m[0][1]);
        scanf("%d%d",&bb.m[1][0],&bb.m[1][1]);
        AA();
        bb=quick_mod(num);
        printf("%d %d\n",bb.m[0][0],bb.m[0][1]);
        printf("%d %d\n",bb.m[1][0],bb.m[1][1]);
    }
    return 0;
}