可以把线段分成无数个点,这些点经过的面积就是答案,显然这些点的轨迹是无数个同心圆,只要把最大的圆包围的面积减去最小的圆包围的面积就是答案了
我是用三分求最小的圆的,这样比较方便,最大的圆肯定是线段的两个端点之一,比较一下就好了。
其他的画画图推推式子就好啦
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
const double eps=1e-9,pi=acos(-1);
int t;
namespace work{
double h,w,x1,x2,y1,y2;
inline double dis(double x,double y){
return x*x+y*y;
}
inline double calc(double r){
if(2*r<=h) return pi*r*r;
if(r*r>dis(h/2.0,w/2.0)) return h*w;
double ret=sqrt(r*r-h*h/4.0)*h/4.0;
if(2*r<w){
double ang=pi/2.0-acos(h/2.0/r);
return 4*(ret+pi*r*r*ang/(2*pi));
}
else{
ret+=sqrt(r*r-w*w/4.0)*w/4.0;
double ang=pi/2.0-acos(h/2.0/r)-acos(w/2.0/r);
return 4*(ret+pi*r*r*ang/(2*pi));
}
}
void main(){
scanf("%lf%lf%lf%lf%lf%lf",&h,&w,&x1,&y1,&x2,&y2);
if(h>w) swap(h,w);
double maxr,minr;
if(dis(x1,y1)-dis(x2,y2)<eps) maxr=sqrt(dis(x2,y2));
else maxr=sqrt(dis(x1,y1));
while(fabs(x1-x2)>eps||fabs(y1-y2)>eps){
double x3=(x2-x1)/3+x1,x4=(x2-x1)/3*2+x1,
y3=(y2-y1)/3+y1,y4=(y2-y1)/3*2+y1;
if(dis(x3,y3)-dis(x4,y4)<eps) x2=x4,y2=y4;
else x1=x3,y1=y3;
}
minr=sqrt(dis(x1,y1));
printf("%.10lf\n",calc(maxr)-calc(minr));
}
}
int main(){
scanf("%d",&t);
while(t--) work::main();
return 0;
}