agc035 E Develop （dp）

题意

黑板上有所有整数，每次你可以选中一个$[1,n]$中存在的整数，把他删去并补写上$x-2与x+K$（如果原来不存在的话）。你可以进行这个操作任意次（可能不进行），求最终黑板上数字的可能状态有多少种.
$n,k\le 150$

思路

• 从x到x-2,x+K连有向边，不难发现一个环是永远也消不掉的。
• 反过来想，一个删去集合S是合法的，当且仅当$S$无环。（否则可以按照拓扑序删）
• 现在问题变成了求这个图无环的子集个数。
• 图很特殊。如果K是偶数，问题就变成了不能连续选k/2+1个奇偶性一样的。简单dp。
• K是奇数的话类似这种形式
  

转自这个博客
上图中K=3.

• 发现简单环只能是一个八字形或者三角形。
• 再发现一下，简单环都是长度为K+2的。
• 只需要限制从左边某个点开始，不存在长为K+2并且包含一条横边的的路径即可。
• 设$f[i][j][k]$表示做到第i层，从当前左边那个点开始的最长包含横边的路径长度为j，右边连续选了k个。逐层转移决策即可。
• $O(n^3)$

#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;
	}
}