意甲冠军:

给一些直。这架飞机被分成了很多这些线性块。每个块的需求面积封闭曲线图。

分析:

①我们应要求交点22的直线;

②每行上的交点的重排序,借此来离散一整行(正反两条边);

③对于连向一个点的几条线段,对它们进行极角排序,相邻的两条线段我们给它们之间连一条边,我们脑补一下应该能够知道如何能够保证逆时针连边;

④找循环,利用叉积求面积。

ps. vector的调试真心不爽…

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <vector>
#define pb push_back
#include <utility>
#define fi first
#define se second
#define mp make_pair
#include <stack>
using namespace std;
const int MAXN = 100;
const int MAXE = MAXN*MAXN<<1;
const double eps = 1e-8;
typedef double DB; int n;
#define sqr(x) ((x)*(x))
struct pot
{
double x, y;
pot(double x = 0, double y = 0):x(x), y(y){}
inline double len() {return sqrt(sqr(x)+sqr(y));}
inline void read() {scanf("%lf%lf", &x, &y);}
inline pot unit()
{
double d = len();
return pot(x/d, y/d);
}
inline double ang() {return atan2(y, x);}
};
typedef pot vec;
struct Line
{
pot p;
vec v;
Line(){}
Line(pot p, vec v):p(p), v(v){}
}lines[MAXN];
vector<pot> points; struct Edge
{
int sz;
int head[MAXE], to[MAXE], ne[MAXE];
Edge()
{
sz = 0;
memset(head, -1, sizeof(head));
}
inline void add(int u, int v)
{
to[sz] = v;ne[sz] = head[u];
head[u] = sz++;
}
}E; int next[MAXE]; bool vis[MAXE]; vector<DB> ans; inline int dcmp(double x)
{
if(fabs(x) <= eps) return 0;
else return x>0?1:-1;
} inline pot operator + (const vec &a, const vec &b) {return vec(a.x+b.x, a.y+b.y);}
inline pot operator - (const vec &a, const vec &b) {return vec(a.x-b.x, a.y-b.y);}
inline double operator * (const vec &a, const vec &b) {return a.x*b.x+a.y*b.y;}
inline double operator ^ (const vec &a, const vec &b) {return a.x*b.y-a.y*b.x;}
inline pot operator * (const vec &a, double k) {return vec(a.x*k, a.y*k);}
inline pot operator / (const vec &a, double k) {return vec(a.x/k, a.y/k);}
inline bool operator < (const vec &a, const vec &b) {return dcmp(a.x-b.x) < 0 || (dcmp(a.x-b.x) == 0 && dcmp(a.y-b.y) < 0);}
inline bool operator == (const vec &a, const vec &b) {return dcmp(a.x-b.x) == 0 && dcmp(a.y-b.y) == 0;} inline pot get_line_intersection(const pot &p, const vec &v, const pot &q, const vec &w)
{
vec u = p-q;
double t = (w^u)/(v^w);
return p+v*t;
} inline bool parallel(const Line &a, const Line &b) {return dcmp(a.v^b.v) == 0;} inline int get_pot_id(const pot &p) {return lower_bound(points.begin(), points.end(), p)-points.begin();} inline bool cmp(DB a, DB b) {return dcmp(a-b) == 0;} inline void init()
{
scanf("%d", &n);
for(int i = 0; i < n; ++i)
{
pot a, b;
a.read();b.read();
lines[i] = Line(a, b-a);
}
for(int i = 0; i < n; ++i)
for(int j = i+1; j < n; ++j)
if(!parallel(lines[i], lines[j]))
points.pb(get_line_intersection(lines[i].p, lines[i].v, lines[j].p, lines[j].v));
sort(points.begin(), points.end());
points.erase(unique(points.begin(), points.end()), points.end());
for(int i = 0; i < n; ++i)
{
vector<DB> nodes;
vec d = lines[i].v.unit();
for(int j = 0; j < n; ++j)
if(!parallel(lines[i], lines[j]))
{
pot inter = get_line_intersection(lines[i].p, lines[i].v, lines[j].p, lines[j].v);
nodes.pb((inter-lines[i].p)*d);
}
sort(nodes.begin(), nodes.end());
nodes.erase(unique(nodes.begin(), nodes.end(), cmp), nodes.end());
for(int j = 1, sz = nodes.size(); j < sz; ++j)
{
int a = get_pot_id(lines[i].p+d*nodes[j]);
int b = get_pot_id(lines[i].p+d*nodes[j-1]);
E.add(a, b);E.add(b, a);
}
}
memset(next, -1, sizeof(next));
for(int i = 0, sz = points.size(); i < sz; ++i)
{
vector<pair<DB, int> > PA;
for(int j = E.head[i]; j != -1; j = E.ne[j])
PA.pb(mp((points[E.to[j]]-points[i]).ang(), j));
sort(PA.begin(), PA.end());
for(int j = 0, sz = PA.size(); j < sz; ++j)
next[PA[(j+1)%sz].se^1] = PA[j].se;
}
} inline void work()
{
for(int i = 0; i < E.sz; ++i)
if(!vis[i])
{
stack<int> s;
int j = i;
do
{
if(!s.empty() && (s.top()^1) == j)
s.pop();
else s.push(j);
vis[j] = true;
j = next[j];
if(j == -1) break;
}while(!vis[j]);
if(i == j)
{
DB area = 0;
while(!s.empty())
{
area += (points[E.to[s.top()^1]]^points[E.to[s.top()]]);
s.pop();
}
area *= 0.5;
if(dcmp(area) > 0)
ans.pb(area);
}
}
} inline void print()
{
printf("%d\n", ans.size());
sort(ans.begin(), ans.end());
for(int i = 0, sz = ans.size(); i < sz; ++i)
printf("%.4lf\n", ans[i]);
} int main()
{
init();
work();
print();
return 0;
}

