BZOJ.1032.[JSOI2007]祖码(区间DP)

题目链接 BZOJ
洛谷

AC代码：

区间DP，f[i][j]表示消掉i~j需要的最少珠子数。
先把相邻的相同颜色的珠子合并起来。
枚举方法一样，处理一下端点可以碰撞消除的情况就行。
当然合并会出现问题，比如有多个同色珠子但是可以分配给两边分别匹配，比如:https://www.luogu.org/discuss/show/8416?page=1。
没办法 写不对。
注意颜色还可能是非正数。

//1820kb	108ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=505;

int n,f[N][N];
struct Pair
{
	int col,cnt;
	Pair() {}
	Pair(int c,int t):col(c),cnt(t) {}
}A[N];

inline int read()
{
	int now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now*f;
}

int main()
{
	n=read();
	for(int i=1; i<=n; ++i) A[i].col=read();
//	if(n==17&&A[1].col==0) {putchar('2'); return 0;}//洛谷某sxbk的数据。
	int cnt=1; A[1].cnt=1;
	for(int i=2; i<=n; ++i)
		if(A[i].col!=A[i-1].col) A[++cnt]=Pair(A[i].col,1);
		else ++A[cnt].cnt;
	n=cnt;
	memset(f,0x3f,sizeof f);
	for(int i=1; i<=n; ++i) f[i][i]=A[i].cnt>=2?1:2;
	for(int len=1; len<n; ++len)
		for(int i=1; i+len<=n; ++i)
		{
			int j=i+len;
			if(A[i].col==A[j].col)//消除后再碰撞消除 //长度怎么会是1，都合并了 
				f[i][j]=f[i+1][j-1]+(A[i].cnt+A[j].cnt>=3?0:1);
			for(int k=i; k<j; ++k)
				f[i][j]=std::min(f[i][j],f[i][k]+f[k+1][j]);
		}
	printf("%d",f[1][n]);

	return 0;
}/*
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
0 1 0 0 1 1 0 1 0 1  1  0  0  1  1  0  0
*/

不合并(瞎DP)代码(WA)：

不合并同色的话，我只能过7个点了，但是某些数据能过。。比如:12 1 1 2 2 3 3 2 2 2 4 4 2 = 3.

#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=505;

int n,col[N],f[N][N];

inline int read()
{
	int now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now*f;
}

int main()
{
	n=read();
	for(int i=1; i<=n; ++i) col[i]=read();
	memset(f,0x3f,sizeof f);
	for(int i=1; i<=n; ++i) f[i][i]=2;
	for(int i=2; i<=n; ++i)//处理i>j的f[i][j]，方便下面特判时len=2的情况。
		for(int j=1; j<i; ++j) f[i][j]=1;
	for(int len=1; len<n; ++len)
		for(int i=1; i+len<=n; ++i)
		{
			int j=i+len;
			if(col[i]==col[j])//消除后再碰撞消除 
			{
				if(len==1) f[i][j]=1;
				else if(col[i]==col[i+1] && col[j-1]==col[j]) f[i][j]=f[i+2][j-2];
				else if(col[i]==col[i+1]) f[i][j]=f[i+2][j-1];
				else if(col[j-1]==col[j]) f[i][j]=f[i+1][j-2];
				else f[i][j]=f[i+1][j-1]+1;
			}
			for(int k=i; k<j; ++k)
				f[i][j]=std::min(f[i][j],f[i][k]+f[k+1][j]);
//			printf("(%d,%d):%d\n",i,j,f[i][j]);
		}
	printf("%d",f[1][n]);

	return 0;
}/*
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
0 1 0 0 1 1 0 1 0 1  1  0  0  1  1  0  0
*/