题目链接 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
*/