poj 1271 || UVA 10117 Nice Milk

最新推荐文章于 2019-09-30 13:41:26 发布
最新推荐文章于 2019-09-30 13:41:26 发布
阅读量1.5k 收藏
点赞数
分类专栏： poj UVa 计算几何 文章标签： parallel hp struct c 优化 测试
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zxy_snow/article/details/6626849
版权
poj 同时被 3 个专栏收录
211 篇文章 0 订阅
订阅专栏
计算几何
134 篇文章 1 订阅
订阅专栏
UVa
41 篇文章 0 订阅
订阅专栏

纠结了一下午啊！！！超时啊！！！我用WS的方法测了好久又测了下我的时间复杂度！！！感觉应该就卡在数据的最后一组啊！！！！

在UVA wa了啊！！！说明UVA 机器比POJ快啊！！！

思路：求蘸牛奶的最大面积，也就是半平面交后中间区域最小，所以只需要求得半平面交中间的最小区域即可。这个实现需要DFS枚举哈。

一直TLE，后来看网上的还是不知道从哪入手。后来想起一种优化方法，因为里面的边是平行的，所以如果选择里面的边，那么外面的那个边就不需要了（因为会在排序后删掉的），那么直接覆盖那条边就行了，还省点时间。这么一弄，居然WA了 = = 后来发现个很2的错误，输入放在输出前了 T T 。。。输出0.00后直接continue了 = = 改掉后AC。

两个代码都是DFS的差异，一个是查找到一条边就加入然后到k条后求面积。另一种是把DFS方案全部列出来，然后再纠结 = =。。。

第一种快点，因为如果找到面积为0的我就直接return了。而且，如果剩下可选的元素比需要选的元素个数少，也return了。不过最少还600+MS.。。。T T.。。

P.S. 我怀疑过unique的删除，因为如果用半平面交的话，第一次删除斜率相等的点的话，应该留下最左边的（我定义最左边为有效方向）。那么我怕unique乱删。不过刚才测试了下，unique删的都是第一个元素，那么我就放心了。

P.S. 刚查了下unique的内部实现，确实是保留第一个元素，呵呵。

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#include <queue>
#include <stack>
#include <climits>

using namespace std;

const int MAX = 25;
const double inf = 1e20;
struct point {double x,y;};
struct line { double ang; point a,b;};
point p[MAX],s[MAX];
line ln[MAX],deq[MAX],ltmp[MAX],ll[MAX],lk[MAX];;
double h,mmin;
const double eps = 1e-6;
int k;
bool dy(double x,double y)	{	return x > y + eps;}	// x > y 
bool xy(double x,double y)	{	return x < y - eps;}	// x < y 
bool dyd(double x,double y)	{ 	return x > y - eps;}	// x >= y 
bool xyd(double x,double y)	{	return x < y + eps;} 	// x <= y 
bool dd(double x,double y) 	{	return fabs( x - y ) < eps;}  // x == y
double disp2p(point a,point b) //  a b 两点之间的距离 
{
	return sqrt( ( a.x - b.x ) * ( a.x - b.x ) + ( a.y - b.y ) * ( a.y - b.y ) );
}
double crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 
{
	return (c.x - a.x)*(b.y - a.y) - (b.x - a.x)*(c.y - a.y);
}
void changepoint(point a,point b,point &c,point &d)
{
	double len = disp2p(a,b);
    double dx = h / len * ( a.y - b.y );
	double dy = h / len * (-a.x + b.x );
	c.x = a.x + dx; c.y = a.y + dy;
	d.x = b.x + dx; d.y = b.y + dy;
}
void makeline_hp(point a,point b,line &l) // 线段(向量ab)右侧区域有效 
{
	l.a = a; l.b = b;
	l.ang = atan2(a.y - b.y,a.x - b.x);	// 如果是左侧区域，改成b.y - a.y,b.x - a.x 
}
bool equal_ang(line a,line b)	// 第一次unique的比较函数 
{
	return dd(a.ang,b.ang);
}
bool cmphp(line a,line b)	// 排序的比较函数 
{
	if( dd(a.ang,b.ang) ) return xy(crossProduct(b.a,b.b,a.a),0.0);
	return xy(a.ang,b.ang);
}
bool equal_p(point a,point b)//第二次unique的比较函数 
{
	return dd(a.x,b.x) && dd(a.y,b.y);
}
bool parallel(line u,line v)
{
	return dd( (u.a.x - u.b.x)*(v.a.y - v.b.y) - (v.a.x - v.b.x)*(u.a.y - u.b.y) , 0.0 );
}
point l2l_inst_p(line l1,line l2)
{
	point ans = l1.a;
	double t = ((l1.a.x - l2.a.x)*(l2.a.y - l2.b.y) - (l1.a.y - l2.a.y)*(l2.a.x - l2.b.x))/
			   ((l1.a.x - l1.b.x)*(l2.a.y - l2.b.y) - (l1.a.y - l1.b.y)*(l2.a.x - l2.b.x));
	ans.x += (l1.b.x - l1.a.x)*t;
	ans.y += (l1.b.y - l1.a.y)*t;
	return ans;
}
double area_polygon(point p[],int n)
{
	double s = 0.0;
	for(int i=0; i<n; i++)
		s += p[(i+1)%n].y * p[i].x - p[(i+1)%n].x * p[i].y;
	return fabs(s)/2.0;
}
void inst_hp_nlogn(line *ln,int n,point *s,int &len)
{
	len = 0;
	sort(ln,ln+n,cmphp);
	n = unique(ln,ln+n,equal_ang) - ln;
	int bot = 0,top = 1;
	deq[0] = ln[0]; deq[1] = ln[1];
	for(int i=2; i<n; i++)
	{
		if( parallel(deq[top],deq[top-1]) || parallel(deq[bot],deq[bot+1]) )
			return ;
		while( bot < top && dy(crossProduct(ln[i].a,ln[i].b,
			l2l_inst_p(deq[top],deq[top-1])),0.0) )
			top--;
		while( bot < top && dy(crossProduct(ln[i].a,ln[i].b,
			l2l_inst_p(deq[bot],deq[bot+1])),0.0) )
			bot++;
		deq[++top] = ln[i];
	}
	while( bot < top && dy(crossProduct(deq[bot].a,deq[bot].b,
		l2l_inst_p(deq[top],deq[top-1])),0.0) )	top--;
	while( bot < top && dy(crossProduct(deq[top].a,deq[top].b,
		l2l_inst_p(deq[bot],deq[bot+1])),0.0) )	bot++;
	
	if( top <= bot + 1 ) return ;
	
	for(int i=bot; i<top; i++)
		s[len++] = l2l_inst_p(deq[i],deq[i+1]);
	if( bot < top + 1 ) s[len++] = l2l_inst_p(deq[bot],deq[top]);
	len = unique(s,s+len,equal_p) - s;
}

