湖南国庆模拟赛day1 分组

/*
排序，设状态为i,j,k，分别表示前i个数还有j组未充满，不和谐程度为k的方案数，分四种情况讨论：这个数自成一组、这个数为最小值开一个2个数以上的组，把这个数填到别的组中去，把这个数作为最大值填到别的组中去，如果单纯计算k，那么复杂度o(n^2*sumk)，显然超时，可以用差分的办法去处理这个极差，这样复杂度o(n^2*k)
*/
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int mod = 1000000007;
const int maxn = 205,maxm = 1005;
ll n,m,s[maxn],f[2][maxn][maxm],cur;
int read(){
    char ch=getchar();
    int x=0,f=1;
    while(!(ch>=‘0‘&&ch<=‘9‘)){if(ch==‘-‘)f=-1;ch=getchar();};
    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+(ch-‘0‘);ch=getchar();};
    return x*f;
}
int main(){
    n = read();
    m = read();
    for(int i = 1;i <= n;i++)s[i] = read();
    sort(s+1,s+1+n);
    f[0][0][0] = f[0][1][0] = 1;
    ll tv = 0,tk = 0;
    for(int i = 2;i <= n;i++){
        for(int j = 0;j <= n;j++){
            for(int k = 0;k <= m;k++){
                tv = f[cur][j][k];
                f[cur][j][k] = 0;
                if(!tv) continue;
                tk = (s[i]-s[i-1])*j + k;
                if(tk > m) continue;
                f[cur^1][j][tk] = (f[cur^1][j][tk] + tv) % mod;
                f[cur^1][j+1][tk] = (f[cur^1][j+1][tk] + tv) % mod;
                if(j) f[cur^1][j-1][tk] = (f[cur^1][j-1][tk] + (tv * j) % mod) % mod;
                f[cur^1][j][tk] = (f[cur^1][j][tk] + (tv * j) % mod) % mod;
            }
        }
        cur ^= 1;
    }
    ll ans = 0;
    for(int i = 0;i <= m;i++) ans = (ans + f[cur][0][i]) % mod;
    cout<<ans;
    return 0;
}