結果

問題 No.1533 Don't be Negative!
ユーザー googol_S0googol_S0
提出日時 2021-05-23 18:13:18
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 721 ms / 8,000 ms
コード長 1,484 bytes
コンパイル時間 5,328 ms
コンパイル使用メモリ 268,508 KB
実行使用メモリ 93,836 KB
最終ジャッジ日時 2024-04-20 03:40:41
合計ジャッジ時間 38,102 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 159 ms
29,960 KB
testcase_01 AC 643 ms
81,508 KB
testcase_02 AC 700 ms
93,580 KB
testcase_03 AC 703 ms
93,744 KB
testcase_04 AC 340 ms
51,464 KB
testcase_05 AC 626 ms
81,380 KB
testcase_06 AC 632 ms
81,380 KB
testcase_07 AC 332 ms
51,592 KB
testcase_08 AC 331 ms
51,588 KB
testcase_09 AC 681 ms
93,584 KB
testcase_10 AC 635 ms
81,380 KB
testcase_11 AC 706 ms
93,704 KB
testcase_12 AC 340 ms
51,464 KB
testcase_13 AC 161 ms
30,084 KB
testcase_14 AC 163 ms
29,952 KB
testcase_15 AC 341 ms
51,592 KB
testcase_16 AC 654 ms
81,508 KB
testcase_17 AC 166 ms
30,084 KB
testcase_18 AC 639 ms
81,380 KB
testcase_19 AC 635 ms
81,508 KB
testcase_20 AC 331 ms
51,464 KB
testcase_21 AC 641 ms
81,380 KB
testcase_22 AC 332 ms
51,592 KB
testcase_23 AC 171 ms
30,084 KB
testcase_24 AC 333 ms
51,716 KB
testcase_25 AC 721 ms
93,708 KB
testcase_26 AC 701 ms
93,580 KB
testcase_27 AC 339 ms
51,464 KB
testcase_28 AC 634 ms
81,380 KB
testcase_29 AC 641 ms
81,380 KB
testcase_30 AC 693 ms
93,708 KB
testcase_31 AC 630 ms
81,328 KB
testcase_32 AC 640 ms
81,380 KB
testcase_33 AC 702 ms
93,712 KB
testcase_34 AC 337 ms
51,420 KB
testcase_35 AC 161 ms
30,088 KB
testcase_36 AC 697 ms
93,712 KB
testcase_37 AC 690 ms
93,708 KB
testcase_38 AC 688 ms
93,836 KB
testcase_39 AC 624 ms
81,380 KB
testcase_40 AC 631 ms
81,512 KB
testcase_41 AC 624 ms
81,504 KB
testcase_42 AC 690 ms
93,580 KB
testcase_43 AC 628 ms
81,508 KB
testcase_44 AC 627 ms
81,508 KB
testcase_45 AC 636 ms
81,368 KB
testcase_46 AC 627 ms
81,504 KB
testcase_47 AC 684 ms
93,708 KB
testcase_48 AC 638 ms
81,636 KB
testcase_49 AC 657 ms
81,636 KB
testcase_50 AC 703 ms
93,708 KB
testcase_51 AC 695 ms
93,708 KB
testcase_52 AC 696 ms
93,580 KB
testcase_53 AC 691 ms
93,708 KB
testcase_54 AC 701 ms
93,708 KB
testcase_55 AC 706 ms
93,708 KB
testcase_56 AC 701 ms
93,708 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder;
typedef long long ll;
using mint=modint998244353;
const int log2NMAX=15;
int main(){
  int N,M,K;
  cin>>N>>M>>K;
  vector<vector<mint>> X(log2NMAX+1,vector<mint>(0));
  vector<vector<vector<mint>>> Y(log2NMAX+1,vector<vector<mint>>(0));
  X[0].assign(2*M+1,1);
  X[0][M-K]=0;
  X[0][M+K]=0;
  for(int i=0;i<log2NMAX;i++){
    X[i+1]=convolution(X[i],X[i]);
  }for(int i=0;i<log2NMAX+1;i++){
    Y[i].assign(1<<(log2NMAX+1-i),vector<mint>(1,1));
  }
  Y[log2NMAX][1]=X[log2NMAX];
  vector<int> t(log2NMAX+3,0);
  t[1]=1;
  for(int i=2;i<log2NMAX+3;i++){
    t[i]=t[i-1]*2;
  }
  for(int i=log2NMAX-1;i>=0;i--){
    for(int j=0;j<(int)(Y[i+1].size());j++){
      Y[i][j<<1]=Y[i+1][j];
      Y[i][(j<<1)+1]=convolution(Y[i+1][j],X[i]);
      vector<mint> Z((t[i+1]-1)*2*M+1);
      if(Y[i][j<<1].size()>=Z.size()){
        for(int k=0;k<(int)(Z.size());k++){
          Z[k]=Y[i][j<<1][(Y[i][j<<1].size()>>1)+k-(Z.size()>>1)];
        }
        Y[i][j<<1]=Z;
      }if(Y[i][(j<<1)+1].size()>=Z.size()){
        for(int k=0;k<(int)(Z.size());k++){
          Z[k]=Y[i][(j<<1)+1][(Y[i][(j<<1)+1].size()>>1)+k-(Z.size()>>1)];
        }
        Y[i][(j<<1)+1]=Z;
      }
    }
  }
  mint ANS=0;
  mint P=2*M-1;
  if(K==0){
    P+=1;
  }
  mint V=P.pow(N);
  mint W=P.inv();
  mint U=1;
  for(int i=1;i<=N;i++){
    V*=W;
    U*=P;
    ANS+=(U-Y[0][i][0])*V*(N-i+1);
  }
  cout<<(ANS/2).val()<<"\n";
}
0