usaco3.4 fence4从一点看多边形，输出能看到的边

最新推荐文章于 2022-02-25 19:38:17 发布
最新推荐文章于 2022-02-25 19:38:17 发布
阅读量509 收藏
点赞数
分类专栏： usaco 文章标签： ini c
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ahero_happy/article/details/6592346
版权

        题意：从不在多边形上的一点看多边形，输出能看到的边（看到边的部分也算，当然只看到一个端点的不算）

        分析：我们可以把那点看作一个光源，然后枚举从光源发出的光线，直接枚举每条光线是不可能的。只能枚举一些特殊的光线，即经过可见端点的光线（枚举时，看作一个很小的光束）。对于每条这样的光线，它的光(因为是光束，它的左右也有光)最多可以到达另外一条边，且这条边一定是顺时针或逆时针寻找到达的第一条边。那么这条边到底可不可以被这束光照到呢？要被照到要满足的条件：光源到这条边的过程中，它的光束左右两边不能同时被挡住。

        实际操作：我们先把光源能直接看到的端点求出来，那么边只能通过下面两种途径照到（1）光源能直接看到的端点的附近(坐标+0.5)的点也能被看到（2）通过光源与直接看到的端点的光束照到。。。对于第二种情况，注意可能光线还经过别的端点，那时就要看光线的左右是否都被挡住了

/*
ID:kgdpgfy1
LANG:C++
PROG:fence4
*/
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<iostream>
using namespace std;

struct point
{
	int x,y;
}p[210];
int n,x,y,q[210],l,r;	
bool flag[210],f[210],tag,l1,r1;

int Cr(int k,int i,int j)//叉积
{
	if((p[i].x-p[k].x)*(p[j].y-p[k].y)-(p[j].x-p[k].x)*(p[i].y-p[k].y)==0) return 0;
	else if((p[i].x-p[k].x)*(p[j].y-p[k].y)-(p[j].x-p[k].x)*(p[i].y-p[k].y)<0) return -1;
	return 1;
}
bool C(int k,int i,int j)//(n,k)与(i,j)相交？
{
//	p[n].x=x,p[n].y=y;
	if(Cr(n,k,i)*Cr(n,k,j)>0) return false;
	if(Cr(i,j,n)*Cr(i,j,k)>0) return false;
	return true;
}
bool Cross(int k)//k点是否被看见
{
	int i;
	for(i=0;i<n;i++)//枚举每条线段，是不是把k这个点挡住了
		if(k!=i&&k!=(i+1)%n&&C(k,i,(i+1)%n))
			return false;
	return true;
}

