結果

問題 No.269 見栄っ張りの募金活動
コンテスト
ユーザー あり
提出日時 2023-02-01 10:25:05
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 11 ms / 5,000 ms
コード長 2,011 bytes
コンパイル時間 733 ms
コンパイル使用メモリ 77,776 KB
実行使用メモリ 20,844 KB
最終ジャッジ日時 2024-07-01 07:26:24
合計ジャッジ時間 1,778 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <string>
#include <set>
#include <map>
#include <algorithm>
#include <numeric>
#include <queue>
#include <cassert>
using namespace std;

const int mod = 1000000007;
struct mint {
  typedef long long ll;
  ll x;
  mint(ll x=0):x((x%mod+mod)%mod){}
  mint operator-() const { return mint(-x);}
  mint& operator+=(const mint a) {
    if ((x += a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator-=(const mint a) {
    if ((x += mod-a.x) >= mod) x -= mod;
    return *this;
  }
  mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;}
  mint operator+(const mint a) const { return mint(*this) += a;}
  mint operator-(const mint a) const { return mint(*this) -= a;}
  mint operator*(const mint a) const { return mint(*this) *= a;}
  mint pow(ll t) const {
    if (!t) return 1;
    mint a = pow(t>>1);
    a *= a;
    if (t&1) a *= *this;
    return a;
  }

  // for prime mod
  mint inv() const { return pow(mod-2);}
  mint& operator/=(const mint a) { return *this *= a.inv();}
  mint operator/(const mint a) const { return mint(*this) /= a;}
};
istream& operator>>(istream& is, mint& a) { return is >> a.x;}
ostream& operator<<(ostream& os, const mint& a) { return os << a.x;}

mint dp[20010][110] = {};
int main() {
  int n, s, k;
  cin >> n >> s >> k;
  s -= k*n*(n-1)/2;
  if (s < 0) {
    cout << 0 << endl;
    return 0;
  }
  
  dp[0][0] = 1;
  for (int i = 0; i <= s; i++) {
    for (int j = 1; j <= n; j++) {
      if (i-j >= 0) {
        assert(i >= 0);
        assert(j >= 0);
        assert(j-1 >= 0);
        assert(i-j >= 0);
        assert(i < 20010);
        assert(j < 110);
        assert(j-1 < 110);
        assert(i-j < 20010);
        dp[i][j] = dp[i][j-1] + dp[i-j][j];
      } 
      else {
        assert(i >= 0);
        assert(j >= 0);
        assert(j-1 >= 0);
        assert(i < 20010);
        assert(j < 110);
        assert(j-1 < 110);

        dp[i][j] = dp[i][j-1];
      }
    }
  }
  cout << dp[s][n] << endl;
}
0