做没有好的网站你懂的,阿里云企业邮箱收费标准,企业管理咨询服务是什么,网页设计作业样例题目链接 UVa12304 2D Geometry 110 in 1! 题意 这是一个拥有6#xff08;二进制是110#xff09;个子问题的2D几何问题集。 1 CircumscribedCircle x1 y1 x2 y2 x3 y3#xff1a;求三角形(x1,y1)-(x2,y2)-(x3,y3)的外接圆。这3点保证不共线。答案应格式化成(x,y,r…题目链接 UVa12304 2D Geometry 110 in 1! 题意 这是一个拥有6二进制是110个子问题的2D几何问题集。         1 CircumscribedCircle x1 y1 x2 y2 x3 y3求三角形(x1,y1)-(x2,y2)-(x3,y3)的外接圆。这3点保证不共线。答案应格式化成(x,y,r)表示圆心为(x,y)半径为r。         2、InscribedCircle x1 y1 x2 y2 x3 y3求三角形(x1,y1)-(x2,y2)-(x3,y3)的内切圆。这3点保证不共线。答案应格式化成(x,y,r)表示圆心为(x,y)半径为r。         3、TangentLineThroughPoint xc yc r xp yp给定一个圆心在(xc,yc)半径为r 的圆求过点(xp,yp)并且和这个圆相切的所有切线。每条切线格式化为angle表示直线的极角角度0≤angle180。整个答案应格式化为列表见后。如果无解应打印空列表。         4、CircleThroughAPointAndTangentToALineWithRadius xp yp x1 y1 x2 y2 r求出所有经过点(xp,yp)并且和直线(x1,y1)-(x2,y2)相切的半径为r 的圆。每个圆格式化为(x,y)因为半径已经给定。整个答案应格式化为列表同上。          5、CircleTangentToTwoLinesWithRadius x1 y1 x2 y2 x3 y3 x4 y4 r给出两条不平行直线(x1,y1)-(x2,y2)和(x3,y3)-(x4,y4)求所有半径为r 并且同时和这两条直线相切的圆。每个圆格式化为(x,y)因为半径已经给定。整个答案应格式化为列表同上。         6、CircleTangentToTwoDisjointCirclesWithRadius x1 y1 r1 x2 y2 r2 r给定两个相离的圆(x1,y1,r1)和(x2,y2,r2)求出所有和这两个圆外切的半径为r 的圆。注意因为是外切求出的圆不能把这两个给定圆包含在内部。每个圆格式化为(x,y)因为半径已经给定。整个答案应格式化为列表同上。         对于上述所有直线输入的两个点保证不重合。当格式化实数列表时所有数应从小到大排列当格式化二元组(x,y)时先按x 从小到大排序当x 相同时按y 从小到大排序。 分析 第一个问题设外接圆的圆心坐标为x、y半径为r列方程组求解。 第二个问题利用内心坐标公式求解设的顶点坐标分别为、、三个角对应的边长分别为a、b、c则内心坐标为。求得内心坐标后再求一下它与三角形某条边的距离即可得到内切圆半径。 第三个子问题切线格式化为极角要注意反三角函数的值域最后一定要确保在。 《训练指南》上讲第四、五、六个问题的解法 第四个问题因为已知半径为r所以要想和直线L相切圆心到直线的距离一定为r。满足这个条件的点的轨迹是两条直线。而要想过定点圆心点该点的距离一定为r满足这个条件的点的轨迹是一个圆。求出圆和这两条直线的交点即可。         第五个问题解法类似根据每条直线得到两条新直线再两两求交点即可。         第六个问题解法也类似根据两个圆得到两个新的圆求交点即可。 这里说一下第五、六两个问题的直接几何解法 第五个问题求出两直线的交点P以及pc与两条直线的夹角α则并且将某直线的方向向量旋转α角就得到pc的方向向量由此算出一个Cpc的方向向量再依次旋转90度三下可以算出另外三个C点要注意旋转90度和270度时pc长度变了。 第六个问题本质是求下图的C点坐标借助余弦定理把∠ A求出然后将AB的方向向量旋转∠ A就得到AC的方向向量进而算出C点坐标。 AC代码 #include iostream #include cmath #include iomanip #include algorithm using namespace std;#define eps 1e-10 struct Point {double x, y;Point(double x 0, double y 0) : x(x), y(y) {}void Normalize() {double l sqrt(x*x y*y); x / l; y / l;} }; typedef Point Vector;Vector operator (Vector A, Vector B) {return Vector(A.x B.x, A.y B.y); }Vector operator- (Point A, Point B) {return Vector(A.x - B.x, A.y - B.y); }Vector operator* (Vector A, double p) {return Vector(A.x * p, A.y * p); }Vector operator/ (Vector A, double p) {return Vector(A.x / p, A.y / p); }bool operator (const Point a, const Point b) {return a.x b.x || (a.x b.x a.y b.y); }int dcmp(double x) {return abs(x) eps ? 0 : (x 0 ? -1 : 1); }double Dot(Vector A, Vector B) {return A.x * B.x A.y * B.y; }double Length(Vector A) {return sqrt(Dot(A, A)); }double Angle(Vector A, Vector B) {return acos(Dot(A, B) / Length(A) / Length(B)); }double Cross(Vector A, Vector B) {return A.x * B.y - A.y * B.x; }Vector Rotate(Vector A, double rad) {return Vector(A.x * cos(rad) - A.y * sin(rad), A.x * sin(rad) A.y * cos(rad)); }Vector Normal(Vector A) {double L Length(A);return Vector(-A.y / L, A.x / L); }Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) {Vector u P - Q;double t Cross(w, u) / Cross(v, w);return P v * t; }double DistanceToLine(Point P, Point A, Point B) {Vector v1 B - A, v2 P - A;return abs(Cross(v1, v2)) / Length(v1); }void solveU() {double x1, y1, x2, y2, x3, y3; cin x1 y1 x2 y2 x3 y3;double x ((x1*x1y1*y1)*(y2-y3) (x2*x2y2*y2)*(y3-y1) (x3*x3y3*y3)*(y1-y2)) /((y1-y3)*(x1-x2) - (y1-y2)*(x1-x3)) / 2;double y ((x1*x1y1*y1)*(x2-x3) (x2*x2y2*y2)*(x3-x1) (x3*x3y3*y3)*(x1-x2)) /((x1-x3)*(y1-y2) - (x1-x2)*(y1-y3)) / 2;cout ( x , y , sqrt((x-x1)*(x-x1) (y-y1)*(y-y1)) ) endl; }void solveI() {Point p1, p2, p3; cin p1.x p1.y p2.x p2.y p3.x p3.y;double a Length(p2-p3), b Length(p3-p1), c Length(p1-p2); Point p (p1*ap2*bp3*c)/(abc);cout ( p.x , p.y , DistanceToLine(p, p1, p2) ) endl; }void solveT() {int x, y, r, xp, yp; cin x y r xp yp;int d (x-xp)*(x-xp) (y-yp)*(y-yp) - r*r;if (d 0.) {cout [] endl; return;} else if (d 0) {double a 180*atan2(yp-y, xp-x)/M_PI; a a0 ? a 450 : a 90;cout [ (dcmp(a-360)0 ? a-360 : (dcmp(a-180)0 ? a-180 : a)) ] endl; return;} else {double a 180*atan2(yp-y, xp-x)/M_PI, b 180*asin(r/sqrt((x-xp)*(x-xp) (y-yp)*(y-yp)))/M_PI;if (a 0) a 360;double x a b, y ab ? 180a-b : a-b;x dcmp(x-360)0 ? x-360 : (dcmp(x-180)0 ? x-180 : x);y dcmp(y-360)0 ? y-360 : (dcmp(y-180)0 ? y-180 : y);cout [ min(x, y) , max(x, y) ] endl;} }void solveA() {int xp, yp, x1, y1, x2, y2, r; cin xp yp x1 y1 x2 y2 r;if ((xp-x1)*(yp-y2) (xp-x2)*(yp-y1)) {double a atan2(x2-x1, y1-y2);Point p1(xp r*cos(a), yp r*sin(a)), p2(xp - r*cos(a), yp - r*sin(a)), t;if (p2 p1) t p1, p1 p2, p2 t;cout [( p1.x , p1.y ),( p2.x , p2.y )] endl;} else {Point a(x1, y1), b(x2, y2); Vector v Normal(Vector(x2-x1, y2-y1));double d DistanceToLine(Point(xp, yp), a, b);Point p(xp v.x * r, yp v.y * r);if (d 2*r) {cout [] endl;} else if (d 2*r) {if (DistanceToLine(p, a, b) 2*r) p.x xp - v.x * r, p.y yp - v.y * r;cout [( p.x , p.y )] endl;} else {double d1 DistanceToLine(p, a, b);if ((dr dcmp(d1d-r) 0) || (dr dcmp(d1r-d) 0)) v.x -v.x, v.y -v.y;double a atan2(v.y, v.x), b acos((r-d)/r);Point p1(xp r*cos(ab), yp r*sin(ab)), p2(xp r*cos(a-b), yp r*sin(a-b)), t;if (p2 p1) t p1, p1 p2, p2 t;cout [( p1.x , p1.y ),( p2.x , p2.y )] endl;}} }void solveL() {Point p1, p2, p3, p4; int r; cin p1.x p1.y p2.x p2.y p3.x p3.y p4.x p4.y r;Vector v p2-p1, w p4-p3; Point p GetLineIntersection(p1, v, p3, w); double a Angle(v, w) / 2;v.Normalize(); v Rotate(v, Cross(v, w) 0. ? -a : a);double r1 r / sin(a), r2 r / cos(a);Point c[] {Point(p.x v.x*r1, p.y v.y*r1), Point(p.x - v.x*r1, p.y - v.y*r1),Point(p.x - v.y*r2, p.y v.x*r2), Point(p.x v.y*r2, p.y - v.x*r2)};sort(c, c4);cout [( c[0].x , c[0].y ),( c[1].x , c[1].y ),( c[2].x , c[2].y ),( c[3].x , c[3].y )] endl; }void solveD() {int x1, y1, r1, x2, y2, r2, r; cin x1 y1 r1 x2 y2 r2 r;double d sqrt((x1-x2)*(x1-x2) (y1-y2)*(y1-y2)); Vector v(x2-x1, y2-y1); int cmp dcmp(d - r1 - r2 - 2*r);if (cmp 0 || (dcmp(dmin(r1,r2)-max(r1,r2)) 0)) {cout [] endl;} else if (cmp 0) {double t (r1r) / d;cout [( x1 v.x * t , y1 v.y * t )] endl;} else {double a acos(((rr1)*(rr1) d*d - (rr2)*(rr2)) / (rr1) / d / 2);v.Normalize(); Vector w Rotate(v, -a); v Rotate(v, a);Point p1(x1 v.x*(rr1), y1 v.y*(rr1)), p2(x1 w.x*(rr1), y1 w.y*(rr1)), t;if (p2 p1) t p1, p1 p2, p2 t;cout [( p1.x , p1.y ),( p2.x , p2.y )] endl;} }char s[48];void solve() {if (s[0] I) solveI();else if (s[0] T) solveT();else if (s[4] u) solveU();else if (s[13] A) solveA();else if (s[18] L) solveL();else solveD(); }int main() {cout fixed setprecision(6);while (cin s) solve();return 0; }