1033 To Fill or Not to Fill (25分)待解待解

这道题做不来，书上说是贪心算法，我感觉这应该是图，，呃呃，是直线，，，不过有点像图中加权路径求最短。。。，，
我现在有感觉和图比简单些，因为，它是一条直线，也就是说局部最优可以构成整体最优。具有最优子结构性质，，，，再做做

明天做！昨天的思路错了：将所有加油站按单价从低到高排序。然后只要能在满油耗尽之前到那里，就在上一站加那么多距离的油。。。反例：但是这样可能中途有比出发站便宜的油，但又比下一个目的地油贵的，这样就错过了，至少在改段行程中不是最优解决方案
正确思路：将加油站按距离排序。把加油站看成点，而油费就是两点之间的路径的权重，每次只有选择最短权重就可以。但如果下一站价格更贵，又要分情况讨论吗》是否加满，改：是的，如果下一站的更贵，那就把油加满，去到最小的更贵的点加油。
改：每次确实选择权重最小的路，但是本题问的是在各个点加还是不加油，所以重点应该放在选择哪个点，选择后加满还是不加满。

3分的代码，思维混乱，以后再看

#include<iostream>
#include<algorithm>
using namespace std;
struct node{
    double p,k;
}s[510];
bool cmp(node a,node b)
{
    return a.k<b.k;
}
int main()
{
    double m,n,k,z,pos=0,init1,sum=0;
    cin>>m>>n>>k>>z;
    for(int i=0;i<z;i++)
        scanf("%lf %lf",&s[i].p,&s[i].k);
    sort(s,s+(int)z,cmp);
    init1=s[0].p;
    if(z>1&&s[0].k==0)
    {
        if(s[1].k>m*k)
        {
            printf("The maximum travel distance = %.2lf\n",m*k);
            return 0;
        }
    }
    else
    {
        if(n-s[0].k>m*k)
        {
             printf("The maximum travel distance = %.2lf\n",m*k);
             return 0;
        } 
    }
    if(s[0].k!=0)
    {
        printf("The maximum travel distance = 0.00\n");
    }
    else
    {
        for(int i=1;;)
        {
            if(n-pos<m*k)
            {
                while(s[i].p>init1&&i!=z)
                    i++;
                if(i==z)
                {
                    sum+=init1*(n-pos);
                    break;
                }
                else
                {
                    init1=s[i].p;
                    sum+=(s[i].k-pos)*s[i].p;
                    pos=s[i].k;
                }
            }
            if(n-pos>m*k&&i==z-1)
            {
                pos=s[i].k+m*k;
                printf("The maximum travel distance = %.2lf",pos);
                return 0;
            }
            if(s[i].p<init1&&(s[i].k-pos)<=m*k&&n-pos>m*k)
            {
                sum+=(s[i].k-pos)*init1;
                init1=s[i].p;
                pos=s[i].k;
            }
            else if(s[i].p>init1&&s[i].k-pos<=m*k&&n-pos>m*k)
            {
                while(s[i].k-pos>m*k)
                    i++;
                init1=s[--i].p;
                sum+=(s[i].k-pos)*s[i].p;
                pos=s[i].k;
            }
        }
    }
    printf("%.2lf",sum);
}

改：
1⃣️应该以点

#include<iostream>
#include<algorithm>
#include<iomanip>
using namespace std;
double cmax,d,davg,n,nowp=0,nowl=0,maxd,sum=0,nowi,rem=0;
struct station{
    double price;
    int location;
}p[501];
bool cmp(station a,station b){
    return a.location<b.location;
}
int findop(int j){
    int temp=0;
    ++j;
    double minp=999999;
    while(p[j].location<=(nowl+maxd))
    {
        if(p[j].price<nowp)     //离起点最近的油价更低点
            return j;
        if(p[j].price<=minp){      //都比起点油价高时，油价最低的点
            minp=p[j].price;
            temp=j;
        }
        ++j;
    }
    return temp;
}
int main(){
    cin>>cmax>>d>>davg>>n;
    maxd=cmax*davg;
    for(int i=0;i<n;i++){
        cin>>p[i].price>>p[i].location;
    }
    p[int(n)].price=-1;
    p[int(n)].location=d;
    sort(p,p+int(n+1),cmp);
    nowp=p[0].price;
    nowi=0;
    if(p[0].location!=0)
        cout<<"The maximum travel distance = "<<"0.00"<<endl;
    else{
        for(int i=1;i<=n;i++){
            int op=findop(nowi);
            if(p[op].price==-1){        //下一个最优点是终点
                sum+=(p[op].location-nowl)/davg*nowp;
                printf("%.2lf\n",sum);
                return 0;
            }
            else if(op&&p[op].price<=nowp){              //这个最优点价格更低，加油到这里
                double t=rem*davg;
                sum+=(p[op].location-nowl-t)/davg*nowp;
                nowl=p[op].location;
                nowi=op;
                nowp=p[op].price;
                rem=0;
            }
            else if(!op)                 //下一个最优点是0
            {
               double t=nowl+maxd;
                cout<<"The maximum travel distance = ";
                printf("%.2lf",t);
                return 0;
            }
            else{
                sum+=nowp*(cmax-rem);
                rem=cmax-(p[op].location-nowl)/davg;
                nowi=op;
                nowl=p[op].location;
                nowp=p[op].price;
            }
            i=op;
        }
    }
    return 0;
}