hdu 4284 Travel （2012 天津网络赛 1007 ）

这题当场没有做出来， 看了芳姐代码后面还弄出了个超时，spfa 加位压缩的bfs过的；

先spfa算出目标点两两之间的最短距离；

bfs搜的是dp[i][j]; 当前所在点为i， 状态是j， 2进制j的第x位为1表示第x个目标点已经购买过了通行证；

复杂度O(15*(2^15))+O（15*15*spfa（））；险过

#include<cstdio>
#include<cstring>
#include<queue>
#define maxn 113
#define maxm 11005
using namespace std;
int vv[maxm], nxt[maxm], ww[maxm], h[maxn], e;
int d[maxn], svis[maxn];

void link( int u, int v, int c )
{
	vv[e] = v; ww[e] = c; nxt[e] = h[u]; h[u] = e++;
}

void spfa( int s )
{
	memset( svis, 0, sizeof(svis) );
	memset( d, -1, sizeof(d) );
	d[s] = 0;
	queue<int>q;
	q.push( s );
	while( !q.empty() ){
		int u = q.front();
		q.pop();
		svis[u] = 0;
		for( int j = h[u]; j+1; j = nxt[j] ){
			int v = vv[j];
			if( d[v] == -1 || d[v] > d[u] + ww[j] ){
				d[v] = d[u] + ww[j];
				if( !svis[v] )
					svis[v] = 1, q.push( v );
			}
		}
	}
}

int C[20], D[20], dd[20][20];
int dp[20][70000], vis[20][70000];
int num[20], H, S;
int bit[20], n, m, money;

int bfs(int is)
{
	for( int i = 1; i <= H; i++ ){
		memset( dp[i], -1, sizeof(int)*(bit[H+1]+5) );
		memset( vis[i], 0, sizeof(int)*(bit[H+1]+5) );
	}
	//memset(dp, -1, sizeof(dp));
	//memset(vis, 0, sizeof(vis));
	queue<int> q;
	if(is){
		q.push(70000*S);
		dp[S][0] = money;
		vis[S][0] = 1;
	}
	else{
		for( int i = 1; i <= H; i++ ){
			if( dd[S][i] >= 0 && money >= dd[S][i] + D[i] ){
				dp[i][bit[i]] = money - dd[S][i] - D[i] + C[i];
				q.push( i * 70000 + bit[i] );
				vis[i][bit[i]] = 1;
			}
		}
	}
	while( !q.empty() ){
		int u = q.front(); q.pop();
		int st = u / 70000, state = u % 70000;
		vis[st][state] = 0;
		if( dp[st][state] < 0 ) continue;
		for( int i = 1; i <= H; i++ ){
			if( dp[st][state] - dd[st][i] - D[i] < 0 )
				continue;
			if( (state&bit[i]) == 0 && dd[st][i] >= 0 ){
				int stv = state | bit[i];
				if( dp[i][stv] < dp[st][state] - dd[st][i] - D[i] + C[i] ){
						dp[i][stv] = dp[st][state] - dd[st][i] - D[i] + C[i];
						int v = i*70000 + stv;
						if( !vis[i][stv] )
							vis[i][stv] = 1, q.push( v );
				}
			}
		}
	}
	for( int i = 1; i <= H; i++ )
		if( dd[i][S] >= 0 && dp[i][bit[H+1]-1] - dd[i][S] >= 0 )
			return 1;
	return 0;
}
int main()
{
	bit[1] = 1;
	for( int i = 2; i <= 18; i++ ) bit[i] = bit[i-1]<<1;
	int T;
	int a, b, c;
	scanf( "%d", &T );
	while( T-- ){
		scanf( "%d%d%d", &n, &m, &money );
		e = 0;
		memset( h, -1, sizeof(h) );
		for( int i = 0; i < m; i++ ){
			scanf( "%d%d%d", &a, &b, &c );
			link( a, b, c );
			link( b, a, c );
		}
		scanf( "%d", &H );
		S = 0;
		for( int i = 1; i <= H; i++ ){
			scanf( "%d%d%d", &num[i], &C[i], &D[i] );
			if( num[i] == 1 )
				S = i;
		}
		int s = S;
		if( !S ){
			num[++H] = 1;
			C[H] = 0; D[H] = 0; 
			S = H;
		}
		for( int i = 1; i <= H; i++ ){
			spfa( num[i] );
			for( int j = 1; j  <= H; j++ )
				dd[i][j] = d[num[j]];
		}
		if( !s ) H--;  // 这一步非常重要， 它把最大复杂度O（2^16）降为O(2^15)，这个跑了3s, 乘2就超时了；
		if( bfs(s) )
			printf( "YES\n" );
		else
			printf( "NO\n" );
	}
}