算法提高 递推求值

问题描述

　　已知递推公式：

　　F(n, 1)=F(n-1, 2) + 2F(n-3, 1) + 5,

　　F(n, 2)=F(n-1, 1) + 3F(n-3, 1) + 2F(n-3, 2) + 3.

　　初始值为：F(1, 1)=2, F(1, 2)=3, F(2, 1)=1, F(2, 2)=4, F(3, 1)=6, F(3, 2)=5。
　　输入n，输出F(n, 1)和F(n, 2)，由于答案可能很大，你只需要输出答案除以99999999的余数。

输入格式

　　输入第一行包含一个整数n。

输出格式

　　输出两行，第一行为F(n, 1)除以99999999的余数，第二行为F(n, 2)除以99999999的余数。

样例输入

4

样例输出

14

21

数据规模和约定

　　1<=n<=10^18。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;

const int mod=99999999;

struct Martix{
	ll mar[7][7];
    Martix operator*(const Martix& b)const{
		Martix ans;
		for(int i=0;i<7;i++){
			for(int j=0;j<7;j++){
				ans.mar[i][j]=0;
				for(int k=0;k<7;k++){
					ans.mar[i][j]=(ans.mar[i][j]+((*this).mar[i][k]*b.mar[k][j])%mod+mod)%mod;
				}
			}
		}
		return ans;
	}
}E,ST;
ll num[]={6,5,1,4,2,3,1};
ll st[7][7]={{0,1,0,0,2,0,5},{1,0,0,0,3,2,3}};

void init(){
	for(int i=0;i<7;i++)E.mar[i][i]=1;
	for(int i=0;i<2;i++){
		for(int j=0;j<7;j++)
			ST.mar[i][j]=st[i][j];
	}
	ST.mar[2][0]=ST.mar[3][1]=ST.mar[4][2]=ST.mar[5][3]=ST.mar[6][6]=1;
}

Martix pow_(ll n){
	Martix ans=E,base=ST;
	while(n){
		if(n&1)
			ans=ans*base;
		n/=2;
		base=base*base;
	}
	return ans;
}

/*
这道题怎么讲，我会，但是不熟，一个点就是我当时赋值的时候直接就错了，应该是大的值在上面，小的值在下面

还有就是我们编程的细节的问题，犯了很多小的错误
他说是答案取余，我们只是把过程中取了余，最后的答案，也就是最重要的一个点，我们没有去管，就得了60
还有就是类型的传递的问题，一个long long 的类型，生生的传成了int类型，自己不是很想说话
还是心态不稳，不能稳如老狗啊
*/


int main(){
	init();
	ll n;
	scanf("%I64d",&n);
	//printf("%I64d**\n",n);
	if(n<4){
		//printf("*****\n");
		int pos=(3-n)*2;
		//printf("pos:%d\n",pos);
		printf("%I64d\n",num[pos]);
		printf("%I64d\n",num[pos+1]);
		return 0;
	}
	ST=pow_(n-3);
	ll ans1=0,ans2=0;
	for(int i=0;i<7;i++){
		ans1=(ans1+num[i]*ST.mar[0][i])%mod;
		ans2=(ans2+num[i]*ST.mar[1][i])%mod;
	}
	printf("%I64d\n",ans1%mod);
	printf("%I64d\n",ans2%mod);
	return 0;
}