【分块】【常数优化】【Orz faebdc】洛谷 P1083 NOIP2012提高组 借教室

分块90分。 By AutSky_JadeK 【重点在下面】

 #include<cstdio>
 #include<cmath>
 using namespace std;
 #define N 1000001
 #define INF 2147483647
 #define min(a,b) (((a)<(b))?(a):(b))
 int n,m,a[N],sum,sz,l[],r[],minv[],delta[],num[N],x,y,v;
 int ans,re2;
 int Res;char C;
 inline int G()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=Res*+(C-'');C=getchar();}
     return Res;
 }
 void makeblock()
 {
     sz=(int)sqrt((double)n*0.1); if(!sz) sz=;
     for(sum=;sum*sz<n;sum++)
       {
         l[sum]=r[sum-]+; r[sum]=sum*sz;
         minv[sum]=INF;
         for(int i=l[sum];i<=r[sum];i++) {minv[sum]=min(minv[sum],a[i]); num[i]=sum;}
       }
     l[sum]=r[sum-]+; r[sum]=n;
     minv[sum]=INF;
     for(int i=l[sum];i<=r[sum];i++) {minv[sum]=min(minv[sum],a[i]); num[i]=sum;}
 }
 inline void work(const int &L,const int &R)
 {
     minv[num[L]]=INF;
     for(int i=L;i<=R;i++) a[i]-=v;
     for(int i=l[num[L]];i<=r[num[L]];i++) minv[num[L]]=min(minv[num[L]],a[i]);
 }
 inline void update()
 {
     if(num[x]==num[y]) work(x,y);
     else
       {
         work(x,r[num[x]]); work(l[num[y]],y);
         for(int i=num[x]+;i<num[y];i++) delta[i]-=v;
       }
 }
 inline int query()
 {
     ans=INF;
     if(num[x]==num[y]) {for(int i=x;i<=y;i++) ans=min(ans,a[i]+delta[num[x]]); }
     else
       {
         for(int i=x;i<=r[num[x]];i++) ans=min(ans,a[i]+delta[num[x]]);
         for(int i=num[x]+;i<num[y];i++) ans=min(ans,minv[i]+delta[i]);
         for(int i=l[num[y]];i<=y;i++) ans=min(ans,a[i]+delta[num[y]]);
       } return ans;
 }
 int main()
 {
     n=G(); m=G();
     for(int i=;i<=n;i++) a[i]=G();
     makeblock();
     for(int i=;i<=m;i++)
       {
         v=G(); x=G(); y=G();
         if(query()<v)
           {
             printf("-1\n%d\n",i);
             return ;
           }
         else update();
       }
     puts("");
     return ;
 }

分块 常数优化 100分。 By faebdc 金牌爷の力Orzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！

 #include<cstdio>
 #include<cmath>
 using namespace std;
 #define N 1000001
 #define INF 2147483647
 int n,m,a[N],sum,sz,hf,l[],r[],minv[],delta[],num[N],x,y,v;
 int ans,re2;
 int Res;char C;
 inline int G()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=(Res<<)+(Res<<)+C-'';C=getchar();}
     return Res;
 }
 void makeblock()
 {
     sz=(int)sqrt((double)n*0.125); if(!sz) sz=;
     hf=sz>>;
     for(sum=;sum*sz<n;sum++)
       {
         l[sum]=r[sum-]+; r[sum]=sum*sz;
         minv[sum]=INF;
         for(int i=l[sum];i<=r[sum];i++) {if(minv[sum]>a[i])minv[sum]=a[i]; num[i]=sum;}
       }
     l[sum]=r[sum-]+; r[sum]=n;
     minv[sum]=INF;
     for(int i=l[sum];i<=r[sum];i++) {if(minv[sum]>a[i])minv[sum]=a[i]; num[i]=sum;}
 }
 inline void work(const int &L,const int &R)
 {
     int *t=&(minv[num[L]]);
     *t=INF;
     int *z;
     for(int i=L,*z=a+i;i<=R;i++,z++) {*z-=v;if(*t>*z)*t=*z;}
     for(int i=r[num[L]],*z=a+i;i>R;i--,z--) if(*t>*z)*t=*z;
     for(int i=l[num[L]],*z=a+i;i<L;i++,z++) if(*t>*z)*t=*z;
 }
 inline int update()
 {
     if(num[x]==num[y]) {work(x,y);return minv[num[x]]+delta[num[x]];}
         work(x,r[num[x]]); work(l[num[y]],y);
         int z;
         if(minv[num[x]]+delta[num[x]]>minv[num[y]]+delta[num[y]])
             z=minv[num[y]]+delta[num[y]];
         else
             z=minv[num[x]]+delta[num[x]];
         for(int i=num[x]+;i<num[y];i++) {delta[i]-=v;if(z>minv[i]+delta[i])z=minv[i]+delta[i];}
       return z;
 }
 int main()
 {
     n=G(); m=G();
     for(int i=;i<=n;i++) a[i]=G();
     makeblock();
     for(int i=;i<=m;i++)
       {
         v=G(); x=G(); y=G();
         int z=update();
         if(z<)
           {
             printf("-1\n%d\n",i);
             return ;
           }
       }
     puts("");
     return ;
 }