题解

#include<iostream>
using namespace std;
const int N=10007;
long long s[100][100];
int main(){
    int n,r;
    cin>>n>>r;
    for(int i=0;i<n;i++)s[0][i]=1;
    for(int i=1;i<r;i++){
        for(int j=i;j<n;j++){
            s[i][j]=s[i][j-1]*(i+1)+s[i-1][j-1];
            s[i][j]%=N;
        }
    }
    cout<<s[r-1][n-1];
    return 0;
}