[ZOJ3213] Beautiful Meadow

　　插头DP。。。网格图，有障碍，格子上有权值，求总权值最大的简单路径。

　　因为路径的起始点不确定。。所以多开一维表示当前已经有多少个独立插头。。

　　只要不合并相同的联通块，并且已经用了2个独立插头，那就是一条简单路径了。。。

　　需要特判路径上只有一个点的情况。

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #define ll long long
 using namespace std;
 const int modd=,maxzt=;
 struct zs1{
     struct zs{
         int pre;ll too;
     }e[maxzt];int tot,last[modd];
     int f[maxzt];ll zt[maxzt];
     inline int get(ll v){
         int i,x=v%modd;
         for(i=last[x];i&&e[i].too!=v;i=e[i].pre);
         if(i)return i;
         e[++tot].too=v,e[tot].pre=last[x],last[x]=tot,
         f[tot]=,zt[tot]=v;
         return tot;
     }
 }hm[][];

 int i,j,k,n,m;
 bool can[][];int val[][];
 int mp[],id[];bool u[];

 inline void clr(bool now,int num){
     memset(hm[now][num].last,,modd<<);
     hm[now][num].tot=;
 }
 inline void upd(int &a,int b){if(b>a)a=b;}
 inline void decode(ll x){
     for(int i=m;i>=;i--)mp[i]=x&,x>>=;
 }
 inline ll encode(){
     int i,tt=;ll x=;
     memset(u,,);
     for(i=;i<=m;mp[i]=id[mp[i]],x=x<<|mp[i],i++)
         if(!u[mp[i]]&&mp[i]>)u[mp[i]]=,id[mp[i]]=++tt;
     return x;
 }
 inline void dp_blank(int x,int y,bool pre){
     int i,left,up,f;ll zt;bool now=pre^;
     for(int num=;num<=;num++){
 //        printf("    %d,%d %d\n",x,y,num);
         clr(now,num);
         for(i=;i<=hm[pre][num].tot;i++){
             zt=hm[pre][num].zt[i],f=hm[pre][num].f[i];
             if(y==)zt>>=;
             decode(zt);
 //            for(int j=0;j<=m;j++)printf(" %d",mp[j]);printf("    %d    zt:%lld\n",f,zt);
             left=mp[y-],up=mp[y];
             if(left&&up&&left!=up){
                 mp[y-]=mp[y]=;
                 for(int j=;j<=m;j++)if(mp[j]==up){mp[j]=left;break;}
                 upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
             }
             if(!left&&!up){
                 upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f );
                 if(can[x+][y]&&can[x][y+])
                     mp[y-]=mp[y]=,
                     upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
             }
             if((!left)^(!up)){
                 int j=left|up;
                 if(can[x+][y])
                     mp[y-]=j,mp[y]=,
                     upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
                 if(can[x][y+])
                     mp[y-]=,mp[y]=j,
                     upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
             }
         }
         if(num>)
             for(i=;i<=hm[pre][num-].tot;i++){
                 zt=hm[pre][num-].zt[i],f=hm[pre][num-].f[i];
                 if(y==)zt>>=;
                 decode(zt);
 //                for(int j=0;j<=m;j++)printf(" %d",mp[j]);printf("    (dl) %d   zt:%lld\n",f,zt);
                 left=mp[y-],up=mp[y];
                 if(!left&&!up){
                     if(can[x+][y])
                         mp[y-]=,mp[y]=,
                         upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
                     if(can[x][y+])
                         mp[y-]=,mp[y]=,
                         upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );//,printf("   %lld %d\n",encode(),y);
                 }
                 if((!left)^(!up))
                     mp[y-]=mp[y]=,
                     upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
             }
     }
 }
 inline void dp_bar(int x,int y,bool pre){
     int i,left,up,f;ll zt;bool now=pre^;
     for(int num=;num<=;num++){
         clr(now,num);
         for(i=;i<=hm[pre][num].tot;i++){
             zt=hm[pre][num].zt[i],f=hm[pre][num].f[i];
             if(y==)zt>>=;
             decode(zt),
             left=mp[y-],up=mp[y];
             if(!left&&!up)
                 upd( hm[now][num].f[ hm[now][num].get(encode()) ] , f+val[x][y] );
         }
     }
 }

 int ra;char rx;
 inline int read(){
     rx=getchar(),ra=;
     while(rx<''||rx>'')rx=getchar();
     while(rx>=''&&rx<='')ra*=,ra+=rx-,rx=getchar();return ra;
 }
 int main(){
     for(int T=read();T;T--){
         n=read(),m=read();
         memset(can,,sizeof(can));
         int ans=;
         for(i=;i<=n;i++)for(j=;j<=m;j++)val[i][j]=read(),can[i][j]=val[i][j]!=,ans=max(ans,val[i][j]);
         bool pre=,now=;
         clr(pre,),clr(pre,),clr(pre,);
         hm[pre][].f[ hm[pre][].get() ]=;
         for(i=;i<=n;i++)for(j=;j<=m;j++,swap(pre,now))
             if(can[i][j])dp_blank(i,j,pre);else dp_bar(i,j,pre);
         for(i=;i<=hm[pre][].tot;i++)
             upd(ans,hm[pre][].f[i]);
         printf("%d\n",ans);
     }
     return ;
 }