版权声明:本文博主原创文章,博客,未经同意不得转载。

sgu209:Areas(计算几何)的更多相关文章

  1. ASP.NET Core 中文文档 第四章 MVC(4.6)Areas(区域)

    原文:Areas 作者:Dhananjay Kumar 和 Rick Anderson 翻译:耿晓亮(Blue) 校对:许登洋(Seay) Areas 是 ASP.NET MVC 用来将相关功能组织成 ...

  2. 【无私分享:ASP.NET CORE 项目实战(第九章)】创建区域Areas,添加TagHelper

    目录索引 [无私分享:ASP.NET CORE 项目实战]目录索引 简介 在Asp.net Core VS2015中,我们发现还有很多不太简便的地方,比如右击添加视图,转到试图页等功能图不见了,虽然我 ...

  3. ASP.NET MVC Routing、Areas、URLs

    webForm页面运行起来url一般是这样的:localhost:****/index.aspx,这个过程就是当你运行页面的时候,vs开发工具自带的微型服务器会打开你存在硬盘上的这个文件然后显示在浏览 ...

  4. ACM/ICPC 之 计算几何入门-叉积-to left test(POJ2318-POJ2398)

    POJ2318 本题需要运用to left test不断判断点处于哪个分区,并统计分区的点个数(保证点不在边界和界外),用来做叉积入门题很合适 //计算几何-叉积入门题 //Time:157Ms Me ...

  5. HDU 2202 计算几何

    最大三角形 Time Limit: 5000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submi ...

  6. [转]【无私分享:ASP.NET CORE 项目实战(第九章)】创建区域Areas,添加TagHelper

    本文转自:http://www.cnblogs.com/zhangxiaolei521/p/5808417.html 目录索引 [无私分享:ASP.NET CORE 项目实战]目录索引 简介 在Asp ...

  7. ACM 计算几何中的精度问题(转)

    http://www.cnblogs.com/acsmile/archive/2011/05/09/2040918.html 计算几何头疼的地方一般在于代码量大和精度问题,代码量问题只要平时注意积累模 ...

  8. Areas on the Cross-Section Diagram

    Areas on the Cross-Section Diagram  Aizu - ALDS1_3_D Areas on the Cross-Section Diagram 地域の治水対策として.洪 ...

  9. Asp.Net Mvc Areas 的用法与好处

    前言 在项目中为什么要使用Areas 进行分离 大家都知道,一般的Web应用都有前台(面向用户)和后台(面向管理员)两部分,我们希望以/localhost/Admin 开始的Url 是用户的后台管理地 ...

随机推荐

  1. Apache的.htaccess项目根文件夹伪静态设置规则

    watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQv/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/ ...

  2. struts2笔记---struts2的执行过程

    1.服务器启动: 加载项目web.xml 创建struts核心过滤器对象,执行filter-->init() struts-default.xml     核心功能的初始化 struts-plu ...

  3. ArcGIS中ObjectID,FID和OID字段区别

    lysc_forever 原文 ArcGIS中ObjectID,FID和OID字段有什么区别 ArcGIS Desktop 独立的表和属性表都有一个ObjectID字段.这个字段中包含一个唯一的,长整 ...

  4. iOS ASIHTTPRequest

    ASIHTTPRequest对CFNetwork API进行了封装,并且使用起来非常简单,用Objective-C编写,可以很好的应用在Mac OS X系统和iOS平台的应用程序中.ASIHTTPRe ...

  5. php中的转义字符(用反斜杠\来输出,和C语言一样)

    php中的转义字符(用反斜杠\来输出,和C语言一样) 一.总结 1.引号中的变量:双引号会替换变量的值,而单引号会把它当做字符串输出. 2.引号中的转义字符:双引号将用变量的值(test)代替它的名称 ...

  6. 25、写一个USB摄像头驱动程序(有ioctrl分析)

    videobuf2-core.h中的vb2_buffer,记录了v4l2_buffer ,驱动可以对vb2_buffer的v4l2_buffer进行操控, vb2_buffer是v4l2框架层的代码, ...

  7. PJSIP开源库详解

    PJSIP是一个包含了SIP.SDP.RTP.RTCP.STUN.ICE等协议实现的开源库.它把基于信令协议SIP的多媒体框架和NAT穿透功能整合成高层次.抽象的多媒体通信API,这套API能够很容易 ...

  8. POJ 3252 Round Numbers(数位dp&amp;记忆化搜索)

    题目链接:[kuangbin带你飞]专题十五 数位DP E - Round Numbers 题意 给定区间.求转化为二进制后当中0比1多或相等的数字的个数. 思路 将数字转化为二进制进行数位dp,由于 ...

  9. C# 关于反射事件

    在frmMain类中的代码 private void StartRun(string tag, string date, bool tipType)        {            var d ...

  10. centos / Linux 服务环境下安装 Redis 5.0.3

    原文:centos / Linux 服务环境下安装 Redis 5.0.3 1.首先进入你要安装的目录 cd /usr/local 2.下载目前最新稳定版本 Redis 5.0.3 wget http ...