[USACO15JAN]踩踏Stampede

题目描述

DJ站在原点上向y轴正半轴看，然后有一群奶牛从他眼前飞过。这些奶牛初始都在第二象限，尾巴在(Xi,Yi)，头在(Xi+1,Yi)，每Ci秒向右走一个单位。 DJ能看见一匹奶牛当且仅当它身体任意某部位x坐标为0时，没有其它y坐标小于此奶牛的奶牛身体某部位x坐标为0。问DJ能看见多少奶牛？

输入格式：

第一行：一个数，给出n，表示n头牛。

第二行：三个数，给出Xi和Yi，描述第i头牛的位置，给出Ci，描述这头牛多少秒走一个单位。

输出格式：

一行，给出ans，即DJ能看到多少头牛。

输入输出样例

输入样例#1：

输出样例#1：

样例说明

DJ能看到第一和第二头牛，看不到第三头牛。

思路

　　首先，每一头牛的Xi，Yi，Ci都可能非常的大，所以用普通的暴力搜索是过不去的，所以就要换一个方法去想问题，题目中给出了速度，路程，那么我们可以很容易就想到时间。以样例为例:第一头牛的头到达y轴的时间记为B1=(-2+1)*(-1)*3=3，第一头牛的尾巴到达y轴的时间记为E1=B1+3=6，就像这样我们能知道在第3秒到第6秒时第一头牛是在y轴上的；同理，我们可以知道第二头牛在第6秒到第9秒时是在y轴上的，第三头牛在第4秒到第5秒时是在y轴上的。有人会问知道这些有什么用呢？由于知道时间和纵坐标的位置，我们可以轻易的推出来任意两头牛是否有重叠，方法是当且仅当第i头牛的时间段与第j头牛的时间段有重叠，并且第i头牛的纵坐标比第j头牛的纵坐标小才能当上第j头牛。

　　由于我只需要看到这头牛的一部分就算看到，那么我们每一次只需要看一下当前牛的时间段是否被覆盖了，如果没有被全被覆盖则ans++，这头牛可以被看到，并且把这个时间段全部覆盖。那么现在有一个问题了，怎么能保证每一头牛能在恰当的时间进行验证呢？我们需要把牛按纵坐标进行排序，在一一进行验证，这样就可以保证纵坐标小的牛不被纵坐标大的牛遮挡住。既然思路已经出来，那么算法十分容易分析出来，区间查询+区间修改=线段树。我们再看一下每一个数是十分大的，如果直接进行线段树覆盖会MLE，那么是不是可以把时间离散化一下？我们把时间点排序，在更形成序号就可以啦，再看样例这些牛中出现了3,6,9,4,5五个时间点，排序过后是3,4,5,6,9，那么第一头牛就可以化为在第1秒到第4秒在y轴上，第二头牛就可以化为在第4秒到第5秒在y轴上，第一头牛就可以化为在第2秒到第3秒在y轴上。这样每一次操作就十分方便了。

　　注意：由于时间点重合可能是像2~3，3~4和4~5这样的时间段，那么如果直接按照时间点进行查询和覆盖就可能出现问题，比如说先覆盖2~3和4~5，那么查找3~4时便会出现问题，所以每一次覆盖和查询的是线段，比如2~6是覆盖和查询2~5的线段。

 1 #include<stdio.h>
 2 #include<algorithm>
 3 using namespace std;
 4 #define mid ((l+r)>>1)
 5 struct Cow{
 6     int high,order_l,order_r;
 7 }cow[50001];
 8 struct Num{
 9     int order,order2;
10     int x;
11 }num[100001];
12 bool cmp1(const Num &a,const Num &b){
13     return a.x<b.x;
14 }
15 bool cmp2(const Num &a,const Num &b){
16     return a.order<b.order;
17 }
18 bool cmp3(const Cow &a,const Cow &b){
19     if(a.high==b.high)return a.order_r-a.order_l<b.order_r-b.order_l;
20     return a.high<b.high;
21 }
22 int n;
23 int idx,idx2;
24 int ans;
25 bool tree[1000000];
26 void push(int p){
27     tree[p<<1]=tree[(p<<1)|1]=tree[p];
28 }
29 bool find(int p,int l,int r,int x,int y){
30     if(tree[p])push(p);
31     if(l>=x&&r<=y)
32         return tree[p];
33     bool is1=false,is2=false,is3=false,is4=false;
34     if(x<=mid)is1=find(p<<1,l,mid,x,y),is3=true;
35     if(y>mid)is2=find((p<<1)|1,mid+1,r,x,y),is4=true;
36     if(is1==is3&&is2==is4)return true;
37     return false;
38 }
39 void add(int p,int l,int r,int x,int y){
40     if(l>=x&&r<=y){
41         tree[p]=true;
42         return;
43     }
44     if(x<=mid)add(p<<1,l,mid,x,y);
45     if(y>mid)add((p<<1)|1,mid+1,r,x,y);
46     if(tree[p<<1]&&tree[(p<<1)|1])tree[p]=true;
47 }
48 int main(){
49     scanf("%d",&n);
50     for(int i=1;i<=n;i++){
51         int a,speed;
52         scanf("%d%d%d",&a,&cow[i].high,&speed);
53         ++idx;
54         num[idx].x=(-1-a)*speed;
55         num[idx].order=idx+1;
56         ++idx;
57         num[idx].x=num[idx-1].x+speed;
58         num[idx].order=idx+1;
59     }
60     sort(num+1,num+idx+1,cmp1);
61     for(int i=1;i<=idx;){
62         int tmp=0;
63         ++idx2;
64         while(num[i].x==num[i+tmp].x)
65             num[i+tmp].order2=idx2,++tmp;
66         i+=tmp;
67     }
68     sort(num+1,num+idx+1,cmp2);
69     for(int i=1;i<=idx;i++){
70         (num[i].order%2)?
71         cow[num[i].order/2].order_r=num[i].order2:
72         cow[num[i].order/2].order_l=num[i].order2;
73     }
74     sort(cow+1,cow+n+1,cmp3);
75     for(int i=1;i<=n;i++){
76         if((!find(1,1,idx2,cow[i].order_l,cow[i].order_r-1))){
77             ans++;
78         }
79         add(1,1,idx2,cow[i].order_l,cow[i].order_r-1);
80     }
81     printf("%d",ans);
82 }

View Code

转载于:https://www.cnblogs.com/yangsongyi/p/8664649.html