bool O(double xx,double yy,int n,int i,int j)//(xx,yy)点是否被i,j挡住
{
	if(((p[i].x-p[j].x)*(yy-p[j].y)-(xx-p[j].x)*(p[i].y-p[j].y))*((p[i].x-p[j].x)*(p[n].y-p[j].y)-(p[n].x-p[j].x)*(p[i].y-p[j].y))>=0) 
		return false;
	if(((p[i].x-xx)*(p[n].y-yy)-(p[i].y-yy)*(p[n].x-xx))*((p[j].x-xx)*(p[n].y-yy)-(p[n].x-xx)*(p[j].y-yy))>=0)
		return false;
	return true;
}
bool Line(double xx,double yy,int i,int j)//(xx,yy)是否在i,j的直线上
{
	if(fabs((yy-p[i].y)*(p[j].x-p[i].x)-(p[j].y-p[i].y)*(xx-p[i].x))<0.000000001)
		return true;
	return false;
}
bool OK(double xx,double yy)//(xx,yy)是否可以被看见
{
	int i;
	for(i=0;i<n;i++)
	{
		if(O(xx,yy,n,i,(i+1)%n))
			return false;
	}
	return true;
}
bool Cal(int n,int k,int i,int j)//从点n发出的光经过点k，射到i,j上。i,j是否可被看见
{
	int a1,b1,c1,a2,b2,c2;
	double x0,y0;
	a1=p[j].y-p[i].y;
	b1=p[i].x-p[j].x;
	c1=p[j].x*p[i].y-p[i].x*p[j].y;
	a2=p[n].y-p[k].y;
	b2=p[k].x-p[n].x;
	c2=p[n].x*p[k].y-p[k].x*p[n].y;
	x0=(double)(c2*b1-c1*b2)/(a1*b2-a2*b1);//交于(x0,y0)点，往两边移一点点，看能不能被看见
	y0=(double)(a1*c2-a2*c1)/(a2*b1-a1*b2);
	
	if(p[k].x!=p[n].x)
	{
		if(!(l1&&r1))//左右有一端没被挡住，那么这条线段就是可被照到的
			f[i]=true;
		else if((fabs(p[n].x-l)<fabs(p[n].x-x0))&&(fabs(p[n].x-r)<fabs(p[n].x-x0)))
			f[i]=false;
		else f[i]=true;
	}
	else
	{
		if(!(l1&&r1))
			f[i]=true;
		else if((fabs(p[n].y-l)<fabs(p[n].y-y0))&&(fabs(p[n].y-r)<fabs(p[n].y-y0)))
			f[i]=false;
		else f[i]=true;
	}
	return f[i];
}
int main()
{
	freopen("fence4.in","r",stdin);
	freopen("fence4.out","w+",stdout);
	int i,j,k,head,tail,ii,jj;
	double xx,yy;

	
	while(scanf("%d",&n)!=EOF)
	{
		scanf("%d%d",&x,&y);
		p[n].x=x,p[n].y=y;
		for(i=0;i<n;i++)
			scanf("%d%d",&p[i].x,&p[i].y);
		memset(flag,false,sizeof(flag));
		memset(f,false,sizeof(f));
		head=tail=0;
		//扫描可以直接看到的端点
		for(i=0;i<n;i++)
		{
			flag[i]=Cross(i);//i点是否可以被看见
			if(flag[i])
				q[tail++]=i;
		}
		for(i=0;i<n-1;i++)
			if(flag[i]&&flag[i+1])//两个端点都可以被看见，那么这条线段一定能被看见
				f[i]=true;
		if(flag[0]&&flag[n-1])
			f[n-1]=true;

		while(head!=tail)//光源能直接看到的点的附近的点是否也能被看到
		{
			k=q[head++];
			if(!f[k])
			{
				if(p[k].x!=p[(k+1)%n].x)
				{
					if(p[(k+1)%n].x>p[k].x)
						xx=p[k].x+0.5;
					else xx=p[k].x-0.5;
					yy=((p[(k+1)%n].y-p[k].y)*(xx-p[k].x))/(p[(k+1)%n].x-p[k].x)+p[k].y;
				}
				else
				{
					xx=p[k].x;
					if(p[(k+1)%n].y>p[k].y)
						yy=p[k].y+0.5;
					else yy=p[k].y-0.5;
				}
				if(!Line(xx,yy,k,n)&&OK(xx,yy)) f[k]=true;
			}
			if(!f[(k-1+n)%n])
			{
				if(p[k].x!=p[(k+n-1)%n].x)
				{
					if(p[(k-1+n)%n].x>p[k].x)
						xx=p[k].x+0.5;
					else xx=p[k].x-0.5;
					yy=((p[(k-1+n)%n].y-p[k].y)*(xx-p[k].x))/(p[(k-1+n)%n].x-p[k].x)+p[k].y;
				}
				else 
				{
					xx=p[k].x;
					if(p[(k-1+n)%n].y>p[k].y)
						yy=p[k].y+0.5;
					else yy=p[k].y-0.5;
				}
				if(!Line(xx,yy,k,n)&&OK(xx,yy)) f[(k-1+n)%n]=true;
			}
		}

		//假设光源与每个能直接看到的端点发出一束光	
		head=0;
		while(head!=tail)
		{
			k=q[head++];

			l1=r1=false;//标记光的两端是不是有紧挨着的线段
			if(p[k].x!=p[n].x)//求左右挡着光线的最接近k的点
			{
				if(Cr(n,k,(k+1)%n)<0||Cr(n,k,(k-1+n)%n)<0)
				{
					r1=true;
					r=p[k].x;
				}
				if(Cr(n,k,(k+1)%n)>0||Cr(n,k,(k-1+n)%n)>0) 
				{
					l1=true;
					l=p[k].x;
				}
				if(l1&&r1) continue;//如果k的两边都有线段挡着，显然光穿不过去

				for(i=0;i<n;i++)
				{
					if(i==k) continue;

					ii=Cr(n,k,i);
					if(!ii)
					{
						if((p[n].x-p[k].x)*(p[k].x-p[i].x)>0)
						{
							if(((Cr(n,i,(i+1)%n)<0||Cr(n,i,(i-1+n)%n)<0))&&(r1==false||fabs(p[k].x-r)>fabs((double)p[k].x-p[i].x)))
							{
								r1=true;
								r=p[i].x;
							}
							if((Cr(n,i,(i+1)%n)>0||Cr(n,i,(i-1+n)%n)>0)&&(l1==false||fabs(p[k].x-l)>fabs((double)p[k].x-p[i].x)))
							{
								l1=true;
								l=p[i].x;
							}
						}
					}
				}
			}
			else
			{
				if(Cr(n,k,(k+1)%n)<0||Cr(n,k,(k-1+n)%n)<0)
				{
					r1=true;
					r=p[k].y;
				}
				if(Cr(n,k,(k+1)%n)>0||Cr(n,k,(k-1+n)%n)>0) 
				{
					l1=true;
					l=p[k].y;
				}
				if(l1&&r1) continue;

				for(i=0;i<n;i++)
				{
					if(i==k) continue;

					ii=Cr(n,k,i);
					if(!ii)
					{
						if((p[n].y-p[k].y)*(p[k].y-p[i].y)>0)
						{
							if(((Cr(n,i,(i+1)%n)<0||Cr(n,i,(i-1+n)%n)<0))&&(r1==false||fabs(p[k].y-r)>fabs((double)p[k].y-p[i].y)))
							{
								r1=true;
								r=p[j].y;
							}
							if((Cr(n,i,(i+1)%n)>0||Cr(n,i,(i-1+n)%n)>0)&&(l1==false||fabs(p[k].y-l)>fabs((double)p[k].y-p[i].y)))
							{
								l1=true;
								l=p[i].y;
							}
						}
					}
				}
			}

			i=(k+1)%n;
			while(i!=k)//逆时针找第一个相交的
			{
				ii=Cr(n,k,i);
				jj=Cr(n,k,(i+1)%n);
				if(ii*jj<0)
					break;
				i=(i+1)%n;
			}
			if(i!=k&&!f[i])
			{
				if(Cal(n,k,i,(i+1)%n)) continue;
			}

			j=(k-1+n)%n;
			while(j!=k)//顺时针找第一个相交的
			{
				ii=Cr(n,k,j);
				jj=Cr(n,k,(j-1+n)%n);
				if(ii*jj<0)
					break;
				j=(j-1+n)%n;
			}
			if(j!=k&&!f[(j-1+n)%n])
			{
				Cal(n,k,(j-1+n)%n,j);
			}
		}
		
		
		//输出
		for(j=0,i=0;i<n;i++)
			if(f[i]) j++;
		if(!j) printf("NOFENCE\n");
		else
		{
			printf("%d\n",j);
			for(i=0;i<n-2;i++)
				if(f[i])
					printf("%d %d %d %d\n",p[i].x,p[i].y,p[i+1].x,p[i+1].y);
			if(f[n-1])
				printf("%d %d %d %d\n",p[0].x,p[0].y,p[n-1].x,p[n-1].y);
			if(f[n-2])
				printf("%d %d %d %d\n",p[n-2].x,p[n-2].y,p[n-1].x,p[n-1].y);
		}
	}
	return 0;
}


 

 

