图的收缩

一个图在连通情况下，每个点都有一个包含它的圈。一种图去掉任意一条边，都能成为前面那种图，则符合要求，求这种图的最小边权值之和。

#include<cstdio>#include<cstring>#include<string.h>#include<algorithm>#include<iostream>#include<math.h>#include<vector>using namespace std;const int maxn=13;const int inf=100000000;struct Element{int p,q;}a[maxn][maxn];int mark[maxn],b[maxn],esta[maxn];int t,n,m,minn;void search(int u,int h,int c,int e){if(c+e>=minn)return;int cc=c;if(h==2&&a[u][1].q!=inf)cc+=a[u][1].q;if(h!=2&&a[u][1].p!=inf)cc+=a[u][1].p;if(cc!=c&&cc<minn){int i,j;Element aa[maxn][maxn];memcpy(aa,a,sizeof(a[0])*(n+1));int bb[maxn];memcpy(bb,b,sizeof(int)*(h+1));int ee[maxn];memcpy(ee,esta,sizeof(int)*(n+1));bool end=true;for(j=1;j<=n;j++)if(!mark[j]){end=false;for(i=1;i<h;i++)if(a[b[i]][j].q<a[1][j].p)a[1][j].p=a[b[i]][j].p,a[1][j].q=a[b[i]][j].q,a[j][1]=a[1][j];else if(a[b[i]][j].q<a[1][j].q)a[1][j].q=a[b[i]][j].p,a[j][1].q=a[1][j].q;else if(a[b[i]][j].p<a[1][j].p)a[1][j].q=a[1][j].p,a[1][j].p=a[b[i]][j].p,a[j][1]=a[1][j];else if(a[b[i]][j].p<a[1][j].q)a[1][j].q=a[b[i]][j].p,a[j][1].q=a[1][j].p;}if(end)minn=cc;else{search(1,1,cc,e);memcpy(a,aa,sizeof(a[0])*(n+1));memcpy(b,bb,sizeof(int)*(h+1));memcpy(esta,ee,sizeof(int)*(n+1));}}for(int i=1;i<=n;i++)if(!mark[i]&&a[u][i].p!=inf){b[h]=i;mark[i]=h;search(i,h+1,c+a[u][i].p,e-esta[i]);mark[i]=0;}}int main(){int i,j,x,y,c,est;scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)a[i][j].p=a[i][j].q=inf;while(m--){scanf("%d%d%d",&x,&y,&c);if(c<a[x][y].p)a[x][y].q=a[x][y].p,a[x][y].p=c;else if(c<a[x][y].q)a[x][y].q=c;a[y][x]=a[x][y];}est=0;for(int i=1;i<=n;i++){esta[i]=inf;for(j=1;j<=n;j++)if(a[i][j].p<esta[i])esta[i]=a[i][j].p;est+=esta[i];}memset(mark,0,sizeof(mark));mark[1]=1;b[0]=1;minn=inf;search(1,1,0,est-esta[1]);if(minn==inf)printf("impossible\n");else printf("%d\n",minn);}return 0;}