POJ 1436 Horizontally Visible Segments(线段树区间修改)
最后更新于:2022-04-01 15:53:34
题目链接:[点击打开链接](http://poj.org/problem?id=1436)
题意:n条竖直线段,如果两条线段之间可见(即可以用一条水平线段连接而不触碰其他线段),则称它们可见。 如果三条线段任意两条都可见, 则称它们为a triangle of segments, 求a triangle of segments的个数
思路: 一开始真没想到n^3的复杂度可以过。。。 如果这样的话, 问题的关键就是怎样判断任意两个线段是否可见。
那么如果对y坐标建立区间, 这就成了线段树的区间覆盖问题。
另外, 由于是用点代替线段, 所以对于”可见“这个问题, 就会出现问题, 比如 有三条线段[1,4]、[1,2]、[3,4], 则在[1,4]区间中可以看见3条线段, 即可以通过[2,3]的空隙看见第一条线段。 解决方法很简单, 把坐标扩大2倍就行了,这样[1,2]变成[2,4],[3,4]变成[6,8],空出了一个5。
细节参见代码:
~~~
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
typedef long long ll;
const double PI = acos(-1.0);
const double eps = 1e-6;
const int mod = 1000000000 + 7;
const int INF = 1000000000;
const int maxn = 8000 + 10;
int T,n,t,q,l,r,c,col[maxn<<3],y[maxn<<1];
bool vis[maxn][maxn];
struct node {
int ly, ry, x;
node(int l=0, int r=0, int xx=0):ly(l),ry(r),x(xx) {}
bool operator < (const node& rhs) const {
return x < rhs.x || (x == rhs.x && ly < rhs.ly);
}
}a[maxn];
void pushdown(int o) {
if(col[o] >= 0) {
col[o<<1] = col[o<<1|1] = col[o];
col[o] = -1;
}
}
void update(int L, int R, int c, int l, int r, int o) {
int m = (l + r) >> 1;
if(L <= l && r <= R) {
col[o] = c;
return ;
}
pushdown(o);
if(L <= m) update(L, R, c, l, m, o<<1);
if(m < R) update(L, R, c, m+1, r, o<<1|1);
}
void query(int L, int R, int id, int l, int r, int o) {
int m = (l + r) >> 1;
if(col[o] >= 0) {
vis[id][col[o]] = true; return ;
}
if(l == r) return ;
if(L <= m) query(L, R, id, l, m, o<<1);
if(m < R) query(L, R, id, m+1, r, o<<1|1);
}
int main() {
scanf("%d",&T);
while(T--) {
scanf("%d",&n);
int m = 0;
for(int i=0;i<n;i++) {
scanf("%d%d%d",&a[i].ly, &a[i].ry, &a[i].x);
m = max(m, max(a[i].ly, a[i].ry));
}
sort(a, a+n);
memset(col, -1, sizeof(col));
memset(vis, false, sizeof(vis));
for(int i=0;i<n;i++) {
int l = a[i].ly;
int r = a[i].ry;
query(2*l, 2*r, i, 0, 2*m, 1);
update(2*l, 2*r, i, 0, 2*m, 1);
}
int ans = 0;
for(int i=n-1;i>=0;i--) {
for(int j=i-1;j>=0;j--) {
if(vis[i][j]) {
for(int k=j-1;k>=0;k--) {
if(vis[i][k] && vis[j][k]) ++ans;
}
}
}
}
printf("%d\n",ans);
}
return 0;
}
~~~