結果

問題 No.2512 Mountain Sequences
ユーザー momoyuumomoyuu
提出日時 2023-10-29 12:36:34
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 879 ms / 3,000 ms
コード長 1,475 bytes
コンパイル時間 3,315 ms
コンパイル使用メモリ 257,448 KB
実行使用メモリ 10,552 KB
最終ジャッジ日時 2023-10-29 12:36:55
合計ジャッジ時間 20,966 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 82 ms
7,784 KB
testcase_01 AC 82 ms
7,784 KB
testcase_02 AC 82 ms
7,784 KB
testcase_03 AC 82 ms
7,784 KB
testcase_04 AC 82 ms
7,784 KB
testcase_05 AC 82 ms
7,784 KB
testcase_06 AC 83 ms
7,784 KB
testcase_07 AC 83 ms
7,784 KB
testcase_08 AC 83 ms
7,784 KB
testcase_09 AC 106 ms
7,784 KB
testcase_10 AC 855 ms
10,552 KB
testcase_11 AC 875 ms
10,552 KB
testcase_12 AC 877 ms
10,552 KB
testcase_13 AC 877 ms
10,552 KB
testcase_14 AC 847 ms
10,552 KB
testcase_15 AC 869 ms
10,552 KB
testcase_16 AC 872 ms
10,552 KB
testcase_17 AC 873 ms
10,552 KB
testcase_18 AC 852 ms
10,552 KB
testcase_19 AC 879 ms
10,552 KB
testcase_20 AC 640 ms
10,552 KB
testcase_21 AC 613 ms
10,552 KB
testcase_22 AC 645 ms
10,552 KB
testcase_23 AC 615 ms
10,552 KB
testcase_24 AC 641 ms
10,552 KB
testcase_25 AC 371 ms
10,552 KB
testcase_26 AC 393 ms
10,552 KB
testcase_27 AC 380 ms
10,552 KB
testcase_28 AC 544 ms
10,288 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
using ll = long long;

#include<atcoder/modint>
using mint = atcoder::modint998244353;


const int B = 500;
mint fac[4<<17],ifac[4<<17];
mint comb(int n,int r){
    if(r>n) return 0;
    return fac[n] * ifac[r] * ifac[n-r];
}
int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    
    fac[0] = 1;
    for(int i = 1;i<4<<17;i++) fac[i] = fac[i-1] * i;
    for(int i = 0;i<4<<17;i++) ifac[i] = fac[i].inv();
    int t;
    cin>>t;
    vector<int> n(t),m(t);
    for(int i = 0;i<t;i++) cin>>n[i]>>m[i];
    vector<int> idx(t);
    for(int i = 0;i<t;i++) idx[i] = i;
    sort(idx.begin(),idx.end(),[&](int i,int j){
        if(n[i]/B!=n[j]/B) return n[i]/B < n[j]/B;
        if(n[i]/B%2) return m[i] < m[j];
        return m[i] > m[j];});

    int nn = 1,mm = 1;
    mint tmp = 1;
    vector<mint> ans(t);
    mint inv2 = (998244353+1)/2;
    for(int i = 0;i<t;i++){
        int ni = idx[i];
        while(m[ni]>mm){
            tmp += comb(2*mm,nn-1);
            mm++;
        }
        while(m[ni]<mm){
            tmp -= comb(2*(mm-1),nn-1);
            mm--;
        }
        while(n[ni]>nn){
            mint now = comb(2*mm,nn+1);
            tmp = (now-tmp) * inv2;
            nn++;
        }
        while(n[ni]<nn){
            mint now = comb(2*mm,nn);
            tmp = now - mint(2) * tmp;
            nn--;
        }
        ans[ni] = tmp;
    }
    for(int i = 0;i<t;i++) cout<<ans[i].val()<<endl;
}
0