P3959-宝藏【模拟退火】
生活随笔
收集整理的這篇文章主要介紹了
P3959-宝藏【模拟退火】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P3959
題目大意
nnn個點,mmm條邊,求一棵有根生成樹權值最小。對于一條邊(fa,x,w)(fa,x,w)(fa,x,w)會產生權值depfa?wdep_{fa}*wdepfa??w。
解題思路
我們模擬退火每次隨機一個序列,然后貪心一個一個加入生成樹中。然后一系列玄學參數即可。
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #include<cstdlib> #define ll long long using namespace std; const ll N=13,inf=5e6+10; const double delit=0.996; ll n,m,a[N][N],p[N],k[N],v[N],ans=inf; ll calc(){memset(v,0,sizeof(v));ll ans=0;v[p[1]]=1;k[p[1]]=0;for(ll i=2;i<=n;i++){ll k=1e18;for(ll j=1;j<i;j++){ll x=p[i],y=p[j];if(a[y][x]<inf&&a[y][x]*v[y]<k)k=a[y][x]*v[y],v[x]=v[y]+1;}ans+=k;}return ans; } void S_A(){double t=10000;srand(rand());while(t>1e-15){ll x=rand()%n+1,y=rand()%n+1;swap(p[x],p[y]);ll now=calc(),k=now-ans;if(k<0)ans=now;else if(exp(-1.0*k/t)*RAND_MAX<rand())swap(p[x],p[y]);t*=delit;}return; } void Solve(){for(ll i=1;i<=n;i++)p[i]=i;ans=calc();for(ll i=1;i<=100;i++)S_A();return; } int main() {srand(19260817+114514);scanf("%lld%lld",&n,&m);memset(a,0x3f,sizeof(a)); for(ll i=1;i<=m;i++){ll x,y,w;scanf("%lld%lld%lld",&x,&y,&w);a[x][y]=a[y][x]=min(a[x][y],w);}Solve();printf("%lld\n",ans); } /* 5 4 1 2 3 1 3 3 1 4 3 1 5 3*/總結
以上是生活随笔為你收集整理的P3959-宝藏【模拟退火】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: P1337-[JSOI2004]平衡点/
- 下一篇: jzoj6824-[2020.10.17