2分法求 2^n 和 矩阵的n次方

#include "stdafx.h"
#include <iostream>
using namespace std;

//a^n
int f(int a, int n)
{
	if(n == 0)  return 1;
	if(n == 1)  return a;
	int ans = f(a,n/2);
	if(n&1)
	{
		return ans*ans*a;
	}
	else
		return ans*ans;
}

int main()
{
	cout<<f(2,5)<<endl;
	return 0;
}

//2分法求矩阵冥次
#include "stdafx.h"
#include <iostream>
#include <algorithm>
using namespace std;
//如果矩阵过大或者冥过多都会造成栈溢出
const int MAXN = 30;

typedef struct Matrix
{
	int s[MAXN][MAXN];
}M;

int n;

//a*b
M Multiply(M a, M b)
{
	M c;
	memset(c.s,0,sizeof(c.s));
	for(int i=1; i<=n; i++)
	{
		for(int j=1; j<=n; j++)
		{
			for(int k=1; k<=n; k++)
			{
				c.s[i][j] += a.s[i][k] * b.s[k][j];
			}
		}
	}
	return c;
}


//a^t
M Pow(M a, int t)
{
	M b;
	if(t == 0)
	{
		memset(b.s,0,sizeof(b.s));
		for(int i=1; i<=n; i++)
		{
			b.s[i][i] = 1;
		}
		return b;
	}
	else
	{
		M k = Pow(a,t/2);
		if(t&1) //t为奇数
		{
			return Multiply(Multiply(k,k),a);
		}
		else
		{
			return Multiply(k,k);
		}
	}
}


int main()
{
	M a;
	int t;
	cin>>n>>t;
	for(int i=1; i<=n; i++)
	{
		for(int j=1; j<=n; j++)
		{
			cin>>a.s[i][j];
		}
	}
	M b = Pow(a,t);
	for(int i=1; i<=n; i++)
	{
		for(int j=1; j<=n; j++)
		{
			cout<<b.s[i][j]<<" ";
		}
		cout<<endl;
	}
	return 0;
}