F(N)---hdu2802（寻找循环节）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2802

f[1] = 1;

f[2] = 7;

f[n] = (f[n-2] - (n-1)*(n-1)*(n-1)+ n*n*n) % 2009, n>2;

给你一个n让你输出f[n];（n<1e9）由于n较大，所以不能直接输出，结果是对2009求余所以一定存在循环节，由于 f(n) 只与 f(n-2) 和 n 有关,

所以我们只需判断 G[n][f[n-2]] 第一次出现和第二次出现的次数即可，差值就是循环节；

#include <stdio.h>
#include <algorithm>
#include<string.h>
#include<queue>
using namespace std;

#define met(a, b) memset(a, b, sizeof(a))
#define MOD 1000000007
#define N 2010
#define INF 0x3f3f3f3f

typedef long long LL;

const int maxn = ;

int G[maxn][maxn], F[maxn*maxn] = {,,,};

int main()
{
    int k = , n;
    G[][] = k++;

    for(int i=; i<maxn*maxn; i++)
    {
        F[i] = (F[i-]-(i-)*(i-)*(i-) + i*i*i)%;

        if(G[i%][F[i-]])///说明n和f[n-2]对应的已经出现过了，说明存在循环节
        {
            /// printf("%d %d\n", G[i%2009][F[i-2]], k);
            break;
        }

        G[i%][F[i-]] = k++;
    }

    while(scanf("%d", &n), n)
    {
        printf("%d\n", F[n%]);
    }
    return ;
}