[HDU1724] Ellipse(심슨 포인트)

2909 단어 문제풀이계산 기하학

제목 설명


전송문

문제풀이


심슨 포인트 누드 문제 답은 2를 곱해야지 타원이 반쪽이니까

코드

#include
#include
#include
#include
#include
using namespace std;

const double eps=1e-10;
int dcmp(double x)
{
    if (x<=eps&&x>=-eps) return 0;
    return (x>0)?1:-1;
}
int T;
double a,b,l,r,ans;

double F(double x)
{
    return b*sqrt(1-((x*x)/(a*a)));
}
double calc(double l,double r)
{
    return (F(l)+F(r)+F((l+r)/2)*4)*(r-l)/6;
}
double simpson(double l,double r,double now)
{
    double mid=(l+r)/2;
    double p=calc(l,mid),q=calc(mid,r);
    if (dcmp(now-p-q)==0) return now;
    else return simpson(l,mid,p)+simpson(mid,r,q);
}
int main()
{
    scanf("%d",&T);
    while (T--)
    {
        scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
        ans=simpson(l,r,calc(l,r))*2;
        printf("%.3lf
",ans);
    }
}
6. 그래프
terraform-provider-mackerel 2초(+)로 URL 외형 모니터링α)개시하다

좋은 웹페이지 즐겨찾기

개발자 우수 사이트 수집

개발자가 알아야 할 필수 사이트 100선 추천 우리는 당신을 위해 100개의 자주 사용하는 개발자 학습 사이트를 정리했습니다
best100-homepage