POJ_3635

? ? 可以把圖上一個點拆成C+1個點然后直接應(yīng)用dij即可，同時總狀態(tài)是一定的，為了減少每一步的決策量，我們可以每次加油都只加1個單位。

? ? 貌似這樣的題目用dij都會比SPFA快一些，因為狀態(tài)比較多，而dij搜到終點就可以退出，這樣相比SPFA減少了對很多無用狀態(tài)的搜索量。

#include<stdio.h>
#include<string.h>
#define MAXD 1010
#define MAXM 20010
#define MAXC 110
#define MAXT 110000
#define INF 0x3f3f3f3f
int N, M, Q, D, S, E, C, e, cost[MAXD], first[MAXD], next[MAXM], v[MAXM], w[MAXM];
int f[MAXT], tree[4 * MAXT];
void update(int i)
{
for(; i ^ 1; i >>= 1)
tree[i >> 1] = f[tree[i]] < f[tree[i ^ 1]] ? tree[i] : tree[i ^ 1];
}
void new_edge(int x, int y, int z)
{
v[e] = y, w[e] = z;
next[e] = first[x], first[x] = e;
++ e;
}
void init()
{
int i, j, k, x, y, z;
for(i = 0; i < N; i ++)
scanf("%d", &cost[i]);
memset(first, -1, sizeof(first[0]) * N);
e = 0;
for(i = 0; i < M; i ++)
{
scanf("%d%d%d", &x, &y, &z);
new_edge(x, y, z), new_edge(y, x, z);
}
}
int get_id(int x, int c)
{
return x * (C + 1) + c + 1;
}
void divide_id(int id, int &x, int &c)
{
c = (id - 1) % (C + 1);
x = (id - 1) / (C + 1);
}
void solve()
{
int i, j, k, x, c, id;
scanf("%d", &Q);
for(k = 0; k < Q; k ++)
{
scanf("%d%d%d", &C, &S, &E);
for(D = 1; D < N * (C + 1); D <<= 1);
memset(f + 1, 0x3f, sizeof(f[0]) * (N * (C + 1)));
f[0] = INF;
memset(tree + 1, 0, sizeof(tree[0]) * (D << 1));
j = get_id(S, 0);
f[j] = 0, tree[j + D] = j, update(j + D);
while(id = tree[1])
{
divide_id(id, x, c);
if(x == E)
break;
tree[D + id] = 0, update(D + id);
if(c < C)
{
j = get_id(x, c + 1);
if(f[id] + cost[x] < f[j])
{
f[j] = f[id] + cost[x];
tree[j + D] = j, update(j + D);
}
}
for(i = first[x]; i != -1; i = next[i])
if(c >= w[i])
{
j = get_id(v[i], c - w[i]);
if(f[id] < f[j])
{
f[j] = f[id];
tree[j + D] = j, update(j + D);
}
}
}
if(tree[1] == 0)
printf("impossible\n");
else
printf("%d\n", f[tree[1]]);
}
}
int main()
{
while(scanf("%d%d", &N, &M) == 2)
{
init();
solve();
}
return 0;
}

總結(jié)

以上是生活随笔為你收集整理的POJ 3635 Full Tank?的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。