我正在解决 uva 在线法官的问题 10263(铁路)(http://uva.onlinejudge.org/index.php?option=onlinejudge&page=show_problem&problem=1204)并为我的代码得到错误的答案。在与朋友讨论后(他也解决了问题并被接受),他查看了我的代码,它与他的非常相似(例如,所有几何函数都是相同的,因为我们正在作为一个团队训练ICPC 和使用相同的代码库)。
因此,在我开始尝试一次更改以查看问题所在之后,我发现了这种奇怪的行为。
我的初始代码(得到错误答案)是这个:
#include <iostream>
#include <cmath>
#include <cstdio>
#include <limits>
#define EPS 1e-9
using namespace std;
struct point {
double x, y;
point(double _x, double _y) { x = _x, y = _y; } //EDITED
point() { x = 0.0, y = 0.0; }
};
double dist(point p1, point p2) {
return hypot(p1.x - p2.x, p1.y - p2.y);
}
double distToLine(point p, point A, point B, point *c) {
double scale = (double) ((p.x - A.x) * (B.x - A.x) + (p.y - A.y) * (B.y - A.y)) / ((B.x - A.x) * (B.x - A.x) + (B.y - A.y) * (B.y - A.y));
c->x = A.x + scale * (B.x - A.x);
c->y = A.y + scale * (B.y - A.y);
return dist(p, *c);
}
double distToLineSegment(point p, point A, point B, point* c) {
if ((B.x-A.x) * (p.x-A.x) + (B.y-A.y) * (p.y-A.y) < EPS) {
c->x = A.x; c->y = A.y;
return dist(p, A);
}
if ((A.x-B.x) * (p.x-B.x) + (A.y-B.y) * (p.y-B.y) < EPS) {
c->x = B.x; c->y = B.y;
return dist(p, B);
}
return distToLine(p, A, B, c);
}
int main() {
int Mx, My;
while (cin >> Mx) {
int N, x, y;
cin >> My >> N;
point A, B, aux, res, M(Mx, My);
double dres = numeric_limits<double>::infinity(), d;
cin >> x >> y;
A = point(x,y);
for (int i=0; i<N; i++) {
cin >> x >> y;
B = point(x,y);
d = distToLineSegment(M, A, B, &aux);
if (d < dres) {
dres = d;
res = aux;
}
A = B;
}
printf("%.4f\n%.4f\n",res.x,res.y);
}
}
现在,在更改此行之后:
int Mx, My;
while (cin >> Mx) {
int N, x, y;
cin >> My >> N;
point A, B, aux, res, M(Mx, My);
至
point M;
while (cin >> M.x) {
int N, x, y;
cin >> M.y >> N;
point A, B, aux, res;
我的解决方案被接受了。谁能帮我理解阅读 int 值然后创建要点的问题是什么?在什么可能的原因下,这段代码会“行为不端”并产生与另一个不同的答案?另外,关于为什么“M”是唯一问题的任何想法?(我一直用和以前一样的方式读“A”和“B”,这并没有影响答案,只是我读“M”的方式必须改变)。
任何帮助我理解正在发生的事情的提示将不胜感激!
编辑:“将其改回错误的版本”时写错了构造函数,对不起。