hdu 4284 Travel(floyd + TSP)

虽然题中有n<=100个点，但实际上你必须走过的点只有H<=15个。而且经过任意点但不消耗C[i]跟D[i]可以为无限次，所以可以floyd预处理出H个点的最短路，之后剩下的。。。就成了裸的TSP了，dp[sta][i]表示已经遍历过了sta集合中的点，现在在i点所需的最少花费。dp[i][j]=-1表示该点不可到达。。但是在最后统计最小解的时候要特判。。。

 

#include<algorithm>

#include<iostream>

#include<cstring>

#include<cstdlib>

#include<fstream>

#include<sstream>

#include<bitset>

#include<vector>

#include<string>

#include<cstdio>

#include<cmath>

#include<stack>

#include<queue>

#include<stack>

#include<map>

#include<set>

#define FF(i, a, b) for(int i=a; i<b; i++)

#define FD(i, a, b) for(int i=a; i>=b; i--)

#define REP(i, n) for(int i=0; i<n; i++)

#define CLR(a, b) memset(a, b, sizeof(a))

#define debug puts("**debug**")

#define LL long long

#define PB push_back

#define MP make_pair

#define eps 1e-10

using namespace std;



const int maxn = 105;

const int INF = 1e9;

int n, m, H, V, T;

int go[20], C[20], D[20], g[maxn][maxn], dp[1<<15][20];



void floyd()

{

    FF(i, 1, n+1)

    {

        FF(j, 1, n+1) g[i][j] = INF;

        g[i][i] = 0;

    }

    int u, v, w;

    while(m--)

    {

        scanf("%d%d%d", &u, &v, &w);

        if(g[u][v] > w) g[u][v] = g[v][u] = w;

    }

    FF(k, 1, n+1) FF(i, 1, n+1) FF(j, 1, n+1) g[i][j] = min(g[i][j], g[i][k] + g[k][j]);

}



int main()

{

    scanf("%d", &T);

    while(T--)

    {

        scanf("%d%d%d", &n, &m, &V);

        floyd();



        scanf("%d", &H);

        REP(i, H) scanf("%d%d%d", &go[i], &C[i], &D[i]);



        int tot = (1<<H) - 1;

        CLR(dp, -1);

        REP(i, H) if(V - g[1][go[i]] >= D[i]) dp[1<<i][i] = g[1][go[i]] - C[i] + D[i];



        FF(i, 1, tot+1) REP(j, H) if(dp[i][j] != -1) //dp[i][j]可到达

        {

            REP(k, H) if((i&(1<<k)) == 0) //k点未到达

            {

                int sta = i|(1<<k);

                if(V - dp[i][j] - g[go[j]][go[k]] >= D[k])

                {

                    if(dp[sta][k] == -1) dp[sta][k] = dp[i][j] + g[go[j]][go[k]] - C[k] + D[k];

                    else dp[sta][k] = min(dp[sta][k], dp[i][j] + g[go[j]][go[k]] - C[k] + D[k]);

                }

            }

        }

        int ans = INF;

        //特判dp[tot][i]可到达，最开始忘了。。。10+wa

        REP(i, H) if(dp[tot][i] != -1) ans = min(ans, dp[tot][i] + g[go[i]][1]);

        printf("%s\n", ans <= V ? "YES" : "NO");

    }

    return 0;

}