【BZOJ】1097: [POI2007]旅游景点atr（spfa+状压dp）

http://www.lydsy.com/JudgeOnline/problem.php?id=1097

首先还是我很sb。。。。想到了分层图想不到怎么串起来，，，以为用拓扑序搞转移，，后来感到不行。。。

QAQ

这种数据那么小，有明确的依赖性为嘛我想不到状压。。。（准确的说是没想到状压和分层图一起做。。。。

还有一个。。。。。。。为什么递推不行。。。（还是我写挫了。。。老wa。。）非得记忆化。。。。。。

（其实记忆化的话能省点不必要的状态吧。。。

。。。。

还有。。。和zyf一样，被卡在了。。!(s&(1<<(v-1)))。。。。。不去掉就tle。。也是喜闻乐见。。。

QAQ

（然后我spfa的vis标记打错了。。。。。。。调了好久。。QAQ

#include <cstdio>

#include <cstring>

#include <cmath>

#include <string>

#include <iostream>

#include <algorithm>

#include <queue>

#include <set>

#include <map>

using namespace std;

typedef long long ll;

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

#define for1(i,a,n) for(int i=(a);i<=(n);++i)

#define for2(i,a,n) for(int i=(a);i<(n);++i)

#define for3(i,a,n) for(int i=(a);i>=(n);--i)

#define for4(i,a,n) for(int i=(a);i>(n);--i)

#define CC(i,a) memset(i,a,sizeof(i))

#define read(a) a=getint()

#define print(a) printf("%d", a)

#define dbg(x) cout << (#x) << " = " << (x) << endl

#define error(x) (!(x)?puts("error"):0)

#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)

#define grdm(g, x, i) for(int i=g.ihead[x]; i; i=g.e[i].next)

inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }



const int N=50005, oo=0x3f3f3f3f;

int d[22][N], vis[N], q[N], ihead[N], cnt, n, m, k;

struct dat { int next, to, w; }e[N*10*2];

void add1(int u, int v, int w) { e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v; e[cnt].w=w; }

void add(int u, int v, int w) { add1(u, v, w); add1(v, u, w); }

void spfa(int s, int *f) {

	memset(vis, 0, sizeof(int)*(n));

	memset(f, 0x3f, sizeof(int)*(n));

	int front=0, tail=0;

	f[s]=0; q[tail++]=s; vis[s]=1;

	while(front!=tail) {

		int u=q[front++], v; if(front==N) front=0; vis[u]=0;

		rdm(u, i) if(f[v=e[i].to]>f[u]+e[i].w) {

			f[v]=f[u]+e[i].w;

			if(!vis[v]) {

				vis[v]=1;

				if(f[q[front]]>f[v]) {

					--front; if(front<0) front+=N;

					q[front]=v;

				}

				else {

					q[tail++]=v; if(tail==N) tail=0;

				}

			}

		}

	}

}

int ans[21][1<<(20)], b[22], all;

int dfs(int u, int s) {

	if(ans[u][s]>=0) return ans[u][s];

	if(s==all) return d[u][n-1];

	int &a=ans[u][s]; a=oo;

	for1(v, 1, k) if((b[v]&s)==b[v]) a=min(a, d[u][v]+dfs(v, s|(1<<(v-1))));

	return a;

}

int main() {

	read(n); read(m); read(k);

	for1(i, 1, m) { int u=getint(), v=getint(), w=getint(); add(u-1, v-1, w); }

	for1(i, 0, k) spfa(i, d[i]);

	read(m);

	for1(i, 1, m) { int u=getint(), v=getint(); b[v-1]|=(1<<(u-2)); }

	all=(1<<k)-1;

	CC(ans, -1);

	printf("%d\n", dfs(0, 0));

	return 0;

}

　　

 

---

 

 

Description

FGD想从成都去上海旅游。在旅途中他希望经过一些城市并在那里欣赏风景，品尝风味小吃或者做其他的有趣的事情。经过这些城市的顺序不是完全随意的，比如说FGD不希望在刚吃过一顿大餐之后立刻去下一个城市登山，而是希望去另外什么地方喝下午茶。幸运的是,FGD的旅程不是既定的，他可以在某些旅行方案之间进行选择。由于FGD非常讨厌乘车的颠簸，他希望在满足他的要求的情况下，旅行的距离尽量短，这样他就有足够的精力来欣赏风景或者是泡MM了^_^. 整个城市交通网络包含N个城市以及城市与城市之间的双向道路M条。城市自1至N依次编号，道路亦然。没有从某个城市直接到它自己的道路，两个城市之间最多只有一条道路直接相连，但可以有多条连接两个城市的路径。任意两条道路如果相遇，则相遇点也必然是这N个城市之一，在中途，由于修建了立交桥和下穿隧道，道路是不会相交的。每条道路都有一个固定长度。在中途，FGD想要经过K(K<=N-2)个城市。成都编号为1，上海编号为N,而FGD想要经过的N个城市编号依次为2,3,…,K+1. 举例来说，假设交通网络如下图。FGD想要经过城市2,3,4,5，并且在2停留的时候在3之前，而在4,5停留的时候在3之后。那么最短的旅行方案是1-2-4-3-4-5-8,总长度为19。注意FGD为了从城市2到城市4可以路过城市3,但不在城市3停留。这样就不违反FGD的要求了。并且由于FGD想要走最短的路径，因此这个方案正是FGD需要的。

Input

第一行包含3个整数N(2<=N<=20000),M(1<=M<=200000),K(0<=K<=20),意义如上所述。以下M行，每行包含3个整数X，Y，Z，(1<=X

Output

只包含一行，包含一个整数，表示最短的旅行距离。

Sample Input

8 15 4
1 2 3
1 3 4
1 4 4
1 6 2
1 7 3
2 3 6
2 4 2
2 5 2
3 4 3
3 6 3
3 8 6
4 5 2
4 8 6
5 7 4
5 8 6
3
2 3
3 4
3 5

Sample Output

19

HINT

 

 上面对应于题目中给出的例子。

 

Source