【BZOJ1845】[Cqoi2005] 三角形面积并 几何+扫描线
【BZOJ1845】[Cqoi2005] 三角形面积并
Description
Input
Output
Sample Input
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 Output
题解:先求出所有直线交点的x坐标(注意eps!),然后从左到右扫描每个区域。发现每个区域都可以看成若干个梯形,所以只需要用一条直线去截这些梯形的中间,得到所有梯形的中位线长乘上高度即为答案。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=110;
const double eps=1e-10;
double now,ans;
int n,m,tot,sum;
struct point
{
double x,y;
point(){}
point(double a,double b){x=a,y=b;}
point operator + (point a) {return point(x+a.x,y+a.y);}
point operator - (point a) {return point(x-a.x,y-a.y);}
double operator * (point a) {return x*a.y-y*a.x;}
point operator * (double a) {return point(x*a,y*a);}
point operator / (double a) {return point(x/a,y/a);}
};
struct line
{
point x,v;
line() {}
line(point a,point b){x=a,v=b;}
}s[310];
struct trangle
{
point a,b,c;
}t[110];
struct node
{
double x;
int v;
}q[210];
double p[100000];
bool cmp(node a,node b)
{
return a.x<b.x;
}
point getpoint(line a,line b)
{
point u=a.x-b.x;
double temp=(b.v*u)/(a.v*b.v);
return a.x+a.v*temp;
}
void calc(line a,line b)
{
if(fabs(a.v*b.v)<eps) return ;
point c=getpoint(a,b);
p[++m]=c.x;
}
int main()
{
scanf("%d",&n);
int i,j;
for(i=1;i<=n;i++)
{
scanf("%lf%lf%lf%lf%lf%lf",&t[i].a.x,&t[i].a.y,&t[i].b.x,&t[i].b.y,&t[i].c.x,&t[i].c.y);
if(t[i].a.x>t[i].b.x) swap(t[i].a,t[i].b);
if(t[i].a.x>t[i].c.x) swap(t[i].a,t[i].c);
if(t[i].b.x>t[i].c.x) swap(t[i].b,t[i].c);
s[i]=line(t[i].a,t[i].b-t[i].a),s[i+n]=line(t[i].a,t[i].c-t[i].a),s[i+n+n]=line(t[i].b,t[i].c-t[i].b);
}
for(i=1;i<=3*n;i++) for(j=i+1;j<=3*n;j++) calc(s[i],s[j]);
sort(p+1,p+m+1);
for(i=1;i<m;i++)
{
if(p[i+1]-p[i]<eps) continue;
now=(p[i]+p[i+1])/2;
for(tot=0,j=1;j<=n;j++)
{
if(now>=t[j].a.x&&now<=t[j].b.x)
{
q[++tot].x=t[j].a.y+(t[j].b.y-t[j].a.y)/(t[j].b.x-t[j].a.x)*(now-t[j].a.x);
q[++tot].x=t[j].a.y+(t[j].c.y-t[j].a.y)/(t[j].c.x-t[j].a.x)*(now-t[j].a.x);
if(q[tot].x<q[tot-1].x) swap(q[tot],q[tot-1]);
q[tot-1].v=1,q[tot].v=-1;
}
if(now>=t[j].b.x&&now<=t[j].c.x)
{
q[++tot].x=t[j].c.y+(t[j].b.y-t[j].c.y)/(t[j].c.x-t[j].b.x)*(t[j].c.x-now);
q[++tot].x=t[j].c.y+(t[j].a.y-t[j].c.y)/(t[j].c.x-t[j].a.x)*(t[j].c.x-now);
if(q[tot].x<q[tot-1].x) swap(q[tot],q[tot-1]);
q[tot-1].v=1,q[tot].v=-1;
}
}
sort(q+1,q+tot+1,cmp);
for(sum=0,j=1;j<=tot;j++)
{
if(sum) ans+=(p[i+1]-p[i])*(q[j].x-q[j-1].x);
sum+=q[j].v;
}
}
printf("%.2lf",ans-eps);
return 0;
}
【BZOJ1845】[Cqoi2005] 三角形面积并 几何+扫描线的更多相关文章
- BZOJ 1845: [Cqoi2005] 三角形面积并 [计算几何 扫描线]
1845: [Cqoi2005] 三角形面积并 Time Limit: 3 Sec Memory Limit: 64 MBSubmit: 1151 Solved: 313[Submit][Stat ...
- BZOJ1845 [Cqoi2005] 三角形面积并 扫描线 计算几何
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ1845 题意概括 给出n个三角形,求其面积并. 题解 有一个很经典的扫描线题目:矩形面积并.那个比较 ...
- BZOJ1845 : [Cqoi2005] 三角形面积并
求出所有交点后从左往右扫描线,用每段的中位线去截所有三角形,算出长度并后乘以该段长度即可,时间复杂度$O(n^3\log n)$. #include<cstdio> #include< ...
- bzoj 1845: [Cqoi2005] 三角形面积并 扫描线
1845: [Cqoi2005] 三角形面积并 Time Limit: 3 Sec Memory Limit: 64 MBSubmit: 848 Solved: 206[Submit][Statu ...
- CQOI2005 三角形面积并 和 POJ1177 Picture
1845: [Cqoi2005] 三角形面积并 Time Limit: 3 Sec Memory Limit: 64 MBSubmit: 1664 Solved: 443[Submit][Stat ...
- [CQOI2005]三角形面积并
[CQOI2005]三角形面积并 题目大意: 求\(n(n\le100)\)个三角形的面积并. 思路: 自适应辛普森法,玄学卡精度可过. 源代码: #include<cmath> #inc ...
- BZOJ 1845: [Cqoi2005] 三角形面积并 (辛普森积分)
大力辛普森积分 精度什么的搞了我好久- 学到了Simpson的一个trick 深度开11,eps开1e-4.跑的比有些扫描线还快- CODE #include <bits/stdc++.h> ...
- [hdu4629 Burning]三角形面积并,扫描线
题意:给n个三角形,分别求覆盖1次~n次的总面积 思路: 对每个y坐标作一条平行于x轴的直线,按直线从下往上处理,每两条直线之间为若干梯形(也可以是三角形)首尾相连的情况,从左扫到右时,用一个变量cn ...
- ytu 1058: 三角形面积(带参的宏 练习)
1058: 三角形面积 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 190 Solved: 128[Submit][Status][Web Boar ...
随机推荐
- Linux System Programming 学习笔记(三) 标准缓冲I/O
1. partial block operations are inefficient. The operating system has to “fix up” your I/O by ensuri ...
- iis 配置 aspnet起始页
起始页配置:在default document上面加上路径就可以了,例如 /home/kickoffjob asp.net 关于起始页的一般设置: 1: 在web.config中,加入form认证: ...
- 我的logging 配置
#encoding=utf-8 import logging.config logging.config.dictConfig({ 'version': 1, 'disable_existing_lo ...
- [Inside HotSpot] Java分代堆
[Inside HotSpot] Java分代堆 1. 宇宙初始化 JVM在启动的时候会初始化各种结构,比如模板解释器,类加载器,当然也包括这篇文章的主题,Java堆.在hotspot源码结构中gc/ ...
- 思科CCIE全新升级,SDN/SD-WAN成重头戏!
CCIE,全称Cisco Certified Internetwork Expert,是美国Cisco公司于1993年开始推出的专家级认证考试.被全球公认为IT业最权威的认证,是全球Internetw ...
- 554 DT:SPM 163 smtp5,D9GowAD3RPYqSvxZjpMaAA--.4817S2 1509706293 坑爹的防垃圾邮件机制
代码如下 package ssmtest; import java.io.File;import java.io.UnsupportedEncodingException;import java.ut ...
- SQL-基础学习3--通配符:LIKE,%,(_); 拼接:+,||,concat;
第六课 用通配符进行过滤 6.1 LIKE操作符 通配符本身实际上是SQL的WHERE子句中有特殊含义的字符,SQL支持几种通配符.为在搜索子句中使用通配符,必须使用LIKE操作符.LIKE指示DB ...
- 记录一下:关于mysql数据误删除恢复的问题
大概看来几篇博客: 1.delete的可以通过回滚(rollback)的方式恢复;但是前提是,你的数据表引擎是InnoDB而不是MyISAM,而且操作不是自动提交的 但是这种方式不可以恢复trunca ...
- Redis 数据结构解析和命令指南
命令參考文档:redis commands - 你或许已经知道Redis并非简单的key-value存储.实际上他是一个数据结构server.支持不同类型的值. 也就是说.你不必只把字符串当作键所指向 ...
- C# 使用xenocode混淆加密【转】
http://www.cnblogs.com/chendaoyin/archive/2013/05/03/3056692.html 之前就了解过混淆加密工具,但这还是第一次使用,选择了xenoco ...