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