HDOJ 1005 Number Sequence

原题连接

斐波拉及数列形题目。以前上课作业有做过。但是明显水很多~几乎不考虑到时间复杂度问题。
 
这道题目有两种解题方法（先忽视时间问题）

首先就是典型的递归算法，易于理解。但是问题很显然，不断的递归必然导致超时。

所以，我采用了较为快速的递推算法。

f(i)=(A*f(i-1)+B*f(i-2))%7

显然知道上式只要利用两个变量就可以实现本题的解答。

但是看看N的大小 100,000,000       光光数零就要好久了。所以，直接循环递推是不可能的。

由数论知识可得，这样的取余必然有一个循环点。这个循环点根据网上某大神给出的是49.但是对于一些特殊的数49也并不是都适用的。所以，必须给出一个更适合的循环数。而我们要做的不是自己拿笔一点点算出来而是让计算去代替我们算。

一开始我们的循环第一个和第二个数必须是1 1显然我们可以利用这两个数来取得循环点：
        for(i=3;i<200;i++){
            s[i]=(A*s[i-1]+B*s[i-2])%7;
            if(s[i-1]==1&&s[i]==1)
                break;}
设置为200以内的原因是既然49不能全部满足那么就尽量设大~200差不多了。
取得循环点后面就很好做了。
代码如下：

<span style="font-size:18px;">#include <stdio.h>
int s[50];
int count(unsigned long A,unsigned long B,unsigned long n)
{
	unsigned long i;
	if(n==1||n==2)
		return 1;
	else
	{
		s[1]=s[2]=1;
		for(i=3;i<200;i++){
			s[i]=(A*s[i-1]+B*s[i-2])%7;
			if(s[i-1]==1&&s[i]==1)
				break;}
		i-=2; //一个周期包含元素个数
		n=n%i; 
	}
	if(n==0)
		return s[i];
	else
	return s[n];
}
int main()
{
	unsigned long A,B,n,res;
	unsigned long t1,t2;
	while(scanf("%ld %ld %ld",&A,&B,&n)!=EOF)
	{
		if(A==0&&B==0&&n==0)
			break;
		res=count(A,B,n);
		printf("%d\n",res);
	}
	return 0;
}</span>