結果

問題 No.3394 Big Binom
コンテスト
ユーザー Mike scarf
提出日時 2025-12-12 15:42:47
言語 C++23
(gcc 13.3.0 + boost 1.89.0)
結果
TLE  
実行時間 -
コード長 1,005 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,986 ms
コンパイル使用メモリ 276,260 KB
実行使用メモリ 35,132 KB
最終ジャッジ日時 2025-12-12 15:43:04
合計ジャッジ時間 17,009 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 19 TLE * 2
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:51:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   51 |     scanf("%lld%lld",&n,&m);
      |     ~~~~~^~~~~~~~~~~~~~~~~~

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll maxn=2000005;
const ll mod=998244353;
ll s[maxn+5],p[maxn+5];
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1)
            ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
ll inv(ll n){
    return qpow(n,mod-2);
}
void init(){
    ll i;
    s[0]=1;
    for(i=1;i<=maxn;i++)
        s[i]=s[i-1]*i%mod;
    p[maxn]=inv(s[maxn]);
    for(i=maxn-1;i>=0;i--)
        p[i]=p[i+1]*(i+1)%mod;
}
ll cal(ll n,ll m){
    if(m<0 or m>n)
        return 0;
    if(m==0 or m==n)
        return 1;
    if(m>n/2)
        m=n-m;
    if(n<=maxn)
        return s[n]*p[m]%mod*p[n-m]%mod;
    ll i,c1=1,c2=1;
    for(i=1;i<=m;i++)
        c1=c1*(n-i+1)%mod,c2=c2*i%mod;
    return c1*inv(c2)%mod;
}
ll luc(ll n,ll m){
    if(m==0)
        return 1;
    return luc(n/mod,m/mod)*cal(n%mod,m%mod)%mod;
}
int main(){
    ll n,m;
    init();
    scanf("%lld%lld",&n,&m);
    printf("%lld\n",luc(n,m));
    return 0;
}
0