这道题也是厉害,刚看题就想到了递归,但一看n的范围就果断放弃了,然后想到用常系数线性齐次递推关系转化为关于r^n的形式求解,但因为涉及根号会有误差,也pass,想了半天没结果,上网一搜发现是用枚举求解的,很显然f[n]=(A*a[n-1]+B*a[n-2])%7中,A和B是常数,而a[n-1]和a[n-2]各只有7种可能,即a[n]只有49种可能结果,最多运算50次就能产生循环。Note:一定要考虑n比较小时没有产生循环的情况啊!代码如下:
#include <iostream>
#include <cstring>
using namespace std;
int main()
{
int a,b,i,j;
int r[60];
long n;
bool mark;
while(cin>>a>>b>>n&&(a!=0||b!=0||n!=0)){
memset(r,-1,sizeof(r));
r[0]=r[1]=r[2]=1;
mark=false;
for(i=3;i<=n;i++){
r[i]=(a*r[i-1]+b*r[i-2])%7;
for(j=2;j<=i-1;j++)
if(r[i]==r[j]&&r[i-1]==r[j-1]) {mark=true;break;}//搜索之前有没有循环,晕,j一定是小于等于i-1的
if(mark) break;
}
if(mark)
cout<<r[j+(n-i)%(i-j)]<<endl;
else
cout<<r[n]<<endl;
}
return 0;
}