POJ 2420 A Star not a Tree?

Description
Luke wants to upgrade his home computer network from 10mbs to 100mbs. His existing network uses 10base2 (coaxial) cables that allow you to connect any number of computers together in a linear arrangement. Luke is particulary proud that he solved a nasty NP-complete problem in order to minimize the total cable length.
Unfortunately, Luke cannot use his existing cabling. The 100mbs system uses 100baseT (twisted pair) cables. Each 100baseT cable connects only two devices: either two network cards or a network card and a hub. (A hub is an electronic device that interconnects several cables.) Luke has a choice: He can buy 2N-2 network cards and connect his N computers together by inserting one or more cards into each computer and connecting them all together. Or he can buy N network cards and a hub and connect each of his N computers to the hub. The first approach would require that Luke configure his operating system to forward network traffic. However, with the installation of Winux 2007.2, Luke discovered that network forwarding no longer worked. He couldn’t figure out how to re-enable forwarding, and he had never heard of Prim or Kruskal, so he settled on the second approach: N network cards and a hub.
Luke lives in a loft and so is prepared to run the cables and place the hub anywhere. But he won’t move his computers. He wants to minimize the total length of cable he must buy.

---

【题目分析】
精度要求不高，考虑一下模拟退火。确定一个初始的速度和退火的速度，然后就是瞎算了，随机性算法，只要能A，怎么都可以。

---

【代码】

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
const int maxn=105;
using namespace std;
int n;
int xi[maxn],yi[maxn];
double sum(double x,double y)
{
    double ans=0;
    for(int i=1;i<=n;i++)
    {
        ans+=sqrt((x-xi[i])*(x-xi[i])+(y-yi[i])*(y-yi[i]));
    }
    return ans;
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=1;i<=n;i++)scanf("%d%d",&xi[i],&yi[i]);
        double step=100000;
        double x=xi[1],y=yi[1];
        while(step>1e-9)
        {
            double rad=rand()%360+1;
            double nx=x+step*cos(rad);
            double ny=y+step*sin(rad);
            if(sum(nx,ny)<sum(x,y))
            {
                x=nx;
                y=ny;
            }
            step*=0.99;
        }
        printf("%.0f\n",sum(x,y));
    }
    return 0;
}