今天计算几何冲锋的第一天,选择了一道看上去颇有难度的扫描线题目。说实在的,这道题我打一开始抱着看题解的心去做的,想不到居然过了。废话不多说,咱们上题(三角形)。
题目描述
给出平面上的n个等腰直角三角形。每个三角形用3个整数描述x, y, m(m>0)。一个三角形的3个顶点分别是(x, y),(x+m, y)and(x, y+m)。你的任务是计算这些三角形覆盖的总面积。
输入格式
输入第一行是整数n(n≤2000)。
接下来n行每行描述一个三角形,包括3个整数xi,yi,mi(1≤i≤n,-10^7 ≤xi≤10^7,
-10^7 ≤yi≤10^7,0<mi≤1000)。
题目描述的意思是,在二维直角坐标系上面,有很多等腰直角三角形,求这些实心三角形覆盖掉的平面面积有多少。
首先,我们发现每个三角形可能会很大,但是三角形的数量却非常的少,因此考虑对坐标进行离散化,先以x为优先级再以y为优先级,最后以m为优先级从小到大把所有的三角形进行排序。
这样做的目的是,当我们处理前面的三角形的时候,后面的三角形不会产生影响。有人可能会问,不对啊,两个三角形重合是相互的,为什么说后面的不会对前面产生影响呢?因为我们这么规定,可以方便处理,也就是说,我们只处理当前加进去的三角形与前面的所有三角形重合的部分而不处理与后面的三角形的重合部分,后面的等后面再说。
那么就非常显然了,我们每一次处理一排x相等的三角形,从底部开始放入,最下面这个直接放进去就可以。注意,这里放进去之后并不直接计算他的贡献,它的贡献要等到下一个x才会求,也就是说我们实际上记录的是两个离散了之后的x之间的那一部分三角形切片的面积。因此每次放进去的是这个三角形的y坐标与其高度(因为是等腰直角三角形,因此只需要这两个值就可以了)。
然后,处理本次之前放进去的三角形的贡献。对于底部的那个三角形,相当于从上一个x到现在这个x要么全部放入,要么切出来一个直角梯形,两者的面积都非常好酸。而后面计算的贡献,只需要剪除与前面的那个三角形的面积交叉值就可以,注意全部重合乃至于包含的情况,此时贡献为0而非负数。
代码如下:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> P;
struct delta
{
ll x, y;
ll m;
} p[2020];
bool cmp(delta a,delta b)
{
if(a.x==b.x)
return a.y < b.y;
else
return a.x < b.x;
}
template<typename T1, typename T2>
bool operator<(const pair<T1, T2> &a, const pair<T1, T2> &b)
{
if(a.first==b.first)
return a.second < b.second;
else
return a.first < b.first;
}
priority_queue<P,vector<P>, greater<P> > que;
queue<P> temp;
int main()
{
ll n;
cin >> n;
for (ll i = 1; i <= n;i++)
{
cin >> p[i].x >> p[i].y >> p[i].m;
}
sort(p + 1, p + 1 + n, cmp);
ll ans = 0;
ll pre = p[1].x;
for (ll i = 1; i <= n||!que.empty(); i++)
{
ll x = p[i].x;
ll del = x - pre;
if(i==n+1)
del = 1e9+10;
ll last_y = -1e9+10;
ll last_m = 0;
//cout << x << " " << del << " :";
while(!que.empty())
{
P top = que.top();
que.pop();
ll add = (top.second + max(0ll, top.second - del)) * min(top.second, del);
if(last_y>top.first)
{
add-=(last_y-top.first + max(0ll, last_y-top.first - del)) * min(last_y-top.first,del);
//cout << "?";
}
last_y = top.first + top.second;
if(del<top.second)
{
temp.push({top.first, top.second - del});
}
last_m = top.second;
ans += max(0ll, add);
//cout << ans << " ";
}
while(!temp.empty())
{
que.push(temp.front());
temp.pop();
}
while(i<=n&&p[i].x==x)
{
que.push({p[i].y, p[i].m});
i++;
}
i--;
//cout << ans << endl;
pre = x;
//cout << endl;
}
//cout << endl;
cout << ans / 2<<".";
if(ans&1)
cout << "5";
else
cout << "0";
cout << endl;
}
本题难点在于扫描线扫面的时候究竟要怎么处理扫描前后这片区域中三角形的贡献值。但是对于这些扫面积的题目,我们不应该看他们的整体长啥样,更多地应该看他当前状态的贡献。