[hdu6580]Milk

考虑定义以下dp数组：
1.g1[i][j]表示第i行从中间出发向左取j瓶牛奶最少要多久
2.g2[i][j]表示第i行从中间出发向右取j瓶牛奶最少要多久
3.g3[i][j]表示在g1[i][j]的基础上回到中间最少要多久
4.g4[i][j]表示在g2[i][j]的基础上回到中间最少要多久
5.f1[i][j]表示第i行从中间出发取到j瓶牛奶最少要多久
6.f2[i][j]表示在f1[i][j]的基础上回到中间最少要多久
7.f3[i][j]表示从(1,1)出发前i行取到j个物品最少要多久
8.f4[i][j]表示在f3[i][j]的基础上回到中间最少要多久(不算向下操作)
可以用dp算出以上八个数组（可以一起算就不用开vector），然后就可以用f3算出ans，时间复杂度$o(k^{2})$

#include<bits/stdc++.h>
using namespace std;
#define N 10005
#define inf 1000000LL*0x3f3f3f3f
#define ll long long
struct ji{
    int id,x,y,z;
    bool operator < (const ji &k)const{
        return (x<k.x)||(x==k.x)&&(y<k.y);
    }
}a[N];
set<int>s;
int T,n,m;
ll sum,g1[N],g2[N],g3[N],g4[N],f1[N],f2[N],f3[N],f4[N];
void calc(ll *x,ll *y,int l,int r,int t,int p){ 
    for(int j=1,k;j!=(r-l)*p+1;j++){
        s.clear();
        sum=0;
        for(k=l;abs(k-l)<j;k+=p){
            s.insert(a[k].z);
            sum+=a[k].z;
        }
        x[j]=min(x[j],sum+p*(a[k-p].y-t));
        y[j]=min(y[j],sum+p*(a[k-p].y-t)+abs(a[k-p].y-n));
        for(;k!=r;k+=p){
            int ma=*(--s.end());
            if (ma>a[k].z){
                sum-=ma-a[k].z;
                s.erase(--s.end());
                s.insert(a[k].z);
            }
            x[j]=min(x[j],sum+p*(a[k].y-t));
            y[j]=min(y[j],sum+p*(a[k].y-t)+abs(a[k].y-n));
        }
    }
}
int main(){
    scanf("%d",&T);
    while (T--){
        scanf("%*d%d%d",&n,&m);
        n=(n+1)/2;
        for(int i=1;i<=m;i++)scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
        sort(a+1,a+m+1);
        for(int i=1;i<=m;i++)f3[i]=f4[i]=inf;
        f3[0]=f4[0]=n-1;
        for(int i=1,j;i<=m;){
            for(j=i;(j<m)&&(a[j].x==a[j+1].x);j++);
            int l=i,r=j,mid;
            for(mid=l;(mid<=r)&&(a[mid].y<n);mid++);
            for(j=1;j<=r-l+1;j++)g1[j]=g2[j]=g3[j]=g4[j]=inf;
            calc(g1,g3,mid-1,l-1,n,-1);
            calc(g2,g4,mid,r+1,n,1);
            for(j=1;j<=r-l+1;j++)f1[j]=min(g1[j],g2[j]);
            for(j=1;j<=r-l+1;j++)f2[j]=min(g3[j],g4[j]);
            for(i=1;i<=mid-l;i++)
                for(j=1;j<=r-mid+1;j++){
                    f1[i+j]=min(f1[i+j],min(g1[i]+g4[j],g3[i]+g2[j]));
                    f2[i+j]=min(f2[i+j],g3[i]+g4[j]);
                }
            if (a[l].x==1){
                calc(f3,f4,1,r-l+2,1,1);
                i=r+1;
                continue;
            }
            for(i=l-1;i>=0;i--)
                for(j=1;j<=r-l+1;j++){
                    f3[i+j]=min(f3[i+j],f4[i]+a[l].x-1+f1[j]);
                    f4[i+j]=min(f4[i+j],f4[i]+f2[j]);
                }
            i=r+1;
        }
        for(int i=1;i<m;i++)printf("%lld ",f3[i]);
        printf("%lld\n",f3[m]);
    }
}