poj 2420 广义费马点

题目概述
给定一个N边形所有顶点座标x,y,求其费马点到所有顶点距离和 
费马点是指到多边形所有顶点距离和最小的点

时限
1000ms/3000ms

输入
第一行正整数N,其后N行,每行两个整数x,y,输入到EOF为止

限制
1<=N<=100;0<=x,y<=10000

输出
每行一个数,为所求距离和,精确到整数

样例输入

0 0 
0 10000 
10000 10000 
10000 0

样例输出
28284

讨论
计算几何?沾点边,主要思想是模拟退火?或者叫随机化变步长贪心?额感觉还是叫瞎猜比较合适,从给定的多边形顶点中随便找一个点开始,沿着上下左右四个方向移动相同距离,并分别计算到所有顶点距离和,取最小的一个(而且要比上一次小),记录移动后点的位置为新的起点,然后减小移动距离,再向四个方向移动,直到移动距离小于精度就可以返回距离和了 
说是模拟退火,但其实只有几个for循环,实际的模拟退火比这个复杂的多,说是随机化变步长贪心,额只看见贪心了,没看见rand(),所以基本上只能叫瞎猜了 
其实这个题用牛顿迭代法一样可解,不过需要一点高数的基础,而且写起来会比较麻烦,于是就没考虑 
看讨论版说数据弱的出奇,看出来了,瞎猜都能猜对

题解状态
156K,0MS,C++,846B

题解代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0x3f3f3f3f
#define MAXN 102
#define memset0(a) memset(a,0,sizeof(a))

int N;//顶点总数
int x[MAXN], y[MAXN];//存放顶点座标
const int dx[4] = { -1,0,1,0 }, dy[4] = { 0,-1,0,1 };//点移动方向向量
double fun()
{
    for (int p = 0; p < N; p++)
        scanf("%d%d", &x[p], &y[p]);//input
    double x0 = x[0], y0 = y[0], least = INF;//取得第一个点 初始化最短距离和
    for (double p = 10000; p > 0.4; p /= 2)//由于座标范围是0到10000 而且没有斜向移动 所以步长从10000开始 终点0.4是为了避免舍入误差 每次除以2有点太快 应该慢一些才是
        for (int i = 0; i < 4; i++) {
            double dis = 0, x1 = x0 + dx[i] * p, y1 = y0 + dy[i] * p;//初始化距离和 并算出移动后的点座标
            for (int u = 0; u < N; u++)
                dis += sqrt((x1 - x[u])*(x1 - x[u]) + (y1 - y[u])*(y1 - y[u]));
            if (dis < least) {//当发现移动后更近时就要更新 这就是贪心
                least = dis;
                x0 = x1, y0 = y1;
            }
        }
    return least;
}
int main(void)
{
    //freopen("vs_cin.txt", "r", stdin);
    //freopen("vs_cout.txt", "w", stdout);

    while (~scanf("%d", &N))//input
        printf("%.0lf\n", fun());//output
}

——————— 
作者:liurui39660 
来源:CSDN 
原文:https://blog.csdn.net/liurui39660/article/details/52127826 
版权声明:本文为博主原创文章,转载请附上博文链接!

 

点赞