[Neerc2013]Green Energy(几何)

题目大意：

在二维世界中，有片地方的地形是折线形的，并且顶点的横坐标严格单调上升，顶点不超过10 000个。 在这么一片地方要修建不超过10 000座太阳能塔。每座塔的塔高均已确定。阳光的角度是 α。 太阳能塔和地形之间会相互遮盖。求一个建设方案，使得每座塔被阳光照射到的长度总和最大。

 

题解：

不太常见的几何技巧。做下面的变换：x′=x；y′=y+xtanα（另外其实不做这个变换本题也可做，只是略显复杂）

这个就是把每条边按比例延长

 

 

……就可以将阳光的角度修正为水平的！

 

 

之后修整地形。在下图所示的坑中放塔是不划算的，也就是塔只放在红线的区域内。

 

之后先将最高的塔放在最高的地方。然后从最左开始，依次将塔放在当前最左的不会被别的塔覆盖的地方。

如果覆盖满了，就将剩余的塔全扔在最后就行了。

 

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define PI acosl(-1)
using namespace std;
const int maxn = 100110;
const double eps = 1e-9;
struct Tpoint{
	long double x, y;
}a[maxn];
long double ans[maxn]; 
int n, m;
struct node{
	int i, h;
}t[maxn];
int k;
void init(){
	memset(ans, 0, sizeof(0));
	memset(a, 0, sizeof(a));
	memset(t, 0, sizeof(t));
	for (int i = 1; i <= n; i ++){
		int a = 0; scanf("%d", &a); 
		t[i].h = a;  t[i].i = i;
	}
	for (int i = 1; i <= m; i ++){
		int x, y; 
		scanf("%d%d", &x, &y);
		a[i].x = x, a[i].y = y + x * tanl(k * PI/ 180);
	}
}
bool cmp(const node &a, const node &b){ return a.h < b.h; }
void work(){
	int newm = m;
	for (int i = 1; i <= m; i ++)
		if (a[i].y > a[newm].y) newm = i;
	m = newm;
	sort(t + 1, t + n + 1, cmp);
	long double last = a[1].y;
	long double la1 = a[1].y; int now = 1;
	for (int i = 1; i < m; i ++){
		while (a[i].x <= a[i + 1].x && a[i + 1].y >= last && now <= n){
			long double x1 = a[i].x, y1 = a[i].y, x2 = a[i + 1].x, y2 = a[i + 1].y;
			long double k = (y2 - y1) / (x2 - x1), b = y2 - k * x2;
			ans[t[now].i] = (last - b) / k;
			last += t[now].h;
			now ++;
		}
	}
	for (int i = now; i <= n; i ++) ans[t[i].i] = a[m].x;
	long double aa = 0; a[1].y = la1; 
	if (now <= n && last <= a[m].y) aa = last - a[1].y + t[n].h;
	if (now <= n && a[m].y < last && last <= a[m].y + t[n].h) aa = a[m].y + t[n].h - a[1].y;
	if (now > n) aa = last - a[1].y;

	printf("%.15f\n", (double)aa); 
	for (int i = 1; i <= n; i ++)
		printf("%.15f\n", (double)ans[i]);
}
int main(){
	while (scanf("%d%d%d", &n, &m, &k) != EOF)
 		init(),work();
	return 0;
}