hdu 4752

计算几何+数值计算的题目；

要用到辛普森积分，没有学过~~~

参考学习了acm_Naruto大神 的代码！

代码：

 #include<cstdio>
 #include<cmath>
 #include<algorithm>
 #define maxn 20005
 using namespace std;
 int n,a,b,c,l,r,x[maxn],y[maxn];
 double sqr(double x)
 {
     return x*x;
 }
 double len(double x)
 {
     double t,tt;
     t=(sqrt((b+*a*x)*(b+*a*x)+)+*a*x+b);
     tt=log(t)/(*a);
     t=sqrt((b+*a*x)*(b+*a*x)+)*(b+*a*x)/(*a);
     tt+=t;
     return tt;
 }

 int main()
 {
     while(scanf("%d%d%d%d%d%d",&n,&a,&b,&c,&l,&r)!=EOF)
     {
         for(int i=; i<n; i++)
             scanf("%d%d",&x[i],&y[i]);
         x[n]=x[],y[n]=y[];
         double ans=;
         for(int i=; i<n; i++)if(x[i]!=x[i+])
             {
                 double x0=x[i],y0=y[i];
                 double x1=x[i+],y1=y[i+];
                 if(x0>x1)swap(x0,x1),swap(y0,y1);
                 double k=(y1-y0)/(x1-x0);
                 double d=y1-k*x1;
                 double u=sqr(b-k)-*a*(c-d);//b*b-4*a*c;
                 double length=;
                 if(u>=)
                 {
                     double t1=(-(b-k)*1.0+sqrt(u))/(2.0*a);//方程的根;
                     double t2=(-(b-k)*1.0-sqrt(u))/(2.0*a);//方程的根;
                     if(t1>t2)swap(t1,t2);
                     t1=max(x0,max(t1,l*1.0));
                     t2=min(x1,min(t2,r*1.0));
                     if(t1<=t2)
                         length=fabs(len(t2)-len(t1));//辛普森积分；
                     if(x[i]<x[i+])ans+=length;
                     else ans-=length;
                 }
             }
         printf("%.2lf\n",fabs(ans));
     }
     return ;
 }