CF-Educational Round102-E-Minimum Path:最短路变形

题目大意:

给你一张无向带权图。询问从1点出发到每个点的最小
${\sum }_{i=1}^k{e}_{k_i}-{\max }_{i=1}^k{e}_{k_i}+{\min }_{i=1}^k{e}_{k_i}(e_1,e_2,...,e_k为从1到x的路径(不一定是简单路径)上的边权)$

$n,m\le 2e5,e_i\le 1e9$

题目思路:

首先，若没有后面两个式子，那么就是一个最短路。所以可以预料到这是一个最短路的变形题.

考虑这样的图:
3 2
1 2 4
2 3 1

那么$dis[2]=3$.因为我们可以这么走:$1->2->3->2$

这是可以来回走的模型，自然就想到要对图分层.(即对状态升维).

将这个问题转化为:跑最短路，你可以且必须进行如下两种操作:
1.选择路径上的一条边，减去它的边权
2.选择路径上的一条边，再加一次它的边权。

那么$dis[i][0/1[0/1]$代表从$1$到$i$点，是否已经使用了操作1，是否已经使用了操作2的最小值.

那么答案就是$dis[i][1][1]$.

正确性：很显然，对于任意一条路径，我们肯定是减去最大边权，加上最小边权。但问题是我们在跑最短路的过程中不知道哪条边会是最终的最小边权/最大边权(有后效性)。所以我们需要对状态升维来消除后效性.

心得：不管式子怎么复杂，其实我们发现它只多了两项，也就是只有两条边最终会受影响。所以可以考虑暴力升维，从而表示出所有状态。

AC代码:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
#define pb push_back
#define mp make_pair
const int maxn = 5e5 + 5;
const int maxm = 5e5 + 5;
const int mod = 1e9 + 7;
int u[maxm] , v[maxm] , nextt[maxm] , first[maxn] , sign;
ll w[maxm];
void addedge(int x , int y , int z){
    u[++sign]=x;
    v[sign]=y;
    w[sign] = z;
    nextt[sign]=first[x];
    first[x]=sign;
}
struct Node{
    ll id , dis;
    int s0 , s1;
    Node (int n = 0, ll d = 1e17 , int s0 = 0 , int s1 = 0)
    {
        id = n;dis = d;
        this->s0 = s0;
        this->s1 = s1;
    }
    bool operator < (const Node & a)const
    {
        return dis > a.dis;
    }
};
int book[maxn][2][2] , n , m;
ll dis[maxn][2][2];
void Dijstra(ll s){
    priority_queue<Node>q;
    while(q.size()) q.pop();
    for (int i = 1 ; i <= n ; i++){
        for (int j = 0 ; j <= 1 ; j++){
            for (int k = 0 ; k <= 1 ; k++){
                book[i][j][k] = 0;
                dis[i][j][k] = 1e16;
            }
        }
    }
    dis[s][0][0] = 0;
    q.push(Node(s , 0 , 0 , 0));
    while(q.size())
    {
        Node g = q.top() ; q.pop();
        int id = g.id;
        ll dist = g.dis;
        int s0 = g.s0 , s1 = g.s1;
        if (book[id][s0][s1]) continue;
        book[id][s0][s1] = 1;
        for (int i = first[id] ; i ; i = nextt[i]){
            if (dis[v[i]][s0][s1] > dist + w[i]){
                dis[v[i]][s0][s1] = dist + w[i];
                q.push(Node(v[i] , dis[v[i]][s0][s1] , s0 , s1));
            }
            if (!s0 && dis[v[i]][1][s1] > dist){
                dis[v[i]][1][s1] = dist;
                q.push(Node(v[i] , dis[v[i]][1][s1] , 1 , s1));
            }
            if (!s1 && dis[v[i]][s0][1] > dist + 2 * w[i]){
                dis[v[i]][s0][1] = dist + 2 * w[i];
                q.push(Node(v[i] , dis[v[i]][s0][1] , s0 , 1));
            }
            if (!s0 && !s1 && dis[v[i]][1][1] > dist + w[i]){
                dis[v[i]][1][1] = dist + w[i];
                q.push(Node(v[i] , dis[v[i]][1][1] , 1 , 1));
            }
        }
    }
    return ;
}
int main()
{
    ios::sync_with_stdio(false);
    cin >> n >> m;
    for (int i = 1 ; i <= m ; i++){
        int x , y , z; cin >> x >> y >> z;
        addedge(x , y , z);
        addedge(y , x , z);
    }
    Dijstra(1);
    for (int i = 2 ; i <= n ; i++){
        cout << dis[i][1][1] << " ";
    }
    cout << endl;
    return 0;
}