[BZOJ4920][Lydsy六月份月赛 .D][数学][三分]薄饼切割

可以把线段分成无数个点,这些点经过的面积就是答案,显然这些点的轨迹是无数个同心圆,只要把最大的圆包围的面积减去最小的圆包围的面积就是答案了

我是用三分求最小的圆的,这样比较方便,最大的圆肯定是线段的两个端点之一,比较一下就好了。
其他的画画图推推式子就好啦

#include 
#include 
#include 
#include 

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*rdouble 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)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)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;
}

你可能感兴趣的:(二分,&,三分,数学)