（2016北京集训十二）【xsy1542】疯狂求导

题解：

这题看起来很难。。。但是实际上并没有想象中的那么难

第一眼看上去不会求导公式怎么办？不要紧，题目背景非常良心的给出了题目中的导数计算公式

求完导合并同类项很恶心怎么办？不要紧，样例解释说明了不需要合并同类项（然后有许多人因为这个爆〇了）

一看这种题目形式明显就是大数据结构，外面的序列明显线段树维护，次数也可以用线段树，但是线段树套线段树容易MLE；

所以用树状数组套线段树实现

具体就是以1~n为下标建线段树，外面用树状数组维护次数，每次在树状数组上查询即可

写完过样例直接1A就是爽

代码：

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 #include<queue>
 #define inf 2147483647
 #define eps 1e-9
 #define lb(x) (x&-x)
 using namespace std;
 typedef long long ll;
 const int N=,bit=;
 struct node{
     int num,ls,rs;
     ll v;
 }t[];
 int n,m,ans=,op,l,r,tot=,rts[];
 int TheseNumbersDoNotMakeAnySense;
 ll s;
 void newn(int &u){
     if(!u)u=++tot;
 }
 void pd(int u,int l,int r){
     if(t[u].num){
         int mid=(l+r)/;
         newn(t[u].ls);
         newn(t[u].rs);
         t[t[u].ls].num+=t[u].num;
         t[t[u].ls].v+=t[u].num*(mid-l+);
         t[t[u].rs].num+=t[u].num;
         t[t[u].rs].v+=t[u].num*(r-mid);
         t[u].num=;
     }
 }
 void updata(int l,int r,int &u,int L,int R,int v){
     newn(u);
     if(L<=l&&r<=R){
         t[u].num+=v;
         t[u].v+=v*(r-l+);
         return;
     }
     int mid=(l+r)/;
     pd(u,l,r);
     if(L<=mid)updata(l,mid,t[u].ls,L,R,v);
     if(mid<R)updata(mid+,r,t[u].rs,L,R,v);
     t[u].v=t[t[u].ls].v+t[t[u].rs].v;
 }
 int query(int l,int r,int u,int L,int R){
     if(!u)return ;
     if(L<=l&&r<=R){
         return t[u].v;
     }
     int mid=(l+r)/,ret=;
     pd(u,l,r);
     if(L<=mid)ret=query(l,mid,t[u].ls,L,R);
     if(mid<R)ret+=query(mid+,r,t[u].rs,L,R);
     return ret;
 }
 void ins(int l,int r,int s){
     for(;s<=N;s+=lb(s)){
         updata(,n,rts[s],l,r,);
     }
 }
 ll ask(int l,int r,int s){
     ll ret=;
     for(;s;s-=lb(s)){
         ret+=query(,n,rts[s],l,r);
     }
     return ret;
 }
 void work(int l,int r,ll s){
     ll ret=ask(l,r,N),tmp=;
     if(ret<=s){
         printf("1 %lld\n",ret);
         ans=;
         return;
     }
     s=ret-s;
     for(int bt=bit;bt;bt/=){
         if(tmp+bt<=N){
             int q=query(,n,rts[tmp+bt],l,r);
             //printf("%lld %lld %lld\n",tmp,bt,q);
             if(s>q)s-=q,tmp+=bt;
         }
     }
     ll rt=ask(l,r,tmp+);
     rt=ret-rt;
     printf("%lld %lld\n",tmp+,rt);
     ans=tmp+;
 }
 int main(){
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++){
         scanf("%d",&TheseNumbersDoNotMakeAnySense);
     }
     for(int i=;i<=m;i++){
         scanf("%d%d%d%lld",&op,&l,&r,&s);
         if(op==){
             s^=ans;
             ins(l,r,s);
         }else{
             work(l,r,s);
         }
     }
     return ;
 }