最大三角形
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;
}