// UVa 438 - The Circumference of the Circle #include <iostream> #include <stdio.h> #include <cmath> using namespace std; struct point { double x, y; }; struct line { double a, b, c; }; line points_to_line(point p, point q) { line l; l.a = p.y - q.y; l.b = q.x - p.x; l.c = -l.a * p.x - l.b * p.y; return l; } line perpendicular(line l, point p) { line per; per.a = -l.b; per.b = l.a; per.c = -per.a * p.x - per.b * p.y; return per; } point intersection(line l1, line l2) { point p; p.y = (l2.c * l1.a - l1.c * l2.a) / (l1.b * l2.a - l2.b * l1.a); p.x = (l2.c * l1.b - l1.c * l2.b) / (l1.a * l2.b - l2.a * l1.b); return p; } double dist(point p, point q) { return sqrt((p.x - q.x) * (p.x - q.x) + (p.y - q.y) * (p.y - q.y)); } point middle(point p, point q) { point m = { (p.x + q.x) / 2, (p.y + q.y) / 2 }; return m; } int main() { point p1, p2, p3; while (cin >> p1.x >> p1.y >> p2.x >> p2.y >> p3.x >> p3.y) { line side1 = points_to_line(p1, p2); line side2 = points_to_line(p2, p3); line per1 = perpendicular(side1, middle(p1, p2)); line per2 = perpendicular(side2, middle(p2, p3)); point c = intersection(per1, per2); double d = dist(c, p1) * 2; double circun = d * M_PI; printf("%.2f\n", circun); } return 0; }
Sunday, June 7, 2015
UVa 438 - The Circumference of the Circle
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment