给出点集,构成曲折线,如下图,每两个线构成一个三角形,三点共线认为面积为零,每次删掉面积最小的三角形和其顶点,余下两点自动连接,如下图,删掉面积S2(假设S2面积最小),p3也被删去,形成下面的情况,一直删到剩余边数为m时停止,把每次删的点输出来;按题意模拟即可;利用向量叉乘求出面积存到set或priority_queue中,记录每个点前驱后驱,删除时不断更新;
坑点:将一个三角形删除后,会产生一个或两个新的三角形;
AC代码:
#include <iostream>
#include <algorithm>
#include <string.h>
#include <set>
#include <math.h>
using namespace std;
const int maxn=3e5+10;
int sto[maxn],pre[maxn],Next[maxn],vis[maxn];
struct point {
double x;
double y;
point(double x=0,double y=0):x(x),y(y){}
}pt[maxn];
struct ppt{
int l;
int mid;
int r;
double area;
bool operator < (const struct ppt &a)const
{
if (this->area == a.area)
return this->mid < a.mid;
else
return this->area < a.area;
}
}pp,p1;
multiset<ppt>se;
multiset<ppt>::iterator it;
point operator - (point a,point b)
{
return point (a.x-b.x,a.y-b.y);
}
double cross(point a,point b)
{
return a.x*b.y-a.y*b.x;
}
double area(point a,point b,point c)
{
return cross(b-a,c-a);
}
int main()
{
int n,m;
memset(vis,0,sizeof(vis));
cin>>n>>m;
for (int i=0;i<=n;i++)
{
cin>>pt[i].x>>pt[i].y;
pre[i]=i-1;
Next[i]=i+1;
}
ppt ps;
for (int i=1;i<=n-1;i++)
{
ps.l=i-1;
ps.mid=i;
ps.r=i+1;
ps.area=fabs(area(pt[i-1],pt[i],pt[i+1]));
//cout<<ps.area<<"%%%%%%"<<endl;
//q.push(ps);
se.insert(ps);
}
int re=n+1,k=0;
while (re>m+1)
{
// ppt pp=q.top();
// q.pop();
it=se.begin();
pp=*it;
se.erase(it);
if (vis[pp.l]==1||vis[pp.r]==1)
continue;
sto[k++]=pp.mid;
vis[pp.mid]=1;
re--;
if (re==2)
break;
if (pre[pp.l]==-1)
{
int nnex=Next[pp.r];
ppt qaq;qaq.l=pp.l;qaq.mid=pp.r;qaq.r=nnex;qaq.area=fabs(area(pt[qaq.l],pt[qaq.mid],pt[qaq.r]));se.insert(qaq);
Next[pp.l]=pp.r;
pre[pp.r]=pp.l;
pre[pp.l]=-1;
}
else if (Next[pp.r]==n+1)
{
int ppre=pre[pp.l];
ppt qaq;qaq.l=ppre;qaq.mid=pp.l;qaq.r=pp.r;qaq.area=fabs(area(pt[qaq.l],pt[qaq.mid],pt[qaq.r]));se.insert(qaq);
pre[pp.r]=pp.l;
Next[pp.l]=pp.r;
Next[pp.r]=n+1;
}
else
{
int nnex=Next[pp.r];
ppt qaq;qaq.l=pp.l;qaq.mid=pp.r;qaq.r=nnex;qaq.area=fabs(area(pt[qaq.l],pt[qaq.mid],pt[qaq.r]));se.insert(qaq);
int ppre=pre[pp.l];
qaq.l=ppre;qaq.mid=pp.l;qaq.r=pp.r;qaq.area=fabs(area(pt[qaq.l],pt[qaq.mid],pt[qaq.r]));se.insert(qaq);
pre[pp.r]=pp.l;
Next[pp.l]=pp.r;
}
}
for (int i=0;i<k;i++)
cout<<sto[i]<<endl;
return 0;
}