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=5000,maxn=15,maxm=1010;

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

int eg[maxm*2][3],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][0]=b;eg[ect][1]=egh[a];eg[ect][2]=v;egh[a]=ect++;
    for(int i=1;i<bin[N+1];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]<=500000) dfs(x+1,y+dis0[s][x],s,e|bin[x]);
        dfs(x+1,y,s,e);
    }
}

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