NOI 2013 书法家

http://uoj.ac/problem/125

我真是日狗了。。。。。。

果然还是没有耐心读题，搞到读题读错了2个地方，结果调试了半天。。。。。。

言归正传。

动态规划。

这种题目很常见。

我们发现竖着做比较麻烦，那么可以横着做。

打竖将"NOI“分成11种类型。

F[i][j][k][l]表示第i列的涂色部分的上边界为j,下边界为k,类型为l。

然后转移。

从2类型->2类型有点麻烦，我们可以借组一个辅助数组b,见程序注释部分。

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<fstream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<utility>
#include<set>
#include<bitset>
#include<vector>
#include<functional>
#include<deque>
#include<cctype>
#include<climits>
#include<complex>
//#include<bits/stdc++.h>适用于CF,UOJ，但不适用于poj

using namespace std;

typedef long long LL;
typedef double DB;
typedef pair<int,int> PII;
typedef complex<DB> CP;

#define mmst(a,v) memset(a,v,sizeof(a))
#define mmcy(a,b) memcpy(a,b,sizeof(a))
#define re(i,a,b)  for(i=a;i<=b;i++)
#define red(i,a,b) for(i=a;i>=b;i--)
#define fi first
#define se second
#define m_p(a,b) make_pair(a,b)
#define SF scanf
#define PF printf
#define two(k) (1<<(k))

template<class T>inline T sqr(T x){return x*x;}
template<class T>inline void upmin(T &t,T tmp){if(t>tmp)t=tmp;}
template<class T>inline void upmax(T &t,T tmp){if(t<tmp)t=tmp;}

const DB EPS=1e-;
inline int sgn(DB x){if(abs(x)<EPS)return ;return(x>)?:-;}
const DB Pi=acos(-1.0);

inline void clear(vector<int> *A,int a,int b){int i,j;A->clear();re(i,,a)re(j,,b)A[i].push_back();}

inline int gint()
  {
        int res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
    }
inline LL gll()
  {
      LL res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
  }

const int maxN=;
const int maxM=;
const int INF=;

int N,M;
int mp[maxN+][maxM+];
int sum[maxM+][maxN+];
int x,y,F[][maxN+][maxN+][];
int ans;

inline int ask(int r,int x,int y){return sum[r][y]-sum[r][x-];}

int b[maxN+];

inline void clear(int x){int j,k,l;re(j,,N+)re(k,,N+)re(l,,+) F[x][j][k][l]=-INF;}

int main()
  {
      freopen("penman.in","r",stdin);
      freopen("penman.out","w",stdout);
      int i,j,k;
      N=gint();M=gint();
      re(i,,N)re(j,,M)mp[i][j]=gint();
      re(i,,M)re(j,,N)sum[i][j]=sum[i][j-]+mp[j][i];
      x=;y=;
      clear(y);
      re(j,,N)re(k,j,N)F[y][j][k][]=ask(,j,k);
      ans=-INF;
      re(i,,M)
        {
            int t;

              swap(x,y);
              clear(y);

            //
            re(j,,N)re(k,j,N)
              {
                upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                    upmax(F[y][j][k][],ask(i,j,k));
                  }

                //
                re(j,,N)
                  {
                      t=F[x][j][N][];
                      red(k,N-,j)
                        {
                            upmax(F[y][j][k][],t+ask(i,j,k));
                            upmax(t,F[x][j][k][]);
                        }
                  }
                /*
                (1)
                l=j
                r=j..k
                */
                re(j,,N)
                  {
                      t=-INF;
                      re(k,j,N)
                        {
                            upmax(t,F[x][j][k][]);
                            upmax(F[y][j][k][],t+ask(i,j,k));
                        }
                  }
                /*
                (2)
                re(l,1,j-1)
                  re(r,j-1,k)
                    F[i-1][l][r][2]
                b[r]表示F[i-1][1..j-1][r][2]的最大值,b[r]随着j的递增的递增
                求b[j-1...k]的最大值
                */
                re(j,,N+)b[j]=-INF;
                re(j,,N)
                  {
                      t=-INF;
                      re(k,j-,N)
                        {
                            upmax(b[k],F[x][j-][k][]);
                            upmax(t,b[k]);
                            if(k>=j) upmax(F[y][j][k][],t+ask(i,j,k));
                        }
                  }
                /*int l,r;
                re(j,1,N)re(k,j,N)re(l,j,j)re(r,j,k)upmax(F[y][j][k][2],F[x][l][r][2]+ask(i,j,k));
                re(j,1,N)re(k,j,N)re(l,1,j-1)re(r,j-1,k)upmax(F[y][j][k][2],F[x][l][r][2]+ask(i,j,k));*/

                //
                re(j,,N)re(k,j,N)upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                re(k,,N)
                  {
                      t=F[x][k][k][];
                      red(j,k-,)
                        {
                            upmax(F[y][j][k][],t+ask(i,j,k));
                            upmax(t,F[x][j][k][]);
                        }
                  }

                //
                re(j,,N)re(k,j,N)upmax(F[y][j][k][],F[x][j][k][]);
                t=-INF;
                re(j,,N)
                  re(k,j,N)
                      upmax(t,F[x][j][k][]);
                re(j,,N)re(k,j,N)upmax(F[y][j][k][],t);

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                    }

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                  }

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                  }

                //
                re(j,,N)re(k,j,N)upmax(F[y][j][k][],F[x][j][k][]);
                t=-INF;
                re(j,,N)re(k,j+,N)upmax(t,F[x][j][k][]);
                re(j,,N)re(k,j,N)upmax(F[y][j][k][],t);

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                  }

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,k));
                  }

                //
                re(j,,N)re(k,j+,N)
                  {
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                      upmax(F[y][j][k][],F[x][j][k][]+ask(i,j,j)+ask(i,k,k));
                      upmax(ans,F[y][j][k][]);
                  }

            }
        cout<<ans<<endl;
        return ;
  }