Submission #1177608

#TimeUsernameProblemLanguageResultExecution timeMemory
1177608mychecksedadSkyscraper (JOI16_skyscraper)C++20
100 / 100
137 ms242788 KiB
/* Author : Mychecksdead  */
#include<bits/stdc++.h>
using namespace std;
#define ll long long int
#define MOD (1000000000+7)
#define MOD1 (998244353)
#define pb push_back
#define all(x) x.begin(), x.end()
#define en cout << '\n'
#define ff first
#define ss second
#define pii pair<int,int>
#define vi vector<int>
const int N = 100+5, M = 1000+10, K = 52, MX = 30;


int n, a[N], L, d[N];
ll dp[N][N][M][3];
void solve(){
  cin >> n >> L;
  for(int i = 1; i <= n; ++i){
    cin >> a[i];
  }
  if(n == 1){
    cout << 1;
    return;
  }
  for(int i = 0; i <= n; ++i){
    for(int j = 0; j <= n; ++j) for(int x = 0; x <= L; ++x) for(int y = 0; y <= 2; ++y) dp[i][j][x][y] = 0;
  }
  sort(a+1, a+1+n);
  for(int i = 1; i <= n; ++i){
    a[i - 1] = a[i];
  }
  for(int i = 1; i <= n; ++i){
    d[i - 1] = a[i] - a[i - 1];
  }
  dp[1][1][2 * d[0]][0] = 1;
  dp[1][1][d[0]][1] = 2;

  for(int i = 1; i < n; ++i){
    for(int j = 1; j <= i + 1; ++j){
      for(int k = 0; k <= L; ++k){
        for(int l = 0; l <= 2; l++) dp[i][j][k][l] %= MOD;
        // l == 0
        if(k + (2*j+2) * d[i] <= L) dp[i + 1][j + 1][k + (2*j+2) * d[i]][0] += dp[i][j][k][0] * (j + 1);
        if(k + (2*j+1) * d[i] <= L) dp[i + 1][j + 1][k + (2*j+1) * d[i]][1] += dp[i][j][k][0] * 2;
        if(k + (2*j-1) * d[i] <= L) dp[i + 1][j][k + (2*j-1) * d[i]][1] += dp[i][j][k][0] * 2;
        // if(i + 1 == 3 && j == 1){
        //   cout << dp[i][j][k][0] << ' ' << k + (2*j-1)*d[i] << ' ' << k << ' ' << " f\n";
        // }
        if(k + (2*j) * d[i] <= L) dp[i + 1][j][k + (2*j) * d[i]][0] += dp[i][j][k][0] * 2 * j;
        if(k + (2*j-2) * d[i] <= L) dp[i + 1][j - 1][k + (2*j-2) * d[i]][0] += dp[i][j][k][0] * (j - 1);

        // l == 1
        if(k + (2*j+1) * d[i] <= L) dp[i + 1][j + 1][k + (2*j+1)*d[i]][1] += dp[i][j][k][1] * j;
        if(k + (2*j) * d[i] <= L) dp[i + 1][j + 1][k + (2*j)*d[i]][2] += dp[i][j][k][1] * 1;
        if(k + (2*j-2) * d[i] <= L && (j > 1 || i == n-1)) dp[i + 1][j][k + (2*j-2)*d[i]][2] += dp[i][j][k][1] * 1;
        if(k + (2*j-1) * d[i] <= L) dp[i + 1][j][k + (2*j - 1)*d[i]][1] += dp[i][j][k][1] * (2*j-1);
        if(k + (2*j-3) * d[i] <= L) dp[i + 1][j - 1][k + (2*j - 3)*d[i]][1] += dp[i][j][k][1] * (j - 1);

        // l == 2
        if(k + (2*j) * d[i] <= L) dp[i + 1][j + 1][k + (2*j) * d[i]][2] += dp[i][j][k][2] * (j - 1);
        if(k + (2*j-2) * d[i] <= L) dp[i + 1][j][k + (2*j-2) * d[i]][2] += dp[i][j][k][2] * (2*j - 2);
        if(k + (2*j-4) * d[i] <= L && (j >= 3 || i == n-1)) dp[i + 1][j - 1][k + (2*j-4) * d[i]][2] += dp[i][j][k][2] * (j - 1);
      }
    }
  }
  ll sum = 0;
  for(int i = 0; i <= L; ++i){
    sum += dp[n][1][i][2];
    sum %= MOD; 
  }
  cout << sum;
}


int main(){
  cin.tie(0); ios::sync_with_stdio(0);
  int tt = 1, aa;
  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
  while(tt--){
    solve();
    en;
  }
  cerr<<"time taken : "<<(float)clock()/CLOCKS_PER_SEC<<" seconds\n";
  return 0;
} 
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...