题目连接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=88584#problem/D
题目大意:给定数列递推式,求第n项mod7余数
题目分析:由于n范围10^8,若用递推必超时。由抽屉原理可知,从某两项开始必会循环,故只需找到循环节即可。注意:最初若干项可能不在循环节内需特判。
代码如下:
#include<stdio.h>
int main()
{
int m,n,i,a,b,t,j,f[1000];
bool x;
scanf("%d%d%d",&b,&a,&n);
while(a!=0 || b!=0 || n!=0)
{
a=a%7; b=b%7;
for (i=1;i<1000;i++) f[i]=0;
f[1]=1; f[2]=1; x=false; i=2;
j=2; m=0;
while (x==false)
{
i=i+1;
f[i]=(a*f[i-2]+b*f[i-1])%7;
for (j=2;j<i;j++)
{
if (f[j]==f[i] && f[j-1]==f[i-1])
{
x=true;
t=i-j; //t为循环节
break;
}
}
}
if (n<j) printf("%d\n",f[n]); //特判循环节前特例
else
{
m=(n-j)%t;
printf("%d\n",f[j+m]);
}
scanf("%d%d%d",&b,&a,&n);
}
return 0;
}