RMRC2017 Polyline Simplification

给出点集，构成曲折线，如下图，每两个线构成一个三角形，三点共线认为面积为零，每次删掉面积最小的三角形和其顶点，余下两点自动连接，如下图，删掉面积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;
}