hdu2802 循环节思维

F(N)

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 5026    Accepted Submission(s): 1788

Problem Description

Giving the N, can you tell me the answer of F(N)?

 

Input

Each test case contains a single integer N(1<=N<=10^9). The input is terminated by a set starting with N = 0. This set should not be processed.

 

Output

For each test case, output on a line the value of the F(N)%2009.

 

Sample Input

  
  

   
   1
2
3
0
  
  

 

Sample Output

  
  

   
   1
7
20
  
  

 

Source

HDU 2009-4 Programming Contest

解题思路：

这道题目算是给我提醒了一个方向，这道题其实不难，就是由于惯性思维，一拿到这个题目就在狂推公式，也用矩阵快速幂试了，发现都不得要得要点。

总的来说，循环节就像是数学中的周期函数，他的值只有一定范围，过了这个范围，就是重复的。例如周期函数：F(0)=1,F(1)=2,F(2)=3,F(3)=4,所以F(n)=F(n-1%4).

具体怎么求出这个范围，你可以打表出来看，写一个循环，自己看规律。以后一般方法做不出来的时候可以往这方面想一想。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
long long f[110000];
void  init()
{
	int i;
	f[1]=1;
	f[2]=7;
	for(i=3;i<=5000;i++)
	f[i]=(f[i-2]+3*i*i-3*i+1)%2009;
}
int main()
{
	long long  n,i;
	init();
	while(scanf("%lld",&n),n)
	{
		//for(i=3;i<=110000;i=i+2)打表法找出循环节
		//if(f[i]==1)
		//break;
		//printf("%I64d\n",i-1);
		//for(i=4;i<=110000;i=i+2)
		//if(f[i]==7)
		//break;
		//printf("%I64d\n",i-2);
	    //printf("%I64d\n",f[4018]);
		printf("%lld\n",f[n%4018]);
	}
	return 0;
}