[状压DP 边双连通分量] BZOJ 3590 [Snoi2013]Quare

陕西出出来的好题
题解

注意一个边双可以拆成一条链和一个小边双 然后就可以了

f[i]:集合状态为i.且使在i中的点双联通的最小权值.
h[i][j][0]:一个端点是j.另一个端点在点集i中的边的最小权值.
h[i][j][1]:一个端点是j.另一个端点在点集i中的边的次小权值.
g[i][j][k].集合状态为i.且使在i中的点构成一条链.两端点分别是j和k的最小权值.

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
using namespace std;

#define read(x) scanf("%d",&(x))

const int N=13;
const int M=45;
struct edge{
  int u,v,w,next;
}G[M<<1];
int head[N],inum;
inline void add(int u,int v,int w,int p){
  G[p].u=u; G[p].v=v; G[p].w=w; G[p].next=head[u]; head[u]=p;
}
#define V G[p].v

int n,m;

int f[1<<12];
int h[N][1<<12][2];
int g[1<<12][N][N];

int main(){
  int T,iu,iv,iw;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(T);
  while (T--){
    read(n); read(m);
    for (int i=1;i<=m;i++) read(iu),read(iv),read(iw),add(iu,iv,iw,++inum),add(iv,iu,iw,++inum);
    for (int i=1;i<=n;i++)
      for (int j=0;j<(1<<n);j++){
    int t1=1<<28,t2=1<<28;
    for (int p=head[i];p;p=G[p].next)
      if ((j>>(V-1)&1))
        if (G[p].w<t1)
          t2=t1,t1=G[p].w;
        else
          t2=min(t2,G[p].w);
    h[i][j][0]=t1; h[i][j][1]=t2;
      }
    memset(g,0x1f,sizeof(g));
    for (int i=1;i<=n;i++) g[1<<(i-1)][i][i]=0;
    for (int s=1;s<(1<<n);s++)
      for (int i=1;i<=n;i++)
    if (s>>(i-1)&1)
      for (int j=1;j<=n;j++)
        if (s>>(j-1)&1){
          for (int p=head[i];p;p=G[p].next)
        if (!(s>>(V-1)&1))
          g[s|(1<<(V-1))][V][j]=min(g[s|(1<<(V-1))][V][j],g[s][i][j]+G[p].w);
          for (int p=head[j];p;p=G[p].next)
        if (!(s>>(V-1)&1))
          g[s|(1<<(V-1))][i][V]=min(g[s|(1<<(V-1))][i][V],g[s][i][j]+G[p].w);
        }
    memset(f,0x1f,sizeof(f));
    for (int i=1;i<=n;i++) f[1<<(i-1)]=0;
    for (int s=1;s<(1<<n);s++){
      int rest=((1<<n)-1)^s;
      for (int t=rest;t;t=(--t)&rest)
    for (int i=1;i<=n;i++)
      if (t>>(i-1)&1)
        for (int j=1;j<=n;j++)
          if (t>>(j-1)&1)
        if (i==j)
          f[s^t]=min(f[s^t],f[s]+g[t][i][j]+h[i][s][0]+h[i][s][1]);
        else
          f[s^t]=min(f[s^t],f[s]+g[t][i][j]+h[i][s][0]+h[j][s][0]);
    }
    if (f[(1<<n)-1]>=1<<28)
      printf("impossible\n");
    else
      printf("%d\n",f[(1<<n)-1]);
    for (int i=1;i<=n;i++) head[i]=0; inum=0;
  }
  return 0;
}