bzoj1084 最大子矩阵【轮廓线dp】

解题思路：

虽说题目数据n$\le$100,k$\le$20，但其实可以用轮廓线O(kn)过。
m=1的情况就是最大M字段和问题，详见：http://blog.csdn.net/cdsszjj/article/details/78270332；
现在说m=2的情况，考虑dp覆盖部分边缘，共有四种情况：

分别记为0,1,2,3。
设f[i][j][k]表示用了i个矩阵，覆盖了j行，边缘情况为0,1,2,3的最优值；
考虑这一行直接并上上一行的所在矩阵：
0可以从0转移；
1可以从1,3转移来；
2可以从2,3转移来；
3可以从3转移来；
由于可以随时新开一个矩阵，所以第j行任何一种情况都可以从1~j-1行的任意一种情况转移来，所以可以维护一个mx[i][j]表示用了i个矩阵覆盖到第j行的最优值，总的dp方程详见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=105;
int n,m,k;
ll ans;
ll a[N],dp[N],mx1[N];
ll b[N][3],f[N][N][4],mx[N][N];

int main()
{
    //freopen("lx.in","r",stdin);
    //freopen("lx.out","w",stdout);
    n=getint(),m=getint(),k=getint();
    if(m==1)
    {
        int i,j;
        for(i=0;i<=n;i++)dp[i]=-1e17;
        for(i=1;i<=n;i++)a[i]=getint();
        for(i=1;i<=k;i++)
        {
            for(j=i;j<=n;j++)
                dp[j]=max(dp[j-1]+a[j],mx1[j-1]+a[j]);
            for(j=0;j<=n;j++)mx1[j]=-1e17;
            for(j=i;j<=n;j++)
                mx1[j]=max(mx1[j-1],dp[j]);
        }
        ans=-1e17;
        for(i=k;i<=n;i++)
            ans=max(ans,dp[i]);
        cout<<ans;
    }
    else 
    {
        int i,j;
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
                b[i][j]=getint();
        for(i=0;i<=k;i++)
            for(j=0;j<=n;j++)
                f[i][j][0]=f[i][j][1]=f[i][j][2]=f[i][j][3]=-1e17;
        for(i=1;i<=k;i++)
        {
            for(j=(i+1)/2;j<=n;j++)
            {
                f[i][j][0]=max(f[i][j][0],f[i][j-1][0]+b[j][1]+b[j][2]);
                f[i][j][0]=max(f[i][j][0],mx[i-1][j-1]+b[j][1]+b[j][2]);
                f[i][j][1]=max(f[i][j][1],mx[i-1][j-1]+b[j][1]);
                f[i][j][1]=max(f[i][j][1],f[i][j-1][1]+b[j][1]);
                f[i][j][1]=max(f[i][j][1],f[i][j-1][3]+b[j][1]);
                f[i][j][2]=max(f[i][j][2],mx[i-1][j-1]+b[j][2]);
                f[i][j][2]=max(f[i][j][2],f[i][j-1][2]+b[j][2]);
                f[i][j][2]=max(f[i][j][2],f[i][j-1][3]+b[j][2]);
                if(i>1)
                {
                    f[i][j][3]=max(f[i][j][3],mx[i-2][j-1]+b[j][1]+b[j][2]);
                    f[i][j][3]=max(f[i][j][3],f[i][j-1][3]+b[j][1]+b[j][2]);
                    f[i][j][3]=max(f[i][j][3],f[i-1][j][1]+b[j][2]);
                    f[i][j][3]=max(f[i][j][3],f[i-1][j][2]+b[j][1]);
                }
            }
            for(j=0;j<=n;j++)mx[i][j]=-1e17;
            for(j=1;j<=n;j++)
            {
                mx[i][j]=max(max(f[i][j][0],f[i][j][1]),max(f[i][j][2],f[i][j][3]));
                mx[i][j]=max(mx[i][j],mx[i][j-1]);
            }
        }
        ans=-1e17;
        for(i=1;i<=n;i++)
            for(j=0;j<4;j++)
                ans=max(ans,f[k][i][j]);
        cout<<ans;
    }
    return 0;
}