CDOJ 1268 Open the lightings

#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <queue>
#include <stack>
#include <map>
#include <vector>
using namespace std;

const int maxn=1000+10;
long long MOD=1e9+7;
long long f[maxn];
long long c[maxn][maxn];

void init()
{
    f[0]=1,f[1]=1;
    for(int i=2;i<=1000;i++) f[i]=(f[i-1]+((i-1)*f[i-2])%MOD)%MOD;

    c[1][1]=1;
    for(int i=1;i<=1000;i++) c[i][0]=1;
    for(int i=2;i<=1000;i++)
    {
        for(int j=1;j<=1000;j++)
        {
            c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;
        }
    }
}

int main()
{
    init();
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        printf("%lld\n",f[m-1]*f[n-m]%MOD*c[n-1][m-1]%MOD);
    }

    return 0;
}