标签:rand P2538 int 55 模拟退火 SCOI2008 include dis
传送门
解题思路
先用Floyd求一遍最短路,
然后对m个点之外的序列模拟退火,取前k个设置城堡,O(N^2)求出当前对小距离,判断并更新ans。
没什么注意的但是我rand()写在while外面了所以调了一个小时把我心态搞崩了草
AC代码
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<ctime>
using namespace std;
const double delta=0.996;
int n,m,k,dis[55][55],r[55],vis[55],ans,cnt,a[55],x;
void floyd(){
for(int k=1;k<=n;k++){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}
}
}
}
int cal(){
if(k==cnt) return 0;
int ans=0;
for(int i=k+1;i<=cnt;i++){
int res=1e9;
for(int j=1;j<=n;j++){
if(vis[j]){
res=min(res,dis[a[i]][j]);
}
}
ans=max(ans,res);
}
return ans;
}
void SA(){
double t=3000;
while(t>1e-15){
int x=rand()%k+1,y=rand()%(cnt-k)+k+1;
swap(a[x],a[y]);
for(int i=1;i<=k;i++) vis[a[i]]=1;
int res=cal();
for(int i=1;i<=k;i++) vis[a[i]]=0;
if(res<ans) ans=res;
else if(exp((ans-res)/t)*RAND_MAX<rand()) swap(a[x],a[y]);
t*=delta;
}
}
int main(){
ios::sync_with_stdio(false);
srand(time(0));
srand(rand());
srand(rand());
memset(dis,0x3f,sizeof(dis));
cin>>n>>m>>k;
if(m+k==n){
cout<<0;
return 0;
}
for(int i=1;i<=n;i++) cin>>r[i],r[i]++;
for(int i=1;i<=n;i++) cin>>x,dis[i][r[i]]=dis[r[i]][i]=min(dis[i][r[i]],x);
for(int i=1;i<=n;i++) dis[i][i]=0;
floyd();
for(int i=1;i<=m;i++) cin>>x,vis[x+1]=1;
for(int i=1;i<=n;i++) if(!vis[i]) a[++cnt]=i;
random_shuffle(a+1,a+cnt+1);
ans=cal();
while((double)clock()/CLOCKS_PER_SEC<=0.7) SA();
cout<<ans;
return 0;
}
标签:rand,P2538,int,55,模拟退火,SCOI2008,include,dis 来源: https://www.cnblogs.com/yinyuqin/p/15027168.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。