题意
黑板上有所有整数,每次你可以选中一个
[
1
,
n
]
[1,n]
[1,n]中存在的整数,把他删去并补写上
x
−
2
与
x
+
K
x-2与x+K
x−2与x+K(如果原来不存在的话)。你可以进行这个操作任意次(可能不进行),求最终黑板上数字的可能状态有多少种.
n
,
k
≤
150
n,k\leq150
n,k≤150
思路
- 从x到x-2,x+K连有向边,不难发现一个环是永远也消不掉的。
- 反过来想,一个删去集合S是合法的,当且仅当 S S S无环。(否则可以按照拓扑序删)
- 现在问题变成了求这个图无环的子集个数。
- 图很特殊。如果K是偶数,问题就变成了不能连续选k/2+1个奇偶性一样的。简单dp。
- K是奇数的话类似这种形式
转自这个博客
上图中K=3.
- 发现简单环只能是一个八字形或者三角形。
- 再发现一下,简单环都是长度为K+2的。
- 只需要限制从左边某个点开始,不存在长为K+2并且包含一条横边的的路径即可。
- 设 f [ i ] [ j ] [ k ] f[i][j][k] f[i][j][k]表示做到第i层,从当前左边那个点开始的最长包含横边的路径长度为j,右边连续选了k个。逐层转移决策即可。
- O ( n 3 ) O(n^3) O(n3)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 160;
int mo, n, k;
int f[N * 2][N][N], o;
#define add(x, y) ((x) = ((x) + (y)) % mo)
ll solve2() {
f[0][0][0] = 1;
for(int i = 1; i <= n; i++) {
for(int x = 0; x < i; x++) {
for(int y = 0; y < i; y++) if(f[i-1][x][y]) {
for(int o = 0; o < 2; o++) {
int _x = x, _y = y;
if (i & 1) {
_x = o ? x + 1 : 0;
if (_x * 2 > k) continue;
} else {
_y = o ? y + 1 : 0;
if (_y * 2 > k) continue;
}
add(f[i][_x][_y], f[i - 1][x][y]);
}
}
}
}
ll ans = 0;
for(int i = 0; i <= n; i++) {
for(int j = 0; j <= n; j++) add(ans, f[n][i][j]);
}
return ans;
}
ll solve1() {
f[o][0][0] = 1;
int h = (k - 1) / 2, mid = n / 2 - h, t = (n + 1) / 2 - mid;
for(int i = 1; i <= h; i++) { //only even
memset(f[1 - o], 0, sizeof f[1 - o]);
for(int j = 0; j < i; j++) if (f[o][0][j])
for(int oo = 0; oo < 2; oo++) {
add(f[1 - o][0][oo ? j + 1 : 0], f[o][0][j]);
}
o = 1 - o;
}
for(int i = 1; i <= mid + t; i++) { //both
memset(f[1 - o], 0, sizeof f[1 - o]);
for(int j = 0; j <= k + 1; j++) {
for(int w = 0; w * 2 <= n; w++) if (f[o][j][w]) {
for(int o1 = 0; o1 < 2; o1++) {
for(int o2 = 0; o2 < (i > mid ? 1 : 2); o2++) {
int _w = o2 ? w + 1 : 0;
int _j = 0;
if (o1 == 0) _j = 0;
else {
if (j) _j = j + 1;
if (o2) _j = max(_j, _w + 1);
}
if (_j <= k + 1) {
add(f[1 - o][_j][_w], f[o][j][w]);
}
}
}
}
}
o = 1 - o;
}
ll ans = 0;
for(int i = 0; i <= k + 1; i++) {
for(int j = 0; j <= n; j++) {
add(ans, f[o][i][j]);
}
}
return ans;
}
int main() {
freopen("e.in", "r", stdin);
cin >> n >> k >> mo;
if (k & 1) {
cout << solve1() << endl;
} else {
cout << solve2() << endl;
}
}