[POJ2288] [2004上海区域赛] Islands and Bridges [状态压缩][dp]

[$\mathfrak{L}\mathfrak{i}\mathfrak{n}\mathfrak{k}$]
给定一个$1\le \mathfrak{n}\le 13$点$1\le \mathfrak{m}\le 100$边无向图，每个点有权值$1\le {\mathfrak{v}}_{\mathfrak{i}}\le 100$，
求权值最大的哈密顿回路（以某个点为起点，经过其他每个点一次的回路）的权值和数量。权值定义如下：
　1.记${\mathfrak{V}}_1$表示$\sum {\mathfrak{v}}_{\mathfrak{i}}$
　2.记${\mathfrak{V}}_2$表示$\sum {\mathfrak{v}}_{\mathfrak{i}}{\mathfrak{v}}_{\mathfrak{i}+1}$当${\mathfrak{v}}_{\mathfrak{i}}$和${\mathfrak{v}}_{\mathfrak{i}+1}$在哈密顿回路中且被边直接相连
　3.记${\mathfrak{V}}_3$表示$\sum {\mathfrak{v}}_{\mathfrak{i}}{\mathfrak{v}}_{\mathfrak{i}+1}{\mathfrak{v}}_{\mathfrak{i}+2}$当${\mathfrak{v}}_{\mathfrak{i}}$、${\mathfrak{v}}_{\mathfrak{i}+1}$和${\mathfrak{v}}_{\mathfrak{i}+2}$在哈密顿回路中且构成一个连通分量
　那么权值$\mathfrak{V}=\sum _{\mathfrak{i}=1}^3{\mathfrak{V}}_{\mathfrak{i}}$。

看见题目的第一反应： ？？？？？？？？？？？？？？？？？？

拆开来考虑每个权值。
在某一点时候的状态由之前的两个点决定。按照访问顺序将这三个点标为$\mathfrak{A},\mathfrak{B},\mathfrak{C}$。
发现一定会有${\mathfrak{V}}_1$和${\mathfrak{V}}_2$。而${\mathfrak{V}}_3$题目也提示了，只要$\mathfrak{A}$和$\mathfrak{C}$有边直接相连就会产生。
这样一来发现其实转移并不难。这道题也确实就只是一个普通的TSP变形。状压即可。

本题细节很多，要注意。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<ctime>
using namespace std;
int n,m;
long long Ans=0,t;
long long BAns;
long long Pt[15]={};
bool lnk[15][15]={};
long long F[8200][15][15]={};
long long C[8200][15][15]={};
int main()
{
	scanf("%d",&n);
	while(~scanf("%d%d",&n,&m))
	{
		Ans=-1; BAns=0;
		memset(lnk,0,sizeof(lnk));
		memset(C,0,sizeof(C));
		memset(F,-1,sizeof(F));
		for(int i=1;i<=n;++i)
		{
			scanf("%lld",&Pt[i]);
		}
		if(n==1){printf("%lld 1\n",Pt[1]);continue;} //***
		for(int u,v,i=1;i<=m;++i)
		{
			scanf("%d%d",&u,&v);
			lnk[u][v]=lnk[v][u]=1;
		}
		for(int a=1;a<=n;++a)
		{
			for(int b=1;b<=n;++b)
			{
				if(a==b||lnk[a][b]==0)continue;
				F[(1<<a-1)|(1<<b-1)][a][b]=Pt[a]*Pt[b]+Pt[a]+Pt[b];
				C[(1<<a-1)|(1<<b-1)][a][b]=1; //*
			}
		}
		for(int i=1;i<(1<<n);++i)
		{
			for(int a=1;a<=n;++a)
			{
				if(!(i&(1<<a-1)))continue;
				for(int b=1;b<=n;++b)
				{
					if(!(i&(1<<b-1)))continue;
					if(lnk[a][b]==0||a==b)continue;
					if(F[i][a][b]==-1)continue;
					for(int v,t,c=1;c<=n;++c)
					{
						if(i&(1<<c-1))continue;
						if(lnk[b][c]==0||b==c||a==c)continue;
						v=i|(1<<c-1);
						t=F[i][a][b]+Pt[b]*Pt[c]+Pt[c]+((lnk[a][c])?(Pt[a]*Pt[b]*Pt[c]):0);
						if(t>F[v][b][c])F[v][b][c]=t,C[v][b][c]=C[i][a][b]; //*
						else if(t==F[v][b][c])C[v][b][c]+=C[i][a][b]; //*
					}
				}
			}
		}
		t=(1<<n)-1;
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=n;++j)
			{
				if(i==j)continue;
				if(F[t][i][j]>Ans)
				{
					Ans=F[t][i][j];
					BAns=C[t][i][j];
				}
				else if(F[t][i][j]==Ans)BAns+=C[t][i][j];
			}
		}
		if(Ans<0)printf("0 0\n");
		else printf("%lld %lld\n",Ans,BAns>>1);
	}
    return 0;
}