OJ Board

Code: Select all

/* @JUDGE_ID: MYID 190 C */

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define NWAY       1
#define ROUND      1
#define DEBUG      0
#define EPS        0.000001
#define ISNULL(x)  (((x) < EPS && (x) > -EPS) ? 1 : 0)

struct vect {
  double x, y;
};

struct line {
  double a, b, c;
};

void
paramtoeq (struct vect *p, struct vect *v, struct line *l)
{
  if (ISNULL(v->y)) {
    l->a = 0;
    l->b = 1;
    l->c = - (p->y);
  } else {
    l->a = 1;
    l->b = - ((v->x) / (v->y));
    l->c = - (l->b * p->y + p->x);
  }
}

void
intersection (struct line *a, struct line *b, struct vect *ip)
{
  ip->x = (b->c * a->b - a->c *  b->b) / (a->a * b->b - a->b * b->a);
  ip->y = - (b->c * a->a - a->c *  b->a) / (a->a * b->b - a->b * b->a);
}

int
main (void)
{
  struct vect p[3], c[2], d[2], cc;
  struct line l[2];
  double r, t, u;
  double g, h, i;

  for (;;)
    {
      if (scanf ("%lf%lf%lf%lf%lf%lf", &(p[0].x), &(p[0].y),
		 &(p[1].x), &(p[1].y), &(p[2].x), &(p[2].y)) < 6)
	break;
      c[0].x = (p[0].x + p[1].x) / 2.0;
      c[0].y = (p[0].y + p[1].y) / 2.0;
      c[1].x = (p[0].x + p[2].x) / 2.0;
      c[1].y = (p[0].y + p[2].y) / 2.0;
      d[0].x = p[1].y - p[0].y;
      d[0].y = - (p[1].x - p[0].x);
      d[1].x = p[2].y - p[0].y;
      d[1].y = - (p[2].x - p[0].x);
#if NWAY
      if (!ISNULL(d[0].y)) {
	u = d[0].x / d[0].y;
	t = (d[0].y * (c[0].x - c[1].x) - d[0].x * (c[0].y - c[1].y)) /
	  (d[1].x * d[0].y - d[0].x * d[1].y);
      } else {
	/* d[1].y will never be null */
	  t = (c[0].y - c[1].y) / d[1].y;
      }
      cc.x = c[1].x + t * d[1].x;
      cc.y = c[1].y + t * d[1].y;
#else
      paramtoeq (c, d, l);
      paramtoeq (c+1, d+1, l+1);
      intersection (l, l+1, &cc);
#endif
      r = hypot (cc.x - p->x, cc.y - p->y);
      g = - 2 * cc.x;
      h = - 2 * cc.y;
      i = cc.x * cc.x + cc.y * cc.y - r * r;
#if ROUND
      cc.x = floor (1000.0*cc.x + 0.5) / 1000.0;
      cc.y = floor (1000.0*cc.y + 0.5) / 1000.0;
      r = floor (1000.0*r + 0.5) / 1000.0;
      g = floor (1000.0*g + 0.5) / 1000.0;
      h = floor (1000.0*h + 0.5) / 1000.0;
      i = floor (1000.0*i + 0.5) / 1000.0;
#endif
      printf ("(x %c %.3lf)^2 + (y %c %.3lf)^2 = %.3lf^2\n", (cc.x < -EPS) ? '+' : '-',
	      fabs(cc.x), (cc.y < -EPS) ? '+' : '-', fabs(cc.y), r);
      printf ("x^2 + y^2 %c %.3lfx %c %.3lfy %c %.3lf = 0\n\n",
	      (g < -EPS) ? '-' : '+', fabs(g), (h < -EPS) ? '-' : '+',
	      fabs(h), (i < -EPS) ? '-' : '+', fabs(i));
    }
  return EXIT_SUCCESS;
}