HDOJ1071 定积分

#include"stdio.h"
int main()
{
    int t;
    double x1,y1,x2,y2,x3,y3,m,n,a,b,c,p,q,area;
    for(scanf("%d", &t); t; --t)
    {
        scanf("%lf%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&x3,&y3);
        m=(y2-y3)/(x2-x3);
        n=y2-m*x2;
        a=(x1*y3-x1*y2+x2*y1-x2*y3+x3*y2-x3*y1)/((x2-x1)*(x3-x1)*(x2-x3));
        b=(y2-y1)/(x2-x1)-a*(x2+ x1);
        c=y1-a*x1*x1-b*x1;
        p=x2>x3?x3:x2;
        q=x2>x3?x2:x3;
        area=a/3.0*(q*q*q-p*p*p)+(b-m)/2.0*(q*q-p*p)+(c-n)*(q-p);
        printf("%.2f\n",area);
    }
    return 0;
}