【jzoj 3290】【luogu P6085】Foodie / 吃货 JYY（数位DP）（欧拉回路）

Foodie / 吃货 JYY

题目链接：jzoj 3290 / luogu P6085

题目大意

有 n 个点，有一些路径一定要走，有一些路径可以走可以不走，都有走的费用。
路径双向，然后问你从 1 号点出发经过所有点再回到 1，所需要的最小费用。

思路

虽然之前做过，但我发现重新写一次（模拟赛碰到）却不太记得了。
——>之前<——
这一次就不讲做法了，主要是讲讲思路。

首先不难看出是要找一个欧拉回路，欧拉回路有一个特点就是每个点的度数都是偶数，而我们就根据这个来 DP。

然后搞三进制的 $f_i$。
主要是你 DP 求出来的 $f_i$ 不直接是答案，你需要通过再 DP 一个二进制的 $g_i$。
对于每个 $f_i$，再判断必连的连了之后，找到要让它全部是偶的需要对上的 $g_j$，然后两个加在一起，取最小的那个才是答案。

代码

#include<cstdio>
#include<vector>
#include<cstring>
#define mi313 15943223

using namespace std;

struct node {
	int to, nxt;
}e[201]; 
int n, k, m, x, y, z, mi3[14];
int dis[14][14], ans, sum, le[14], KK;
int f[2000001], num[14], g[10001];

void add(int x, int y) {
	e[++KK] = (node){y, le[x]}; le[x] = KK;
	e[++KK] = (node){x, le[y]}; le[y] = KK;
}

int main() {
	mi3[0] = 1;
	for (int i = 1; i <= 13; i++)
		mi3[i] = mi3[i - 1] * 3;
	for (int i = 0; i <= 13; i++)
		for (int j = 0; j <= 13; j++)
			dis[i][j] = 1e9;
	ans = 1e9;
	
	scanf("%d %d", &n, &k);
	for (int i = 1; i <= k; i++) {
		scanf("%d %d %d", &x, &y, &z);
		x--; y--;
		sum += z;
		add(x, y);
		num[x]++; num[y]++;
		dis[x][y] = min(dis[x][y], z);
		dis[y][x] = min(dis[y][x], z);
	}
	scanf("%d", &m);
	for (int j = 1; j <= m; j++) {
		scanf("%d %d %d", &x, &y, &z);
		x--; y--;
		dis[x][y] = min(dis[x][y], z);
		dis[y][x] = min(dis[y][x], z);
	}
	
	//Floyed
	for (int i = 0; i < n; i++) dis[i][i] = 0;
	for (int k = 0; k < n; k++)
		for (int i = 0; i < n; i++)
			for (int j = 0; j < n; j++)
				dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
	
	//求 g
	for (int i = 0; i < (1 << n); i++)
		g[i] = 1e9;
	g[0] = 0;
	for (int now = 0; now < (1 << n); now++) {
		for (int i = 0; i < n; i++)
			if (!((now >> i) & 1))
				for (int j = i + 1; j < n; j++)
					if (!((now >> j) & 1))
						g[now | (1 << i) | (1 << j)] = min(g[now | (1 << i) | (1 << j)], g[now] + dis[i][j]);
	}
	
	//求 f
	for (int i = 0; i < mi3[n]; i++)
		f[i] = 1e9;
	f[2] = 0;
	for (int now = 2; now < mi3[n]; now++)
		if (f[now] != 1e9) {
			for (int i = 0; i < n; i++)
				if (now / mi3[i] % 3 == 0) {
					for (int j = le[i]; j; j = e[j].nxt)
						if (now / mi3[e[j].to] % 3) {//一定要选
							int to = now + 2 * mi3[i];
							f[to] = min(f[to], f[now]);
						}
					for (int j = 0; j < n; j++)
						if (now / mi3[j] % 3) {//可以选可以不选
							int to = now + mi3[i];
							if (now / mi3[j] % 3 == 1) to += mi3[j];
								else to -= mi3[j];
							f[to] = min(f[to], f[now] + dis[i][j]);
						}
				}
		}
	
	for (int now = 0; now < mi3[n]; now++) {//枚举作为答案的 f
		bool can = 1;
		for (int i = 0; i < n; i++)//判断一定要选的点有没有选
			if (num[i] && now / mi3[i] % 3 == 0) {
				can = 0;
				break;
			}
		if (can) {
			int newn = now;
			for (int i = 0; i < n; i++)
				if (num[i] & 1) {//一定要选的度数是奇数，使得你总的度数奇偶性变化
					if (newn / mi3[i] % 3 == 1) newn += mi3[i];
						else if (newn / mi3[i] % 3 == 2) newn -= mi3[i];
				}
			int to = 0;
			for (int i = 0; i < n; i++)//然后看把所有奇的变成偶的要的 g
				if (newn / mi3[i] % 3 == 1)
					to |= (1 << i);
			ans = min(ans, f[now] + g[to]);
		}
	}
	
	printf("%d", ans + sum);
	
	return 0;
}