結果

問題 No.3190 Scoring
ユーザー 蜜蜂
提出日時 2025-06-17 11:07:07
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,062 bytes
コンパイル時間 3,148 ms
コンパイル使用メモリ 322,624 KB
実行使用メモリ 22,800 KB
最終ジャッジ日時 2025-06-20 20:54:35
合計ジャッジ時間 14,643 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3 TLE * 1
other -- * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

// g++-14 1.cpp -std=c++20 -O2 -I .
#include <bits/stdc++.h>
using namespace std;

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/tree_policy.hpp>
#include<ext/pb_ds/tag_and_trait.hpp>
using namespace __gnu_pbds;

#include <atcoder/modint>
using namespace atcoder;

using ll = long long;
using ld = long double;
 
using vi = vector<int>;
using vvi = vector<vi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vld = vector<ld>;
using vvld = vector<vld>;
using vst = vector<string>;
using vvst = vector<vst>;

#define fi first
#define se second
#define pb push_back
#define eb emplace_back
#define pq_big(T) priority_queue<T,vector<T>,less<T>>
#define pq_small(T) priority_queue<T,vector<T>,greater<T>>
#define all(a) a.begin(),a.end()
#define rep(i,start,end) for(ll i=start;i<(ll)(end);i++)
#define per(i,start,end) for(ll i=start;i>=(ll)(end);i--)
#define uniq(a) sort(all(a));a.erase(unique(all(a)),a.end())

random_device seed;
mt19937_64 randint(seed());

ll grr(ll mi, ll ma) { // [mi, ma)
    return mi + randint() % (ma - mi);
}

using mint = modint998244353;

const int MAX=510000;
const long long MOD=998244353;

long long fac[MAX],finv[MAX],inv[MAX];

void COMinit(){
  fac[0]=fac[1]=1;
  finv[0]=finv[1]=1;
  inv[1]=1;
  for(int i=2;i<MAX;i++){
    fac[i]=fac[i-1]*i%MOD;
    inv[i]=MOD-inv[MOD%i]*(MOD/i)%MOD;
    finv[i]=finv[i-1]*inv[i]%MOD;
  }
}

long long COM(int n,int k){
  if(n<k) return 0;
  if(n<0||k<0) return 0;
  return fac[n]*(finv[k]*finv[n-k]%MOD)%MOD;
}

int main(){
  std::ios::sync_with_stdio(false);
  std::cin.tie(nullptr);

  COMinit();

  int n,s,m;cin>>n>>s>>m;
  if(n==1){
    cout<<s<<endl;
    return 0;
  }
  mint ans=0;
  rep(x,1,s+1){
    mint q_x=fac[s]*fac[n+s-x-1];
    q_x*=finv[s-x]*finv[n+s-1];
    mint r_x=1-q_x;
    mint iv_r_x=r_x.inv();
    mint q=q_x.pow((m+1)/2);
    mint r=r_x.pow(m-(m+1)/2);
    rep(i,(m+1)/2,m+1){
      ans+=COM(m,i)*q*r;
      q*=q_x;
      r*=iv_r_x;
    }
  }
  ans*=n;
  cout<<ans.val()<<endl;
}
0