Codeforces 803 G. Periodic RMQ Problem

题目链接：http://codeforces.com/problemset/problem/803/G

---

大致就是线段树动态开节点。

然后考虑到如果一个点还没有出现过，那么这个点显然未被修改，就将这个点所代表的区间定位到原序列中，利用ST表查一下区间最小值就可以了。

定位：

 llg minn(llg l,llg r)
 {
     if(r-l+>=n) return mt;
     l%=n;if(!l) l=n;
     r%=n;if(!r) r=n;
     if(l>r) return min(gw(l,n),gw(,r));
     else return gw(l,r);
 }

其中${gw(l,r)}$表示原序列${[l,r]}$中元素的最小值。

 #include<iostream>
 #include<algorithm>
 #include<cstdio>
 #include<cstdlib>
 #include<cmath>
 #include<cstring>
 #include<queue>
 #include<vector>
 #include<map>
 using namespace std;
 #define llg int
 #define inf 0x7fffffff
 #define maxn 500010
 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
 llg ST[maxn][],mt=inf,n,k,Q,pre[maxn],cnt;

 struct point
 {
     llg lc,rc,l,r,set,val;
 }po[maxn*];

 inline llg getint()
 {
     llg w=,q=; char c=getchar();
     while((c<'' || c>'') && c!='-') c=getchar();
     if (c=='-')  q=, c=getchar(); while (c>='' && c<='') w=w*+c-'', c=getchar();
     return q ? -w : w;
 }

 llg gw(llg l,llg r)
 {
     int k=pre[r-l+];
     return min(ST[l][k],ST[r-(<<k)+][k]);
 }

 llg minn(llg l,llg r)
 {
     if(r-l+>=n) return mt;
     l%=n;if(!l) l=n;
     r%=n;if(!r) r=n;
     if(l>r) return min(gw(l,n),gw(,r));
     else return gw(l,r);
 }

 void make_st()
 {
     for(llg i=;i<=n;i++) pre[i]=pre[i>>]+;
     for(llg i=;i<=n;i++) ST[i][]=getint();
     for(llg i=;i<=n;i++) mt=min(mt,ST[i][]);
     for(llg i=;i<=n;i++) pre[i]=pre[i>>]+;
     for(llg j=;j<=;j++)
         for(llg i=,u=(<<j),p=(u>>);i+u-<=n;i++)
             ST[i][j]=min(ST[i][j-],ST[i+p][j-]);
     cnt=,po[].l=,po[].r=n*k,po[].val=mt;
 }

 llg new_po(llg l,llg r)
 {
     llg o=++cnt;
     po[o].val=minn(l,r);
     po[o].l=l,po[o].r=r;
     return o;
 }

 void update(llg o)
 {
     if (po[o].l==po[o].r) return ;
     po[o].val=inf;
     if (po[o].lc) po[o].val=min(po[o].val,po[po[o].lc].val);
     else
     {
         po[o].lc=cnt+;
         new_po(po[o].l,(po[o].l+po[o].r)/);
         po[o].val=min(po[o].val,po[po[o].lc].val);
     }
     if (po[o].rc) po[o].val=min(po[o].val,po[po[o].rc].val);
     else
     {
         po[o].rc=cnt+;
         new_po((po[o].l+po[o].r)/+,po[o].r);
         po[o].val=min(po[o].val,po[po[o].rc].val);
     }
 }

 void pushdown(llg o)
 {
     llg l=po[o].l,r=po[o].r,lc=po[o].lc,rc=po[o].rc,mid=(l+r)>>;
     if (!po[o].set || l==r) return ;
     if (!lc)
     {
         lc=cnt+;
         new_po(l,mid);
     }
     if (!rc)
     {
         rc=cnt+;
         new_po(mid+,r);
     }
     po[lc].set=po[rc].set=po[lc].val=po[rc].val=po[o].set;
     po[o].set=;
     po[o].lc=lc;
     po[o].rc=rc;
     return ;
 }

 void modify(llg o,llg l,llg r,llg L,llg R,llg v)
 {
     pushdown(o);
     if (o==) o=new_po(l,r);
     llg mid=(l+r)>>;
     if (l>=L && r<=R)
     {
         po[o].val=v;
         po[o].set=v;
         return ;
     }
     if (mid>=L)
     {
         llg son=po[o].lc;
         if (po[o].lc==) po[o].lc=cnt+;
         modify(son,l,mid,L,R,v);
     }
     if (mid<R)
     {
         llg son=po[o].rc;
         if (po[o].rc==) po[o].rc=cnt+;
         modify(son,mid+,r,L,R,v);
     }
     update(o);
 }

 llg query(llg o,llg l,llg r,llg L,llg R)
 {
     pushdown(o);
     if (o==) o=new_po(l,r);
     if (l>=L && r<=R) return po[o].val;
     llg mid=(l+r)>>,ans=inf;
     if (mid>=L)
     {
         llg son=po[o].lc;
         if (po[o].lc==) po[o].lc=cnt+;
         ans=min(ans,query(son,l,mid,L,R));
     }
     if (mid<R)
     {
         llg son=po[o].rc;
         if (po[o].rc==) po[o].rc=cnt+;
         ans=min(ans,query(son,mid+,r,L,R));
     }
     return ans;
 }

 int main()
 {
     yyj("seg");
     cin>>n>>k;
     make_st();
     llg T=n*k;
     cin>>Q;
     while (Q--)
     {
         llg ty=getint(),l=getint(),r=getint();
         if (ty==)
             modify(,,T,l,r,getint());
         else
             printf("%d\n",query(,,T,l,r));
     }
     return ;
 }