結果
問題 | No.834 Random Walk Trip |
ユーザー |
![]() |
提出日時 | 2019-02-25 21:15:35 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 19 ms / 2,000 ms |
コード長 | 1,246 bytes |
コンパイル時間 | 486 ms |
コンパイル使用メモリ | 55,644 KB |
実行使用メモリ | 19,072 KB |
最終ジャッジ日時 | 2024-06-24 17:55:34 |
合計ジャッジ時間 | 1,919 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
ソースコード
#include <iostream>#define llint long long#define mod 1000000007using namespace std;llint fact[1000005], fact_inv[1000005];llint modpow(llint a, llint n){if(n == 0) return 1;if(n % 2){return ((a%mod) * (modpow(a, n-1)%mod)) % mod;}else{return modpow((a*a)%mod, n/2) % mod;}}void make_fact(){llint val = 1;fact[0] = 1;for(int i = 1; i < 1000005; i++){val *= i;val %= mod;fact[i] = val;}fact_inv[1000004] = modpow(fact[1000004], mod-2);for(int i = 1000003; i >= 0; i--){fact_inv[i] = fact_inv[i+1] * (i+1) % mod;}}llint comb(llint n, llint k){llint ret = 1;ret *= fact[n];ret *= fact_inv[k], ret %= mod;ret *= fact_inv[n-k], ret %= mod;return ret;}llint n, m;int main(void){cin >> n >> m;make_fact();if(n == 1){cout << 1 << endl;return 0;}llint ans = 0;if(m%2){for(int i = 1; ; i++){if((m-1)/2+i*n > m) break;ans += comb(m, (m-1)/2+i*n), ans %= mod;}for(int i = 0; ; i++){if((m+1)/2+i*n > m) break;ans += comb(m, (m+1)/2+i*n), ans %= mod;}}else{for(int i = 1; ; i++){if(m/2+i*n > m) break;ans += 2 * comb(m, m/2+i*n), ans %= mod;}ans += comb(m, m/2), ans %= mod;}cout << ans << endl;return 0;}