#include <iostream>
#include <cstdio>
#define N 5005
#define mod 666013

using namespace std;

int d[3][N], n, rez;

int main()
{
    freopen("100m.in", "r", stdin);
    freopen("100m.out", "w", stdout);

    scanf("%d\n", &n);
    d[1][1]=d[2][1]=1;
    for(int i=2;i<=n;i++)
    {
        for(int j=2;j<=i;j++)
            d[2][j]=(((d[1][j]+d[1][j-1])%mod)*1LL*j)%mod;
        for(int j=2;j<=i;j++)
            d[1][j]=d[2][j];
    }
    rez=0;
    for(int i=1;i<=n;i++)
        rez=(1LL*rez+d[1][i])%mod;
    printf("%d", rez);
    return 0;
}