UVA-10518 How Many Calls?

问用递归函数求斐波那契数f(n)时，会调用多少次函数

设调用次数为C(n)，C(n)=C(n-1)+C(n-2)+1，C(1)=C(0)=1

可得矩阵：用矩阵快速幂求解

C(n)      | 1   1   1   |    C(n-1)

C(n-1)   |  1   0   0  |    C(n-2)

  1         |  0   0   1  |        1

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
ll MOD;
const int N=3;
struct node
{
	ll a[10][10];
};
node shu,ans,mp;
//shu是输入的矩阵，ans是所求答案
node matrix(node x,node y)
{
    for(int i=1;i<=N;i++)
        for(int j=1;j<=N;j++){
            mp.a[i][j]=0;
            for(int p=1;p<=N;p++)
                mp.a[i][j]=(mp.a[i][j]+x.a[i][p]*y.a[p][j]+MOD)%MOD;
            //矩阵乘法 
        }
    return mp;
}
void work(ll k)
{//矩阵快速幂 
	for(int i=1;i<=N;i++)
			for(int j=1;j<=N;j++)
				ans.a[i][j]=0;
	for(int i=1;i<=N;i++) ans.a[i][i]=1;
	node t=shu;
    while(k){
        if(k&1) 
            ans=matrix(ans,t);
        k>>=1;
        t=matrix(t,t);
    }
}
int main()
{
	ll n,b;
	int kase=0;
	while(~scanf("%lld%lld",&n,&b))
	{
		if(n==0&&b==0) break; 
		MOD=b;
		ll ret;
		if(n<2)
		{
			ret=1%MOD;
		}
		else
		{
			memset(shu.a,0,sizeof(shu.a));
			shu.a[1][1]=shu.a[1][2]=shu.a[1][3]=1;
			shu.a[2][1]=1;
			shu.a[3][3]=1;
			work(n-1);
			ret=(ans.a[1][1]+ans.a[1][2]+ans.a[1][3])%MOD;
		}
		printf("Case %d: %lld %lld %lld\n",++kase,n,b,ret);
	}
	return 0;
}