BZOJ 1845 Simpson积分

思路：

Simpson积分直接上  限制一下递归深度+精度就好了

（难以理解为什么这么多人写扫描线）

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
#define pr pair<double,double>
const int N=;
int n;
double inf=,eps=1e-,l=inf,r=-inf;
pr p[N];
struct Triangle{double x1,y1,x2,y2,x3,y3;}tr[N];
double jd(double x1,double y1,double x2,double y2,double x){
    if(abs(x1-x2)<eps)return y2-y1;
    return y1+(y2-y1)/(x2-x1)*(x-x1);
}
double Cut(double x){
    int top=;
    double t1,t2,last=-inf,tt=;
    for(int i=;i<=n;i++){
        if(tr[i].x1<x&&tr[i].x2<x&&tr[i].x3<x)continue;
        if(tr[i].x1>x&&tr[i].x2>x&&tr[i].x3>x)continue;
        if(tr[i].x2<x&&tr[i].x1<x)
            t1=jd(tr[i].x1,tr[i].y1,tr[i].x3,tr[i].y3,x),t2=jd(tr[i].x2,tr[i].y2,tr[i].x3,tr[i].y3,x);
        else t1=jd(tr[i].x1,tr[i].y1,tr[i].x3,tr[i].y3,x),t2=jd(tr[i].x1,tr[i].y1,tr[i].x2,tr[i].y2,x);
        p[++top]=make_pair(min(t1,t2),max(t1,t2));
    }
    sort(p+,p++top);
    for(int i=;i<=top;i++){
        if(p[i].first>last)last=p[i].second,tt+=p[i].second-p[i].first;
        else if(p[i].second>last)tt+=p[i].second-last,last=p[i].second;
    }return tt;
}
double Simpson(double l,double mid,double r,double Cl,double Cmid,double Cr,int deep){
    double Clmid=Cut((l+mid)/),Crmid=Cut((mid+r)/);
    double Ansl=(mid-l)*(Cl+*Clmid+Cmid)/,Ansr=(r-mid)*(Cr+*Crmid+Cmid)/,Ans=(r-l)*(Cl+*Cmid+Cr)/;
    if(deep>&&abs(Ansl+Ansr-Ans)<1e-)return Ans;
    else return Simpson(l,(l+mid)/,mid,Cl,Clmid,Cmid,deep+)+Simpson(mid,(mid+r)/,r,Cmid,Crmid,Cr,deep+);
}
int main(){
    scanf("%d",&n);
    for(int i=;i<=n;i++){
        scanf("%lf%lf%lf%lf%lf%lf",&tr[i].x1,&tr[i].y1,&tr[i].x2,&tr[i].y2,&tr[i].x3,&tr[i].y3);
        l=min(l,min(tr[i].x1,min(tr[i].x2,tr[i].x3)));
        r=max(r,max(tr[i].x1,max(tr[i].x2,tr[i].x3)));
    }
    for(int i=;i<=n;i++){
        if(tr[i].x3<tr[i].x2)swap(tr[i].x2,tr[i].x3),swap(tr[i].y2,tr[i].y3);
        if(tr[i].x2<tr[i].x1)swap(tr[i].x2,tr[i].x1),swap(tr[i].y2,tr[i].y1);
        if(tr[i].x3<tr[i].x2)swap(tr[i].x2,tr[i].x3),swap(tr[i].y2,tr[i].y3);
    }
    printf("%.2f\n",Simpson(l,(l+r)/,r,,Cut((l+r)/),,));

}