Recursive sequence

1 矩阵快速幂,好记性不如烂笔头,这是我第二次写这一题了,我先写下博客,方便自己第三次看……啊……

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<cmath>

using namespace std;
typedef long long ll;
const ll Mod=2147493647;
const ll maxn=10;
ll N,A,B,T;
struct node
{
	ll jz[maxn][maxn];
};
node Mul(node a,node b)
{
	node res;
	memset(res.jz,0,sizeof(res.jz));
	for(int i=0;i<7;i++)
	{
		for(int j=0;j<7;j++)
		{
			for(int k=0;k<7;k++)
			{
				res.jz[i][j]=(res.jz[i][j]+a.jz[i][k]*b.jz[k][j]%Mod+Mod)%Mod;
			}
		 } 
	}
	return res;
}
node Pow_mod(node a,int b)
{
	node res;
	memset(res.jz,0,sizeof(res.jz));
	for(int i=0;i<7;i++)
	{
		for(int j=0;j<7;j++)
		{
			if(i==j)
			res.jz[i][j]=1;
		}
	 } 
	 while(b)
	 {
	 	if(b&1)
	 	res=Mul(res,a);
		a=Mul(a,a);
		b>>=1; 
	 }
	 return res;
 }
int main()
{
	scanf("%lld",&T);
	while(T--)
	{
		scanf("%lld%lld%lld",&N,&A,&B);
		if(N==1)
		printf("%lld\n",A);
		else if(N==2)
		printf("%lld\n",B);
		else
		{
			node jie,jie1;
			memset(jie.jz,0,sizeof(jie.jz));
			memset(jie1.jz,0,sizeof(jie1.jz));
			jie.jz[0][0]=1;jie.jz[3][3]=1;
			jie.jz[0][1]=2;jie.jz[3][4]=3;
			jie.jz[0][2]=1;jie.jz[3][5]=3;
			jie.jz[1][0]=1;jie.jz[3][6]=1;
			jie.jz[2][2]=1;jie.jz[4][4]=1;
			jie.jz[2][3]=4;jie.jz[4][5]=2;
			jie.jz[2][4]=6;jie.jz[4][6]=1;
			jie.jz[2][5]=4;jie.jz[5][5]=1;
			jie.jz[2][6]=1;jie.jz[5][6]=1;
			jie1.jz[0][0]=B;jie.jz[6][6]=1;
			jie1.jz[1][0]=A;
			jie1.jz[2][0]=pow(3,4);
			jie1.jz[3][0]=pow(3,3);
			jie1.jz[4][0]=pow(3,2);
			jie1.jz[5][0]=3;
			jie1.jz[6][0]=1;
			node res=Pow_mod(jie,N-2);
			res=Mul(res,jie1);
			printf("%lld\n",res.jz[0][0]);	
		}
	}
	return 0;
}
全部评论

相关推荐

感觉这一周太梦幻了,就像一个梦,很不真实~~~感觉这个暑期,我的运气占了99成,实力只有百分之一4.15上午&nbsp;腾讯csig&nbsp;腾讯云部门,面完秒进入复试状态4.16下午&nbsp;美团优选供应链部门,4.18上午发二面4.17晚上&nbsp;阿里国际一面,纯拷打,面完我都玉玉了4.18下午&nbsp;阿里国际二面,是我们leader面的我,很轻松~~4.18晚上&nbsp;约了hr面4.19上午&nbsp;hr面,下午两点口头oc4.19晚上&nbsp;意向书说起来我的暑期好像一次都没挂过~~~~~难道我是天生面试圣体?----------------------------------------------------------------------六个月前,我还是0项目0刷题,当时想的是先把论文发出来再去找实习。结果一次组会,老师打破了我的幻想(不让投B会,只让投刊或者A)我拿头投啊!!!然后就开始物色着找实习,顺便做完了mit的6.s081,但是基本上还是没刷过题目-----------------------------------------------------------------------11月&nbsp;&nbsp;一次偶然的机会,面进了某个耳机厂的手环部门,大概是做嵌入式的,用的是CPP。12月&nbsp;莫名其妙拿到了国创的面试机会,0基础四天速成java基础!居然也给我面过了hhhhh,可能是面试没写题吧入职国创后的几个月,一直没活,天天搁那看剧,都快忘了还有暑期实习这回事了~~~~命运的齿轮在2.26开始转动,因为这一天美团开了,我开始慌了,因为那时的我什么都不会。lc,八股,sql全部是0进度。然后就开始了女娲补天,上班刷题,下班继续做之前的开源,顺便学一学八股。3月到现在,lc也刷到快200了,一天最多提交了47次~~~~~~~~~~八股根据别人的面经总结和博客,写了快十万字的笔记~~~~~~~~~~简历上的实习经历和开源,也努力去深挖了,写了几万字的记录~~~~~~所以面试的时候,基本上都能cover了,面试官问到的基础基本都会,不基础的我就把他往我会的地方引。结果好像还不错,基本上每个面试官评价都挺好的emmmmmmmm
投递阿里巴巴等公司10个岗位
点赞 评论 收藏
转发
1 收藏 评论
分享
牛客网
牛客企业服务