HDU 4967 Handling the Past

题意：

你有一个栈  一些操作发生在栈上  包括进栈、出栈、询问栈顶  每个操作有一个独一无二的时间  当操作读进来时  要把之前处理的本该在本操作之后的操作全撤销  接着完成现在的操作  再把撤销的操作重做一遍  每次询问操作输出栈顶元素

思路：

由于时间唯一  那么可以用时间来对应插入的元素以及操作的类型  所以首先离散化时间（这里不去掉重复也行）

然后我们把push当作+1  pop当作-1  就可以维护栈内元素个数了  再考虑peak询问  即是询问“最靠后的后缀和>0的位置”  这个位置对应的是一个时间（不过刚才我们也说了用时间可以对应到栈内元素）

考虑到peak是一个区间的操作  所以想到线段树维护  我们可以将问题转化成“对于区间[u,v]  它被询问的区间[1,x]包含  那么它里面是不是有一个位置i  可以使sum[i,v]>sum[v+1,x]”  为了计算这个信息  我们需要维护每个区间的和sum  还有内部最大的后缀和rmax  利用这些就可以求出所有的i位置  由于我们要最靠后的  所以只需先处理右边的区间  一旦找到满足的i就不再找了

代码里我还用了一个flag  它是来记录这个区间里是不是被加过值的 （感觉这样代码好写）

哦还有  我在时间序列前后各虚拟了一个点  前面的点来应付peak时栈内无元素的情况  后面的大概是脑洞大- -b

代码：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;#define N 50005#define L(x) (x<<1)#define R(x) ((x<<1)|1)struct op{    int x,y,t;}g[N];int h[N],ans[N];int n,cas=1,sum;struct tree{    int l,r,sum,rmax,flag;}f[N*4];void up(int i){    f[i].sum=f[L(i)].sum+f[R(i)].sum;    if(f[i].flag)    {        if(f[L(i)].flag&&f[R(i)].flag) f[i].rmax=max(f[R(i)].rmax,f[L(i)].rmax+f[R(i)].sum);        else if(f[L(i)].flag) f[i].rmax=f[L(i)].rmax;        else f[i].rmax=f[R(i)].rmax;    }}void init(int l,int r,int i){    f[i].l=l;    f[i].r=r;    if(l) f[i].flag=0;    else f[i].flag=1;    if(l==r)    {        if(!l) f[i].sum=f[i].rmax=N*10;        else f[i].sum=f[i].rmax=0;        return ;    }    int mid=(l+r)>>1;    init(l,mid,L(i));    init(mid+1,r,R(i));    up(i);}void add(int pos,int val,int i){    if(f[i].l==f[i].r)    {        f[i].sum=f[i].rmax=val;        f[i].flag=1;        return ;    }    f[i].flag=1;    int mid=(f[i].l+f[i].r)>>1;    if(pos<=mid) add(pos,val,L(i));    else add(pos,val,R(i));    up(i);}int dfs(int val,int i){    if(f[i].l==f[i].r) return f[i].l;    if(f[R(i)].flag&&f[R(i)].rmax>val) return dfs(val,R(i));    else return dfs(val-f[R(i)].sum,L(i));}int query(int l,int r,int i){    if(f[i].l==l&&r==f[i].r)    {        int tmp;        if(!f[i].flag||f[i].rmax<=-sum) tmp=-1;        else tmp=dfs(-sum,i);        sum+=f[i].sum;        return tmp;    }    int mid=(f[i].l+f[i].r)>>1;    if(r<=mid) return query(l,r,L(i));    else if(l>mid) return query(l,r,R(i));    else    {        int res=-1;        res=query(mid+1,r,R(i));        if(res==-1) res=query(l,mid,L(i));        return res;    }}int main(){    int i,t;    char x[10];    while(~scanf("%d",&n))    {        if(!n) break;        for(i=1;i<=n;i++)        {            scanf("%s",x);            if(!strcmp(x,"push"))            {                g[i].x=1;                scanf("%d%d",&g[i].y,&g[i].t);            }            else if(!strcmp(x,"pop"))            {                g[i].x=-1;                scanf("%d",&g[i].t);            }            else            {                g[i].x=0;                scanf("%d",&g[i].t);            }            h[i]=g[i].t;        }        sort(h+1,h+n+1);        for(i=1;i<=n;i++)        {            t=lower_bound(h+1,h+n+1,g[i].t)-h;            g[i].t=t;            ans[t]=g[i].y;        }        printf("Case #%d:\n",cas++);        ans[0]=-1;        init(0,n+1,1);        for(i=1;i<=n;i++)        {            if(g[i].x) add(g[i].t,g[i].x,1);            else            {                sum=0;                printf("%d\n",ans[query(0,g[i].t,1)]);            }        }    }    return 0;}