海绵的博客

期望少一点,做的多一点

算法提高 递推求值
问题描述
  已知递推公式:

  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;
}

阅读更多
版权声明:小白一个,欢迎各位指错。 https://blog.csdn.net/qq_36424540/article/details/80343902
个人分类: .....快速幂
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

算法提高 递推求值

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