HDU 1005 number sequence

一道数学题，纠结了好长时间，10^8 的话，用O（N）的算法一定会超时，考虑到模7，是一个循环问题，重要的是找出循环节，注意不一定是从f[1]出开始循环。。。

#include<iostream>
#include<cstdio>
using namespace std;
int f[100];
int A,B;
long long N;
int find(int left,int right)
{
for(int i=1;i<left;i++)
{
if(f[i]==f[left]&&f[i+1]==f[right])return i;
}
return 0;
}
int main()
{
int i,j,k;
while(scanf("%d %d",&A,&B))
{
cin>>N;

if(!A&&!B&&!N)break;
f[1]=1;
f[2]=1;
if(N==1||N==2)
{
cout<<"1"<<endl;
continue;
}
int cycle=100;
for(i=3;i<=N;i++)
{
f[i]=(A*f[i-1]+B*f[i-2])%7;
j=find(i-1,i);
if(j){cycle=i-1-j;break;}
}
if(j==0)cout<<f[N]<<endl;
else
{
k=(N-j)%cycle;
if(k==0)cout<<f[j+cycle]<<endl;
else cout<<f[k+j]<<endl;
}
}
return 0;
}