codeforces 1499E DP

       给定两个串 $x$ 和 $y$。定义 $z$ 为从 $x$ 串和 $y$ 串的归并而成的字符串，归并顺序自定义，要求得到的 $z$ 串相邻字符不相同。问你 $x$ 和 $y$ 的所有子串能构成 $z$ 的方案数有多少种。

       定义 $dp[i][j][a][b][c]$，$dp[i][j]$ 表示前面的长度为 $i$ 的 $x$ 串和长度为 $j$ 的 $y$ 串的所有子串的方案数。

       $a$ 代表最后一个字符结尾为 $x$ 串字符还是 $y$ 串字符（$0$ 表示 $x$，$1$ 表示 $y$）。

       $b$ 代表是否存在 $x$ 串字符，（$0$ 表示不存在，$1$ 表示存在）。

       $c$ 代表是否存在 $y$ 串字符，（$0$ 表示不存在，$1$ 表示存在）。

       于是我们存在四种转移情况：

       1.上一个字符为 $x$ 串字符，当前字符为 $y$ 串字符。

       转移后的情况一定同时含有 $x$ 串字符和 $y$ 串字符，且以 $y$ 串字符结尾，所以当前状态为 $dp[i][j][1][1][1]$，而前一个状态可能同时存在 $x$ 串字符和 $y$ 串字符，也可能只有 $x$ 串字符。所以是 $dp[i][j-1][0][1][1]+dp[i][j-1][0][1][0]$。其他情况类似分析即可，不再赘述。

       2.上一个字符为 $y$ 串字符，当前字符为 $x$ 串字符。

       3.上一个字符为 $x$ 串字符，当前字符为 $x$ 串字符。

       4.上一个字符为 $y$ 串字符，当前字符为 $y$ 串字符。

       几个 $tips$：

       1.每加入一个字符时都要额外加入统计，子串方案数加一，也就是把当前字符作为子串来考虑。（对应就是代码的 $dp[i][j][0][1][0]++$ 和 $dp[i][j][1][0][1]++$）

       2.每次统计只统计同时存在 $x$ 串和 $y$ 串字符的子串的方案数，因为题目要求两个子串长度都要大于 $0$。

       3.为什么输入的字符串从 $1$ 开始，统计子串方案时却从 $0$ 开始，因为需要统计只有一个串的字符的情况用于辅助计算，所以需要从 $0$ 开始。

#include<bits/stdc++.h>
using namespace std;
#define IOS ios::sync_with_stdio(false); cin.tie(0); cout.tie(0)
#define endl '\n'
typedef long long LL;
const int maxn = 2e3 + 1;
const LL mod = 998244353;

char x[maxn], y[maxn];
LL dp[maxn][maxn][2][2][2]; // 0 x, 1 y,

int main(){
    IOS;
    LL ans = 0;
    cin >> (x + 1) >> (y + 1);
    x[0] = '?';
    y[0] = '?';
    for(int i = 0; x[i]; i++){
        for(int j = 0; y[j]; j++){
            if(i && x[i] != x[i - 1]){
                for(int k = 0; k < 2; k++){
                    dp[i][j][0][1][k] += dp[i - 1][j][0][1][k];
                }
            }
            if(j && y[j] != y[j - 1]){
                for(int k = 0; k < 2; k++){
                    dp[i][j][1][k][1] += dp[i][j - 1][1][k][1];
                }
            }
            if(i) dp[i][j][0][1][0]++;
            if(j) dp[i][j][1][0][1]++;
            if(i && j && x[i] != y[j]){
                dp[i][j][0][1][1] += dp[i - 1][j][1][1][1] + dp[i - 1][j][1][0][1];
                dp[i][j][1][1][1] += dp[i][j - 1][0][1][1] + dp[i][j - 1][0][1][0];
            }
            for(int a = 0; a < 2; a++)
                for(int b = 0; b < 2; b++)
                    for(int c = 0; c < 2; c++)
                        dp[i][j][a][b][c] %= mod;
            for(int a = 0; a < 2; a++) ans += dp[i][j][a][1][1];
            ans %= mod;
        }
    }
    cout << ans << endl;
}
/*
 2
5 2
1000 2
 */