hdu 4571 Travel in time

         Floyd + Spfa + dp。先用Floyd预处理出所有点对之间的最短路,然后考虑到始点和终点是必须经过的,并且可以不游玩,于是就不能用原来的始点和终点,我们考虑新加入原点和终点,并且原点花费时间和满意值都为0(这样保证不会对结果做出影响),终点花费时间为0,满意值大于满意值的最大值(保证如果可以到达原来终点一定可以到达新终点)。接下来就是从新建图了,新原点到旧原点的路程为0,新原点到别的点的路程为旧原点到该点的距离。终点建边同理。然后别的边正常处理即可。

注意有重边。

        对于dp[ i ][ j ]表示当前点为 i , 使用时间 j 获得的最大满意度。则就变成很明显的Spfa了。 

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#include<cmath>
#define LL long long
#define CLR(a, b) memset(a, b, sizeof(a))

using namespace std;

const int N = 110;
const int M = N * N;
const int INF = 0x3f3f3f3f;

struct Edge
{
    int u, v, c;
}E[M];

struct Node
{
    int u, tim;
};

int fir[N], next[M], tot;
int G[N][N], n, m, t, s, e, cas = 1;
int C[N], S[N], dp[N][N * 3];

void Add_Edge(int u, int v, int c)
{
    E[tot].u = u, E[tot].v = v, E[tot].c = c;
    next[tot] = fir[u], fir[u] = tot ++;
}

void input()
{
    int u, v, c, i, j;
    scanf("%d%d%d%d%d", &n, &m, &t, &s, &e);
    CLR(G, INF);
    for(i = 0; i < n; i ++)
        scanf("%d", &C[i]);
    for(i = 0; i < n; i ++)
    {
        scanf("%d", &S[i]);
        G[i][i] = 0;
    }
    for(i = 0; i < m; i ++)
    {
        scanf("%d%d%d", &u, &v, &c);
        G[u][v] = G[v][u] = min(c, G[u][v]);
    }
}

void floyd()
{
    int i, j, k;
    for(k = 0; k < n; k ++)
        for(i = 0; i < n; i ++)
            for(j = 0 ; j < n; j ++)
                G[i][j] = min(G[i][j], G[i][k] + G[k][j]);
}

void init()
{
    int i, j;tot = 0;CLR(fir, -1);
    S[n] = C[n] = 0;
    S[n + 1] = N;C[n + 1] = 0;
    for(i = 0; i < n; i ++)
    {
        if(i != s) Add_Edge(n, i, G[s][i]);
        else Add_Edge(n, s, 0);
        if(i != e) Add_Edge(i, n + 1, G[e][i]);
        else Add_Edge(e, n + 1, 0);
    }
    for(i = 0; i < n; i ++)
        for(j = 0; j < n; j ++)
            if(S[j] > S[i])
                Add_Edge(i, j, G[i][j]);
    s = n; e = n + 1;n += 2;
}

void Spfa()
{
    int v, i, j, c, ans;
    queue<Node> Q;CLR(dp, 0);
    dp[s][0] = 0;Node A, B;
    A.u = s;A.tim = 0;Q.push(A);
    while(!Q.empty())
    {
        A = Q.front();Q.pop();
        for(i = fir[A.u]; ~i; i = next[i])
        {
            v = E[i].v;c = E[i].c;
            int tmp = A.tim + c + C[v];
            if(tmp <= t && dp[v][tmp] < dp[A.u][A.tim] + S[v])
            {
                dp[v][tmp] = dp[A.u][A.tim] + S[v];
                B.u = v;
                B.tim = tmp;
                Q.push(B);
            }
        }
    }
    ans = 0;
    for(i = 0; i <= t; i ++)
    {
        ans = max(ans, dp[e][i]);
    }
    printf("Case #%d:\n%d\n", cas ++, max(ans - N, 0));
}

int main()
{
    //freopen("input.txt", "r", stdin);
    int w;
    scanf("%d", &w);
    while(w --)
    {
        input();
        floyd();
        init();
        Spfa();
    }
}


  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值