Noip2017 宝藏

题目传送门

感觉这个题被玩坏了……
翻看了一下题解区，发现了\(N\)种神仙做法：

1. 随机化搜索(伪模拟退火)
2. 搜索+剪枝
   
3. 模拟退火
4. 不对但是能A的状压\(DP\)
5. 正确的状压\(DP\)
   其中\(DP\)又分用\(DFS\)的\(DP\)，不用\(DFS\)的\(DP\)，而且各种\(DP\)的时间复杂度也各不相同……

我还有什么话可说呢？

---

学习科学，使用玄学
——zhx

既然标签上都写了搜索和随机化了，那我们肯定要用第一种神仙做法(其实是因为我看不懂状压DP)

这里的随机化用到了一点模拟退火的思想，以一定的概率去接受新的连边方式，这个概率是随着\(n\)的上升而增大的
为了寻找最优解，我们可以让它多跑几遍，取最小值

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define inf 1000000
using namespace std;
int g[15][15],deth[15];
struct zzz{
    int f,t;
}pass[1010]; int top;
bool operator < (struct zzz a,struct zzz b){
    return deth[a.f]*g[a.f][a.t] > deth[b.f]*g[b.f][b.t];
}
bool vis[15];
int n,m;
int bfs(int s){
    priority_queue <zzz> q;
    int cost=0;
    //将一开始的边入堆 
    for(int i=1;i<=n;i++)
      if(g[s][i]<inf)
        q.push(zzz{s,i});
    //因为要多加(n-1)条边,所以循环到n-1 
    for(int i=1;i<n;i++){
        zzz e=q.top(); q.pop();
        while(!q.empty()&&(vis[e.t]||(rand()%n<1))){ //随机找边 
            if(!vis[e.t]) pass[++top]=e;  //跳过的边仍然可能用到,先将它们存起来 
            e=q.top(); q.pop();  //随机一条新边 
        }
        vis[e.t]=1; deth[e.t]=deth[e.f]+1; 
        for(int i=1;i<=top;i++)  //将之前跳过的边在怼进堆里 
          q.push(pass[i]);
        top=0;
        for(int i=1;i<=n;i++)
          if(g[e.t][i]<inf&&!vis[i])  //扩展新的边 
            q.push(zzz{e.t,i});
        cost+=g[e.f][e.t]*deth[e.f];
    }
    return cost;
}
int read(){
    int k=0; char c=getchar();
    while(c<'0'||c>'9') c=getchar();
    while(c>='0'&&c<='9')
      k=k*10+c-48, c=getchar();
    return k;
}
int main(){
    srand(19****17); //打码保平安 
    memset(g,127/3,sizeof(g));
    n=read(),m=read();
    for(int i=1;i<=m;i++){  //因为点少边多,使用邻接矩阵存图更为方便 
        int x=read(),y=read();
        g[x][y]=g[y][x]=min(g[x][y],read());
    }
    int ans=0x7fffffff;
    for(int i=1;i<=500;i++)  //为寻找最优解,多跑几遍 
      for(int i=1;i<=n;i++){  //枚举根节点 
          //初始化 
          memset(deth,127/3,sizeof(deth));
          memset(vis,0,sizeof(vis));
          deth[i]=vis[i]=1; top=0;
          ans=min(ans,bfs(i))
      }
    cout<<ans;
    return 0;
}

• 参考资料
  题解 P3959【宝藏】——EternalAlexander

转载于:https://www.cnblogs.com/wxl-Ezio/p/9929250.html