`

hdu3867

阅读更多

/*
点光源能够射到的棍子的个数,线段树。
线段树的每个节点记录覆盖这段区间的棍子的下标,以及这根 棍子上的 某点到点光源的距离 
*/
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int N = 10005;
const double eps = 1e-7;
const double pi = acos(-1.0);
int sign(double a){
	return a < -eps ? -1 : (a > eps);
}
struct point{
	double x, y;
	void read(){
		scanf("%lf%lf", &x, &y);
	};
};
struct seg{
	int l, r;
	int id;
	double len;
	int mid(){
		return (l + r) >> 1;
	}
	int getSum(){
		return r - l + 1;
	}
}segs[N*16];
int n;
double ks[N*4];
int  kn;
point ss[N][2], s;
double angs[N][2];
bool mark[N];
double getmid(double r, double l){
	if(r - l > pi) return (r + l) * 0.5 + pi;
	return (r + l) * 0.5;
}
bool input(){
	if(scanf("%d", &n) == EOF) return false;
	point a, b;
	s.read();
	int i;
	for(kn = i = 0; i < n; i++){
		a.read(); ss[i][0] = a;
		b.read(); ss[i][1] = b;
		angs[i][0] = ks[kn++] = atan2(a.y-s.y, a.x-s.x);
		angs[i][1] = ks[kn++] = atan2(b.y-s.y, b.x-s.x);
		if(angs[i][0] > angs[i][1]){
			swap(angs[i][0], angs[i][1]);
			swap(ss[i][0], ss[i][1]);
		}

	}
	return true;
}
int hash(double k){
	return lower_bound(ks, ks+kn, k) - ks;
}
inline int ls(int id){
	return (id<<1)+1;
}
inline int rs(int id){
	return (id<<1)+2;
}
void build(int id, int l, int r){
	segs[id].l = l;
	segs[id].r = r;
	segs[id].id = -1;
	if(l < r){
		int mid = (l + r) >> 1;
		build(ls(id), l, mid);
		build(rs(id), mid+1, r);
	}
}
void change(int id, int k){
	double len;
	point a = ss[k][0], b = ss[k][1];
	double ang = getmid(ks[segs[id].r], ks[segs[id].l]);
	len = (b.y-a.y)*(s.x-a.x) - (b.x-a.x)*(s.y-a.y);
	len /= (sin(ang)*(b.x-a.x) - cos(ang)*(b.y-a.y));
	if(segs[id].id < 0 || sign(len - segs[id].len) < 0){
		segs[id].id = k;
		segs[id].len = len;
	}
}
void down(int id){
	if(segs[id].id >= 0){
		change(ls(id), segs[id].id);
		change(rs(id), segs[id].id);
	}
}
void ins(int id, int l, int r, int k){
	if(l <= segs[id].l && segs[id].r <= r){
		change(id, k);
		return;
	}
	if(segs[id].l == segs[id].r) return ;
	down(id);
	int mid = segs[id].mid();
	if(l <= mid) ins(ls(id), l, r, k);
	if(r > mid) ins(rs(id), l, r, k);
}
void qry(int id){
	if(segs[id].l == segs[id].r){
		if(segs[id].id >= 0){
			mark[segs[id].id] = true;
		}
		return;
	}
	down(id);
	qry(ls(id));
	qry(rs(id));
}
void solve(){
	sort(ks, ks+kn);
	kn = unique(ks, ks+kn) - ks;
	int i, ans, l, r, j;
	j = kn;
	for(i = 1; i < j; i++){
		ks[kn++] = getmid(ks[i], ks[i-1]);
	}
	ks[kn++] = getmid(ks[j-1], ks[0]);
	sort(ks, ks+kn);
	kn = unique(ks, ks+kn) - ks;
	build(0, 0, kn - 1);
	for(i = 0; i < n; i++){
		l = hash(angs[i][0]);
		r = hash(angs[i][1]);
		mark[i] = false;
		if(angs[i][1] - angs[i][0] > pi){
			ins(0, 0, l, i);
			ins(0, r, kn-1, i);
		}else{
			ins(0, l, r, i);
		}
	}
	qry(0);
	for(ans = i = 0; i < n; i++){
		if(mark[i]) ans++;
	}
	printf("%d\n", ans);
}
int main(){
	//freopen("in.txt", "r", stdin);
	while(input()) solve();
	return 0;
}
 
分享到:
评论

相关推荐

    hdu.rar_hdu

    HDU(杭州电子科技大学在线评测系统)是一个深受程序员喜爱的在线编程练习平台,它提供了丰富的算法题目供用户挑战,帮助他们提升编程技能和算法理解能力。"hdu.rar_hdu"这个压缩包文件很可能是某位程序员整理的他在...

    HDU_2010.rar_hdu 2010_hdu 20_hdu acm20

    【标题】"HDU_2010.rar"是一个压缩包文件,其中包含了与"HDU 2010"相关的资源,特别是针对"HDU ACM20"比赛的编程题目。"hdu 2010"和"hdu 20"可能是该比赛的不同简称或分类,而"hdu acm20"可能指的是该赛事的第20届...

    HDU题目java实现

    【标题】"HDU题目java实现"所涉及的知识点主要集中在使用Java编程语言解决杭州电子科技大学(HDU)在线评测系统中的算法问题。HDU是一个知名的在线编程竞赛平台,它提供了大量的算法题目供参赛者练习和提交解决方案...

    ACM HDU题目分类

    ACM HDU 题目分类 ACM HDU 题目分类是指对 HDU 在线判题系统中题目的分类,总结了大约十来个分类。这些分类将有助于编程选手更好地理解和解决问题。 DP 问题 DP(Dynamic Programming,动态规划)是一种非常重要...

    hdu1250高精度加法

    ### hdu1250高精度加法 #### 背景介绍 在计算机科学与编程竞赛中,处理大整数运算(特别是加法、减法、乘法等)是常见的需求之一。当数字的位数超过了标准数据类型(如`int`、`long`等)所能表示的最大值时,就需要...

    hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj

    【标题】"hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj" 提供的信息是关于一个压缩文件,其中包含了一个名为 "HDU 1089.cpp" 的源代码文件,这个文件是为了解决杭州电子科技大学(Hangzhou Dianzi ...

    HDU DP动态规划

    【标题】"HDU DP动态规划"涉及到的是在算法领域中的动态规划(Dynamic Programming,简称DP)技术,这是解决复杂问题的一种高效方法,尤其适用于有重叠子问题和最优子结构的问题。动态规划通常用于优化多阶段决策...

    HDU1059的代码

    HDU1059的代码

    hdu1001解题报告

    hdu1001解题报告

    hdu 1574 passed sorce

    hdu 1574 passed sorce

    ACM HDU

    【ACM HDU】指的是在ACM(国际大学生程序设计竞赛,International Collegiate Programming Contest)中,参赛者在杭州电子科技大学(Hangzhou Dianzi University,简称HDU)的在线评测系统上完成并已解决的题目集合...

    hdu2101解决方案

    hdu2101AC代码

    杭电ACMhdu1163

    【标题】:杭电ACMhdu1163 【描述】:这是一道源自杭州电子科技大学(Hangzhou Dianzi University,简称HDU)的ACM编程竞赛题目,编号为1163。这类问题通常需要参赛者利用计算机编程解决数学、逻辑或算法上的挑战,...

    HDU.rar_hdu_hdu07_com_shownv9b_www.563hdu.

    【标题】"HDU.rar_hdu_hdu07_com_shownv9b_www.563hdu." 暗示这是一个与HDU(杭州电子科技大学在线编程平台)相关的压缩包,其中可能包含了该平台上的编程竞赛题目或练习题目的源代码。"hdu07"可能是某个特定题目的...

    Hdu1000—2169部分代码

    HDU是杭州电子科技大学(Hangzhou Dianzi University)举办的一个在线编程竞赛平台,全称为HDU Online Judge。ACM是国际大学生程序设计竞赛(International Collegiate Programming Contest)的缩写,是一个全球性的...

    hdu 5007 Post Robot

    hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。

    HDU acm-PPT课件

    【ACM入门与提高:HDU ACM竞赛课程详解】 ACM(国际大学生程序设计竞赛,International Collegiate Programming Contest,简称ICPC或ACM/ICPC)是一项全球性的竞赛,旨在激发大学生对计算机科学的兴趣,提升他们的...

    hdu_acm_1084.rar_ACM_HDU10_acm10_hdu_hdu 1084

    【标题】"hdu_acm_1084.rar_ACM_HDU10_acm10_hdu_hdu 1084" 提供的是一个关于杭电(HDU)ACM竞赛第1084题的解决方案。该题目可能是在编程竞赛中常见的算法问题,而ACM(国际大学生程序设计竞赛)是全球知名的编程...

    HDU最全ac代码

    HDU(Hangzhou Dianzi University)是国内外知名的在线编程竞赛平台,主要服务于ACM/ICPC(国际大学生程序设计竞赛)以及相关的算法训练。"HDU最全ac代码"这个压缩包很可能是包含了在HDU平台上解题通过的完整源代码...

    hdu动态规划算法集锦

    根据提供的信息,我们可以总结出以下关于“hdu动态规划算法集锦”的知识点: ### 动态规划基础概念 动态规划是一种解决多阶段决策问题的方法,它通过将原问题分解为互相重叠的子问题,利用子问题的解来构建原问题...

Global site tag (gtag.js) - Google Analytics