BOI 2002 双调路径

题目：

城市的道路是双向的，每条道路有固定的旅行时间以及需要支付的费用。路径由连续的道路组成。总时间是各条道路旅行时间的和，总费用是各条道路所支付费用的总和。同样的出发地和目的地，如果路径A比路径B所需时间少且费用低，那么我们说路径A比路径B好。对于某条路径，如果没有其他路径比它好，那么该路径被称为最优双调路径。这样的路径可能不止一条，或者说根本不存在。
给出城市交通网的描述信息，起始点和终点城市，求最优双条路径的条数。城市不超过100个，边数不超过300，每条边上的费用和时间都不超过100。

输入描述

第一行给出有多少个点,多少条边,开始点及结束点.

输出描述

有多少条双调路径

 

样例输入

4 5 1 4
2 1 2 1
3 4 3 1
2 3 1 2
3 1 1 4
2 4 2 4

样例输出 Sample Output

2 

分析：

本题用的一种巧妙的做法。显然这题是一道单源最短路问题，用SPFA或是Dijkstra都是可以的。但是我们要进行SPFA的是对于某个固定的花费下的所用的时间。也就是程序当中的dist[city][cost] = time这个在注释里有写。因为费用最大是100，所以cost开到二维数组开出来不会爆内存的！

本题有几个要注意的小问题：（有些是我愚蠢的错误）

1. 因为是无向边Edge数组记得开2*MAXM
2. 有的变量名字注意一下哈：time叫t，target也叫t，那不是很愚蠢啊……

程序：

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 100 + 1;
const int MAXM = 300 + 1;
int n, m, s, t, EdgeCount = 0, MinTime = 0x3f3f3f3f, Head[MAXN];
int dist[MAXN][MAXN*MAXN];
bool vis[MAXN][MAXN*MAXN];
struct edge
{
    int Next, Aim, Weight, Time;
}Edge[MAXM*2];
struct node
{
    int pos, cost;
};
void insert(int u, int v, int w, int t)
{
    Edge[++EdgeCount] = (edge){Head[u],v,w,t};
    Head[u] = EdgeCount;
}
void SPFA()
{
    memset(vis,false, sizeof(vis));
    memset(dist, 0x3f, sizeof(dist));
    //dist[city][cost] = time
    queue<node> Q;
    Q.push((node){s,0});
    dist[s][0] = 0;
    vis[s][0] = true;
    while (!Q.empty())
    {
        node p = Q.front();
        Q.pop();
        //vis[p.pos][p.cost] = false;
        if(p.cost > n*MAXN)
            continue;
        for (int i = Head[p.pos]; i; i = Edge[i].Next)
        {
            int v = Edge[i].Aim;
            if(dist[v][p.cost + Edge[i].Weight] > dist[p.pos][p.cost] + Edge[i].Time)
            {
                dist[v][p.cost+Edge[i].Weight] = dist[p.pos][p.cost] + Edge[i].Time;
                if(!vis[v][p.cost + Edge[i].Weight])
                {
                    vis[v][p.cost + Edge[i].Weight] = true;
                    Q.push((node){v,p.cost + Edge[i].Weight});
                }
            }
        }
    }
}
int main()
{
    memset(Head,0,sizeof(Head));
    cin >> n >> m >> s >> t;
    for (int i = 1; i <= m; i++)
    {
        int u, v, w, t;
        cin >> u >> v >> w >> t;
        insert(u,v,w,t);
        insert(v,u,w,t);
    }
    SPFA();
    int ans = 0;
    for (int i = 0; i <= MAXN*n; i++)
    {
        //cout << dist[t][i] << endl;
        if ((!(vis[t][i])) || dist[t][i] >= MinTime)
            continue;
        MinTime = dist[t][i];
        ans++;
    }
    cout << ans << endl;
    return 0;
}