結果
問題 |
No.1011 Infinite Stairs
|
ユーザー |
![]() |
提出日時 | 2020-03-20 21:59:37 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 421 ms / 2,000 ms |
コード長 | 2,238 bytes |
コンパイル時間 | 1,726 ms |
コンパイル使用メモリ | 173,080 KB |
実行使用メモリ | 428,416 KB |
最終ジャッジ日時 | 2024-07-17 11:49:47 |
合計ジャッジ時間 | 6,763 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 24 |
ソースコード
#include <bits/stdc++.h> #define rep(X,N) for(ll X = 0; X < (N); X++) #define PI (acos(-1.0)) #define pback push_back #define mpair make_pair #define MODN 1000000007 #define ALL(V) (V).begin(),(V).end() #define INT_MAX_HALF (INT_MAX / 2) #define EPS (1e-10) using namespace std; typedef long long ll; class ModInt{ public: long long n; int p; ModInt(){ } //値と割る値を受け取って、オブジェクトを返す ModInt(long long _n, int _p){ p = _p; n = _n % _p; //負の値を受け取ったときには正の値に戻す if(n < 0){ n += p; } } ModInt operator+(ModInt a){ ModInt tmp = ModInt(*this); tmp.n += a.n; if(tmp.n > tmp.p){ tmp.n -= tmp.p; } return tmp; } ModInt operator-(ModInt a){ ModInt tmp = ModInt(*this); tmp.n -= a.n; if(tmp.n < 0){ tmp.n += tmp.p; } return tmp; } ModInt operator*(ModInt a){ ModInt tmp = ModInt(*this); tmp.n = (tmp.n * a.n) % tmp.p; return tmp; } int modpow(int b, int k){ if(k == 0){ return 1; }else if(k % 2 == 1){ long long tmp = b; tmp = tmp * modpow(b, k - 1) % p; return tmp; }else{ long long tmp = modpow(b, k / 2); return tmp * tmp % p; } } int modinv(int b){ return modpow(b, p - 2); } ModInt operator/(ModInt a){ ModInt tmp = ModInt(*this); tmp.n = tmp.n * modinv(a.n) % tmp.p; return tmp; } }; int main(){ int n,d,k; cin >> n >> d >> k; vector<vector<ModInt>> dp(n + 1, vector<ModInt>(300 * n + 2, ModInt(0, MODN))); dp[0][0] = ModInt(1, MODN); for(int i = 0; i < n; i++){ for(int j = 0; j < 300 * i + 1; j++){ dp[i + 1][j + 1] = dp[i + 1][j + 1] + dp[i][j]; dp[i + 1][j + d + 1] = dp[i + 1][j + d + 1] - dp[i][j]; } for(int j = 1; j <= 300 * n + 1; j++){ dp[i + 1][j] = dp[i + 1][j] + dp[i + 1][j - 1]; } } cout << dp[n][k].n << endl; }