ahero_happy
关注
专栏目录
Binary Tree Right Side View 从左边视角输出二叉树的值
小人物的草稿本
06-01 1545
Binary Tree Right Side View   Given a binary tree, imagine yourself standing on the right side of it, return the values of the nodes you can see ordered from top to bottom. For example: Give
USACO 3.4 Closed Fences (fence4)
Hello World！
03-01 1721
/* Main idea: 1.判断多边形是否合法 任两条边都不相交即合法，注意这里的相交是严格相交，顶点相交不算相交。 2.二分法判断当前线段 seg_a 是否可见 假设观察点为 eye，seg_a 的两个端点分别为st 和 end。判断视线（eye，st）和（eye，end）是否 与其他线段（即 fence）相交。 如果都不相交，则seg_a 可见。 如果两视线均与某一 fence 相交，则s
USACO 3.4
lemondinosaur的博客
11-06 385
解题报告洛谷 1827 美国血统题目代码洛谷 2735 电网题目分析代码洛谷 2736 破锣摇滚乐队代码 洛谷 1827 美国血统 题目 给出一棵树的前序遍历和中序遍历，求它的后序遍历 代码 /* ID:lemondi1 LANG:C++ TASK:heritage */ #include <cstdio> #include <cstring> #define rr re...
USACO 3.4 Electric Fence
weixin_30636089的博客
01-02 74
Electric FenceDon Piele In this problem, `lattice points' in the plane are points with integer coordinates. In order to contain his cows, Farmer John constructs a triangular electric fence by s...
usaco 3.4 Electric Fence 2010.8.4
Batman_的博客
02-13 313
Electric Fence Don Piele In this problem, `lattice points' in theplane are points with integer coordinates.   In order to contain his cows, Farmer Johnconstructs a triangular electric fence by
USACO 3.4 Electric Fence 皮克定理
weixin_30784945的博客
12-11 104
题意：在方格纸上画出一个三角形，求三角形里面包含的格点的数目 因为其中一条边就是X轴，一开始想的是算出两条边对应的数学函数，然后枚举x坐标值求解。但其实不用那么麻烦。 皮克定理：给定顶点坐标均是整点（或正方形格点）的简单多边形，皮克定理说明了其面积A和内部格点数目i、边上格点数目b的关系：A = i + b/2 - 1。 有了这条定理就好办了。 三角形面积直接用公式就能算出来。...
Electric Fences_usaco3.4_皮克定理
olahiuj的博客
10-16 598
Description给定两个整点(m,n)(m,n)和(p,0)(p,0)，与坐标原点(0,0)(0,0)构成三角形，求三角形内的整点数量Analysis找到这么一道神奇的题目以及神奇的定理 皮克定理 一个计算点阵中顶点在格点上的多边形面积公式：S=a+b2−1S=a+\frac{b}{2}-1，其中a表示多边形内部的点数，b表示多边形边界上的点数，s表示多边形的面积
usaco 3.4 Closed Fences 计算几何
Ben的专栏
08-20 1168
题目给定一个多边形,要求1)判断是否简单多边形方法如下 对于临边：假设p1,p2和p2,p3为临边,如果p1,p2,p3重合并且p1,p3在p2的同侧,在存在交叉边,可以通过叉积cross(p1,p2,p3)==0 和点击dot(p1,p2,p3)  对于不相邻的边：如果规范相交，则存在交叉边. 2)判断观察者是否能看到，题目的意思是只要能看到一点就满足条件，但是如
USACO SECTION 3.4 Electric Fence
08-27 249
Electric Fence Don Piele In this problem, `lattice points' in the plane are points with integer coordinates. In order to contain his cows, Farmer John constructs a triangular electric fence by s
usaco3.4解题报告1
08-03
在编程竞赛USACO中，这个问题被称为"Fence4 Closed Fences"，需要检查给定的顶点列表是否形成一个合法的闭合栅栏，并找出从特定点(x, y)可以看到的所有栅栏。 【描述】中提到，闭合栅栏是由不相交的线段首尾相连...
P1827 [USACO3.4]美国血统 American Heritage (数据结构之树的前序中序和后序遍历
yxrzibinanhai的博客
07-11 382
C 　　　　　　 / \ 　　　　　　 /　　\ 　　　　　　 B　 G 　　　　　　/ \　　 / 　　　　 A D H 　　　　　　 / \ 　　　　　　 E F
（dfs）[USACO3.4]“破锣摇滚”乐队 Raucous Rockers
numb_ac的博客
07-03 662
一、算法分析 自己初步分析是一个dfs问题，（据说还可以DP？但是我加上一点点剪枝和回溯之后比好多DP都快，下次再写个DP的）。容易想到如果当前歌曲不选，那么以后永远不会再选了，且对于任何一个新的歌曲，只有三种选择：不选、放到最后一张CD里（此时要求最后一张CD容量足够）或者新开一张CD。但是本题一些细节还是要注意。一个是可以最优性剪枝（详见代码），还有就是在新开一张的时候要处理一些细节（详见代码）。 可能是半夜打代码的原因，写完之后几个点死活不对，然后发现自己把n和m搞混了两次，还有一次把n写成n-1 二
P1827 [USACO3.4]美国血统 American Heritage 【二叉树】
ln2037的博客
10-02 898
题目描述 农夫约翰非常认真地对待他的奶牛们的血统。然而他不是一个真正优秀的记帐员。他把他的奶牛 们的家谱作成二叉树，并且把二叉树以更线性的“树的中序遍历”和“树的前序遍历”的符号加以记录而 不是用图形的方法。 你的任务是在被给予奶牛家谱的“树中序遍历”和“树前序遍历”的符号后，创建奶牛家谱的“树的 后序遍历”的符号。每一头奶牛的姓名被译为一个唯一的字母。（你可能已经知道你可以在知道树的两 种遍历以后可以经常地重建这棵树。）显然，这里的树不会有多于 26 个的顶点。 这是在样例输入和 样例输出中的树的图形表达
K8s 很难么？带你从头到尾捋一遍，不信你学不会
民工哥的博客
02-23 9092
点击关注公众号，回复“1024”获取2TB学习资源！为什么要学习 Kubernetes?虽然 Docker 已经很强大了，但是在实际使用上还是有诸多不便，比如集群管理、资源调度、文件管理等...
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）
北豼不太皮的博客
02-21 3625
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）一、学习目标二、了解单片机STM32F401RET6三、C语言基础 一、学习目标 二、了解单片机STM32F401RET6 4、STM32F401RE特征 三、C语言基础 1.数据类型 常用2的次方： 2^7 = 128 2^8 = 256 2^15 = 32768 2^16= 65536 51单片机常见的数据类型： char: 占内存1字节 取值范围：-128~127 -2^7 ~ 2
Android开发（1） | Fragment 的应用——新闻应用
Jormungand_V的博客
02-21 7119
文章目录 news_item.xml： <TextView xmlns:android="http://schemas.android.com/apk/res/android" android:id="@+id/news_title" android:layout_width="match_parent" android:layout_height="wrap_content" android:lines="1" android:ellipsize="
系统学习 TypeScript（四）——变量声明的初步学习
编程三昧
02-25 1758
认识了 TypeScript 中的基础类型，接下来当然是变量声明的相关学习了，我在这里记录一下我学习过程中的一些总结。
清除浮动的五种方法
weixin_52212950的博客
02-22 2732
为什么要清除浮动？因为往往浮动后的子元素因为脱离了标准流，不能自动撑起父元素的高度，所以会对后续布局产生影响，对此清除浮动不如理解为清除浮动产生的影响更为合理。 例如：我们设置了两个盒子如图所示，粉色为父盒子，只有宽度没有高度，蓝色盒子有宽有高，粉色盒子的高由蓝盒子的高度撑开。 但是给蓝色的子盒子设置了左浮动后，脱离了标准流，无法再撑开粉盒子，可以看到粉盒子高度变成了0； 清除浮动有五种方法： 直接设置父元素的高度 额外标签法 单伪元素法 双伪元素法 ove.
USACO Chapter 4 题解：动态规划与多重背包问题
"usaco chap4 题解" 在USACO（美国计算机奥林匹克）的章节4中，我们主要探讨了两个题目及其解决方案。首先是"BeefMcNuggets(nuggets)"，这是一个典型的背包问题，它涉及到动态规划的运用。 在"BeefMcNuggets"题目...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值