結果

問題 No.3213 depth max K
コンテスト
ユーザー alcea
提出日時 2025-07-25 22:56:25
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 11 ms / 2,000 ms
コード長 2,204 bytes
コンパイル時間 1,543 ms
コンパイル使用メモリ 195,488 KB
実行使用メモリ 7,720 KB
最終ジャッジ日時 2025-07-25 22:56:29
合計ジャッジ時間 3,400 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
template<long long mod>
struct modint{
  long long v;
  modint(long long x=0):v((x%=mod)<0?x+mod:x){}
  constexpr long long val()noexcept{return v;}
  constexpr modint& operator=(const modint o){v=o.v;return *this;}
  constexpr modint& operator+=(const modint o)noexcept{if((v+=o.v)>=mod)v-=mod;return *this;}
  constexpr modint& operator-=(const modint o)noexcept{if((v-=o.v)<0)v+=mod;return *this;}
  constexpr modint& operator*=(const modint o)noexcept{v=v*o.v%mod;return *this;}
  constexpr modint& operator/=(modint o)noexcept{return *this*=o.inv();}
  constexpr modint pow(long long n)noexcept{
    modint ret=1,base=*this;
    while(n>0){if(n&1)ret*=base;base*=base;n/=2;}
    return ret;
  }
  constexpr modint inv()noexcept{return pow(mod-2);}
  constexpr modint operator+(const modint o)noexcept{return modint(*this)+=o;}
  constexpr modint operator+()noexcept{return *this;}
  constexpr modint operator-(const modint o)noexcept{return modint(*this)-=o;}
  constexpr modint operator-()noexcept{return modint()-*this;}
  constexpr modint operator*(const modint o)noexcept{return modint(*this)*=o;}
  constexpr modint operator/(const modint o)noexcept{return modint(*this)/=o;}
  constexpr modint& operator++()noexcept{return *this+=1;}
  constexpr modint operator++(int)noexcept{return (*this+=1)-1;}
  constexpr modint& operator--()noexcept{return *this-=1;}
  constexpr modint operator--(int)noexcept{return (*this-=1)+1;}
  constexpr bool operator==(const modint o)noexcept{return v==o.v;}
  constexpr bool operator!=(const modint o)noexcept{return v!=o.v;}
  friend ostream& operator<<(ostream& os,const modint o)noexcept{os<<o.v;return os;}
};
using mint=modint<998244353>;
const int m=50000;
mint fact[m],ifact[m];
mint comb(int n,int r){
  if(n<0||r<0||n<r) return 0;
  else return fact[n]*ifact[r]*ifact[n-r];
}
mint f(int n,int k){
  mint ans=0;
  for(int i=-n;i<=n;i++){
    ans+=comb(2*n,n+i*(k+2))-comb(2*n,n-1+i*(k+2));
  }
  return ans;
}
int main(){
  fact[0]=ifact[0]=1;
  for(int i=1;i<m;i++){
    fact[i]=fact[i-1]*i;
    ifact[i]=mint(1)/fact[i];
  }
  int n,k;
  cin>>n>>k;
  //cerr<<f(n,0)<<endl;
  cout<<f(n,k)-f(n,k-1)<<endl;
}
0