B - A Star not a Tree?

Time Limit: 20 Sec

Memory Limit: 256 MB

题目连接

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=88808#problem/B

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.

Input

The first line of input contains a positive integer N <= 100, the number of computers. N lines follow; each gives the (x,y) coordinates (in mm.) of a computer within the room. All coordinates are integers between 0 and 10,000.

Output

Output consists of one number, the total length of the cable segments, rounded to the nearest mm.

Sample Input

4
0 0
0 10000
10000 10000
10000 0

Sample Output

28284

HINT

题意

平面上有100个点,让你找到一个点,使得这些点到这个点的距离和最小

就是让你找广义费马点

题解

爬山算法

有一个蠢萌兔子,喝醉了,四处走,然后清醒了,就滚回去了,就是这样……

代码:

//qscqesze
#include <cstdio>
#include <cmath>
#include <cstring>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <set>
#include <bitset>
#include <vector>
#include <sstream>
#include <queue>
#include <typeinfo>
#include <fstream>
#include <map>
#include <stack>
typedef long long ll;
using namespace std;
//freopen("D.in","r",stdin);
//freopen("D.out","w",stdout);
#define sspeed ios_base::sync_with_stdio(0);cin.tie(0)
#define maxn 200051
#define mod 10007
#define eps 1e-9
int Num;
//const int inf=0x7fffffff; //нчоч╢С
const int inf=0x3f3f3f3f;
inline ll read()
{
ll x=,f=;char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
//************************************************************************************** int n;
struct node
{
double x,y;
};
double sqr(double x)
{
return x*x;
}
double dis(double x,double y,node p)
{
return sqrt(sqr(x-p.x)+sqr(y-p.y));
}
node p[];
double get_sum(double x,double y)
{
double ans=;
for(int i=;i<n;i++)
{
ans+=dis(x,y,p[i]);
}
return ans;
} int main()
{
n=read();
for(int i=;i<n;i++)
scanf("%lf%lf",&p[i].x,&p[i].y);
node ans;
ans.x=,ans.y=;
for(int i=;i<n;i++)
ans.x+=p[i].x,ans.y+=p[i].y;
ans.x=(ans.x/(n*1.0)),ans.y=(ans.y/(n*1.0));
double Ans = get_sum(ans.x,ans.y);
double t = ;
double x,y,tmp;
while(t>0.02)
{
x=,y=;
for(int i=;i<n;i++)
{
x+=(p[i].x-ans.x)/dis(ans.x,ans.y,p[i]);
y+=(p[i].y-ans.y)/dis(ans.x,ans.y,p[i]);
}
tmp = get_sum(ans.x+x*t,ans.y+y*t);
if(tmp<Ans)
{
Ans = tmp;
ans.x+=x*t;
ans.y+=y*t;
}
t*=0.9;
}
printf("%.0lf\n",Ans);
}

最新文章

  1. 通过NPOI操作Excel
  2. 用hibernate tools生成对应的sql应用代码
  3. java 导出数据为word文档(保持模板格式)
  4. [已解决] 快速理解RSA算法
  5. C++设计模式-Command命令模式
  6. html5 如何进行自定义属性的定义和查询
  7. Delphi中的函数指针判断是否为空
  8. Cookie工具类
  9. Navicat数据存放位置和备份数据库路径设置
  10. linux查看磁盘使用情况
  11. 九度OJ 1209 最小邮票数 -- 动态规划
  12. 提高你的Java代码质量吧:使用valueof前必须进行校验
  13. Beego学习笔记——Config
  14. ckplayer 项目实战
  15. 微信公众号开发总结(Node.js + express + winston)
  16. java远程备份mysql数据库关键问题(限windows环境,亲测解决)
  17. Datatable get请求传参应用
  18. python 中增加css样式的三种方式
  19. 彻头彻尾理解 LinkedHashMap
  20. timeStamp(时间戳) 事件属性

热门文章

  1. 动态加载so文件
  2. 给table中某一列的文字右对齐
  3. Android-监听sdcard状态
  4. 查询MySQL锁等待的语句
  5. HDU-4882 ZCC Loves Codefires
  6. 【原】Mongodb相关资料
  7. C++11 之 &quot; = delete &quot;
  8. makefile、gdb使用记录
  9. SRM 508 DIV1 500pt(DP)
  10. 最近的bug列表总结(C++)