void DFS(int x,int n,int sum)
{
	if( mmin == 0.0 ) return ;
	if( n - x < k - sum ) return ;
	if( sum == k )
	{
		int len;
		memcpy(ltmp,ln,sizeof(ln));
		inst_hp_nlogn(ltmp,n,s,len);
		mmin = min(mmin,area_polygon(s,len));
		if( mmin == 0.0 ) return ;
		return ;
	}
	for(int i=x; i<n; i++)
	{
		ln[i] = ll[i];
		DFS(i+1,n,sum+1);
		if( mmin == 0.0 ) return ;
		ln[i] = lk[i];
	}
}

int main()
{
	int n;
	
	while( ~scanf("%d%d%lf",&n,&k,&h) )
	{
		if( !(n || k || h) ) break;
		for(int i=0; i<n; i++)
			scanf("%lf%lf",&p[i].x,&p[i].y);
		if( h == 0 || k == 0 )
		{
			printf("0.00\n");
			continue;
		}
		p[n] = p[0];
		double area = area_polygon(p,n);
		for(int i=0; i<n; i++)
		{
			point c,d;
			changepoint(p[i],p[i+1],c,d);
			makeline_hp(c,d,ll[i]);
			makeline_hp(p[i],p[i+1],ln[i]);
		}
		memcpy(lk,ln,sizeof(ln));
		if( k > n ) k = n;
		mmin = inf;
		DFS(0,n,0);
		double ans = area - mmin;
		printf("%.2lf\n",ans);
	}
return 0;
}


#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#include <queue>
#include <stack>
#include <climits>

using namespace std;

