雷达安装UVA

时间:2015-05-29 02:41:52

标签: c++ algorithm greedy

我正在尝试了解UVA问题1193的示例解决方案:

问题陈述:

UVA 1193 - Radar Installation

解决方案:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <cstdlib>
#include <cmath>
#include <bitset>

using namespace std;

#define Max 10000
const double eps = 1e-10;

struct Interval {
    double st, en;
    Interval() {}
    Interval(double s, double e) {
        st = s, en = e;
    }
    bool operator < (const Interval &i) const {
        return (i.en == en) ? (st < i.st) : (en < i.en);
    }
};
long double x[Max], y[Max];
Interval inter[Max];

//bujhlam na baal
int main(void) {
    int n, d, testcase = 0;
    while(scanf("%d %d", &n, &d) == 2 && !(n == 0 && d == 0)) {
        for(int i = 0; i < n; i++)
            scanf("%Lf %Lf", &x[i], &y[i]);
        int count = 0, ok = true;
        for(int i = 0; i < n; i++) {
            if(d < y[i]) { // if island is out of radar radious
                ok = false; // that means at least one of the islands is not reachable results in -1
                break;
            } else {
                long double sqrtd = sqrt( d * d - y[i] * y[i] );
                inter[i] = Interval(x[i] - sqrtd, x[i] + sqrtd);
            }
        }
        if(!ok) {
            printf("Case %d: %d\n", ++testcase, -1);
            continue;
        }
        sort(inter, inter + n);

        for(int i = 0; i < n;) {
            int j;
            for(j = 0; j < n; j++) {
                if(inter[j].st > inter[i].en)
                    break;
            }
            i = j;
            count++;
        }
        printf("Case %d: %d\n", ++testcase, count);
    }
    return 0;
}

我似乎无法按照作者的方法来解决这个问题。让我卡住的部分如下所示:

long double sqrtd = sqrt( d * d - y[i] * y[i] );
inter[i] = Interval(x[i] - sqrtd, x[i] + sqrtd);

似乎作者正在使用毕达哥拉斯定理?我看不出这个的目的。

另外,为什么要使用排序?

sort(inter, inter + n);
有人可以赐教我吗?感谢。

1 个答案:

答案 0 :(得分:4)

关于你的第一个问题:

long double sqrtd = sqrt( d * d - y[i] * y[i] );
inter[i] = Interval(x[i] - sqrtd, x[i] + sqrtd);

这是计算可以覆盖岛屿i的雷达的范围。

            . (x,y)
           /|\
          / | \ d
 0 ______/__|__\________  
        A   x   B

所以,看看上图,覆盖位置(x,y)的岛屿,我们可以将雷达放在x - (d^2 - y^2)x + (d^2 - y^2)

的范围内

一些解释:

调用点A和B在Ox轴上与点(x,y)的距离等于d的两个点,所以我们有一个正方形三角形(A,(x,y),(x,0) )),使用毕达哥拉斯定理,我们可以很容易地计算出A和B的位置

A = x - (d^2 - y^2)

B = x + (d^2 + y^2)

关于你的第二个问题:

Also, why is sorting being used?

sort(inter, inter + n);

为了覆盖所有岛屿,我们需要开始将雷达从最左边的岛​​屿放到第二个最远的岛屿......直到我们覆盖所有岛屿。因此,我们可以贪婪地完成这个过程,将雷达放在第一岛的右端x + (d^2 - y^2),这有助于覆盖该岛右侧最大数量的岛屿,并继续这个过程与下一个岛屿发现岛屿直到最后。