快速幂+矩阵快速幂

---

快速幂：

计算a^b
复杂度：log（n）
例如计算a的10次方，10用二进制表示就是1010（重点！！！）。用二进制表示之后我们就可以明显的看出来1010=1000+0010，也就是10=8+2。那么a^10= a^2 * a^8，这也就是快速幂的思想。

#include <iostream>
using namespace std;

int qPow(int a,int b)
{
	int ans=1;//开始为单位元
	while( b!=0 )
	{
		if( b&1 ) ans*=a;//如果指数在这一位是1，就乘起来
		a*=a;//随着每一位的移动，权值也要增加（a->a^2->a^4->a^8...）
		b>>=1;
	}
	return ans;
}

int main()
{
	int a,b;
	cin >> a >> b; 
	cout << qPow(a,b) << endl;
	return 0;
 } 

---

总的来说，快速幂没什么用，但是思想很关键。然后从它的原理就可以推出矩阵快速幂。

矩阵快速幂

操作一样的，只是把数改成了矩阵，乘法改成矩阵乘法
洛谷3390

#include <iostream>
using namespace std;
#define mod 1000000007
#define ll long long
int n; 

struct Mat{
	ll m[105][105];
};

Mat e;

Mat mul(Mat x,Mat y)
{
	Mat c;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
			c.m[i][j]=0;
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			for(int k=1;k<=n;k++)
			{
				c.m[i][j]+=x.m[i][k]*y.m[k][j];
				c.m[i][j]%=mod;
			}
		}
	}
	return c;
}

Mat qPow(Mat x,ll k)
{
	Mat ans=e;//单位元
	while( k )
	{
		if( k&1 )
			ans=mul(ans,x);//矩阵乘法
		x=mul(x,x);
		k>>=1;
	}
	return ans;
}

int main()
{
	ll k;
	Mat a;
	cin >> n >> k;
	for(int i=1;i<=n;i++)
	{
		e.m[i][i]=1;
		for(int j=1;j<=n;j++)
		{
			cin >> a.m[i][j];
		}
	}
	a=qPow(a,k);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if( j==n ) cout << a.m[i][j] << endl;
			else cout << a.m[i][j] << " ";
		}
	}
	return 0;
}