循环节4018

1 2 3 0
1 7 20
/* ***********************************************
Author :CKboss
Created Time :2014年12月14日 星期日 12时29分33秒
File Name :HDOJ2802.cpp
************************************************ */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cstdlib>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <cmath>
using namespace std;
const int mod=2009;
int f[5000];
int Get(int x)
{
int part1=((x*x)%mod*x)%mod;
int y=x-1;
int part2=((y*y)%mod*y)%mod;
return ((f[x-2]+part1)%mod+mod-part2)%mod;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int n;
f[1]=1; f[2]=7;
for(int i=3;i<=4018;i++) f[i]=Get(i);
while(scanf("%d",&n)!=EOF&&n)
{
n=n%4018;
if(n==0) puts("0");
else printf("%d\n",f[n]);
}
return 0;
}
原文:http://blog.csdn.net/ck_boss/article/details/41924343