比赛 201712练习 评测结果 AAAAAAAAAA
题目名称 数学序列 最终得分 100
用户昵称 Ceres 运行时间 0.004 s
代码语言 C++ 内存使用 0.29 MiB
提交时间 2018-01-08 15:49:35
显示代码纯文本
#include <cstdio>
#include <cstring>
#define mod 7
using namespace std;
struct matrix
{
	int E[2][2];
	void init(int a,int b)
	{
		memset(E,0,sizeof(E));
		E[0][1]=1;
		E[1][1]=a;
		E[1][0]=b;
	}
	matrix operator*(matrix& a)const
	{
		matrix ans;
		memset(ans.E,0,sizeof(ans.E));
		for(int i=0;i<2;i++)
			for(int j=0;j<2;j++)
				for(int k=0;k<2;k++)
					ans.E[i][j]+=E[i][k]*a.E[k][j];
		return ans;
	}
	matrix operator%(int a)const
	{
		matrix ans;
		memset(ans.E,0,sizeof(ans.E));
		for(int i=0;i<2;i++)
			for(int j=0;j<2;j++)
				ans.E[i][j]=E[i][j]%a;
		return ans;
	}
}f,base;
void quick_pow(int n)
{
	memset(f.E,0,sizeof(f.E));
	f.E[0][0]=f.E[1][0]=1;
	while(n)
	{
		if(n&1)
			f=(base*f)%mod;//f=(f*base)%mod;注意这里,矩阵乘法不具有交换律 
		base=(base*base)%mod;
		n>>=1;
	}
}
int a=0,b=0,n=0;
int main()
{
	freopen("number1.in","r",stdin);
	freopen("number1.out","w",stdout);
	while(scanf("%d%d%d",&a,&b,&n)!=EOF)
	{
		a%=mod;
		b%=mod;
		if(n<3)
			printf("1\n");
		else
		{
			base.init(a,b);
			quick_pow(n-2);
			printf("%d\n",f.E[1][0]);
		}
	}
	return 0;
}