[CF1178F1]Short Colorful Strip

最新推荐文章于 2021-07-26 23:34:42 发布

OneInDark

最新推荐文章于 2021-07-26 23:34:42 发布

阅读量190

点赞数

分类专栏：数学 # 区间 dp 动态规划

本文链接：https://blog.csdn.net/qq_42101694/article/details/107641326

版权

数学同时被 3 个专栏收录

280 篇文章 1 订阅

订阅专栏

动态规划

220 篇文章 2 订阅

订阅专栏

区间 dp

14 篇文章 0 订阅

订阅专栏

题目

传送门 to luogu

思路

看看题目里的条件……

$\tt ...\; such\; that\; the\; segment$ $a_i,b_i]$ $\tt is\; currently\; painted\; with\; a\; single\; colour\; ...$

为啥一定要求，接下来的涂色全是在纯色格子上进行？翻译一下这句话，不就是 新的染色区间不能包含已有染色区间的边界 吗？

既然已有染色区间是不可逾越的，它就会把整个数轴切开。然后就变成了很多个 独立子问题。然后就可以进行区间 $\tt dp$ 了。

用 $f (l, r)$ 表示，目前 $l$ 到 $r$ 是纯色区间（且没有已经被染好的格子），染色区间不超过该范围，有多少种染色方案。那么，对于最小的颜色，也就是最先染上的颜色，其位置为 $p$ ，我们可以枚举其染色区间，然后进一步将区间切割。写成方程式就是

$f(l,r)=\sum_{l\le i\le p}\;\sum_{p\le j\le r}f(l,i{\tt -}1)\cdot f(i,p{\tt -}1)\cdot f(p{\text +}1,j)\cdot f(j{\text +}1,r)$

好像是 $\mathcal O(n^4)$ 的，实际上只需要一个很小的数学优化就变成了 $\mathcal O(n^3)$ 。

$f(l,r)=\bigg(\sum_{l\le i\le p}f(l,i{\tt -}1)f(i,p{\tt -}1)\bigg)\bigg(\sum_{p\le j\le r}f(p{\text +}1,j)f(j{\text +}1,r)\bigg)$

边界自然是 $f (x, x - 1) = 1$ ，即什么都不用涂，也不能涂。

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
using namespace std;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
template < class T >
void getMax(T&a,T b){ if(a < b) a = b; }
template < class T >
void getMin(T&a,T b){ if(b < a) a = b; }

const int MaxN = 502;
const int Mod = 998244353;
int dp[MaxN][MaxN], c[MaxN];

int main(){
	int n = readint();
	readint(); // 要它何用？
	for(int i=1; i<=n; ++i)
		c[i] = readint();
	for(int i=1; i<=n+1; ++i)
		dp[i][i-1] = dp[i][i] = 1;
	c[0] = Mod;
	for(int len=1; len<n; ++len)
		for(int i=1; i+len<=n; ++i){
			int id = 0, L = 0, R = 0;
			for(int j=i; j<=i+len; ++j)
				if(c[j] < c[id]) id = j;
			for(int j=i; j<=id; ++j)
				L = (L+1ll*dp[i][j-1]*
					dp[j][id-1])%Mod;
			for(int j=id; j<=i+len; ++j)
				R = (R+1ll*dp[id+1][j]*
					dp[j+1][i+len])%Mod;
			dp[i][i+len] = 1ll*L*R%Mod;
		}
	printf("%d\n",dp[1][n]);
	return 0;
}