const int MAX = 25;
const double inf = 1e20;
const double pi = acos(-1.0);
struct point {double x,y;};
struct line { double ang; point a,b;};
point p[MAX],s[MAX];
line ln[MAX],deq[MAX],ltmp[MAX],ll[MAX];
double h,mmin;
const double eps = 1e-6;
int k,cnt;
bool dy(double x,double y)	{	return x > y + eps;}	// x > y 
bool xy(double x,double y)	{	return x < y - eps;}	// x < y 
bool dyd(double x,double y)	{ 	return x > y - eps;}	// x >= y 
bool xyd(double x,double y)	{	return x < y + eps;} 	// x <= y 
bool dd(double x,double y) 	{	return fabs( x - y ) < eps;}  // x == y
double disp2p(point a,point b) //  a b 两点之间的距离 
{
	return sqrt( ( a.x - b.x ) * ( a.x - b.x ) + ( a.y - b.y ) * ( a.y - b.y ) );
}
double crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 
{
	return (c.x - a.x)*(b.y - a.y) - (b.x - a.x)*(c.y - a.y);
}
void changepoint(point a,point b,point &c,point &d)
{
	double len = disp2p(a,b);
    double dx = h / len * ( a.y - b.y );
	double dy = h / len * (-a.x + b.x );
	c.x = a.x + dx; c.y = a.y + dy;
	d.x = b.x + dx; d.y = b.y + dy;
}
void makeline_hp(point a,point b,line &l)
{
	l.a = a; l.b = b;
	l.ang = atan2(a.y - b.y,a.x - b.x);
}
bool equal_ang(line a,line b)	// 第一次unique的比较函数 
{
	return dd(a.ang,b.ang);
}
bool cmphp(line a,line b)	// 排序的比较函数 
{
	if( dd(a.ang,b.ang) ) return xy(crossProduct(b.a,b.b,a.a),0.0);
	return xy(a.ang,b.ang);
}
bool equal_p(point a,point b)//第二次unique的比较函数 
{
	return dd(a.x,b.x) && dd(a.y,b.y);
}
bool parallel(line u,line v)
{
	return dd( (u.a.x - u.b.x)*(v.a.y - v.b.y) - (v.a.x - v.b.x)*(u.a.y - u.b.y) , 0.0 );
}
point l2l_inst_p(line l1,line l2)
{
	point ans = l1.a;
	double t = ((l1.a.x - l2.a.x)*(l2.a.y - l2.b.y) - (l1.a.y - l2.a.y)*(l2.a.x - l2.b.x))/
			   ((l1.a.x - l1.b.x)*(l2.a.y - l2.b.y) - (l1.a.y - l1.b.y)*(l2.a.x - l2.b.x));
	ans.x += (l1.b.x - l1.a.x)*t;
	ans.y += (l1.b.y - l1.a.y)*t;
	return ans;
}
double area_polygon(point p[],int n)
{
	if( n < 3 ) return 0.0;
	double s = 0.0;
	for(int i=0; i<n; i++)
		s += p[(i+1)%n].y * p[i].x - p[(i+1)%n].x * p[i].y;
	return fabs(s)/2.0;
}
line lk[MAX];
void inst_hp_nlogn(line *ln,int n,point *s,int &len)
{
	len = 0;
	sort(ln,ln+n,cmphp);
	n = unique(ln,ln+n,equal_ang) - ln;
	int bot = 0,top = 1;
	deq[0] = ln[0]; deq[1] = ln[1];
	for(int i=2; i<n; i++)
	{
		if( parallel(deq[top],deq[top-1]) || parallel(deq[bot],deq[bot+1]) )
			return ;
		while( bot < top && dy(crossProduct(ln[i].a,ln[i].b,
			l2l_inst_p(deq[top],deq[top-1])),0.0) )
			top--;
		while( bot < top && dy(crossProduct(ln[i].a,ln[i].b,
			l2l_inst_p(deq[bot],deq[bot+1])),0.0) )
			bot++;
		deq[++top] = ln[i];
	}
	while( bot < top && dy(crossProduct(deq[bot].a,deq[bot].b,
		l2l_inst_p(deq[top],deq[top-1])),0.0) )	top--;
	while( bot < top && dy(crossProduct(deq[top].a,deq[top].b,
		l2l_inst_p(deq[bot],deq[bot+1])),0.0) )	bot++;
	
	if( top <= bot + 1 ) return ;
	
	for(int i=bot; i<top; i++)
		s[len++] = l2l_inst_p(deq[i],deq[i+1]);
	if( bot < top + 1 ) s[len++] = l2l_inst_p(deq[bot],deq[top]);
	len = unique(s,s+len,equal_p) - s;
}
int t[150000][10];
int tt[10];
void DFS(int x,int n,int sum)
{
	if( n - x < k - sum ) return ;
	if( sum == k )
	{
		memcpy(t[cnt++],tt,sizeof(tt));
		return ;
	}
	for(int i=x; i<n; i++)
	{
		tt[sum] = i;
		DFS(i+1,n,sum+1);
	}
}

