第十二周 项目四 点、圆的关系（3）

问题及代码：

/*
*Copyright (c)2014,烟台大学计算机与控制工程学院
*All rights reserved.
*文件名称：Project.cpp
*作    者：chenqin.
*完成日期：2015年5月29日
*版 本 号：v1.0
*
*问题描述:与圆心相连的直线：给定一点p，其与圆心相连成的直线，会和圆有两个交点，如图。在上面定义的Point(点)类和Circle(圆)类基础上，设计一种方案，输出这两点的坐标。
*程序输入：略
*程序输出：略
*/
#include <iostream>
#include <cmath>
using namespace std;
class Circle;<span style="font-family: Consolas, 'Courier New', Courier, mono, serif; line-height: 18px;">//由于在Point中声明友元函数crossover_point中参数中用了Circle，需要提前声明</span>
class Point
{
public:
    Point(double xx=0,double yy=0);
    friend ostream& operator<<(ostream &output,Point &p);
    friend void crossover_point(Point &p,Circle &c,Point &p1,Point &p2);
protected:
    double x,y;
};
Point::Point(double xx,double yy):x(xx),y(yy) {}
ostream &operator<<(ostream &output,Point &p)
{
    output<<"("<<p.x<<","<<p.y<<")";
    return output;
}
class Circle:public Point
{
public:
    Circle(double xx=0,double yy=0,double r=0);
    friend ostream& operator<<(ostream &output,Circle &c);
    friend void crossover_point(Point &p,Circle&c,Point &p1,Point &p2);
protected:
    double rad;
};
Circle::Circle(double xx,double yy,double r):Point(xx,yy),rad(r) {}
ostream &operator<<(ostream &output,Circle &c)
{
    output<<"("<<c.x<<","<<c.y<<")"<<'\t'<<"rad="<<c.rad;
    return output;
}
void crossover_point(Point &p,Circle &c,Point &p1,Point &p2)
{
    p1.x=(c.x+sqrt(c.rad*c.rad/(1+((c.y-p.y)/(c.x-p.x))*((c.y-p.y)/(c.x-p.x)))));
    p2.x = (c.x - sqrt(c.rad*c.rad/(1+((c.y-p.y)/(c.x-p.x))*((c.y-p.y)/(c.x-p.x)))));
    p1.y = (p.y + (p1.x -p.x)*(c.y-p.y)/(c.x-p.x));
    p2.y = (p.y + (p2.x -p.x)*(c.y-p.y)/(c.x-p.x));
}
int main()
{
    Circle c1(3,2,4);
    Point p1(1,1),p2,p3;
    crossover_point(p1,c1,p2,p3);
    cout<<"点p1："<<p1<<endl;
    cout<<"圆c1："<<c1<<endl;
    cout<<"交点有："<<endl;
    cout<<p2<<"\t"<<p3<<endl;
    return 0;
}

运行结果：

 知识点总结：学会利用引用类型参数获得结果。