*51nod 1409

https://blog.csdn.net/stay_accept/article/details/81476358

不懂啊

#include <map>
#include <queue>
#include <string>
#include <math.h>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string.h>
#include <algorithm>
using namespace std;
inline bool scan_d(int &num)  {
    char in;
    bool IsN=false;
    in=getchar();
    if(in==EOF) return false;
    while(in!='-'&&(in<''||in>'')) in=getchar();
    if(in=='-') {
        IsN=true;
        num=;
    } else num=in-'';
    while(in=getchar(),in>=''&&in<='') {
        num*=,num+=in-'';
    }
    if(IsN) num=-num;
    return true;
}
const int siz=;
int a[siz][siz];
long long g[siz],dp[][siz][];
int main() {
    long long ans;
    int n,m,i,j,k,x,y,op;
    while(scanf("%d%d",&n,&m)!=EOF) {           //dp[i][j][k]表示到dii列第j行向上走还是向下走
        for(i=; i<=n; i++)
            for(j=; j<=m; j++)
                scan_d(a[i][j]);
        for(i=; i<=n; i++)
            dp[][i][]=dp[][i][]=;
        for(i=; i<=m; i++) {
            x=i&,y=-(i&);
            for(j=; j<=n; j++)
                dp[x][j][]=dp[x][j][]=-;
            dp[x][][]=dp[x][][]=;
            dp[x][n+][]=dp[x][n+][]=;
            for(j=; j<=n; j++) {
                if(a[j][i]==-)
                    continue;
                for(k=; k<=; k++) {
                    if(dp[y][j][k]!=-)
                        dp[x][j][]=max(dp[x][j][],dp[y][j][k]+a[j][i]);
                    if(dp[y][j][k]!=-)
                        dp[x][j][]=max(dp[x][j][],dp[y][j][k]+a[j][i]);
                }
            }                                   //先算出从前一列过来的得分
            for(j=; j<=n; j++) {
                if(dp[x][j-][]!=-&&a[j][i]!=-)
                    dp[x][j][]=max(dp[x][j][],dp[x][j-][]+a[j][i]);
            }                                   //然后本列进行转移
            for(j=; j<=n; j++)
                g[j]=-;
            g[]=;
            if(dp[x][n][]!=-) {               //算出经过传送门的得分
                op=n;                           //想下走就找尽可能往下更新的点
                while(op>) {
                    if(dp[x][op][]==dp[x][op-][]+a[op][i]&&dp[x][op-][]!=-)
                        op--;
                    else
                        break;
                }
                for(j=; j<op; j++) {
                    if(a[j][i]==-)
                        break;
                    g[j]=g[j-]+a[j][i];
                }
            }                                   //往上走跟往下走是一样的思路
            for(j=; j<n; j++)
                dp[x][j][]=max(dp[x][j][],g[j]);
            for(j=n-; j>=; j--) {
                if(dp[x][j+][]!=-&&a[j][i]!=-)
                    dp[x][j][]=max(dp[x][j][],dp[x][j+][]+a[j][i]);
            }
            for(j=; j<=n; j++)
                g[j]=-;
            g[n+]=;
            if(dp[x][][]!=-) {
                op=;
                while(op<n) {
                    if(dp[x][op][]==dp[x][op+][]+a[op][i]&&dp[x][op+][]!=-)
                        op++;
                    else
                        break;
                }
                for(j=n; j>op; j--) {
                    if(a[j][i]==-)
                        break;
                    g[j]=g[j+]+a[j][i];
                }
            }
            for(j=n; j>; j--)
                dp[x][j][]=max(dp[x][j][],g[j]);
        }
        ans=-;
        for(i=; i<=n; i++)
            ans=max(ans,max(dp[m&][i][],dp[m&][i][]));
        printf("%I64d\n",ans);
    }
    return ;
}