int main()
{
	int n,len;
	point c,d;
	while( ~scanf("%d%d%lf",&n,&k,&h) )
	{
		if( !(n || k || h) ) break;
		cnt = 0;

		for(int i=0; i<n; i++)
			scanf("%lf%lf",&p[i].x,&p[i].y);
		if( h == 0 || k == 0 )
		{
			printf("0.00\n");
			continue;
		}
		p[n] = p[0];
		double area = area_polygon(p,n);
		for(int i=0; i<n; i++)
		{
			changepoint(p[i],p[i+1],c,d);
			makeline_hp(c,d,ll[i]);
			makeline_hp(p[i],p[i+1],ln[i]);
		}
		if( k > n )	k = n;
		mmin = inf;
		DFS(0,n,0);
		for(int i=0; i<cnt; i++)
		{
			memcpy(ltmp,ln,sizeof(ln));
			for(int j=0; j<k; j++)
				ltmp[t[i][j]] = ll[t[i][j]];
			inst_hp_nlogn(ltmp,n,s,len);
			mmin = min(mmin,area_polygon(s,len));
			if( mmin == 0 )
				break;
		}
		double ans = area - mmin;
		printf("%.2lf\n",ans);
	}
return 0;
}


zxy_snow
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
专栏目录
POJ 1265 pick 定理
MR CODER
08-28 1044
首先来脑补一下几个几何知识 非常奇妙的证明：图形必在格点之外 从cut-the-knot上看到的。 问题：     设想一个平面上布满间距为1的横纵直线，形成由一个个1×1正方形组成的网格。任意给一个面积小于1个单位的图形，证明这个图形总能放在网格中而不包含任何一个格点。         证明：     我们可以这样考虑这个问题：把图形随意放在网格中，如何移动网格使每个格
Nice Milk POJ - 1271 UVA 10117 半平面相交
wood
09-01 575
题目链接: POJ1271 , UVa 10117 Little Tomy likes to cover his bread with some milk. He does this by putting it in the cup so that one of its sides (called bottom side) touches the bottom of the cup, just as the pi
POJ 1271 && uva 10117 Nice Milk
05-06 439
题意：给你一个凸多边形的面包，然后
UVa Problem 10117 Nice Milk （美味的牛奶）
寂静山林
11-16 2855
// Nice Milk （美味的牛奶） // PC/UVa IDs: 111408/10117, Popularity: C, Success rate: low Level: 4 // Verdict: Accepted // Submission Date: 2011-11-15 // UVa Run Time: 0.456s // // 版权所有（C）2011，邱秋。metaphysis
poj1271 Nice Milk （dfs+半平面交）
java开发指南博客 【转载】
08-03 127
//题目链接： http://poj.org/problem?id=1271 //题意：给一个凸多边形的面包，和一瓶深度为h的牛奶，要沾牛奶k次，求可以沾到牛奶的面包的面积 //思路：保存各个边被牛奶沾后的位置，用dfs找出要沾的边，然后对多边形执行k个半面交操作（这里我用O(n^2)的联机算法），求面积。。。 //曾经纠结了半天，现在重写了一下，没特判wa了一次。。。。。。//代码如下: ...
Vijos 1217-P1217乒乓球
小I的博客
03-01 1015
华华通过以下方式进行分析，首先将比赛每个球的胜负列成一张表，然后分别计算在11分制和21分制下，双方的比赛结果（截至记录末尾）。 比赛规则和实际比赛规则相同, 在11分制下比分为10:11时还要继续比赛，直到分数差距达到2分；同理21分制下比分为20:21的时候也还要继续比赛。 比如现在有这么一份记录，（其中W表示华华获得一分，L表示华华对手获得一分）： WWWWWWWWWWWWWWWWWWWWWWLW 在11分制下，此时比赛的结果是华华第一局11比0获胜，第二局11比0获胜，正在进行第三局，当前比分1比1
POJ.rar_poj java_poj1048
最新发布
09-24
【标题】"POJ.rar_poj java_poj1048" 涉及的知识点主要围绕编程竞赛中的“约瑟夫环”问题，这里是一个加强版，使用Java语言进行解决。 【描述】"POJ1048，加强版的约瑟夫问题 难度中等" 提示我们，这个问题是编程...
poj.rar_poj
09-23
标题中的"poj.rar_poj"暗示了这是一个与POJ（Programming Online Judge）相关的压缩文件，POJ是一个在线编程竞赛平台，主要供程序员们练习和提交算法解决方案。在这个压缩包中，很可能包含了用户在POJ上参与编程挑战...
POJ.rar_poj
09-19
【标题】"POJ.rar_poj" 指的是一份与编程竞赛相关的压缩文件，主要涉及POJ（Problem Overloading Judge）平台上的问题解答。POJ是一个在线的编程竞赛平台，它为参赛者提供了大量的算法题目进行练习和比赛，以提升...
poj-1009.rar_poj
09-23
【标题】"POJ-1009"是北大在线编程竞赛平台上的一个题目，它涉及到计算机科学中的算法和编程技巧。POJ（Problem Set of Peking University）是中国北京大学主办的一个在线编程竞赛平台，旨在提高学生的算法设计和...
POJ 1271 Nice Milk 笔记
woniupengpeng的专栏
05-21 231
n多边形面包（给出顶点坐标）沾k次高为h的牛奶，求沾牛奶面包的最大面积。
poj1271——半平面交+dfs枚举
mumei的博客
09-30 176
题目链接：poj.org/problem?id=1271 Little Tomy likes to cover his bread with some milk. He does this by putting it in the cup so that one of its sides (called bottom side) touches the bottom of the cup, ju...
HDU-1271 整数对
飘过的小牛
11-25 4170
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1271 解题思路： 这道题看之后没一点思路，感觉自己跟外星来的差不多。所以就搜下解题报告，发现网上的代码也几乎一样。。估计思路太BT，没几个人做出来吧。。。。 下面就稍微分析一下吧。。。。。 首先假设X的第k位拿走，然后加上加上X的和正好等于N！ 这样的话 我们可以把X 分解成：X= a+
UVA 10173 Smallest Bounding Rectangle(旋转卡壳求最小面积外接矩形)
小媛在努力~
06-15 3965
神呐。。。终于过了 T T 。。。 一个卡壳卡两天啊啊啊啊 ！！！至少两天啊啊啊啊 ！！！T T 。。。 这两天都是，不想复习了，就在纸上画画卡壳。。。画画凸包。。。然后想这个题怎么实现比较好。。。 我开始还是枚举旋转最小角度，代码十分之繁琐。。。其实我写差不多了，但是应该是有问题的，如果凸包是三角形，怎么都枚举不对T T 。。 丫丫的，索性今天下午想了想，改成每次
poj 1264 || UVA 109 SCUD Busters
小媛在努力~
06-09 2516
UVA这题过得蛮多人的，觉得应该是个水题吧。 一读题，也不麻烦，给你N个王国，求下凸包，再求面积。给你一些炮弹，问炮弹炸掉的面积。（一个炮弹炸的话，整个王国都被炸了）。 直接求凸包后，求出各个王国的面积，然后判断炮弹在哪个王国里，这个直接用判断点是否在多边形内。 UVA WA了一次，如果多个炮弹打到同一个王国里，面积只算一次的，改了后，A掉了。 poj上一搜，居然也
poj 1410 || UVA 191 Intersection
小媛在努力~
06-23 2133
判断线段是否和矩形相交，如果在矩形里也算相交。 初始化的坐标赋值写错了一个数，WA了半天T T 。。还有就是，输入的矩形坐标不一定是左上和右下。。真黑啊  #include #include #include #include #include #include #include #include #include #include us
UVA 356 Square Pegs And Round Holes
小媛在努力~
06-25 2112
给你n，让你求2n*2n的半径为2n-1的圆里面有多少个格子还有多少个格子被圆的线盖着。 开始以为有公式，结果自己画图来看，园内格子数不符合公式T T 。。不过最多就150，套两层循环求格子数就好了。 #include #include #include #include #include #include #include #include #in
UVA 10678 The Grazing Cow
小媛在努力~
06-25 1882
水题啊水题。。 一头牛被栓在两个柱子之间，牛可以在绳子上来回移动，问它能吃到草的最大面积。 高中学椭圆的定义，椭圆是平面上到两定点的距离之和为常值的点之轨迹， 也可定义为到定点距离与到定直线间距离之比为常值的点之轨迹。 所以可知，它的轨迹一定是个椭圆，然后计算出长半轴和短半轴长度求面积（PI*a*b）即可。  #include #include #inc
UVa 11045 My T-shirt suits me
小媛在努力~
03-14 1782
第二道UVA的题目，1A，嘻嘻。类似之前做的那个衣服分配问题，不过这个比那个还简单了。每个人只有两种衣服适合，而且每种衣服的数量是相等的。建图很随意哈，把每种衣服拆成一件一件的，然后跟队员连线求最大二分匹配即可。#include #include #include #define MAX 650 using namespace std; int map[MAX][MAX]; int used[MAX],mat[MAX]; char Size[7][5] = {"","XXL","XL",
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值