HDU 4752 Polygon（抛物线长度积分）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4752

题意：给出一个抛物线和一个简单多边形。求抛物线在多边形内部的长度。

思路：首先求出多边形所有边和抛物线的交点。这里要注意：

（1）边和抛物线相切不算，因为整个抛物线在边的另一侧；

（2）如下，上面两种相交都是不能算的，下面一种要算。

算出所有交点之后，那么每相邻两个交点与抛物线的关系必然是抛物线交替进入-离开。下面对于一段[x1,x2]计算抛物线的长度。

struct point

{

    double x,y;



    void get()

    {

        cin>>x>>y;

    }



    point(){}

    point(double _x,double _y)

    {

        x=_x;

        y=_y;

    }



    point operator+(point a)

    {

        return point(x+a.x,y+a.y);

    }



    point operator-(point a)

    {

        return point(x-a.x,y-a.y);

    }



    point operator*(double t)

    {

        return point(x*t,y*t);

    }

};



point p[N];

int n;

double L,R,a,b,c;





int sgn(double x)

{

    if(x>1e-8) return 1;

    if(x<-1e-8) return -1;

    return 0;

}



vector<double> V;



double func(double x)

{

    return a*x*x+b*x+c;

}





void deal(point p,point q)

{



    if(sgn(p.x-q.x)==0)

    {

        double y=func(p.x);

        if(p.y>q.y) swap(p,q);

        if(sgn(y-p.y)>0&&sgn(y-q.y)<0) V.pb(p.x);

        return;

    }



    double k=(q.y-p.y)/(q.x-p.x);

    double A=a;

    double B=b-k;

    double C=c+k*p.x-p.y;

    double temp=B*B-4*A*C;

    if(sgn(temp)<=0) return;

    temp=sqrt(temp);

    double x1=(-B-temp)/(2*A);

    double x2=(-B+temp)/(2*A);

    if(p.x>q.x) swap(p,q);

    if(sgn(x1-p.x)>0&&sgn(x1-q.x)<0) V.pb(x1);

    if(sgn(x2-p.x)>0&&sgn(x2-q.x)<0) V.pb(x2);

}





void deal(point p1,point p,point p2)

{

    if(sgn(func(p.x)-p.y)) return;

    p1=p+(p1-p)*1e-3;

    p2=p+(p2-p)*1e-3;

    if(sgn(func(p1.x)-p1.y)*sgn(func(p2.x)-p2.y)==-1) V.pb(p.x);

}



double cal(double x)

{

    double A=4*a*a;

    double B=4*a*b;

    double C=b*b+1;



    double temp=A*x*x+B*x+C;

    return (2*A*x+B)/(4*A)*sqrt(temp)+(4*A*C-B*B)/(8*A*sqrt(A))*log(fabs(2*A*x+B+2*sqrt(A)*sqrt(temp)));

}



double cal(double l,double r)

{

    if(l<L) l=L;

    if(r>R) r=R;

    if(l>=r) return 0;

    return cal(r)-cal(l);

}





int main()

{

    Rush(n)

    {

        scanf("%lf%lf%lf%lf%lf",&a,&b,&c,&L,&R);

        int i;

        FOR0(i,n) p[i].get();

        p[n]=p[0]; p[n+1]=p[1];

        V.clear();

        FOR0(i,n) deal(p[i],p[i+1]);

        FOR1(i,n) deal(p[i-1],p[i],p[i+1]);

        sort(V.begin(),V.end());

        double ans=0;

        for(i=0;i+1<SZ(V);i+=2) ans+=cal(V[i],V[i+1]);

        PR(ans);

    }

}