【Ybt OJ】[动态规划 第2章] 区间DP [后半章]

$「$动态规划$」$第$2$章 区间$DP$$($后$2$题$)$
目录：

C.消除木块
D.棋盘分割

$C.$ 例题$3$ 消除木块

洛谷$link$

分析：

区间$dp$
$co{l}_k$ 表示$k$区间的颜色 $le{n}_k$表示$k$区间的长度
先预处理这些
$f_{l,r,k}$表示 消除$[l,r]$区间 在$r$后有$k$个与$r$颜色相同方块的最大分数值
首先考虑直接消除 $[l,r]$区间 那就是$f_{l,r-1,k}+$当前$r$的分数

其次就是区间$dp$ 区间分别为$f_{l,i,len[r]+k},f_{i+1,r-1,0}$
但是要颜色相同 这样就可以消除$[i+1,r-1]$ 再计算$[l,i]$的
$l=r$ 就直接计算当前的分数
记搜实现

CODE：

#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
const int N=255;
int T,n,col[N],len[N],f[N][N][N<<1],r;
int dfs(int l,int r,int k)
{
	if(f[l][r][k]) return f[l][r][k];
	if(l==r) return (len[r]+k)*(len[r]+k);
	f[l][r][k]=dfs(l,r-1,0)+(len[r]+k)*(len[r]+k);
	for(int i=l;i<r-1;i++)
		if(col[i]==col[r])
			f[l][r][k]=max(f[l][r][k],dfs(l,i,len[r]+k)+dfs(i+1,r-1,0));
	return f[l][r][k];
} 
int main(){
	scanf("%d",&T);
	int E=T;
	while(T--)
	{
		scanf("%d",&n);
		int ans=0;
		memset(len,0,sizeof(len));
		memset(col,0,sizeof(col));
		memset(f,0,sizeof(f));
		int k=0;r=-1;
		for(int i=1;i<=n;i++)
		{
			int x;
			scanf("%d",&x);
			if(x==r) len[k]++;
			else
			{
				col[++k]=x;
				len[k]++;
			}
			r=x;
		}
		ans=dfs(1,k,0);
		printf("Case %d: %d\n",E-T,ans);
	}
	
	return 0;
}

$D.$ 例题$4$ 棋盘分割

$ybtoj$上的题面有问题 先放$luogu$的


洛谷$link$

分析：

$\sqrt{\frac{{\sum }_{i=1}^n(x_i-\bar{x})}{n}}$最小 那就是$\frac{{\sum }_{i=1}^n(x_i-\bar{x})}{n}$最小
把它展开 就是$\frac{{\sum }_{i=1}^n(x_i^2-2x_i^2\bar{x}+{\bar{x}}^2)}{n}=\frac{{\sum }_{i=1}^n{x}_i^2}{n}-{\bar{x}}^2$
$\bar{x}=\frac{{\sum }_{i=1}^8{\sum }_{j=1}^8{a}_{i,j}}{n}$
所以${\bar{x}}^2$是定值 那满足${\sum }_{i=1}^n{x}_i^2$最小就好了

后面的就是$2$维区间$dp$ 分别划分$[x1,x2]$ $[y1,y2]$

$f_{x1,y1,x2,y2,k}$表示左上角为$(x1,y1)$ 右下角为$(x2,y2)$的矩阵划分成$k$个子矩阵的最小${\sigma }^2$
需要$2$维前缀和 $s_{i,j}={\sum }_{k=1}^i{\sum }_{l=1}^j{a}_{k,l}$
每次转移需要计算子矩阵和 令函数$val(x1,y1,x2,y2)$表示左上角$(x1,y1)$ 右下角$(x2,y2)$的矩阵和

因为 a n s = m i n { ∑ i = 1 n x i 2 n − x ˉ 2 } ans=min\{\frac{\sum_{i=1}^nx_i^2}{n}-\bar{x}^2\} ans=min{n∑i=1n​xi2​​−xˉ2}
做完 m i n { f x 1 , y 1 , x 2 , y 2 , k } min\{f_{x1,y1,x2,y2,k}\} min{fx1,y1,x2,y2,k​}后 $ans$就等于$\sqrt{f_{1,1,8,8,n}-{\bar{x}}^2}$

数据只有$8\times 8,n<=15$ 就直接$5$循环$dp$了

CODE：

#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
const int N=1e1-1;
int n,a[N][N];
double f[N][N][N][N][N<<1];
double val(int x1,int y1,int x2,int y2)
{
	double res=a[x2][y2]-a[x2][y1-1]-a[x1-1][y2]+a[x1-1][y1-1];
	return res*res/n;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=8;i++)
		for(int j=1;j<=8;j++)
		{
			scanf("%d",&a[i][j]);
			a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
		}
	double X=(double)a[8][8]/n;
	for(int k=1;k<=n;k++)
	for(int x1=1;x1<=8;x1++)
	for(int y1=1;y1<=8;y1++)
	for(int x2=1;x2<=8;x2++)
		for(int y2=1;y2<=8;y2++)
		{
			f[x1][y1][x2][y2][k]=1e9;
			if(k==1)
			{
				f[x1][y1][x2][y2][k]=val(x1,y1,x2,y2);
				continue;
			}
			for(int i=x1;i<x2;i++)
			{
				f[x1][y1][x2][y2][k]=min(f[x1][y1][x2][y2][k],f[x1][y1][i][y2][k-1]+val(i+1,y1,x2,y2));
				f[x1][y1][x2][y2][k]=min(f[x1][y1][x2][y2][k],f[i+1][y1][x2][y2][k-1]+val(x1,y1,i,y2));
			}
			for(int i=y1;i<y2;i++)
			{
				f[x1][y1][x2][y2][k]=min(f[x1][y1][x2][y2][k],f[x1][y1][x2][i][k-1]+val(x1,i+1,x2,y2));
				f[x1][y1][x2][y2][k]=min(f[x1][y1][x2][y2][k],f[x1][i+1][x2][y2][k-1]+val(x1,y1,x2,i));
			}
		}	
	printf("%.3lf",sqrt(f[1][1][8][8][n]-X*X));
	
	return 0;
}