codevs 1344 模拟退火

1344 线型网络

 

 时间限制: 1 s

 空间限制: 128000 KB

 题目等级 : 钻石 Diamo

 

题目描述 Description

有 N ( <=20 ) 台 PC 放在机房内，现在要求由你选定一台 PC，用共 N-1 条网线从这台机器开始一台接一台地依次连接他们，最后接到哪个以及连接的顺序也是由你选定的，为了节省材料，网线都拉直。求最少需要一次性购买多长的网线。（说白了，就是找出 N 的一个排列 P1 P2 P3 ..PN 然后 P1 -> P2 -> P3 -> ... -> PN 找出 |P1P2|+|P2P3|+...+|PN-1PN| 长度的最小值)

输入描述 Input Description

第一行 N ，下面 N 行，每行分别为机器的坐标(x,y) ( 实数 -100<=x,y<=100 )

输出描述 Output Description

最小的长度，保留两位小数。

样例输入 Sample Input

3
0 0
1 1
1 -1

样例输出 Sample Output

2.83

______________________________________________________________________________________________________________________________

题目数据范围看似很小，实际细算，肯定超时。常规方法不行。

随机了一下，结果能够过5个点。

没办法了，只好模拟退火。临时学的，更确切的说是临时抄的。

 

模拟退火，一种启发式搜索方法，据说是模拟物理上的退火过程而发明的。

个人理解：

爬山算法就是用贪心的方法求出局部最优解。就是说你只要只走上坡路就一定可以走到一个山顶，但不一定是珠穆朗玛峰。

而模拟退火就是从一个点试着向四周蹦，如果蹦到点比这里高，那就上去（这个好理解，总忘高处走，这就是爬山算法嘛）；如果蹦到的点比较低，则有定的概率上去。

为什么一定要接受更低的点呢？只有这样才可以从局部最优的情况中跳出来，从而试着去爬更高的山。

 

关键为题在于概率如何计算。也就是随着时间的推移，接受的概率要逐渐降低。

代码如下

bool accept(double delta,double cur)
{
    if(delta<=0)return 1;　　　　//若果跳到的点更符合要求
4     return rand()<exp((-delta)/cur)*RAND_MAX;
5　　　//若果跳到的点更差....,cur为当前温度，它会随着时间降低（退火），RAND_MAX为最大随机整数，exp(负数)产生0-1间的实数，明白了吗？
}

 

还不明白就看大神的吧！

https://zhuanlan.zhihu.com/p/23968011

______________________________________________________________________________________________________________________________

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
int cx[21],cxf[21];
double x[21],y[21],dis[21][21];
double ans;
int n;
double getdis(int * s)
{
    double tp=0;
    for(int i=1;i<n;i++)tp+=dis[s[i-1]][s[i]];
    return tp;
}
bool accept(double delta,double cur)
{
    if(delta<=0)return 1;
    return rand()<exp((-delta)/cur)*RAND_MAX;
}
double solve()
{
    double dec=0.999;
    double tp,ttp=ans;
    double cur=10000;
    while(cur>0.01)
    {
        memcpy(cxf,cx,sizeof(cx));
        int a=rand()%n,b=rand()%n;
        swap(cxf[a],cxf[b]);
        tp=getdis(cxf);
        if(accept(tp-ttp,cur))
        {
            memcpy(cx,cxf,sizeof(cx));
            ttp=tp;
        }
        cur*=dec;
    }
    return ttp;
}
int main()
{
    srand(2323);
    scanf("%d",&n);
    for(int i=0;i<n;++i)
        scanf("%lf%lf",x+i,y+i);
    for(int i=0;i<n;++i)
        for(int j=0;j<i;++j)
        {
            dis[i][j]=dis[j][i]=hypot(x[i]-x[j],y[i]-y[j]);
        }
    for(int i=0;i<n;i++)cx[i]=i;
    ans=getdis(cx);
    int t=150;
    while(t--)
    {
        ans=min(ans,solve());
    }
    printf("%.2lf\n",ans);
    return 0;
}