1. 程式人生 > >bzoj 1187: [HNOI2007]神奇遊樂園【插頭dp】

bzoj 1187: [HNOI2007]神奇遊樂園【插頭dp】

要判邊界!!要判邊界!!要判邊界!!if(j!=m)!!!
我真是zz橫著轉移要判斷到底能不能向右邊出邊……
然後剩下的和1814差不多,九十因為不要求經過所有格子,所以左右括號隨時可以合併,但是注意合併的時候輪廓線上不能有別的括號,然後還有是(0,0)可以轉移到(0,0),相當一不經過這個格子

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=300005,mod=299989;
int n,m,nw,la,h[N],b[20],c[2];
long long a[105][105],ans=-1e18;
struct qwe
{
    int ne,to[2];
    long long va[2];
}f[N];
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>'9'||p<'0')
    {
        if(p=='-')
            f=-1;
        p=getchar();
    }
    while(p>='0'&&p<='9')
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
void add(int x,long long v)
{//cerr<<v<<endl;
    int u=x%mod+1;
    for(int i=h[u];i;i=f[i].ne)
        if(f[i].to[nw]==x)
        {
            f[i].va[nw]=max(f[i].va[nw],v);
            return;
        }
    c[nw]++;
    f[c[nw]].ne=h[u];
    f[c[nw]].to[nw]=x;
    f[c[nw]].va[nw]=v;
    h[u]=c[nw];
}
int main()
{
    n=read(),m=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]=read();
    b[0]=1;
    for(int i=1;i<=10;i++)
        b[i]=b[i-1]<<2;
    c[0]=1,f[1].va[0]=0,f[1].to[0]=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=c[nw];j++)
            f[j].to[nw]<<=2;
        for(int j=1;j<=m;j++)
        {//cerr<<i<<" "<<j<<endl;
            la=nw,nw^=1;
            memset(h,0,sizeof(h));
            c[nw]=0;
            for(int k=1;k<=c[la];k++)
            {
                int x=f[k].to[la],b1=(x>>(j*2-2))%4,b2=(x>>(j*2))%4;
                long long v=f[k].va[la];
                if(!b1&&!b2)
                {
                    add(x,v);
                    if(j!=m)
                        add(x+b[j-1]+2*b[j],v+a[i][j]);
                }
                else if(!b1&&b2)
                {
                    add(x+b2*b[j-1]-b2*b[j],v+a[i][j]);
                    if(j!=m)
                        add(x,v+a[i][j]);
                }
                else if(b1&&!b2)
                {
                    add(x,v+a[i][j]);
                    if(j!=m)
                        add(x-b1*b[j-1]+b1*b[j],v+a[i][j]);
                }
                else if(b1==1&&b2==1)
                {
                    for(int t=1,l=j+1;l<=m;l++)
                    {
                        if((x>>(l*2))%4==1)
                            t++;
                        if((x>>(l*2))%4==2)
                            t--;
                        if(!t)
                        {
                            add(x-b[j-1]-b[j]-b[l],v+a[i][j]);
                            break;
                        }
                    }
                }
                else if(b1==2&&b2==2)
                {
                    for(int t=1,l=j-2;l>=0;l--)
                    {
                        if((x>>(l*2))%4==1)
                            t--;
                        if((x>>(l*2))%4==2)
                            t++;
                        if(!t)
                        {
                            add(x-2*b[j-1]-2*b[j]+b[l],v+a[i][j]);
                            break;
                        }
                    }
                }
                else if(b1==2&&b2==1)
                    add(x-2*b[j-1]-b[j],v+a[i][j]);
                else if(x==b[j-1]+b[j]*2)
                        ans=max(ans,v+a[i][j]);
            }
        }
    }
    printf("%lld\n",ans);
    return 0;
}