圈水池
-
描述
-
有一个牧场,牧场上有很多个供水装置,现在牧场的主人想要用篱笆把这些供水装置圈起来,以防止不是自己的牲畜来喝水,各个水池都标有各自的坐标,现在要你写一个程序利用最短的篱笆将这些供水装置圈起来!(篱笆足够多,并且长度可变)
-
输入
-
第一行输入的是N,代表用N组测试数据(1<=N<=10)
第二行输入的是m,代表本组测试数据共有m个供水装置(3<=m<=100)
接下来m行代表的是各个供水装置的横纵坐标
输出
- 输出各个篱笆经过各个供水装置的坐标点,并且按照x轴坐标值从小到大输出,如果x轴坐标值相同,再安照y轴坐标值从小到大输出 样例输入
-
1 4 0 0 1 1 2 3 3 0
样例输出
-
0 0 2 3
3 0
-
代码实现:
-
//Graham扫描(凸背包问题)
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
struct point
{
int x,y;
}pnt[105],res[105];
bool cmp(const point &p1,const point &p2)
{
return p1.x<p2.x||(p1.x==p2.x&&p1.y<p2.y);
}
//叉积函数(判断法向量的方向,即三者位置的顺序关系)
bool mult(point sp,point ep,point op)
{
return (sp.x-op.x)*(ep.y-op.y)>(sp.y-op.y)*(ep.x-op.x);
}
//凸背包算法
int graham(point pnt[],int n,point res[])
{
int i,len,top;
top=1;
sort(pnt,pnt+n,cmp);
if(n==0) return 0;
res[0]=pnt[0];
if(n==1) return 1;
res[1]=pnt[1];
if(n==2) return 2;
res[2]=pnt[2];
for(i=2;i<n;i++)
{
while(top&&mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
len=top;
res[++top]=pnt[n-2];
for(i=n-3;i>=0;i--)
{
while(top!=len&& mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
return top;
}
int main()
{
int N,m;
cin>>N;
while(N--)
{
cin>>m;
int i,j;
for(i=0;i<m;i++)
cin>>pnt[i].x>>pnt[i].y;
j=graham(pnt,m,res);
sort(res,res+j,cmp);
for(i=0;i<j;i++)
cout<<res[i].x<<" "<<res[i].y<<endl;
}
return 0;
}
-
本题实质是利用凸背包算法来解决实际问题,并且其中涉及到向量的叉乘,属于计算几何问题,对于我来说还是比较有难度的,特别是对于Graham扫描(凸背包问题)的具体代码理解。以下代码是转载别人的代码:
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
struct point
{
int x,y;
}pnt[105],res[105];
bool cmp(const point &p1,const point &p2)
{
return p1.x<p2.x||(p1.x==p2.x&&p1.y<p2.y);
}
//叉积函数(判断法向量的方向,即三者位置的顺序关系)
bool mult(point sp,point ep,point op)
{
return (sp.x-op.x)*(ep.y-op.y)>(sp.y-op.y)*(ep.x-op.x);
}
//凸背包算法
int graham(point pnt[],int n,point res[])
{
int i,len,top;
top=1;
sort(pnt,pnt+n,cmp);
if(n==0) return 0;
res[0]=pnt[0];
if(n==1) return 1;
res[1]=pnt[1];
if(n==2) return 2;
res[2]=pnt[2];
for(i=2;i<n;i++)
{
while(top&&mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
len=top;
res[++top]=pnt[n-2];
for(i=n-3;i>=0;i--)
{
while(top!=len&& mult(pnt[i],res[top],res[top-1]))
top--;
res[++top]=pnt[i];
}
return top;
}
int main()
{
int N,m;
cin>>N;
while(N--)
{
cin>>m;
int i,j;
for(i=0;i<m;i++)
cin>>pnt[i].x>>pnt[i].y;
j=graham(pnt,m,res);
sort(res,res+j,cmp);
for(i=0;i<j;i++)
cout<<res[i].x<<" "<<res[i].y<<endl;
}
return 0;
}
-
第一行输入的是N,代表用N组测试数据(1<=N<=10)