ural 1170 desert

http://acm.timus.ru/problem.aspx?space=1&num=1170

 

大体解释在《算法艺术与信息学竞赛（第二版）》205页。

这个题目需要推算出一次导数（求导数零点）还有二次导数（判断极大极小），数学演算能力有限啊。。。公式推导了好久。。。

 

设斜率为k，若穿过改矩形，则可以按照205页四种情况算出a，b，给a，b加上标号编程ai，bi表示第i个矩形的a，b值，若没有穿过这个矩形，则ai = bi = 0；

 

用ci表示第i个矩形的时间参数，在矩形中走过的总时间为 T1 = ∑ci * (ai + bi / k) * √(1 + k ^ 2)

 

L表示总路程长度，cc表示不在矩形中的时间参数，不在矩形中的总时间为 T2 = cc * (L - ∑(ai + bi / k)  * √(1 + k ^ 2))

 

总时间 y = T1 + T2;

 

可以抽取出关于ai的常数表达式，和关于bi的常数表达式

 

A = ai * (ci - cc)    B = bi * (ci - cc)

 

y = T1 + T2

   = (A + B / k) * √(1 + k ^ 2)

 

导数

dy = (A * k ^ 3 - B) / (√(1 + k ^ 2) * k ^ 3)

 

二次导数

ddy = (A * k ^ 3 + B * (3 * k ^ 2 + 2)) / (1 + k ^ 2) ^ (3 / 2)

 

代码中用 x1，x2，slp， slope_min 或 slope_max 代替斜率 k， 这些 x 和点坐标中的 x 需要加以区分，除了点结构体里面的x是表示点的 x 坐标，其余的 x 全部表示和上述函数 y 相关的自变量。

 

有一对全局变量 A, B

 

每次在斜率最相近的两个点间求出 A， B，然后在这两个点构成的区间 [x1, x2] 中求最小值

 

所有区间的最小值中最小的一个就是所求最小时间

 

 这道题网上搜不到可以AC的代码，自己debug了很久才成功了，编码能力有待提高。呵呵。。。

 

#include <cstdio> #include <cstdlib> #include <cmath> #include <cassert> #include <algorithm> using namespace std; #define RECT_MAX 500 struct point_t { double x, y; }; struct rect_t { point_t lower_left; point_t upper_right; point_t lower_right; point_t upper_left; }; rect_t rect[RECT_MAX]; point_t point[RECT_MAX * 4]; double c[RECT_MAX]; double cc; int point_sz; int rect_sz; double cbrt(double d) {//立方根 if(d < 0) return -pow(-d, 1.0 / 3.0); return pow(d, 1.0 / 3.0); } double slope(const point_t &a) {//斜率 return a.y / a.x; } struct point_slope_cmp { bool operator()(const point_t &a, const point_t &b) { return slope(a) < slope(b); } }; double A, B, L; #ifdef _DEBUG void print_point(const point_t p) { fprintf(stderr, "(%lf, %lf)/n", p.x, p.y); } void print_rect(const rect_t r) { fprintf(stderr, "lower_left : "); print_point(r.lower_left); fprintf(stderr, "lower_right : "); print_point(r.lower_right); fprintf(stderr, "upper_left : "); print_point(r.upper_left); fprintf(stderr, "upper_right : "); print_point(r.upper_right); } #endif void calc_A_B(double slp_min, double slp_max) { A = B = 0.0; double a, b; double slp = slp_min / 2 + slp_max / 2; //穷举所有矩形，计算全局 A, B for(int i=0; i<rect_sz; ++i) { if(slope(rect[i].upper_left) < slp || slope(rect[i].lower_right) > slp) { a = 0; b = 0; } else { if(slope(rect[i].lower_left) < slp) { if(slope(rect[i].upper_right) > slp) {//书 p.205 第三种情况 a = rect[i].upper_right.x - rect[i].lower_left.x; b = 0; } else {//书 p.205 第四种情况 a = -rect[i].upper_left.x; b = rect[i].upper_right.y; } } else { if(slope(rect[i].upper_right) < slp) {//书 p.205 第二种情况 a = 0; b = rect[i].upper_right.y - rect[i].lower_left.y; } else {//书 p.205 第一种情况 a = rect[i].upper_right.x; b = -rect[i].lower_left.y; } } } B += b * (c[i] - cc); A += a * (c[i] - cc); } } //y值 double y(double x) { return A * sqrt(x * x + 1.0) + B * sqrt(1.0 / x / x + 1.0) + cc * L; } //y一次导数 double dy(double x) { return A * x / sqrt(x * x + 1.0) - B / x / x / sqrt(x * x + 1.0); } //y二次导数 double ddy(double x) { return (A * x * x * x + B * (3.0 * x * x + 2)) / pow(x * x + 1, 1.5); } int solve() { sort(point, point + point_sz, point_slope_cmp()); double min_time = 1e18;//最小时间，即答案 double min_slp = 1e18;//最小时间对应斜率 for(int i=1; i<point_sz; ++i) { #ifdef _DEBUG print_point(point[i - 1]); print_point(point[i]); #endif double x1 = slope(point[i - 1]); double x2 = slope(point[i]); calc_A_B(x1, x2); double dy_zero_x = cbrt(B / A);//导数为 0 的点 x 值 double min_x; //下面求 y(x) 在区间 [x1, x2] 上的最小值 if(dy_zero_x < x1 || x2 < dy_zero_x) { if(dy(x1) < 0) min_x = x2; else min_x = x1; } else if(::ddy(dy_zero_x) > 0) { min_x = dy_zero_x; } else { if(y(x1) < y(x2)) min_x = x1; else min_x = x2; } //跟新最小时间及对应斜率 if(min_time > y(min_x)) { min_time = y(min_x); min_slp = min_x; } #ifdef _DEBUG fprintf(stderr, "min_x = %lf/ndif_zero_x = %lf/nddy(dif_zero_x) = %lf/nx1 = %lf, x2 = %lf/ny(x1) = %lf/ny(x2) = %lf/ny(min_x) = %lf/nA = %lf, B = %lf/ndy(x1) = %lf/n/n", min_x, dy_zero_x, ::ddy(dy_zero_x), x1, x2, y(x1), y(x2), y(min_x), A, B, dy(x1)); #endif } printf("%.6lf/n%.6lf %.6lf/n", min_time, L / sqrt(min_slp * min_slp + 1.0), L * min_slp / sqrt(min_slp * min_slp + 1.0)); return 0; } int main() { scanf("%d", &rect_sz); point_sz = rect_sz * 4; for(int i=0; i<rect_sz; ++i) { int k = i * 4; scanf("%lf%lf%lf%lf%lf", &point[k].x, &point[k].y, &point[k + 1].x, &point[k + 1].y, &c[i]); point[k + 2].x = point[k].x; point[k + 2].y = point[k + 1].y; point[k + 3].x = point[k + 1].x; point[k + 3].y = point[k].y; rect[i].lower_left = point[k]; rect[i].upper_right = point[k + 1]; rect[i].lower_right = point[k + 3]; rect[i].upper_left = point[k + 2]; } scanf("%lf%lf", &cc, &L); solve(); return 0; }