网站开发与维护能做什么职业,wordpress在哪编辑模板,网站开发的研究计划书,qq电脑版HDU5120 - Intersection 做法#xff1a;答案就是\(A1 \cap B1 - A2 \cap B1 - A1 \cap B2 A2 \cap B2\)#xff0c;圆形面积交码了好久#xff0c;先把两个圆转到x轴上#xff0c;大的放到原点#xff0c;小的放在大圆的左边#xff0c;计算一下交点坐标#xff0c;然… HDU5120 - Intersection 做法答案就是\(A1 \cap B1 - A2 \cap B1 - A1 \cap B2 A2 \cap B2\)圆形面积交码了好久先把两个圆转到x轴上大的放到原点小的放在大圆的左边计算一下交点坐标然后讨论即可。 #include bits/stdc.h
#define pb push_back
typedef long long ll;
const double eps 1e-10;
const double PI acos(-1.0);
using namespace std;
int sgn(double x) {if(fabs(x) eps) return 0;if(x 0) return -1;return 1;
}
int n;
struct cir{double x,y,r;cir(){}cir(double _x,double _y,double _r) {x_x; y_y; r_r;}void output() {cout x y r endl;}
}a1, a2, b1, b2;
double dis(double ax,double ay,double bx,double by) {return sqrt((ax-bx)*(ax-bx)(ay-by)*(ay-by));
}
double area_t(double a,double b,double c) {if(ab c) return 0;if(ac b) return 0;if(cb a) return 0;if(abs(a-b) c) return 0;if(abs(c-b) b) return 0;if(abs(a-c) a) return 0;double p (abc)*0.5;return sqrt((p-a)*(p-b)*(p-c)*p);
}
void chg(cir a, cir b) {if(a.r b.r) swap(a,b);double D dis(a.x,a.y,b.x,b.y);a.x a.y 0;b.x D; b.y 0;
}
double cal(cir a, cir b) {chg(a, b);double D dis(a.x,a.y,b.x,b.y), ans 0;if(sgn(D- (a.rb.r)) 0) return 0;if(sgn(a.r-b.r-D)0) return PI*b.r*b.r;double cx (a.r*a.r D*D - b.r*b.r)*0.5/D;double cy sqrt(a.r*a.r - cx*cx);if(sgn(cx - b.x) 0) {double S1 D*cy;double tha asin(cy/a.r)*2;double thb asin(cy/b.r)*2;double Sa a.r*a.r*tha*0.5;double Sb b.r*b.r*thb*0.5;ans Sa Sb - S1;}else {double tha asin(cy/a.r)*2;double thb 2*PI-asin(cy/b.r)*2;double Sa a.r*a.r*tha*0.5 - area_t(a.r,a.r,cy*2);double Sb b.r*b.r*thb*0.5 area_t(b.r,b.r,cy*2);ans Sa Sb;}return ans;
}
int T, CC 0;
int main() {scanf(%d,T);while(T--) {double r,R,tx1,tx2,ty1,ty2;scanf(%lf%lf,r,R);scanf(%lf%lf,tx1,ty1);scanf(%lf%lf,tx2,ty2);a1 cir(tx1,ty1,R);a2 cir(tx1,ty1,r);b1 cir(tx2,ty2,R);b2 cir(tx2,ty2,r);double ans cal(a1,b1) - cal(a2,b1) - cal(a1,b2) cal(a2,b2);printf(Case #%d: %.6f\n,CC,ans);}
}转载于:https://www.cnblogs.com/RRRR-wys/p/9710574.html
