HDU 1724 Ellipse积分再换成三角函数(周赛数论J题）

这题比较郁闷……在实验室的时候忘了积分了，大一的积分竟然忘了，唉，刚才找了大一的课本又把积分给复习了……终于巩固了一下……

用积分求面积：
y = b√1-x^2/a^2
S =∫[l 到 r]（2*b√(1-x^2/a^2)）dx    //这里刚开始没看课本的时候看别人的答案就不太懂，然后看了课本又问了宿舍的人，才知道别人的解释是错的，应该是S =∫[l 到 r]（b√(1-x^2/a^2)-(-b√(1-x^2/a^2)）dx   ,所以别人忘了乘2.
用三角换元：y=b*sinθ，x=a*cosθ
x=l对应θ1，x=r对应θ2，则变换为：
S =∫[θ1 到 θ2]{[bsinθ]}d(a*cosθ) = ∫[θ1 到 θ2]{-ab(sinθ)^2}dθ
利用cos2θ = 1-2*(sinθ)^2就很容易积分出来，积得到的答案是：∫[θ1 到 θ2]（a*b*(cos2θ-1)）dθ 

因为给的是L和r，而这两个值对应的是x，而x=a*cosθ,由课本中可知：θ1=acos(l/a),θ2=acos(r/a)

积分出来面积就是：a*b*(sin(2*acos(r/a))/2-sin(2*acos(l/a))/2-acos(r/a)+acos(l/a))

#include <iostream>
#include <map>
#include <deque>
#include <queue>
#include <stack>
#include <string>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <algorithm>
#include <map>
#include <set>
using namespace std;
int main()
{
    double a,b,l,r;
    int t;
    cin>>t;
    while(t--)
    {
        cin>>a>>b>>l>>r;
        printf("%.3f\n",a*b*(sin(2*acos(r/a))/2-sin(2*acos(l/a))/2-acos(r/a)+acos(l/a)));
    }
    return 0;
}