UVALive 3181 Fixing the Great Wall

题目大意：现在长城上有n个点要修，每个点有三个属性，x、c、d，x表示位置，其他两个表示如果在 t 时刻修它，那么费用就是 c + t*d，修理点不占时间。有一个人要去修这些点，现在告诉你n、它的速度v、起始位置s，问你总的最小费用。

思路：这道题的关键地方有两个：（1）由于修理不占时间，那么时间全花在走路上，对于某一个区间，最小费用情况下，他修理完这些点的最终的位置肯定是两个端点其中之一。（2）光是想到上面那个还不行，因为有 t 的存在，这里计算费用需要变换一下，如果他现在从 i 点走到 j 点，设这个时间是 tt，本来话费计算方法应该是 原来的话费 + （cj + （前面用的时间 + tt）*d），现在我们用这个 原来的话费 + （cj + （所有没走过的点的 d 之和，包括 j ）*tt），相当于时间移动了。然后就是进行区间DP了，先排序，设 d[ i ][ j ][ dir ] 表示 i 号点 ~ j 号点的最小话费，dir == 0 表示最终在左端点，也就是i，dir == 1表示最终停在右端点，也就是 j ，状态转移方程为 d[ i ][ j ][ 0 ] = min(d[ i + 1 ][ j ][ 0 ] +( sum_d[ i ] + sum_d[ n ] - sum_d[ j ]) * (x(i+1) -x[ i ] )/v  + ci , d[ i+1 ][ j ][ 1 ] + (sum_d[ i ] + sum_d[ n ] - sum_d[ j ]) * ( x[ j ]  - x[ i ])/v + ci)，dir == 1 时也是两个取小，就不写了，具体看代码吧，我设的是d[ len ][ i ][ dir ] 表示长度为 len，以 i 为左端点，一样。。

自己想的时候，第一想法是前 i 个，然后 再多一维表示最终位置，可就是 t 这里搞不来，时间空间都不行。。 想不出其他的了，然后就去看了别人的代码，终于发现了他的思路。。。 原来 t 是这么处理的，再区间这样来搞，厉害啊，又学习了。。 然后敲完又 wa 了几遍，原因是因为都说是 int ，然后我也设的都是 int ，/v 的时候都是 int ，然后就悲剧了。。。所以简单起见，直接设 v 为 double，然后为了防止 sum_d 和 xi - xj 相乘溢出，就先 /v 。。 = =

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;

const double INF = 1e22;

const int MAXN = 1111;

struct Sec
{
    int x,c,d;
    void read()
    {
        scanf("%d%d%d",&x,&c,&d);
    }
    bool operator < (const  Sec &tmp) const
    {
        return x < tmp.x;
    }
} sec[MAXN];

double d[MAXN][MAXN][2];

int sum[MAXN];

int main()
{
    int n,s;
    double v;
    while(~scanf("%d%lf%d",&n,&v,&s)&&n)
    {
        for(int i = 1;i<=n;i++)
        {
            sec[i].read();
        }
        sort(sec+1,sec+1+n);
        sum[0] = 0;
        for(int i = 1;i<=n;i++)
            sum[i] = sum[i-1] + sec[i].d;
        for(int len = 1;len<=n;len++)
            for(int i = 1;i<=n;i++)
                d[len][i][0] = d[len][i][1] = INF;
        for(int i = 1;i<=n;i++)
        {
            d[1][i][0] = d[1][i][1] = sum[n]*(abs(s - sec[i].x)/v) + sec[i].c;
        }
        //for(int i = 1;i<=n;i++)
            //printf("len = 1,d0 = %f,d1 = %f\n",d[1][i][0],d[1][i][1]);
        for(int len = 2;len <= n;len++)
        {
            for(int i = 1;i<=n;i++)
            {
                int j = i+len-1;
                if(j>n) break;
                d[len][i][0] = min(d[len-1][i+1][0] + (sum[i] + sum[n] - sum[j])*((sec[i+1].x - sec[i].x)/v) + sec[i].c,
                                   d[len-1][i+1][1] + (sum[i] + sum[n] - sum[j])*((sec[j].x - sec[i].x)/v) + sec[i].c);
                d[len][i][1] = min(d[len-1][i][0] + (sum[i-1] + sum[n] - sum[j-1])*((sec[j].x - sec[i].x)/v) + sec[j].c,
                                   d[len-1][i][1] + (sum[i-1] + sum[n] - sum[j-1])*((sec[j].x - sec[j-1].x)/v) + sec[j].c);
                //printf("len = %d,i = %d,d0 = %f,d1 = %f\n",len,i,d[len][i][0],d[len][i][1]);
            }
        }
        printf("%.0f\n",min(d[n][1][0],d[n][1][1]));
    }
    return 0;
}