递推式

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998'244'353;
long long n, k, S[1005][1005];
int main()
{
    cin >> n >> k;
    S[0][0] = 1;
    for (int i = 1; i <= n; i++)
    {
        S[i][0] = 0;
        for (int j = 1; j <= i; j++)
            S[i][j] = (S[i - 1][j - 1] + j * S[i - 1][j] % MOD) % MOD;
    }
    cout << S[n][k] << "\n";
    return 0;
}