HDU-6578 Blank

HDU-6578 Blank

联想（口胡）（可跳过）

限制为几种不一样的数—>存每个数最后出现的位置

又是一道神仙DP题。但如果思路正确还是挺简单的。

由于限制是限制一个区间有几个不同的数，所以我们不必知道每个位置具体放的是什么数。于是我们定义$dp[i][j][k][q]$为前i个位置中四种数的最后出现位置为i,j,k,q的方案数这里排序一下使得$i\ge j\ge k\ge q$(因为0的存在，所以是大于等于）。然后枚举第i+1个位置放这四种数中的哪一个数即可。

对于一个状态，我们枚举以i为右端点的限制，对于不合法的状态将其方案数清零处理。(所以我们一开始就把所有的限制加在它的右端点上)

由于内存限制，把i这维滚动。

复杂度$O(n^4)$

#include<cstdio>
#include<cstring>
#define M 105
#define P 998244353
struct E{//正向表存限制 
	int L,d,nx;
}Lim[M];
int head[M],tot;
void Addlim(int a,int b,int d){
	Lim[++tot].L=a;
	Lim[tot].d=d;
	Lim[tot].nx=head[b];
	head[b]=tot;
}
int dp[2][M][M][M];
void Add(int &x,int y){
	x=1ll*(x+y)%P;
}
bool check(int a4,int a3,int a2,int a1){//此处传参将逆序转为正序（a1<=a2<=a3<=a4)（为了好看） 
	for(int i=head[a4];i;i=Lim[i].nx){
		int d=Lim[i].d,L=Lim[i].L;
		if(d==1&&L<=a3)return false;
		if(d==2&&(L<=a2||L>a3))return false;
		if(d==3&&(L<=a1||L>a2))return false;
		if(d==4&&L>a1)return false;
	}
	return true;
}
void Init(){
	memset(head,0,sizeof(head));
	tot=0;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		int n,m;
		scanf("%d%d",&n,&m);
		Init();
		for(int i=1;i<=m;i++){
			int L,R,d;
			scanf("%d%d%d",&L,&R,&d);
			Addlim(L,R,d);
		}
		dp[0][0][0][0]=1;
		int cur=0;
		int ans=0;
		for(int i=0;i<=n;i++){
			for(int j=0;j<=i;j++)
				for(int k=0;k<=j;k++)
					for(int q=0;q<=k;q++)dp[!cur][j][k][q]=0;//手动清空，不然会T(反正我memset就T了) 
			for(int j=0;!j||j<i;j++){
				for(int k=0;!k||k<j;k++){
					for(int q=0;!q||q<k;q++){
						if(!check(i,j,k,q)){dp[cur][j][k][q]=0;continue;}//不合法状态 
						if(!dp[cur][j][k][q])continue;
						if(i==n){Add(ans,dp[cur][j][k][q]);dp[cur][j][k][q]=0;continue;}
						//注意上面的清空可以使下次使用dp数组时不用清空 
						Add(dp[!cur][j][k][q],dp[cur][j][k][q]);
						Add(dp[!cur][i][k][q],dp[cur][j][k][q]);
						Add(dp[!cur][i][j][q],dp[cur][j][k][q]);
						Add(dp[!cur][i][j][k],dp[cur][j][k][q]);
					}
				}
			}
			cur=!cur;
		}
		printf("%d\n",ans);
	}
	return 0;
}