【NOIP模拟】艾比所特

Description

阿良良木历将要迎来人生（不，是吸血鬼生涯）的第二次战斗——与身为人类和吸血鬼混血儿的艾比所特在直江津高中的操场solo，以取回Heartunderblade的左脚。
艾比所特个子不大，却单手扛着一副形状比自己大三倍，重量是他体重乘三次方的巨大十字架。吸血鬼是怕十字架的，所以历要远远地躲着。更麻烦的是艾比所特能变成雾气，历攻击不到他，他却因为是混血儿，吸血鬼的弱点降到了mininum，可以碰十字架。直到历用沙子使艾比所特显形，历才获得胜利。
可以把直江津高中的操场看成排成一条直线的n个点，其中第i个点位于数轴上的位置x[i]。如果从i跳到j，首先要花费时间|x[j]-x[i]|。假如i

Solution

这道题，题意很奇怪，字句间暗藏玄机。
我开始还以为是只用找一个环，那么就是值最大的欧拉回路。
然后其实是可以有很多个简单环的。

那么就有一个很显然的性质

每个点的入度和出度只有1。
我们设f[i][j][k]表示做到第i个点，在1~i之间有j个点的还没有出度，有k个点还没有入度。
我们考虑从第i个点推向第i+1个点。
第一种情况：前i个点有一个点连向第i+1个点，第i+1个点又连回来——>$f[i+1][j-1][k-1]=max(f[i+1][j+1][k+1],f[i][j][k]+x[i+1]*2+a[i+1]+c[i+1]);$
第二种情况：前i个点有一个点连向第i+1个点，第i+1个点往后面连——>$f[i+1][j][k]=max(f[i+1][j][k],f[i][j][k]+a[i+1]+d[i+1]);$为什么j和k没有-1呢？因为推到i+1是又会多一个点的，那么影响j和k的点原本是又会多一个，但是这里会抵消一个，上面的那种情况会抵消两个。
第三种情况：第i+1个点连向前i个点的一个点，i+1后面的点连向第i+1个点——>$f[i+1][j][k]=max(f[i+1][j][k],f[i][j][k]+b[i+1]+c[i+1]);$
第四种情况：i+1后面的点连向第i+1个点，第i+1个点又连向i+1后面的点——>$f[i+1][j+1][k+1]=max(f[i+1][j+1][k+1],f[i][j][k]-x[i+1]*2+b[i+1]+d[i+1]);$这种情况不会抵消。
注意边界条件，前三种情况是在j>0||k>0的时候才能向后推的，因为都与i前面的点有关，前面的点用完了，自然不能连。

缩小空间

我们发现j和k都是同时加减的，所以三维可以合并为二维。
其实还可以继续缩小，当然是滚动啦。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a;i>=b;i--)
#define rep(i,a) for(i=first[a];i;i=next[i])
using namespace std;
typedef long long ll;
const int maxn=10007;
ll i,j,k,l,t,n,m,ans,da;
ll x[maxn],a[maxn],b[maxn],c[maxn],d[maxn],p,q;
ll f[2][maxn];
int main(){
    scanf("%lld",&n);
    fo(i,1,n)scanf("%lld",&x[i]);
    fo(i,1,n)scanf("%lld",&a[i]);
    fo(i,1,n)scanf("%lld",&b[i]);
    fo(i,1,n)scanf("%lld",&c[i]);
    fo(i,1,n)scanf("%lld",&d[i]);
    memset(f,128,sizeof(f));da=f[0][0];
    f[0][0]=0;p=0;
    fo(i,0,n-1){
        q=p^1;
        fo(j,0,i){
            if(f[p][j]==da)continue;
            if(j){
                f[q][j-1]=max(f[q][j-1],f[p][j]+x[i+1]*2+a[i+1]+c[i+1]);
                f[q][j]=max(f[q][j],f[p][j]+a[i+1]+d[i+1]);
                f[q][j]=max(f[q][j],f[p][j]+b[i+1]+c[i+1]);  
            }
            f[q][j+1]=max(f[q][j+1],f[p][j]-x[i+1]*2+b[i+1]+d[i+1]);  
        }
        p=p^1;
    }
    printf("%lld\n",f[p][0]);
}