HDOJ 2202 最大三角形 凸包旋转卡壳求最大三角形面积
2014-12-28 20:17
375 查看
凸包旋转卡壳求最大三角形面积
Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3316 Accepted Submission(s): 1119
Problem Description
老师在计算几何这门课上给Eddy布置了一道题目,题目是这样的:给定二维的平面上n个不同的点,要求在这些点里寻找三个点,使他们构成的三角形拥有的面积最大。
Eddy对这道题目百思不得其解,想不通用什么方法来解决,因此他找到了聪明的你,请你帮他解决这个题目。
Input
输入数据包含多组测试用例,每个测试用例的第一行包含一个整数n,表示一共有n个互不相同的点,接下来的n行每行包含2个整数xi,yi,表示平面上第i个点的x与y坐标。你可以认为:3 <= n <= 50000 而且 -10000 <= xi, yi <= 10000.
Output
对于每一组测试数据,请输出构成的最大的三角形的面积,结果保留两位小数。
每组输出占一行。
Sample Input
3
3 4
2 6
3 7
6
2 6
3 9
2 0
8 0
6 6
7 7
Sample Output
1.50
27.00
Author
Eddy
/* ***********************************************
Author :CKboss
Created Time :2014年12月28日 星期日 19时28分15秒
File Name :HDOJ2202.cpp
************************************************ */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <queue>
#include <set>
#include <map>
using namespace std;
struct Point
{
double x,y;
Point operator - (Point a) const
{
return (Point){x-a.x,y-a.y};
}
bool operator < (Point a) const
{
if(x!=a.x) return x<a.x;
return y<a.y;
}
bool operator == (Point a) const
{
return (x==a.x)&&(y==a.y);
}
};
int n,m;
double Cross(Point A,Point B) { return A.x*B.y-A.y*B.x; }
double Area2(const Point A,const Point B,const Point C)
{
return fabs(Cross(B-A,C-A));
}
vector<Point> vp;
vector<Point> ch;
void CovexHull(vector<Point> &p)
{
sort(p.begin(),p.end());
p.erase(unique(p.begin(),p.end()),p.end());
int n=p.size();
m=0;
ch=vector<Point>(n+1);
for(int i=0;i<n;i++)
{
while(m>1&&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<0) m--;
ch[m++]=p[i];
}
int k=m;
for(int i=n-2;i>=0;i--)
{
while(m>k&&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<0) m--;
ch[m++]=p[i];
}
if(n>1) m--;
ch.resize(m);
}
double Rotating_Calipers(vector<Point>& p,int n)
{
double ans=0;
for(int i=0;i<n;i++)
{
int j=(i+1)%n;
int k=(j+1)%n;
while(j!=i&&k!=i)
{
ans=max(ans,Area2(p[i],p[j],p[k]));
while(Cross(p[i]-p[j],p[(k+1)%n]-p[k])<0)
k=(k+1)%n;
j=(j+1)%n;
}
}
return ans;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(scanf("%d",&n)!=EOF)
{
vp.clear();
for(int i=0;i<n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
vp.push_back((Point){x,y});
}
CovexHull(vp);
double area=Rotating_Calipers(ch,ch.size())/2.;
printf("%.2lf\n",area);
}
return 0;
}
最大三角形
Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3316 Accepted Submission(s): 1119
Problem Description
老师在计算几何这门课上给Eddy布置了一道题目,题目是这样的:给定二维的平面上n个不同的点,要求在这些点里寻找三个点,使他们构成的三角形拥有的面积最大。
Eddy对这道题目百思不得其解,想不通用什么方法来解决,因此他找到了聪明的你,请你帮他解决这个题目。
Input
输入数据包含多组测试用例,每个测试用例的第一行包含一个整数n,表示一共有n个互不相同的点,接下来的n行每行包含2个整数xi,yi,表示平面上第i个点的x与y坐标。你可以认为:3 <= n <= 50000 而且 -10000 <= xi, yi <= 10000.
Output
对于每一组测试数据,请输出构成的最大的三角形的面积,结果保留两位小数。
每组输出占一行。
Sample Input
3
3 4
2 6
3 7
6
2 6
3 9
2 0
8 0
6 6
7 7
Sample Output
1.50
27.00
Author
Eddy
/* ***********************************************
Author :CKboss
Created Time :2014年12月28日 星期日 19时28分15秒
File Name :HDOJ2202.cpp
************************************************ */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstdlib>
#include <vector>
#include <queue>
#include <set>
#include <map>
using namespace std;
struct Point
{
double x,y;
Point operator - (Point a) const
{
return (Point){x-a.x,y-a.y};
}
bool operator < (Point a) const
{
if(x!=a.x) return x<a.x;
return y<a.y;
}
bool operator == (Point a) const
{
return (x==a.x)&&(y==a.y);
}
};
int n,m;
double Cross(Point A,Point B) { return A.x*B.y-A.y*B.x; }
double Area2(const Point A,const Point B,const Point C)
{
return fabs(Cross(B-A,C-A));
}
vector<Point> vp;
vector<Point> ch;
void CovexHull(vector<Point> &p)
{
sort(p.begin(),p.end());
p.erase(unique(p.begin(),p.end()),p.end());
int n=p.size();
m=0;
ch=vector<Point>(n+1);
for(int i=0;i<n;i++)
{
while(m>1&&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<0) m--;
ch[m++]=p[i];
}
int k=m;
for(int i=n-2;i>=0;i--)
{
while(m>k&&Cross(ch[m-1]-ch[m-2],p[i]-ch[m-2])<0) m--;
ch[m++]=p[i];
}
if(n>1) m--;
ch.resize(m);
}
double Rotating_Calipers(vector<Point>& p,int n)
{
double ans=0;
for(int i=0;i<n;i++)
{
int j=(i+1)%n;
int k=(j+1)%n;
while(j!=i&&k!=i)
{
ans=max(ans,Area2(p[i],p[j],p[k]));
while(Cross(p[i]-p[j],p[(k+1)%n]-p[k])<0)
k=(k+1)%n;
j=(j+1)%n;
}
}
return ans;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(scanf("%d",&n)!=EOF)
{
vp.clear();
for(int i=0;i<n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
vp.push_back((Point){x,y});
}
CovexHull(vp);
double area=Rotating_Calipers(ch,ch.size())/2.;
printf("%.2lf\n",area);
}
return 0;
}
相关文章推荐
- POJ 2079 Triangle(凸包+旋转卡壳,求最大三角形面积)
- poj 2079 Triangle 凸包+旋转卡壳 求最大三角形面积
- ZOJ 2419-- Triangle-凸包+旋转卡壳求最大面积三角形(计算几何)
- hdu 3934&&poj 2079 (凸包+旋转卡壳+求最大三角形面积)
- HDU 2202 最大三角形 (凸包+旋转卡壳)
- poj 2079(旋转卡壳求解凸包内最大三角形面积)
- POJ 2079 Triangle(凸包_旋转卡壳之最大三角形面积)
- POJ 2079 求最大三角形面积 (凸包+旋转卡壳)
- POJ 2079 Triangle(凸包+旋转卡壳,求最大三角形面积)
- POJ 2079 Triangle(凸包+旋转卡壳求最大三角形面积)
- BZOJ 1069: [SCOI2007]最大土地面积 凸包,旋转卡壳
- bzoj1069 [SCOI2007]最大土地面积(凸包+旋转卡壳)
- SCOI-2007 最大土地面积(凸包+旋转卡壳)
- hdu 最大三角形(凸包+旋转卡壳)
- [BZOJ1069][SCOI2007]最大土地面积 凸包+旋转卡壳
- 计算几何 ( 求凸包,计算三角形面积 )——最大三角形 ( HDU 2202 )
- 强大的旋转卡壳 POJ 2187 最远点对 POJ 2079点集中面积最大的三角形
- 最大三角形 HDU - 2202 旋转卡壳模板
- POJ 2079 Triangle (凸包中的最大三角形&旋转卡壳)
- [BZOJ1069][SCOI2007][凸包][旋转卡壳]最大土地面积