[bzoj1910] [Ctsc2002] Award 颁奖典礼

　　应该是第一次写这种图形类的DP。。

　　一个“I”可以分成三个矩形。。令f[1..3][i][j][k]表示第几个矩形，下边界为第i行的j~k列，的最大面积。

　　然后就是各种优化啊什么的。。。时间复杂度O(nm²)

　　一开始一个辅助的区间DP写挂然后调了半天TAT

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int maxn=205;
const int inf=1002333;
int f[2][3][maxn][maxn];
int sm[maxn][maxn],mx1[maxn][maxn],mx2[maxn][maxn];
int n,m,pre,now,ans,tmpmx;
 
int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0'||rx>'9')rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra;
}
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a<b?a:b;}
int main(){
    register int i,j,k;
    n=read(),m=read();
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)sm[i][j]=sm[i][j-1]+read();
    for(i=1;i<=m;i++)sm[0][i]=i;
     
    now=1,pre=0;
    for(i=1;i<=2;i++)memset(f[0][i],200,sizeof(f[0][i]));
    for(i=1;i<=n;i++,swap(now,pre)){
        for(j=1;j<m;j++)for(k=j+1;k<=m;k++)
            if(sm[i][k]==sm[i][j-1])
                f[now][0][j][k]=(sm[i-1][k]==sm[i-1][j-1]?f[pre][0][j][k]:0)+k-j+1;
            else f[now][0][j][k]=-inf;
         
         
        for(j=1;j<m;j++)for(k=m,tmpmx=-inf;k>j;k--)
            tmpmx=max(tmpmx,f[pre][0][j][k]),
            mx1[j][k]=max(mx1[j-1][k],tmpmx);
        for(j=2;j<m;j++)for(k=j;k<m;k++)
            if(sm[i][k]==sm[i][j-1]){
                f[now][1][j][k]=f[pre][1][j][k]+k-j+1;
                if(sm[i-1][k+1]==sm[i-1][j-2])
                    f[now][1][j][k]=max(f[now][1][j][k],mx1[j-1][k+1]+k-j+1);
            }
            else f[now][1][j][k]=-inf;
         
         
        for(j=2;j<m;j++)mx2[j][j]=f[pre][1][j][j];
        for(j=1;j<m;j++)
            for(k=2;k<m-j;k++)
            mx2[k][k+j]=max(mx2[k][k+j-1],max(mx2[k+1][k+j],f[pre][1][k][k+j]));
 
        for(j=1;j<m;j++)for(k=j+2;k<=m;k++){
            if(sm[i][k]==sm[i][j-1])
                f[now][2][j][k]=max(f[pre][2][j][k],mx2[j+1][k-1])+k-j+1;
            else f[now][2][j][k]=-inf;
            if(f[now][2][j][k]>ans)ans=f[now][2][j][k];
        }
    }
    printf("%d
",ans);
    return 0;
}