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
*/

转载于:https://www.cnblogs.com/SovietPower/p/9011339.html