結果

問題 No.1696 Nonnil
コンテスト
ユーザー eSeF
提出日時 2021-03-17 00:01:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,144 bytes
コンパイル時間 2,109 ms
コンパイル使用メモリ 196,660 KB
最終ジャッジ日時 2025-01-19 17:41:08
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample WA * 4
other AC * 2 WA * 37
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define rep(i,n) for(int i=0;i<n;i++)
using namespace std;
typedef long long ll;
const ll inf = 4611686018427387903LL;
const int MOD=998244353;
//ll dp1[2001][2001][2];//縦
//ll dp2[2001][2001][2];//斜め
ll dp[1501][1501][2];
ll mp(ll x,ll n){
  if(n==0)return 1;
  if(n%2==0)return mp((x*x)%MOD,n/2);
  return (x*mp(x,n-1))%MOD;
}
int main()
{
  int N,K,M;
  cin >> N >> K;
  cin >> M;
  vector<int>R(K+1,-1);
  rep(i,M){
    int l,r;
    cin >> l >> r;
    R[r]=max(R[r],l);
  }
  for(int i=K;i>=1;i--){
    bool ok=true;
    for(int j=i-1;j>=0;j--){
      ok&=(R[j]<R[i]);
    }
    if(!ok)R[i]=-1;
  }
  dp[0][0][0]=1;
  for(int j=1;j<=K;j++){
    for(int k=0;k<=K;k++){
      if(R[j]!=-1&&k>0){
        int L=R[j];
        for(int l=0;l<2;l++){
          if(dp[j][k][l]==0)continue;
          for(int m=j-1;m>=L&&k-(j-m)>=0;m--){
            dp[j][k][l]+=dp[m][k-(j-m)][l^1];
            dp[j][k][l]%=MOD;
          }
          for(int m=L-1;m>=0&&k-(j-L)-1>=0;m--){
            dp[j][k][l]+=dp[m][k-(j-L)-1][l^1];
            dp[j][k][l]%=MOD;
          }
        }
      }
    }
  }
  ll ans=0;
  for(int i=0;i<=K;i++){
    for(int j=0;j<=K;j++){
      //printf("i=%d,j=%d,dp=%lld,mp=%lld\n",i,j,(dp[i][j][0]-dp[i][j][1]+MOD)%MOD,mp(K-j,N));
      ans+=(dp[i][j][0]-dp[i][j][1]+MOD)*mp(K-j,N);
      ans%=MOD;
    }
  }
  cout << ans << endl;
  /*dp1[0][0][0]=dp2[0][0][0]=1;
  for(int j=1;j<=K;j++){
    for(int k=0;k<=K;k++){
      ll nx[2]={0,0};
      if(R[j]!=-1&&k>0){
        int L=R[j];
        for(int l=0;l<2;l++){
          nx[l]+=dp2[j-1][k-1][l^1];
          if(k-(j-L)-1>=0)nx[l]-=dp2[L-1][k-(j-L)-1][l^1];

          if(k-(j-L)-1>=0)nx[l]+=dp1[L-1][k-(j-L)-1][l^1];
          nx[l]%=MOD;
          if(nx[l]<0)nx[l]+=MOD;
        }
      }

      for(int l=0;l<2;l++){
        dp1[j][k][l]=(dp1[j-1][k][l]+nx[l])%MOD;
        dp2[j][k][l]=nx[l];
        if(j-1>=0&&k-1>=0)dp2[j][k][l]+=dp2[j-1][k-1][l];
        dp2[j][k][l]%=MOD;
      }
    }
  }
  ll ans=0;
  for(int i=0;i<=K;i++){
    ans+=(dp1[K][i][0]-dp1[K][i][1]+MOD)*mp(K-i,N);
    ans%=MOD;
  }
  cout << ans << endl;*/
  return 0;
}
0