最大三角形
Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3932 Accepted Submission(s): 1387
Eddy对这道题目百思不得其解,想不通用什么方法来解决,因此他找到了聪明的你,请你帮他解决这个题目。
每组输出占一行。
思路:利用三角形外接圆,即找到一个最小的圆能够包围所有的点(利用凸包),那么在这个圆上至少存在三个点以上,所以,此时就可以使用暴力求出最大三角形面积了;
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<iomanip>
using namespace std;
const int Max=50005;
struct Node
{
int x;
int y;
}a[Max],b[Max];
int cmp(Node m,Node n)
{
if(m.x==n.x)
return m.y<n.y;
else
return m.x<n.x;
}
int Count(Node i,Node j,Node k) //就算三角形的面积;
{
int x1=i.x-j.x;
int y1=i.y-j.y;
int x2=k.x-j.x;
int y2=k.y-j.y;
return(x1*y2-x2*y1);
}
int convex(int n) //凸包计算;求出在圆上的点得数量m;和座标。
{
sort(a,a+n,cmp);
int m=0,i,j,k;
for(i=0;i<n;i++)
{
while(m>1&&Count(b[m-1],a[i],b[m-2])<=0)
m--;
b[m++]=a[i];
}
k=m;
for(i=n-2;i>=0;i--)
{
while(m>k&&Count(b[m-1],a[i],b[m-2])<=0)
m--;
b[m++]=a[i];
}
if(n>1)
m--;
return m;
}
int main()
{
int i,j,k,n,m;
int sum=0;
while(cin>>n)
{
for(i=0;i<n;i++)
cin>>a[i].x>>a[i].y;
sum=0;
m=convex(n);
//cout<<m<<endl;
for(i=0;i<m;i++)
for(j=i+1;j<m;j++)
for(k=j+1;k<m;k++)
{
sum=max(sum,Count(b[i],b[j],b[k]));
}
//sum=sum/2.0;
printf("%.2lf\n",0.5*sum);
// cout<<fixed<<setprecision(2)<<sum<<endl;
}
return 0;
}
最大三角形
Time Limit: 5000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3932 Accepted Submission(s): 1387
Eddy对这道题目百思不得其解,想不通用什么方法来解决,因此他找到了聪明的你,请你帮他解决这个题目。
每组输出占一行。
思路:利用三角形外接圆,即找到一个最小的圆能够包围所有的点(利用凸包),那么在这个圆上至少存在三个点以上,所以,此时就可以使用暴力求出最大三角形面积了;
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<iomanip>
using namespace std;
const int Max=50005;
struct Node
{
int x;
int y;
}a[Max],b[Max];
int cmp(Node m,Node n)
{
if(m.x==n.x)
return m.y<n.y;
else
return m.x<n.x;
}
int Count(Node i,Node j,Node k) //就算三角形的面积;
{
int x1=i.x-j.x;
int y1=i.y-j.y;
int x2=k.x-j.x;
int y2=k.y-j.y;
return(x1*y2-x2*y1);
}
int convex(int n) //凸包计算;求出在圆上的点得数量m;和座标。
{
sort(a,a+n,cmp);
int m=0,i,j,k;
for(i=0;i<n;i++)
{
while(m>1&&Count(b[m-1],a[i],b[m-2])<=0)
m--;
b[m++]=a[i];
}
k=m;
for(i=n-2;i>=0;i--)
{
while(m>k&&Count(b[m-1],a[i],b[m-2])<=0)
m--;
b[m++]=a[i];
}
if(n>1)
m--;
return m;
}
int main()
{
int i,j,k,n,m;
int sum=0;
while(cin>>n)
{
for(i=0;i<n;i++)
cin>>a[i].x>>a[i].y;
sum=0;
m=convex(n);
//cout<<m<<endl;
for(i=0;i<m;i++)
for(j=i+1;j<m;j++)
for(k=j+1;k<m;k++)
{
sum=max(sum,Count(b[i],b[j],b[k]));
}
//sum=sum/2.0;
printf("%.2lf\n",0.5*sum);
// cout<<fixed<<setprecision(2)<<sum<<endl;
}
return 0;
}