結果

問題 No.834 Random Walk Trip
ユーザー Yumeno UkihashiYumeno Ukihashi
提出日時 2019-05-25 23:55:35
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 162 ms / 2,000 ms
コード長 1,063 bytes
コンパイル時間 2,118 ms
コンパイル使用メモリ 193,068 KB
最終ジャッジ日時 2025-01-07 04:05:13
ジャッジサーバーID
(参考情報)
judge3 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
#define REP(x,y,z) for(int x=y;x<=z;x++)
#define MSET(x,y) memset(x,y,sizeof(x))
#define M 1000005
#define MOD 1000000007
using namespace std;
using LL = long long;
int n,m;
LL fact[M];
LL pw(LL x,LL y) {
    LL res=1, tmp=x;
    while (y) {
        if (y&1) res = res * tmp % MOD;
        tmp = tmp * tmp % MOD;
        y /= 2;
    }
    return res;
}
LL C(int x,int y) {
    if (y>x) return 0;
    return fact[x] * pw(fact[y], MOD-2) % MOD * pw(fact[x-y], MOD-2) % MOD;
}
int main()
{
    fact[0] = 1;
    REP(i,1,M-1) fact[i] = fact[i-1] * i % MOD;
    
    while (~scanf("%d %d",&n,&m)) {
        if (n==1) {
            puts("1");
            continue;
        }

        LL ans = 0;
        int mo = m%2;
        for (int i=mo; abs(i)<=m; i+=2*n) {
            int b = (m-i)/2;
            ans += C(m, b);
            ans %= MOD;
        }
        for (int i=-2*n+mo; abs(i)<=m; i-=2*n) {
            int b = (m-i)/2;
            ans += C(m, b);
            ans %= MOD;
        }
        printf("%lld\n", ans);
    }
    return 0;
}
0