【基础 DP，但注意写法上的细节】

AcWing 80 场周赛 4727.摆放棋子

算法：动态规划

1. $f(i,j,k,u)$ —— 表示选 $i$ 个黑子，$j$ 个白子，末尾有连续 $k$ 个棋子，$u=0$ 表示末位是黑子，$u=1$ 表示末位是白子。

2. 状态转移方程：若末位只有一个连续棋子需要特判

$$f(i,j,k,0)=f(i-1,j,k-1,0)$$

$$f(i,j,k,1)=f(i,j-1,k-1,1)$$

$$f(i,j,1,0)+=f(i-1,j,k,1)$$

$$f(i,j,1,1)+=f(i,j-1,k,0)$$

时间复杂度 $O(n_1{n}_2\times (k_1+k_2))$

#include <iostream>
#include <algorithm>
#include <cstring>

#define endl '\n'

#define fup(i, a, b) for (int i = a; i <= b; i ++ )
#define fdn(i, a, b) for (int i = a; i >= b; i -- )

using namespace std;

typedef long long LL;

const int N = 110, M = 15, MOD = 1e8;

int n, m, x, y;
int f[N][N][M][2];

void add(int &a, int b)
{
    a = (a + b) % MOD;
}

int main()
{
    cin >> n >> m >> x >> y;
    
    f[1][0][1][0] = f[0][1][1][1] = 1;
    fup(i, 0, n) fup(j, 0, m)
    {
        if (i) 
        {
            fup(k, 2, x) f[i][j][k][0] = f[i - 1][j][k - 1][0];
            fup(k, 1, y) add(f[i][j][1][0], f[i - 1][j][k][1]);
        }
        if (j)
        {
            fup(k, 2, y) f[i][j][k][1] = f[i][j - 1][k - 1][1];
            fup(k, 1, x) add(f[i][j][1][1], f[i][j - 1][k][0]);
        }
    }
    
    int res = 0;
    fup(i, 1, x) add(res, f[n][m][i][0]);
    fup(i, 1, y) add(res, f[n][m][i][1]);
    
    cout << res << endl;
    
    return 0;
}

---

简洁写法

1. $f(i,j,k)$ 表示摆好前 $i$ 个棋子，其中有 $j$ 个白棋，$k=0$ 表示最后一位为黑棋，$k=1$ 表示最后一位为白棋。

2. 状态转移方程：若最后一位是白棋，则枚举砍掉最后多少位后剩下的序列的最后一位是黑棋，这里枚举砍掉 $k$ 位的范围是 $1⩽k⩽min(j,x)$，表示最后砍掉的部分不能超过当前枚举的连续段的上限以及题目所给的连续长度的上限，下面同理；若最后一位是黑棋，则枚举砍掉最后多少位后剩下的序列的最后一位是白棋，这里枚举砍掉 $k$ 位的范围是 $1⩽k⩽min(i-j,y)$。

$$f(i,j,1)+=f(i-k,j-k,0)$$

$$f(i,j,0)+=f(i-k,j,1)$$

时间复杂度 $O((n_1+n_2)\times k_2)$

#include <iostream>
#include <algorithm>
#include <cstring>

#define endl '\n'

#define fup(i, a, b) for (int i = a; i <= b; i ++ )
#define fdn(i, a, b) for (int i = a; i >= b; i -- )

using namespace std;

const int N = 110, M = N << 1, S = 15, MOD = 1e8;

int n, m, x, y;
int f[M][N][S];

int add(int &a, int b) 
{
    return a = (a + b) % MOD;
}

int main()
{
    cin >> n >> m >> x >> y;
    n += m;
    
    f[0][0][0] = f[0][0][1] = 1;
    fup(i, 1, n) fup(j, 0, min(i, m))
    {
        fup(k, 1, min(j, y)) add(f[i][j][1], f[i - k][j - k][0]);
        fup(k, 1, min(i - j, x)) add(f[i][j][0], f[i - k][j][1]);
    }
    
    cout << add(f[n][m][0], f[n][m][1]) << endl;
    
    return 0;
}