bzoj 1845: [Cqoi2005] 三角形面积并

https://www.lydsy.com/JudgeOnline/problem.php?id=1845

将所有三角形的端点、交点按x坐标排序,从左往右扫描线

每相邻两根扫描线a和b之间的形状是若干个不相交的梯形或三角形

用公式:中位线长度*高 计算面积

高就是两条扫描线之间的距离

中位线,计算x=(xa+xb)/2 被覆盖了多长即可

因为已将三角形端点考虑在内,所以线段相交只需要考虑规范相交即可

直接套lrj的板子

#include<cmath>
#include<cstdio>
#include<algorithm>

using namespace std;

const double eps=1e-10;
const int inf=2e6;

int n;
struct Point
{
    double x,y;

    Point(double x_=0,double y_=0) : x(x_),y(y_) {}
}tri[101][4],seg[101];

typedef Point Vector;

double line[100001];

Vector operator - (Point A,Point B) { return Vector(A.x-B.x,A.y-B.y); }
Vector operator + (Vector A,Vector B) { return Vector(A.x+B.x,A.y+B.y); }
Vector operator * (Vector A,double p) { return Vector(A.x*p,A.y*p); } 

int dcmp(double x)
{
    if(fabs(x)<eps) return 0;
    else return x<0 ? -1 : 1;
}

double Cross(Vector A,Vector B)
{
    return A.x*B.y-A.y*B.x;
}

bool SegmentProperIntersection(Point a1,Point a2,Point b1,Point b2)
{
    double c1=Cross(a2-a1,b1-a1),c2=Cross(a2-a1,b2-a1),c3=Cross(b2-b1,a1-b1),c4=Cross(b2-b1,a2-b1);
    return dcmp(c1)*dcmp(c2)<0 && dcmp(c3)*dcmp(c4)<0;
}

Point GetLineIntersection(Point P,Vector v,Point Q,Vector w)
{
    Vector u=P-Q;
    double t=Cross(w,u)/Cross(v,w);
    return P+v*t;
}

bool cmp(Point a,Point b)
{
    return a.x<b.x;
}

double cal(double x)
{
    Point up(x,-inf),down(x,inf);
    int m=0,k;
    double py[3];
    for(int i=1;i<=n;++i)
    {
         k=0;
         for(int j=0;j<3;++j)
             if(SegmentProperIntersection(tri[i][j],tri[i][j+1],up,down))
                 py[++k]=GetLineIntersection(tri[i][j],tri[i][j+1]-tri[i][j],up,down-up).y;
         if(k) seg[++m]=Point(min(py[1],py[2]),max(py[1],py[2]));
    }
    if(!m) return 0;
    sort(seg+1,seg+m+1,cmp);
    double sum=seg[1].y-seg[1].x,r=seg[1].y;
    for(int i=2;i<=m;++i)
        if(dcmp(seg[i].y-r)>0)
        {
            sum+=seg[i].y-max(r,seg[i].x);
            r=seg[i].y;
        }
    return sum;
}

int main()
{
    freopen("data.txt","r",stdin);
    freopen("my.txt","w",stdout);
    int m=0;
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
    {
        for(int j=0;j<=2;++j)
        {
            scanf("%lf%lf",&tri[i][j].x,&tri[i][j].y);
            line[++m]=tri[i][j].x;
        }
        tri[i][3]=tri[i][0];
    }
    for(int i=1;i<n;++i)
        for(int j=i+1;j<=n;++j)
            for(int k=0;k<3;++k)
                for(int l=0;l<3;++l)
                    if(SegmentProperIntersection(tri[i][k],tri[i][k+1],tri[j][l],tri[j][l+1]))
                        line[++m]=GetLineIntersection(tri[i][k],tri[i][k+1]-tri[i][k],tri[j][l],tri[j][l+1]-tri[j][l]).x;
    double ans=0;
    sort(line+1,line+m+1);
    for(int i=2;i<=m;++i)
        if(dcmp(line[i]-line[i-1]))
             ans+=(line[i]-line[i-1])*cal((line[i]+line[i-1])/2);
    printf("%.2lf",ans-eps);
    return 0;
}

原文地址:https://www.cnblogs.com/TheRoadToTheGold/p/12221288.html

时间: 2024-10-07 00:59:14

bzoj 1845: [Cqoi2005] 三角形面积并的相关文章

BZOJ 1845: [Cqoi2005] 三角形面积并 [计算几何 扫描线]

