计算f(n)=1^2+2^2+3^2+.......+n^2.

/*
f(n)=1^2+2^2+3^2+.......+n^2.
f(n) = 1/6*n*(n+1)*(2n+1)
1/6可以看做上面的数除以6.而对于除法，我们有以下的规定：
(a+b)%c = (a%c+b%c)%c;
(a*b)%c=((a%c)*(b%c))%c;
(a^b)%c=(a%c)^b;
我们可以将6分成2 3.其中（n）*(n+1)这两个数一定可以被2整除，所以我们只需要考虑剩下的3便可以。

*/
using namespace std;
int main()
{
    long long n,ans = 1,num;
    while(cin>>n)
    {
        ans = n*(n+1);
        ans = ans/2;
        if(ans%3==0)
        {
            ans = ans/3;
            num = (ans%1007)*((2*n+1)%1007);
            num =num %1007;
        }
        else
        {
            num = (2*n+1)/3;
            num = num%1007;
            num = (ans%1007)*num;
            num = num%1007;
        }
        cout<<num<<endl;
    }
    return 0;
}