题目:http://acm.hdu.edu.cn/showproblem.php?pid=1005
/*
本题的n很大不用快速幂肯定超时,根据递推公式f(n) =A*f(n-1)+Bf(n-2)构造矩阵
*/
#include<iostream>
#include <stdio.h>
#include<algorithm>
using namespace std;
int a,b,n;
struct Mat{
int num[3][3];
};
Mat mult(Mat a,Mat b){
Mat c;
for(int i=1;i<=2;i++)
for(int j=1;j<=2;j++){
c.num[i][j]=0;
for(int k=1;k<=2;k++)
c.num[i][j]=(c.num[i][j]+a.num[i][k]*b.num[k][j])%7;
}
return c;
}
Mat my_pow(Mat x,int y){
Mat ans;
for(int i=1;i<=2;i++)
for(int j=1;j<=2;j++){
if(i==j)
ans.num[i][j]=1;
else
ans.num[i][j]=0;
}
for(;y;y>>=1){
if(y&1)
ans=mult(ans,x);
x=mult(x,x);
} return ans;
}
int main(){
while(scanf("%d%d%d",&a,&b,&n),(a||b||n)){
Mat A;
A.num[1][1]=a;
A.num[1][2]=b;
A.num[2][1]=1;
A.num[2][2]=0;
if(n ==1){
printf("%d\n",a);//不判会超时
continue;
}
if(n ==2){
printf("%d\n",b);
continue;
}
A=my_pow(A,n-2);
printf("%d\n",(A.num[1][1]+A.num[1][2])%7);//虽然在上面mod7了,结果竟然忘了,WA了几次
}
return 0;
}
/*
还有一个更好的办法,代码简洁,但是很不容易想到,利用数论的知识,结果mod7之后是有周
期性的(7*7=49),结果是在0--6之间的数,直接打表就可以了。
*/
#include<iostream>
#include <stdio.h>
#include<algorithm>
using namespace std;
int a,b,n,i;
int arry[60];
int main(){
while(scanf("%d%d%d",&a,&b,&n),(a||b||n)){
arry[1]= arry[2] = 1;
for(i =3;i < 60;i++)
arry[i] = (a*arry[i-1] +b*arry[i-2])%7;
n%=49;
printf("%d\n",arry[n]);
}
return 0;
}