1845: [Cqoi2005] 三角形面积并 Time Limit: 3 Sec  Memory Limit: 64 MBSubmit: 1151  Solved: 313[Submit][Status][Discuss] Description 给出n个三角形,求它们并的面积. Input 第一行为n(N < = 100), 即三角形的个数 以下n行,每行6个整数x1, y1, x2, y2, x3, y3,代表三角形的顶点坐标.坐标均为不超过10 ^ 6的实数,输入数据保留1位小数 Out

bzoj 1845: [Cqoi2005] 三角形面积并 扫描线

1845: [Cqoi2005] 三角形面积并 Time Limit: 3 Sec  Memory Limit: 64 MBSubmit: 848  Solved: 206[Submit][Status][Discuss] Description 给出n个三角形,求它们并的面积. Input 第一行为n(N < = 100), 即三角形的个数 以下n行,每行6个整数x1, y1, x2, y2, x3, y3,代表三角形的顶点坐标.坐标均为不超过10 ^ 6的实数,输入数据保留1位小数 Outp

BZOJ 1845 [Cqoi2005] 三角形面积并 计算几何扫描线

题意:链接 方法:扫描线 解析: 扫描线裸题. 先搞出所有边. 之后求一下所有边的交点. 按照x坐标排序. 之后每两个不同的坐标间的面积其实就是中间截出来的所有梯形以及三角形的中位线的和乘以高度. 至于中位线,再用这个与y轴平行的线截一下所有的三角形即可. 统计长度.把所有的加起来即可. 卡精! 代码: #include <cmath> #include <cstdio> #include <iomanip> #include <cstring> #incl

BZOJ1845CQOI2005三角形面积并

1845: [Cqoi2005] 三角形面积并 Description 给出n个三角形,求它们并的面积. Input 第一行为n(N < = 100), 即三角形的个数 以下n行,每行6个整数x1, y1, x2, y2, x3, y3,代表三角形的顶点坐标.坐标均为不超过10 ^ 6的实数,输入数据保留1位小数 Output 输出并的面积u, 保留两位小数 Sample Input 2 0.0 0.0 2.0 0.0 1.0 1.0 1.0 0.0 3.0 0.0 2.0 1.0 Sample

BZOJ 1845 CQOI 2005 三角形面积并 扫描线

题目大意:给出一些三角形,求这些三角形面积的并. 思路:应该可以辛普森积分,但是应该会很麻烦.. 以前扫描线就写过矩形的用数据结构维护的那种,和计算几何不占边,这次才是好好写了一次正宗的扫描线.不得不说这个算法还是很靠谱的. 其实这个思路不仅限于三角形面积的并,所有凸多边形的面积并应该都可以解决. 对于任意由线段组成的图形,对这个图形进行多次的划分,总可以将这个图形划分成梯形,面积也很好计算.那么按照什么划分呢?将所有三角形的边都求交点,不难发现,以这些点为划分依据的话,相邻两点之间一定是一个或

叉乘、快速排斥与跨立实验及求取三角形面积

<pre name="code" class="cpp">叉乘 (一)判断方向 (二)判断线段相交 (三)求三角形面积 (一)判断方向 叉乘的性质如下: (1). P x Q > 0; 表示P在Q的顺时针方向; (2). p x Q < 0; 表示P在Q的逆时针方向; (3). P x Q = 0; 表示P和Q是共线的 P(x1,y1),Q(x2,y2), P*Q=x1y2-x2y1 判断结果三种状态 模版为: struct point {

三角形面积

如图1所示.图中的所有小方格面积都是1. 那么,图中的三角形面积应该是多少呢? 请填写三角形的面积.不要填写任何多余内容或说明性文字.

三角形面积(海伦公式)

来自:http://www.oschina.net/code/snippet_149523_14180 描述 给你三个点,表示一个三角形的三个顶点,现你的任务是求出该三角形的面积 输入 每行是一组测试数据,有6个整数x1,y1,x2,y2,x3,y3分别表示三个点的横纵坐标.(坐标值都在0到10000之间)输入0 0 0 0 0 0表示输入结束测试数据不超过10000组 输出 输出这三个点所代表的三角形的面积,结果精确到小数点后1位(即使是整数也要输出一位小数位) 样例输入 0 0 1 1 1

rwkj 1363 正方形 长方形 三角形面积

C++:重载函数2(计算面积)时间限制(普通/Java):1000MS/3000MS 运行内存限制:65536KByte总提交:370 测试通过:241 描述 定义重载函数area(),分别计算正方形.长方形和三角形的面积. 输入 共计有3行. 第1行有一个实数,为正方形的边长: 第2行有二个实数,为长方形的两条边长: 第3行有三个实数,为三角形的三边长. 输出 正方形.长方形和三角形的面积(保留3位小数). 样例输入 3.54.2 5.623.0 4.0 5.0 样例输出 12.25023.6