poj3046

dp,可以再优化。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 using namespace std;
 int t,a,s,b,tmp;
 int x[];
 int dp[][];
 const int mod = 1e6;
 void solve(int res)
 {
     for(int i = ; i <= x[]; i++)
     {
         dp[][i] = ;
     }
     for(int i = ; i < t+; i++)
     {
         memset(dp[i & ], , sizeof(dp[i & ]));
         for(int j = ; j <= res; j++)
         {
             for(int k = ; k <= x[i] && k <= j; k++)
             {
                 dp[i & ][j] = (dp[i & ][j] + dp[(i-) & ][j-k]) % mod;
             }
         }
     }
 }

 int main()
 {
     cin >> t >> a >> s >> b;
     for(int i = ; i < a; i++)
     {
         scanf("%d", &tmp);
         x[tmp]++;
     }
     int cnt = ;
     solve(a);
     for(int i = s; i <= b; i++)
     {
         cnt = (cnt + dp[t & ][i]) % mod;
     }
     cout << cnt << endl;
     return ;
 }