首页 > 代码库 > UVA - 10674-Tangents
UVA - 10674-Tangents
题意:给出两个圆,求它们的公切线,并按照一定格式输出
做法:模拟
代码:
#include<iostream> #include<map> #include<string> #include<cstring> #include<cstdio> #include<cstdlib> #include<cmath> #include<queue> #include<vector> #include<algorithm> using namespace std; const double eps=2e-5; const double pi=acos(-1.0); int dcmp(double x){return fabs(x)<eps?0:x<0?-1:1;} struct dot { double x,y; dot(){} dot(double a,double b){x=a;y=b;} dot operator +(dot a){return dot(x+a.x,y+a.y);} dot operator -(dot a){return dot(x-a.x,y-a.y);} dot operator *(double a){return dot(x*a,y*a);} double operator *(dot a){return x*a.y-y*a.x;} dot operator /(double a){return dot(x/a,y/a);} double operator /(dot a){return x*a.x+y*a.y;} bool operator ==(dot a){return x==a.x&&y==a.y;} void in(){scanf("%f%f",&x,&y);} void out(){printf("%f %f\n",x,y);} dot norv(){return dot(-y,x);} dot univ(){double a=mod();return dot(x/a,y/a);} dot ro(double a){return dot(x*cos(a)-y*sin(a),x*sin(a)+y*cos(a));} double mod(){return sqrt(x*x+y*y);} double dis(dot a){return sqrt(pow(x-a.x,2)+pow(y-a.y,2));} }; typedef pair<dot,dot> LL; LL ans[10]; bool operator <(LL a,LL b) { return dcmp(a.first.x-b.first.x)!=0?dcmp(a.first.x-b.first.x)<0: dcmp(a.first.y-b.first.y)<0; } int work(dot a,double r1,dot b,double r2) { int cnt; double ang; if(dcmp(a.x-b.x)==0&&dcmp(a.y-b.y)==0&&dcmp(r1-r2)==0) return -1; if(dcmp(a.dis(b)+min(r1,r2)-max(r1,r2))<0) return 0; if(dcmp(a.dis(b)+min(r1,r2)-max(r1,r2))==0) { if(dcmp(r1-r2)>0) ans[0].first=ans[0].second=a+(b-a).univ()*r1; else ans[0].first=ans[0].second=b+(a-b).univ()*r2; return 1; } ang=acos((r1-r2)/a.dis(b)); ans[0].first=a+(b-a).ro(ang).univ()*r1; ans[1].first=a+(b-a).ro(-ang).univ()*r1; ang=pi-ang; ans[0].second=b+(a-b).ro(-ang).univ()*r2; ans[1].second=b+(a-b).ro(ang).univ()*r2; cnt=2; if(dcmp(a.dis(b)-r1-r2)==0) { ans[2].first=ans[2].second=a+(b-a).univ()*r1; cnt=3; } if(dcmp(a.dis(b)-r1-r2)>0) { ang=acos((r1+r2)/a.dis(b)); ans[2].first=a+(b-a).ro(ang).univ()*r1; ans[3].first=a+(b-a).ro(-ang).univ()*r1; ans[2].second=b+(a-b).ro(ang).univ()*r2; ans[3].second=b+(a-b).ro(-ang).univ()*r2; cnt=4; } sort(ans,ans+cnt); return cnt; } int main() { int i,cnt; dot a,b; double r1,r2; while(cin>>a.x>>a.y>>r1>>b.x>>b.y>>r2) { if(a.x==0&&a.y==0&&r1==0&&b.x==0&&b.y==0&&r2==0) return 0; cnt=work(a,r1,b,r2); cout<<cnt<<endl; for(i=0;i<cnt;i++) printf("%.5f %.5f %.5f %.5f %.5f\n",ans[i].first.x,ans[i].first.y, ans[i].second.x,ans[i].second.y,ans[i].first.dis(ans[i].second)); } }
Description
Problem B
Tangents
Input: standard input
Output: standard output
Time Limit: 2 seconds
You can see in the pictures below that two different circles can have at most four common tangents. Given the center and radius of two circles your job is to find the length of their common tangents and also the points where they touch the two circles.
Input
The input file contains several lines of inputs.
Each line contains six integers x1 (-100<=x1<=100), y1 (-100<=y1<=100),r1 (0<r1<=200), x2 (-100<=x2<=100),y2 (-100<=y2<=100), r2 (0<r2<=200). Here(x1, y1) and (x2, y2) are the coordinates of the center of the first circle and second circle respectively,r1 is the radius of the first circle and r2 is the radius of the second circle.
Input is terminated by a line containing six zeroes.
Output
For each line of input you should produce one of more lines of output. The description of this output is given below.
First line of the output for each line of input contains an integern, which denotes the number of different tangents between the two circles. If there is infinite number of tangents between the two circles then the value ofn should be –1. If n is positive then nextn lines contains the description of each tangent. The description of the tangent contains five floating-point numbersSx, Sy, Tx,Ty, L in a single line. Here (Sx, Sy) is the point at which the tangent touches the first circle and(Tx, Ty) is the point where the tangent touches the second circle andL is the length of the tangent. All the floating-point numbers have five digits after the decimal point. Errors less than2*10-5 will be ignored. The tangents should be printed in ascending order ofSx and in case of a tie they should be printed in ascending order ofSy.
Sample Input Output for Sample Input
10 10 5 20 20 5 10 10 10 20 20 10 10 10 5 20 10 5 0 0 0 0 0 0 | 4 6.46447 13.53553 16.46447 23.53553 14.14214 10.00000 15.00000 20.00000 15.00000 10.00000 13.53553 6.46447 23.53553 16.46447 14.14214 15.00000 10.00000 15.00000 20.00000 10.00000 2 2.92893 17.07107 12.92893 27.07107 14.14214 17.07107 2.92893 27.07107 12.92893 14.14214 3 10.00000 5.00000 20.00000 5.00000 10.00000 10.00000 15.00000 20.00000 15.00000 10.00000 15.00000 10.00000 15.00000 10.00000 0.00000
|
Problem setter: Shahriar Manzoor. Member of Elite Problem Setters’ Panel
Special Thanks: Derek Kisman, Member of Elite Problem Setters’ Panel
Input
Output
Sample Input
Sample Output
Hint
Source
UVA - 10674-Tangents