結果
問題 | No.834 Random Walk Trip |
ユーザー |
|
提出日時 | 2019-05-25 16:16:20 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 52 ms / 2,000 ms |
コード長 | 1,663 bytes |
コンパイル時間 | 623 ms |
コンパイル使用メモリ | 75,680 KB |
実行使用メモリ | 15,172 KB |
最終ジャッジ日時 | 2024-09-17 15:05:19 |
合計ジャッジ時間 | 1,549 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
ソースコード
#include <iostream>#include <algorithm>#include <vector>#include <string>using namespace std;const int MOD = 1000000007;struct mint {int n;mint(int n_ = 0) : n(n_) {}};mint operator+(mint a, mint b) { a.n += b.n; if (a.n >= MOD) a.n -= MOD; return a; }mint operator-(mint a, mint b) { a.n -= b.n; if (a.n < 0) a.n += MOD; return a; }mint operator*(mint a, mint b) { return a.n == 0 || b.n == 0 ? 0 : (long long)a.n * b.n % MOD; }mint &operator+=(mint &a, mint b) { return a = a + b; }mint &operator-=(mint &a, mint b) { return a = a - b; }mint &operator*=(mint &a, mint b) { return a = a * b; }ostream &operator<<(ostream &o, mint a) { return o << a.n; }vector<mint> F_{1, 1}, R_{1, 1}, I_{0, 1};void check_fact(int n) {for (int i = I_.size(); i <= n; i++) {I_.push_back(I_[MOD % i] * (MOD - MOD / i));F_.push_back(F_[i - 1] * i);R_.push_back(R_[i - 1] * I_[i]);}}mint I(int n) { check_fact(n); return I_[n]; }mint F(int n) { check_fact(n); return n < 0 ? 0 : F_[n]; }mint R(int n) { check_fact(n); return n < 0 ? 0 : R_[n]; }mint C(int n, int r) { return F(n) * R(n - r) * R(r); }mint P(int n, int r) { return F(n) * R(n - r); }mint H(int n, int r) { return n == 0 ? (r == 0) : C(n + r - 1, r); }int modulo(int x, int y) {x %= y;if (x < 0) x += y;return x;}int main() {int N, M;cin >> N >> M;if (N == 1) {cout << 1 << '\n';return 0;}mint ans;for (int i = -M; i <= M; i++) {if (modulo(i, N*2) == 0 || modulo(i, N*2) == N*2-1) {int d = abs(i);if ((M + d) % 2 == 0) {ans += C(M, (M + d) / 2);}}}cout << ans.n << '\n';}