#include using namespace std; double calc(double r1,double r2,double d){ const double pi=3.14159265358979323846; if(r1+d<=r2){ return r1*r1*pi; } if(r2+d<=r1){ return r2*r2*pi; } if(r2+r1<=d){ return 0; } double d1,d2,h; d1=(d+(r1*r1-r2*r2)/d)/2; d2=(d-(r1*r1-r2*r2)/d)/2; h=sqrt(r2*r2-d2*d2); return r1*r1*asin(h/r1)+r2*r2*asin(h/r2)-d2*h-d1*h; } int main(){ int N; double d; cin>>N>>d; vector dist(N+1); for(int i=0;i>dist[i]; } vector rightsum(N),leftsum(N),rightmax(N),leftmax(N); rightsum[0]=dist[0]; rightmax[0]=dist[0]; for(int i=1;i=0;i--){ leftsum[i]=leftsum[i+1]+dist[i+1]; leftmax[i]=max(leftmax[i+1],dist[i+1]); } double a,b,am,bm; for(int i=0;i