(2016北京集训十二)【xsy1542】疯狂求导-LMLPHP

题解:

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

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

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

一看这种题目形式明显就是大数据结构,外面的序列明显线段树维护,次数也可以用线段树,但是线段树套线段树容易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 ;
}
05-20 00:46