結果

問題 No.1474 かさまＪ
ユーザー 👑 NachiaNachia
提出日時 2021-04-10 00:25:05
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 1,291 bytes
コンパイル時間 2,882 ms
コンパイル使用メモリ 213,380 KB
実行使用メモリ 9,856 KB
最終ジャッジ日時 2024-06-25 13:10:08
合計ジャッジ時間 5,436 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 3 ms
5,248 KB
testcase_01 AC 3 ms
5,376 KB
testcase_02 AC 3 ms
5,376 KB
testcase_03 AC 3 ms
5,376 KB
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 AC 23 ms
6,808 KB
testcase_08 AC 44 ms
8,448 KB
testcase_09 RE -
testcase_10 RE -
testcase_11 AC 4 ms
5,376 KB
testcase_12 AC 4 ms
5,376 KB
testcase_13 RE -
testcase_14 AC 57 ms
5,820 KB
testcase_15 AC 39 ms
5,544 KB
testcase_16 AC 7 ms
5,376 KB
testcase_17 AC 54 ms
6,452 KB
testcase_18 RE -
testcase_19 WA -
testcase_20 AC 288 ms
9,732 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/modint>
using namespace std;
using ll = long long;
using ull = unsigned long long;
#define rep(i,n) for(int i=0; i<(n); i++)

const ll MOD = 1000000007;
using mll = atcoder::static_modint<1000000007>;

int N,P,Q,L;
int S[40];
vector<vector<mll>> dp;

vector<mll> F,I,iF;
mll C(int n,int r){ return F[n]*iF[r]*iF[n-r]; }

void initC(int Z){
  F.resize(Z+1);
  F[0]=1; for(int i=1; i<=Z; i++) F[i]=F[i-1]*i;
  I.resize(Z+1);
  I[1]=1; for(int i=2; i<=Z; i++) I[i]=-(MOD/i*I[MOD%i]);
  iF.resize(Z+1);
  iF[0]=1; for(int i=1; i<=Z; i++) iF[i]=iF[i-1]*I[i];
}

int main(){
  scanf("%d%d%d%d",&N,&P,&Q,&L);
  rep(i,N) scanf("%d",&S[i]);
  dp.assign(N+1,vector<mll>(P+1,0));
  dp[0][0]=1;
  rep(i,N){
    vector<vector<mll>> buf = dp;
    rep(j,N+1) rep(q,Q) dp[j][q+1] += dp[j][q];
    rep(j,N) for(int q=1; q<=Q; q++){
      int l = max(0,q-S[i]);
      int r = q-1;
      if(l==0) buf[j+1][q] += dp[j][r];
      else buf[j+1][q] += dp[j][r] - dp[j][l-1];
    }
    dp = move(buf);
  }

  initC(40000);

  mll ans=0;

  for(int i=0; i<=N; i++){
    if(i*L>P+Q) break;
    for(int q=0; q<=Q; q++){
      int freeP = P - (i*L-q);
      if(freeP<0) continue;
      ans += dp[i][q] * C(freeP+N-1,N-1);
    }
  }

  printf("%u\n",ans.val());

  return 0;
}
0