【最短路】CheapestRoute

题目描述

有一列$n$个单元格组成的信息通道，从左至右标号为$0$至$n-1$。有一个信息要从最左边传送到最右边的单元格中。每一次你可以从一个单元格向相邻的左或右单元格传送。每一个单元格有一个价格$p[i]$，如果要传入单元格$i$，则要花费$p[i]$，另外，如果$p[i]$为$-1$，则表示不可以进入第i格。
还有$m$部电话。第$i$部电话的开始位置在$enter[i]$单元格，另一头在$exit[i]$单元格，你可以直接把信息从$enter[i]$传递到$exit[i]$，但你要花费的代价为：$K+x$，$K$是一个给定的定值，$x$则为之前你使用电话的次数。如果你使用电话进入一个单元格$j$，则不须花费$p[j]$。另外类似上面，价格为$-1$的单元格也是不可以进入使用电话的。
现要你找一个花费$C$为最小的到达最右边单元格的方案。花费$C$为最少的前提下，再求移动次数$M$最少的方案。进入相邻格或使用电话都算一次移动。

输入

第一行一个整数$Z\,(1\le Z\le 20)$，表示有多少组测试数据。
对于每一组数据
第一行：第一个数为$N$，后面有$N$个整数，对应$p[i]$；
第二行：第一个数为$M$，后面有$M$个整数，相应表示$enter[i]$；
第三行：第一个数为$M$，后面有$M$个整数，相应表示$exit[i]$；
第四行：一个整数 $K$。

输出

对应每组测试数据，有一行两个整数：$C \; M\;$（行末有一个空格） 。
如果无解，输出两个”$-1$”。 两个“$-1$”之间两个空格。

样例输入

2
1 100
1 0
1 0
1000
4 0 -1 0 0
1 0
1 2
1000

样例输出

0 0
1000 2

---

这啥玩意儿？

算法一：DFS

维护一个$x$一个dist来个可行性剪枝大概就好了？
实测得分100【11个测试点】【跑的还很快】
卧槽数据太水了。

算法二：图论有关

第二眼看，这是一道【边权可变】的最短路问题。
建图：相邻的非$-1$单元格连边，标记非电话，代价从$p$数组找。
电话边标记为电话，权值随$x$改变。

想法一：普通spfa

$dist$数组普通维护肯定是不行的，因为如果$dist$一样松弛的条件还与$x$和$step$有关。
那优先级是$x$还是$step$高呢？
本地测试发现都会WA掉一个点。优先级还是不能确定。
难道这题的后效性导致真的不能用最短路算法？

想法二：加dp

$dist$数组开两维，$dist[i][j]$表示用$j$次电话到达$i$单元格所需要的最小代价。
每次更新$dist$的时候像dp一样，松弛条件和转移方程就是：

• 对于从i到j的电话边，使用$x$次电话：
  $$dist[j][x+1]=min(dist[j][x+1], dist[i][x]+x+K).$$
• 对于从i到j的普通边，使用$x$次电话：
  $$dist[j][x]=min(dist[j][x], dist[i][x]+p[j]).$$

$x$可循环求解。
还可在$dist$数组中再开一维，表示到达当前状态所需要的最少行动数。

上代码【一贯糟烂的代码风格】：

#include <iostream>
#include <cstdio>
#include <queue>
#include <climits>
using namespace std;
struct node{int to; bool phone; node* next; node(int t, int p):to(t), phone(p), next(NULL){}};
struct li{
    node *pre, *last;
    void push(int t, int p){if(pre)last=last->next=new node(t, p);else pre=last=new node(t, p);}
    li():pre(NULL), last(NULL){}
}graph[50];                                                 //邻接表
int z, n, m ,k, p[50], phone[50][2], dist[50][101][2];
//phone:电话起止点 dist:点，使用电话数，{距离，行动数}
bool in[50]={0};
void spfa(){
    for(int j=0;j<n;j++){
        for(int k=0;k<=2*m;k++){
            dist[j][k][0]=INT_MAX;
            dist[j][k][1]=0;
        }
        in[j]=0;
    }
    dist[0][0][0]=0;
    queue<int> q;
    q.push(0);
    in[0]=1;
    int w;
    while(!q.empty()){
        w=q.front();
        in[w]=0;
        q.pop();
        for(node *t=graph[w].pre;t;t=t->next){
            if(t->phone){                                    //电话边
                for(int i=0;i<=2*m;i++)                      //电话使用循环到2*m【或者更小】即可
                    if(dist[w][i][0]==INT_MAX)
                        continue;
                    else if(dist[t->to][i+1][0]>dist[w][i][0]+i+k){
                        dist[t->to][i+1][0]=dist[w][i][0]+i+k;
                        dist[t->to][i+1][1]=dist[w][i][1]+1;
                        if(!in[t->to]){
                            in[t->to]=1;
                            q.push(t->to);
                        }
                    }
                    else if(dist[t->to][i+1][0]==dist[w][i][0]&&dist[t->to][i+1][1]>dist[w][i][1]+1)
                        dist[t->to][i+1][1]=dist[w][i][1]+1;
            }
            else                                             //非电话边
                for(int i=0;i<=2*m;i++)
                    if(dist[w][i][0]==INT_MAX)
                        continue;
                    else if(dist[t->to][i][0]>dist[w][i][0]+p[t->to]){
                        dist[t->to][i][0]=dist[w][i][0]+p[t->to];
                        dist[t->to][i][1]=dist[w][i][1]+1;
                        if(!in[t->to]){
                            in[t->to]=1;
                            q.push(t->to);
                        }
                    }
                    else if(dist[t->to][i][0]==dist[w][i][0]+p[t->to]&&dist[t->to][i][1]>dist[w][i][1]+1)
                        dist[t->to][i][1]=dist[w][i][1]+1;
        }
    }
}
int main()
{
    scanf("%d", &z);
    for(int o=0;o<z;o++){
        scanf("%d", &n);
        for(int i=0;i<n;i++)
            scanf("%d", &p[i]);
        for(int i=0;i<n;i++){                          //单元格建图
            if(p[i]==-1)
                continue;
            if(i>0&&p[i-1]!=-1)
                graph[i].push(i-1, 0);
            if(i<n-1&&p[i+1]!=-1)
                graph[i].push(i+1, 0);
        }
        for(int i=0;i<2;i++){
            scanf("%d", &m);
            for(int j=0;j<m;j++)
                scanf("%d", &phone[j][i]);
        }
        for(int i=0;i<m;i++)                           //电话边建图
            if(p[phone[i][1]]!=-1)
                graph[phone[i][0]].push(phone[i][1], 1);
        scanf("%d", &k);
        spfa();
        int res=INT_MAX, q;
        for(int i=0;i<=m*2;i++)                        //取答案
            if(dist[n-1][i][0]==INT_MAX)
                continue;
            else if(dist[n-1][i][0]<res){
                res=dist[n-1][i][0];
                q=dist[n-1][i][1];
            }
            else if(dist[n-1][i][0]==res&&dist[n-1][i][1]<q)
                q=dist[n-1][i][1];
        if(res!=INT_MAX)
            printf("%d %d \n", res, q);                 //本地数据有问题，若评测对格式要求较严格请自行更改
        else
            printf("-1  -1\n");
    }
    return 0;
}

---

spfa真是神奇，套上各种东西乱跑还真挺好使。