hdu 3001 状态压缩~最优值

一开始的想法是将n个点扩充成2n个节点，复杂有2k万，感觉有点高，静下心想了一下，感觉可以用三进制表示状态，发现确实可以，初始化没写好，郁闷了半小时。

ACcode：

#include<cstdio>
#include<cstring>

const int NS=11;
const int LIM=600000;
const int INF=(1<<30);

int n,m;
int g[NS][NS],p3[NS];
int dp[LIM][NS],bit[LIM][NS];

int Min(int a1,int b1)
{
return a1<b1?a1:b1;
}

void prepare()
{
p3[0]=1;
for (int i=1;i<NS;i++)
p3[i]=p3[i-1]*3;
for (int i=0;i<LIM;i++)
for (int x=i,j=0;x>0&&j<NS;j++)
bit[i][j]=x%3,x/=3;
}

int judge(int x)
{
for (int i=0;i<n;i++)
if (!bit[x][i]) return 0;
return 1;
}

int main()
{
prepare();
int x,y,z,lim,res;
while (~scanf("%d%d",&n,&m))
{
lim=p3
,res=INF;
for (int i=0;i<lim;i++)
for (int j=0;j<n;j++)
dp[i][j]=INF;
for (int i=0;i<n;i++)
dp[p3[i]][i]=0;
for (int i=0;i<n;i++)
for (int j=0;j<n;j++)
g[i][j]=i==j?0:INF;

for (int i=0;i<m;i++)
{
scanf("%d%d%d",&x,&y,&z);
--x,--y;
g[y][x]=g[x][y]=Min(g[x][y],z);
}
for (int i=1;i<lim;i++)
for (int flag=judge(i),j=0;j<n;j++)
{
if (!bit[i][j]) continue;
y=i-p3[j];
for (int k=0;k<n;k++)
if (bit[y][k]&&g[j][k]<INF)
dp[i][j]=Min(dp[i][j],dp[y][k]+g[j][k]);
if (flag) res=Min(res,dp[i][j]);
}
if (res>=INF) res=-1;
printf("%d\n",res);
}
return 0;
}