BZOJ4565 HAOI2016字符合并（区间dp+状压dp）

　　设f[i][j][k]为将i~j的字符最终合并成k的答案。转移时只考虑最后一个字符是由哪段后缀合成的。如果最后合成为一个字符特殊转移一下。

　　复杂度看起来是O(n³2^k)，实际常数极小达到O(玄学)。

　　upd：突然发现根本没在bzoj上交。bzoj的数据输入中没有空格。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
    while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
    return x*f;
}
#define N 310
int n,m,a[N],c[N],w[N],len[N];
long long f[N][N][N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4565.in","r",stdin);
    freopen("bzoj4565.out","w",stdout);
    const char LL[]="%I64d
";
#else
    const char LL[]="%lld
";
#endif
    n=read(),m=read();
    memset(f,200,sizeof(f));
    for (int i=1;i<=n;i++) a[i]=read(),f[i][i][a[i]]=0,len[i]=(i-1)%(m-1)+1;
    for (int i=0;i<(1<<m);i++) c[i]=read(),w[i]=read();
    for (int k=2;k<=n;k++)
        for (int i=1;i<=n-k+1;i++)
        {
            int j=i+k-1;
            for (int d=j-1;d>=i;d-=m-1)
                if (len[k]!=1)
                    for (int v=0;v<(1<<len[k]);v++)
                        f[i][j][v]=max(f[i][j][v],f[i][d][v>>1]+f[d+1][j][v&1]);
                else 
                    for (int u=0;u<(1<<m-1);u++)
                    f[i][j][c[u<<1]]=max(f[i][j][c[u<<1]],f[i][d][u]+f[d+1][j][0]+w[u<<1]),
                    f[i][j][c[u<<1|1]]=max(f[i][j][c[u<<1|1]],f[i][d][u]+f[d+1][j][1]+w[u<<1|1]);
        }
    for (int i=1;i<(1<<m);i++) f[1][n][0]=max(f[1][n][0],f[1][n][i]);
    cout<<f[1][n][0];
    return 0;
}