bzoj4137 [FJOI2015]火星商店问题

  1 #include<cstdio>
  2 #include<set>
  3 #include<vector>
  4 #define N 1000010
  5 #define M 10000010
  6 using namespace std;
  7 int l[N],r[N],n,m,i,a,typ,b,L,R,ti,ans[N],tot,p[N],cnt,flag;
  8 int s[M][2],t[M];
  9 vector<pair<int,int> > vec[N];
 10 vector<int> id[N];
 11 void build(int x,int a,int b)
 12 {
 13     int m;
 14     l[x]=a;r[x]=b;
 15     if (b-a>1)
 16     {
 17         m=(a+b)>>1;
 18         build(2*x,a,m);
 19         build(2*x+1,m,b);
 20     }
 21 }
 22 void insert(int x,int a,int b,pair<int,int> c,int typ)
 23 {
 24     int m;
 25     if (typ==0)
 26     {
 27         vec[x].push_back(c);
 28         id[x].push_back(typ);
 29     }
 30     if ((a<=l[x])&&(r[x]<=b))
 31     {
 32         if (typ)
 33         {
 34             vec[x].push_back(c);
 35             id[x].push_back(typ);
 36         }
 37         return; 
 38     }
 39     m=(l[x]+r[x])>>1;
 40     if (a<m) insert(2*x,a,b,c,typ);
 41     if (m<b) insert(2*x+1,a,b,c,typ);
 42 } 
 43 void cl(int x)
 44 {
 45     int i;
 46     for (i=1;i<=17;i++)
 47     {
 48         p[17-i+1]=x%2;
 49         x=x/2;
 50     }
 51 }
 52 void change(int x,int y)
 53 {
 54     int now,i;
 55     cl(x);
 56     now=1;
 57     for (i=1;i<=17;i++)
 58     {
 59         if (s[now][p[i]]==0) 
 60             s[now][p[i]]=++tot;
 61         now=s[now][p[i]];
 62         t[now]=max(t[now],y);
 63     }
 64 }
 65 int query(int x,int y)
 66 {
 67     int now,i,ans=0;
 68     cl(x);
 69     now=1;
 70     for (i=1;i<=17;i++)
 71     if ((s[now][1-p[i]])&&(t[s[now][1-p[i]]]>=y))
 72     {
 73         now=s[now][1-p[i]];
 74         ans=ans+(1<<(17-i));
 75     }
 76     else
 77     now=s[now][p[i]];
 78     return ans;
 79 }
 80 void gao(int x)
 81 {
 82     int i,a,b;
 83     for (i=1;i<=tot;i++)
 84     s[i][0]=s[i][1]=t[i]=0;
 85     tot=1;
 86     for (i=0;i<vec[x].size();i++)
 87     {
 88         a=vec[x][i].first;
 89         b=vec[x][i].second;
 90         if (id[x][i]==0)
 91         change(a,b);
 92         else
 93         ans[id[x][i]]=max(ans[id[x][i]],query(a,b));
 94     }
 95     if (r[x]-l[x]==1)
 96     return;
 97     gao(2*x);
 98     gao(2*x+1);
 99 } 
100 int main()
101 {
102     scanf("%d%d",&n,&m);
103     build(1,0,n);
104     for (i=1;i<=n;i++)
105     {
106         scanf("%d",&a); 
107         insert(1,i-1,i,make_pair(a,m+1),0);
108     }
109     ti=1;
110     for (i=1;i<=m;i++)
111     {
112         scanf("%d",&typ);
113         if (typ==1)
114         {
115             scanf("%d%d%d%d",&L,&R,&a,&b);
116             insert(1,L-1,R,make_pair(a,ti-b+1),++cnt);
117         }
118         else
119         {
120             ti++;
121             scanf("%d%d",&a,&b);
122             insert(1,a-1,a,make_pair(b,ti),0);
123         }
124     }
125     gao(1);
126     for (i=1;i<=cnt;i++)
127     printf("%d\n",ans[i]);
128 }