C-Bombs

#include <bits/stdc++.h>
#pragma warning (disable : 26451)
#pragma warning (disable : 4996)
#define mod 666013
using namespace std;
ifstream fin("bombs.in");
ofstream fout("bombs.out");
int main() {
    long long n;
    fin >> n;
    vector<vector<long long>> dp(n + 1, vector<long long>(n + 1));
    dp[0][0] = 1;
    for (long long i = 1; i <= n; ++i) {
        for (long long j = 1; j <= i; ++j) {
            dp[i][j] = dp[i - 1][j - 1] + j * dp[i - 1][j];
            dp[i][j] %= mod;
        }
    }
    long long ans = 0;
    for (long long i = 1; i <= n; ++i) {
        ans += dp[n][i];
        if (ans >= mod) {
            ans -= mod;
        }
    }
    fout << ans;
}