codevs 1391 伊吹萃香

*ps：老实说我真不想调试这个程序(T ^ T)

go to the problem

题目描述 Description

在幻想乡，伊吹萃香是能够控制物体密度的鬼王。因为能够控制密度，所以萃香能够制造白洞和黑洞，并可以随时改变它们。某一天萃香闲着无聊，在妖怪之山上设置了一些白洞或黑洞，由于引力的影响，给妖怪们带来了很大的麻烦。于是他们决定找出一条消耗体力最少的路，来方便进出。已知妖怪之山上有N个路口(编号1..N)，每个路口都被萃香设置了一定质量白洞或者黑洞。原本在各个路口之间有M条单向路，走过每一条路需要消耗一定量的体力以及1个单位的时间。由于白洞和黑洞的存在，走过每条路需要消耗的体力也就产生了变化，假设一条道路两端路口黑白洞的质量差为delta：

1. 从有白洞的路口走向有黑洞的路口，消耗的体力值减少delta，若该条路径消耗的体力值变为负数的话，取为0。

2. 从有黑洞的路口走向有白洞的路口，消耗的体力值增加delta。

3. 如果路口两端均为白洞或黑洞，消耗的体力值无变化。

由于光是放置黑洞白洞不足以体现萃香的强大，所以她决定每过1个单位时间，就把所有路口的白洞改成黑洞，黑洞改成白洞。当然在走的过程中你可以选择在一个路口上停留1个单位的时间，如果当前路口为白洞，则不消耗体力，否则消耗s[i]的体力。现在请你计算从路口1走到路口N最小的体力消耗。保证一定存在道路从路口1到路口N。

输入描述 Input Description

第1行：2个正整数N, M

第2行：N个整数，第i个数为0表示第i个路口开始时为白洞，1表示黑洞

第3行：N个整数，第i个数表示第i个路口设置的白洞或黑洞的质量w[i]

第4行：N个整数，第i个数表示在第i个路口停留消耗的体力s[i]

第5..M+4行：每行3个整数，u, v, k，表示在没有影响的情况下，从路口u走到路口v需要消耗k的体力。

输出描述 Output Description

第1行：1个整数，表示消耗的最小体力

样例输入 Sample Input
4 5

1 0 1 0

10 10 100 10

5 20 15 10

1 2 30

2 3 40

1 3 20

1 4 200

3 4 200

样例输出 Sample Output
130

数据范围及提示 Data Size & Hint

对于30%的数据：1 <= N <= 100, 1 <= M <= 500

对于60%的数据：1 <= N <= 1,000, 1 <= M <= 5,000

对于100%的数据：1 <= N <= 5,000, 1 <= M <= 30,000

其中20%的数据为1 <= N <= 3000的链

1 <= u,v <= N, 1 <= k,w[i],s[i] <= 200

按照1 -> 3 -> 4的路线。

---

可以看出是一个变化版的最短路，需要在跑SPFA时进行状态的维护。
开一个结构体记录一个点的编号和当前为白洞还是黑洞，更新时也分白洞和黑洞两种。

需要注意的地方有很多：

1. 走过每一条路需要消耗一定量的体力以及1个单位的时间，因此如果一条路的终点是白洞，这条路消耗的体力值用白洞计算，而更新的是终点为黑洞时的最短路（走过去就变为了黑洞）。
2. 不需要记走到i所花费的时间，每次只需要判断当前点的状态（黑洞还是白洞）与初始状态，如果不变，那终点也不变，如果变了，那终点也变了（每个时刻图中各点的变化是同步的）。
3. 从有白洞的路口走向有黑洞的路口，消耗的体力值减少delta，若该条路径消耗的体力值变为负数的话，取为0。
4. 每次更新时所用的队头元素要分停留与不停留两种情况进行讨论。

那么。。。
上代码吧？(o°ω°o)

代码_(:з」∠)_

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

int N,M,cnt;
int de[5010][2],used[5010][2],k[5010],w[5010],s[5010];
int first[5010],next[30010];
struct edge{
    int i,k;
}; 
struct maple{
    int f,t,d;
}Rode[30010];

queue<edge> q;
void SPFA()
{
    memset(de,63,sizeof(de));
    q.push((edge){1,1});
    q.push((edge){1,0});
    used[1][0]=1,used[1][1]=1;
    de[1][0]=0,de[1][1]=0;
    if(k[1]) de[1][0]=s[1];
    while(!q.empty())
    {
        edge a=q.front();
        q.pop();
        used[a.i][a.k]=0;
        for(int i=first[a.i];i;i=next[i])
        {
            int delta=w[Rode[i].t]-w[a.i],d=0,K=k[Rode[i].t];
            if(a.k!=k[a.i]) K=K^1; // 去反，不能用~   ~0=-1
            if(delta<0) delta*=-1;
            if(K!=a.k)
               if(!a.k) d=-1*delta;
               else d=delta;
            //不停留
            int D=Rode[i].d+d;
            if(D<0) D=0;
            if(de[Rode[i].t][K^1]>de[a.i][a.k]+D)
            {
                de[Rode[i].t][K^1]=de[a.i][a.k]+D;
                if(!used[Rode[i].t][K^1])
                {
                    used[Rode[i].t][K^1]=1;
                    q.push((edge){Rode[i].t,K^1});
                }
            }  
            //停留 
            D=Rode[i].d-d;
            if(D<0) D=0;
            if(a.k) D+=s[a.i];
            if(de[Rode[i].t][K]>de[a.i][a.k]+D)
            {
                de[Rode[i].t][K]=de[a.i][a.k]+D;
                if(!used[Rode[i].t][K])
                {
                    used[Rode[i].t][K]=1;
                    q.push((edge){Rode[i].t,K});
                }
            } 
        }
    }  
}
int main()
{
    scanf("%d%d",&N,&M);
    for(int i=1;i<=N;++i)  scanf("%d",&k[i]);
    for(int i=1;i<=N;++i)  scanf("%d",&w[i]);
    for(int i=1;i<=N;++i)  scanf("%d",&s[i]);
    for(int i=1;i<=M;++i)
    {
        scanf("%d%d%d",&Rode[i].f,&Rode[i].t,&Rode[i].d);
        next[i]=first[Rode[i].f];
        first[Rode[i].f]=i;
    }
    SPFA();
    cout<<min(de[N][0],de[N][1]);
    return 0;
}