循环节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