标签:
#include<iostream> #include<cmath> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; typedef double db; const db PI=acos(-1.0); struct point { db x,y; }p[2]; struct circle { point O; db r; }C[4]; db distcal(point a,point b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } db areacal(circle a,circle b) { if(a.r>b.r)swap(a,b); db d=distcal(a.O,b.O); if(d==0.0||d+a.r<=b.r){db t=a.r;return PI*t*t;} if(a.r+b.r<=d)return 0.0; db th1=acos((a.r*a.r+d*d-b.r*b.r)/(2*a.r*d)); db th2=acos((b.r*b.r+d*d-a.r*a.r)/(2*b.r*d)); db are1=a.r*a.r*th1; db are2=b.r*b.r*th2; db are3=a.r*d*sin(th1); return are1+are2-are3; } int main() { int T; scanf("%d",&T); for(int kase=1;kase<=T;kase++) { db r,R; scanf("%lf%lf",&r,&R); scanf("%lf%lf%lf%lf",&p[0].x,&p[0].y,&p[1].x,&p[1].y); C[0].O=p[0],C[0].r=R; C[1].O=p[0],C[1].r=r; C[2].O=p[1],C[2].r=R; C[3].O=p[1],C[3].r=r; db ans=areacal(C[0],C[2])+areacal(C[1],C[3])-2.0*areacal(C[0],C[3]); printf("Case #%d: %.6f\n",kase,ans); } return 0; } /* 1 1 10 0 0 1 0 */
标签:
原文地址:http://www.cnblogs.com/homura/p/5424228.html