BZOJ4350: 括号序列再战猪猪侠【区间DP】

Description

括号序列与猪猪侠又大战了起来。

众所周知,括号序列是一个只有(和)组成的序列,我们称一个括号序列S合法,当且仅当:

1.( )是一个合法的括号序列。

2.若A是合法的括号序列,则(A)是合法的括号序列。

3.若A,B是合法的括号序列,则AB是合法的括号序列。

我们考虑match[i]表示从左往右数第i个左括号所对应的是第几个右括号,现在他得到了一个长度为2n的括号序列,给了你m个信息,第i个信息形如ai,bi,表示match[ai]<match[bi],要你还原这个序列。

但是你发现这个猪猪侠告诉你的信息,可能有多个括号序列合法;甚至有可能告诉你一个不存在合法括号序列的信息!

你最近学了取模运算,你想知道答案对998244353(7172^23+1)取模的结果,这个模数是一个质数。

Input

第一行一个正整数T,T< = 5,表示数据组数。

对于每组数据,第一行一个n,m,n表示有几个左括号,m表示信息数。

接下来m行,每行两个数ai,bi,1< = ai,bi< = n。

Output

对于每组数据,输出一个数表示答案。

Sample Input

5
1 0
5 0
3 2
1 2
2 3
3 2
2 1
2 3
3 3
1 2
2 3
3 1

Sample Output

1
42
1
2
0

HINT

对于前两个点,是卡特兰数的情况。

对于第三个点,合法的情况只可能是 ()()()。

对于第四个点,合法情况可能是 (()()) 或者 (())()

对于第五个点,由于拓扑关系形成了环,显然无解。

对于 100% 的数据,保证 n < = 300


思路

考虑区间DP

(dp_{l,r})表示满足([l,r])的左区间满足条件的方案数

然后你每次考虑在([l+1,r])的个区间中加入l这个括号

有三种情况:

  1. 全部包含后面
  2. 和后面相离
  3. 把后面分成两半

然后发现我们要处理出两个区间分离没有任何冲突的方案数

这个东西可以对match的二维矩阵做一个前缀和sum

然后([l_1,r_1])([l_2,r_2])的冲突个数就是(l_1,r_1)~(l_2, r_2)子矩阵的和


#include<bits/stdc++.h>
using namespace std;
const int N = 310;
const int Mod = 998244353;
int f[N][N], p[N][N], q[N][N], sum[N][N];
int a[N][N], b[N][N];
int n, m;
 
int add(int a, int b) {
  return (a += b) >= Mod ? a - Mod : a;
}
 
int mul(int a, int b) {
  return 1ll * a * b % Mod;
}
 
int calc(int x1, int y1, int x2, int y2) {
  return sum[x2][y2] - sum[x1 - 1][y2] - sum[x2][y1 - 1] + sum[x1 - 1][y1 - 1];
}
 
void solve() {
  scanf("%d %d", &n, &m);
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= n; j++) 
      f[i][j] = sum[i][j] = 0;
  bool bk = 0;
  for (int i = 1; i <= m; i++) {
    int x, y; scanf("%d %d", &x, &y);
    sum[x][y] = 1;
    if (x == y) bk = 1;
  }
  if (bk) {
    printf("0
");
    return;
  }
  for (int i = 1; i <= n; i++) {
    f[i][i] = 1;
    for (int j = 1; j <= n; j++) {
      sum[i][j] += sum[i][j - 1] + sum[i - 1][j] - sum[i - 1][j - 1];
    }
  }
  for (int len = 2; len <= n; len++) {
    for (int l = 1; l + len - 1 <= n; l++) {
      int r = l + len - 1;
      if (!calc(l, l + 1, l, r)) f[l][r] = add(f[l][r], f[l + 1][r]);
      if (!calc(l + 1, l, r, l)) f[l][r] = add(f[l][r], f[l + 1][r]);
      for (int k = l + 1; k <= r - 1; k++) {
        if (!calc(l, l + 1, l, k) && !calc(k + 1, l, r, k))
          f[l][r] = add(f[l][r], mul(f[l + 1][k], f[k + 1][r])); 
      } 
    } 
  }
  printf("%d
", f[1][n]);
}
int main() {
#ifdef dream_maker
  freopen("input.txt", "r", stdin);
  freopen("output.txt", "w", stdout);
#endif
  int T; scanf("%d", &T);
  while (T--) solve();
  return 0;
}
原文地址:https://www.cnblogs.com/dream-maker-yk/p/9922012.html