DP专题（一）

写在前面

从接触到DP的第一天起就一直想知道很多转移方程为何这么写状态，每次看到都觉得很奇妙的话那一定是自己太垃圾了。我想是因为自己理论知识不到位所以不能站在一定的高度看待问题，但每次找论文有时竟然深入到最优化理论这种程度，又实在是静不下心来研究那些数学公式理论，所以渐渐明白做题和搞理论不一样，理论是搞那些总的规律，而做题是从具体问题中提取特殊的信息，具体问题具体分析，虽然也有规律，但绝不会是一劳永逸的规律。
这几篇DP试图探讨一些DP的规律和思想

hdu6578 Blank【2019多校赛第一场】

传送带

【思考】

这部分是我的思考过程可以跳过
分析问题，总结一句话就是：区间[l,r]中不同数恰好为x个的方案总数
直接归约问题的话这题就很像一个区间相关的题【什么线段树、莫队，总之自己学的不精就会想一大堆】，但这题有很多附加信息可以使题目变得特殊化，梳理一下：
1）这题并不是给定序列再进行区间上的增删改查，而是一个计数问题
2）这个区间很小，最多是100
3）能填的数很少，只能是1，2，3，4
4）要满足m个区间的条件交集，m最多为100
第一个信息就很显然说明解法不是线段树之类的，而且虽然看到第二个信息第一反应是区间dp，但这题不是求最优解类型的问题，而是计数问题，有点要构造的意思，第三第四两个信息为构造又营造了良好的时空复杂度，所以可以通过高维（n小，就像在对你说：你看我这么小你怎么能不多加几维呢）的计数dp来求解。
怎么构造呢，首先由于dp有很多种优化方式，所以我们能接受的空间范围大概在(空间上限*某个题目数据)之内。本题可以从能填的这4个数入手，也能从m个条件入手，还能从区间的长度入手，如果从m个条件入手，也就是满足前i个条件的数量，会用到容斥原理之类的，但这个条件不好表示，不像概率期望那样加加减减就行了，不好想，遂放弃；再看从4个数入手，有区间的限制条件，但也有些信息是冗余的并不需要知道，举个栗子，1 1 2 1和1 3 1 1，对于区间[2,3]都有贡献，对于区间[3,4]只有前者有贡献，所以数的排列和数的种类不重要，重要的是随便哪个数要至少出现一次，将这个问题归约为随便哪个数最后一次出现的位置在l到r之间，而这“随便哪个数”最多也只有四种情况，所以愉快地开个四位数组吧！
我才不会告诉你我是看了题解之后才凹的思路，不过也大差不差吧

【题解】

记录最后一次出现位置i j k cur，i<j<k<cur（限制顺序可以优化时间），并且cur可以通过滚动数组实现
考虑第当前cur位置放哪一个数f[i][j][k][cur-1]会转移到哪
放第一个数可以转移到f[j][k][cur-1][cur]
放第二个数可以转移到f[i][k][cur-1][cur]
放第三个数可以转移到f[i][j][cur-1][cur]
放第四个数可以转移到f[i][j][k][cur]
由于本题的区间操作没有修改操作，所以对于区间的限制可以将循环变量作为右端点，只对左端点进行判断

【代码】

注意取模的写法，容易写T

#include <bits/stdc++.h>
using namespace std; 
#define mp make_pair
const int MAXN=101;
const int P=998244353;
int T,n,m,l,r,x;
int f[MAXN][MAXN][MAXN][3],ans;
vector<pair<int,int> >d[MAXN]; 

void init(){
	ans=0;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		d[i].clear();
		d[i].push_back(mp(i,1));
	}
	for(int i=0;i<m;i++){
		scanf("%d%d%d",&l,&r,&x);
		d[r].push_back(mp(l,x));
	}
	memset(f,0,sizeof(f));
	f[0][0][0][0]=1;
	for(int cur=1;cur<=n;cur++){
		int c=cur&1;
		for(int i=0;i<=cur;i++)
			for(int j=i;j<=cur;j++)
				for(int k=j;k<=cur;k++)
					f[i][j][k][c]=0;
		for(int i=0;i<=cur;i++)
			for(int j=i;j<=cur;j++)
				for(int k=j;k<=cur;k++){
					(f[j][k][cur-1][c]+=f[i][j][k][c^1])%=P;
	                (f[i][k][cur-1][c]+=f[i][j][k][c^1])%=P;
	                (f[i][j][cur-1][c]+=f[i][j][k][c^1])%=P;
	                (f[i][j][k][c]+=f[i][j][k][c^1])%=P;
				}
		for(int i=0;i<=cur;i++)
			for(int j=i;j<=cur;j++)
				for(int k=j;k<=cur;k++)
					for(auto pi:d[cur]){
						l=pi.first;
						r=cur;
						x=pi.second;
						if((i>=l)+(j>=l)+(k>=l)+(cur>=l)!=x) f[i][j][k][c]=0;
					}
	}
	for(int i=0;i<=n;i++)
		for(int j=i;j<=n;j++)
			for(int k=j;k<=n;k++)
				(ans+=f[i][j][k][n&1])%=P;
	printf("%d\n",ans);
}

int main(){
	//freopen("hdu6578.in","r",stdin);
	//freopen("hdu6578.out","w",stdout);
	scanf("%d",&T);
	while(T--){
		init(); 
	}
	return 0;
}