【POJ 2288】Islands and Bridges（状压DP）

Islands and Bridges

题目链接：POJ 2288

题目大意

给你一个图，然后要你找一个哈密顿路径（就是一条路径经过所有点不重不漏）。
然后一条路径的分数是每个点的权值加上相邻两个点权值积加上两两有边的相邻三个点的积。
然后要你求出最大的分数以及对于路径的方案数。
相反的路径看做同一条。

思路

考虑直接状压 DP 转移。
$f_{i,j,k}$ 为当前点是否到的状态时 $i$，最后到的点是 $j$，倒数第二个到的点是 $k$。
然后枚举下一个到的点然后转移。

然后具体的实现细节可以看代码。

代码

#include<cstdio>
#include<cstring>
#define ll long long

using namespace std;

const int N = 13 + 1;
const int B = 1 << 13; 
int q, n, m;
ll f[B][N][N], a[N], g[B][N][N];
bool in[N][N];

int main() {
	scanf("%d", &q);
	while (q--) {
		scanf("%d %d", &n, &m);
		for (int i = 1; i <= n; i++) scanf("%lld", &a[i]);
		for (int i = 1; i <= m; i++) {
			int x, y; scanf("%d %d", &x, &y);
			in[x][y] = in[y][x] = 1;
			
			//一开始的出发边
			f[(1 << (x - 1)) | (1 << (y - 1))][x][y] = f[(1 << (x - 1)) | (1 << (y - 1))][y][x] = a[x] + a[y] + a[x] * a[y];
			g[(1 << (x - 1)) | (1 << (y - 1))][x][y] = g[(1 << (x - 1)) | (1 << (y - 1))][y][x] = 1;
		}
		for (int i = 0; i < (1 << n); i++)//枚举现在状态
			for (int j = 1; j <= n; j++)//枚举倒数第二个到的点
				if ((i >> (j - 1)) & 1)
					for (int k = 1; k <= n; k++)//枚举最后到的点
						if (j != k && (i >> (k - 1)) & 1 && in[j][k] && g[i][k][j])//要有方案有边先
							for (int l = 1; l <= n; l++)//枚举新到的点
								if (!((i >> (l - 1)) & 1) && in[k][l]) {//至少要有边能走
									ll tov = f[i][k][j] + a[l] + a[l] * a[k] + (in[l][j] ? a[j] * a[k] * a[l] : 0);//算放上这个新的点之后的贡献
									if (tov > f[i | (1 << (l - 1))][l][k]) {
										f[i | (1 << (l - 1))][l][k] = tov;
										g[i | (1 << (l - 1))][l][k] = g[i][k][j];
									}
									else if (tov == f[i | (1 << (l - 1))][l][k]) {
										g[i | (1 << (l - 1))][l][k] += g[i][k][j];
									}
								}
		
		if (n == 1) {//特判只有一个点
			printf("%lld 1\n", a[1]);
		}
		else {
			ll ans = 0, an = 0;
			for (int i = 1; i <= n; i++)
				for (int j = 1; j <= n; j++) {
					if (ans < f[(1 << n) - 1][i][j]) {
						ans = f[(1 << n) - 1][i][j];
						an = g[(1 << n) - 1][i][j];
					}
					else if (ans == f[(1 << n) - 1][i][j]) {
						an += g[(1 << n) - 1][i][j];
					}
				}
			if (!ans || !an) printf("0 0\n");//无解
			else printf("%lld %lld\n", ans, an / 2);//因为反着走是同一条，所以方案数要除二
		}
		
		memset(in, 0, sizeof(in));//清空数组
		memset(f, 0, sizeof(f)); memset(g, 0, sizeof(g));
	}
	
	return 0;
}