NOIp2017D2T2(luogu3959) 宝藏 (状压dp)

时隔多年终于把这道题锅过了

数据范围显然用搜索剪枝状压dp。

可以记还有哪些点没到（或者已到了哪些点）、我们最深已到的是哪些点、这些点的深度是多少，然后一层一层地往下推。

但其实是没必要记最深的那一层的，只要强行装作每次更新都是用最深的深度更新就可以。这样的话，虽然会有很多情况偏大，但是能正确更新的情况其实是都已经包括了。

因为你如果想以当前状态去更新，但用的还不是最深一层的点的话，干脆就可以在之前你想用那个点处于最后一层的时候去更新。

代码写的很捉急..最后常数也很捉急...

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<vector>
 #include<queue>
 #include<cmath>
 #include<ctime>
 #define LL long long int
 #define inf 0x3f3f3f3f
 using namespace std;
 const int maxs=,maxn=,maxm=;

 LL rd(){
    LL x=;char c=getchar();int neg=;
    while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
    while(c>=''&&c<='') x=x*+c-'',c=getchar();
    return x*neg;
 }

 int eg[maxm*][],egh[maxn],ect;
 int f[maxn][maxs],dis0[maxs][maxn],dis[maxs][maxs],bin[maxn];
 int N,M;

 inline void adeg(int a,int b,int v){
     eg[ect][]=b;eg[ect][]=egh[a];eg[ect][]=v;egh[a]=ect++;
     for(int i=;i<bin[N+];i++) if((bin[a]&i) && !(bin[b]&i)) dis0[i][b]=min(dis0[i][b],v);
 }

 void dfs(int x,int y,int s,int e){
     if(x>N) dis[s][e]=y;
     else{
         if(dis0[s][x]<=) dfs(x+,y+dis0[s][x],s,e|bin[x]);
         dfs(x+,y,s,e);
     }
 }

 int main(){
     int i,j,k;
     N=rd(),M=rd();
     memset(egh,-,sizeof(egh));memset(dis,,sizeof(dis));
     memset(dis0,,sizeof(dis0));memset(f,,sizeof(f));
     for(i=,j=;i<=N+;i++,j<<=) bin[i]=j;
     for(i=;i<=M;i++){
         int a=rd(),b=rd(),c=rd();
         adeg(a,b,c);adeg(b,a,c);
     }
     for(i=;i<bin[N+];i++) dfs(,,i,);
     for(i=;i<=N;i++) f[][bin[i]]=;
     for(i=;i<N;i++){
         for(j=;j<bin[N+];j++){
             if(f[i][j]>=9e7) continue;
             for(k=;k<bin[N+];k++){
                 if(dis[j][k]>) continue;
                 f[i+][j|k]=min(f[i+][j|k],f[i][j]+i*dis[j][k]);
             }
         }
     }int ans=0x3f3f3f3f;
     for(i=;i<=N;i++) ans=min(f[i][(<<N)-],ans);
     printf("%d\n",ans);
     return ;
 }