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
解题思路:
这道题目算是给我提醒了一个方向,这道题其实不难,就是由于惯性思维,一拿到这个题目就在狂推公式,也用矩阵快速幂试了,发现都不得要得要点。
总的来说,循环节就像是数学中的周期函数,他的值只有一定范围,过了这个范围,就是重复的。例如周期函数: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;
}