題目鏈接 Problem Description You have got a cylindrical cup. Its bottom diameter is 2 units and its height is 2 units as well.The height of liquid level i ...
Problem Description You have got a cylindrical cup. Its bottom diameter is 2 units and its height is 2 units as well.
The height of liquid level in the cup is d (0 ≤ d ≤ 2). When you incline the cup to the maximal angle such that the liquid inside has not been poured out, what is the area of the surface of the liquid?
Input The first line is the number of test cases. For each test case, a line contains a float-point number d.
Output For each test case, output a line containing the area of the surface rounded to 5 decimal places.
Sample Input 4 0 1 2 0.424413182
Sample Output 0.00000 4.44288 3.14159 3.51241 題意：現在有一個圓柱形的水杯，底面直徑為2，水杯高為2，現在水杯中裝有高為d的水，求將水杯傾斜到最大角度（即水到杯口時，水不能流出）時，水面的面積？ 思路：由高中學過的幾何知識可以知道：水面相當於對於一個很長的圓柱體傾斜的用刀切開，那麼這個切麵就是一個完整的橢圓，當然如果不傾斜則得到特殊的橢圓——圓，如果水面經過杯底，那麼水面就是一個缺少一部分的橢圓，所以我們需要分開討論水面經過杯底和不經過杯底兩種情況。 那麼這兩種情況的d的臨界值是多少呢？ 可以發現對於水剛到杯底的時候，有水和無水的部分各占一半，所以分界點d=1； 1、水面不經過杯底（d>=1）
這種情況如上圖所示，h+(2-h)/2=d, 所以h=2*d-2,那麼可以求出水面這個完整橢圓的長半徑a=sqrt(2*2+(2-h)*(2-h))/2,而橢圓的短半經是b=1，所以水面面積為S=PI*a*b.
2、水面經過杯底（d<1）
上圖即是水面經過杯底的樣子（我畫的不是水杯傾斜，直接讓水傾斜了，湊合看吧）。
對於上圖中mid越大則水的體積越大，那麼我們可以根據體積二分mid求出mid真實長，最後根據真實的mid求出水面的面積，可以知道二分範圍為（0，2）。每次我們需要根據當前mid求出水的體積，因為水體不規則所以必須積分求水的體積。
積分：我們根據水的高度積分，如上圖所示利用 y 積分求體積，那麼我們需要根據 y 求出每個水體截面的長 t （類似於杯底的mid），求相似三角形 t/mid=(2-y)/2 得 t=(2-y)*mid/2 ; 然後根據水截面長 t 求出當前水體截面的面積 S；
可以知道水體截面為一個扇形減去一個三角形組成，面積及體積如下圖所示：
求出真實的mid以後，那麼就可以求出水面的面積了。
如上圖所示利用二分求出的mid得 len=sqrt(2*2+mid*mid)，設水面與杯底的一個交點為（x，h），h=sqrt(1-(1-mid)*(1-mid)) (由杯底水面交線所在圓很容易求出 h ),那麼len=a-x 且x^2/a^a+h^2/b^2=1 ,解上面兩個方程 求得a= a=len/(1+（flag表示正負）flag*sqrt(1-h*h))，x=a-len 可以發現mid<1時水面為小半個橢圓，這時x>0,所以flag<0,反之mid>1時，水面為大半個橢圓，這時x<0,flag>0 。
現在橢圓方程已求出，X範圍（x，a）（不是完整的橢圓），所以需要積分，如下圖所述：
面積終於求出來了，太驚喜了！
那麼現在就該上代碼了：
#include <iostream> #include <algorithm> #include <cstdio> #include <cmath> using namespace std; const double PI=acos(-1); const double eps=1e-12; double cal2(double x) { double sum=x+sin(2*x)/2; return sum; } double area(double a,double x) { double sum=cal2(0.5*PI); sum=sum-cal2(asin(x/a)); return sum*a; } double cal(double x) { double ans=sin(x)-sin(x)*sin(x)*sin(x)/3-x*cos(x); return ans; } double getV(double mid) { double V=cal(acos(1))-cal(acos(1-mid)); V=V*(-2)/mid; return V; } int main() { int T; cin>>T; while(T--) { double d; scanf("%lf",&d); if(d>1) { double h=2*d-2; double a=sqrt(4+(2-h)*(2-h))/2; printf("%.5f\n",PI*a); continue; } double l=0.0,r=2.0; for(int i=0;i<50;i++) { double mid=(l+r)/2; //double du=acos(1-mid); //double S=du-sin(du)*(1-mid); double V=getV(mid); if(fabs(V-PI*d)<eps) break; if(V>PI*d) r=mid; else l=mid; } double mid=l; if(mid==0.0) mid=eps; int flag=1; if(mid<1) flag=-1; double len=sqrt(mid*mid+4); double h=sqrt(1-(1-mid)*(1-mid)); double a=len/(1+flag*sqrt(1-h*h)); double x=a-len; double res=area(a,x); printf("%.5f\n",res); } return 0; }
