(2016北京集训十二)【xsy1542】疯狂求导
生活随笔
收集整理的這篇文章主要介紹了
(2016北京集训十二)【xsy1542】疯狂求导
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題解:
這題看起來很難。。。但是實際上并沒有想象中的那么難
第一眼看上去不會求導公式怎么辦?不要緊,題目背景非常良心的給出了題目中的導數(shù)計算公式
求完導合并同類項很惡心怎么辦?不要緊,樣例解釋說明了不需要合并同類項(然后有許多人因為這個爆〇了)
一看這種題目形式明顯就是大數(shù)據(jù)結構,外面的序列明顯線段樹維護,次數(shù)也可以用線段樹,但是線段樹套線段樹容易MLE;
所以用樹狀數(shù)組套線段樹實現(xiàn)
具體就是以1~n為下標建線段樹,外面用樹狀數(shù)組維護次數(shù),每次在樹狀數(shù)組上查詢即可
寫完過樣例直接1A就是爽
?代碼:
1 #include<algorithm> 2 #include<iostream> 3 #include<cstring> 4 #include<cstdio> 5 #include<cmath> 6 #include<queue> 7 #define inf 2147483647 8 #define eps 1e-9 9 #define lb(x) (x&-x) 10 using namespace std; 11 typedef long long ll; 12 const int N=1000000,bit=524288; 13 struct node{ 14 int num,ls,rs; 15 ll v; 16 }t[10000001]; 17 int n,m,ans=0,op,l,r,tot=0,rts[1000001]; 18 int TheseNumbersDoNotMakeAnySense; 19 ll s; 20 void newn(int &u){ 21 if(!u)u=++tot; 22 } 23 void pd(int u,int l,int r){ 24 if(t[u].num){ 25 int mid=(l+r)/2; 26 newn(t[u].ls); 27 newn(t[u].rs); 28 t[t[u].ls].num+=t[u].num; 29 t[t[u].ls].v+=t[u].num*(mid-l+1); 30 t[t[u].rs].num+=t[u].num; 31 t[t[u].rs].v+=t[u].num*(r-mid); 32 t[u].num=0; 33 } 34 } 35 void updata(int l,int r,int &u,int L,int R,int v){ 36 newn(u); 37 if(L<=l&&r<=R){ 38 t[u].num+=v; 39 t[u].v+=v*(r-l+1); 40 return; 41 } 42 int mid=(l+r)/2; 43 pd(u,l,r); 44 if(L<=mid)updata(l,mid,t[u].ls,L,R,v); 45 if(mid<R)updata(mid+1,r,t[u].rs,L,R,v); 46 t[u].v=t[t[u].ls].v+t[t[u].rs].v; 47 } 48 int query(int l,int r,int u,int L,int R){ 49 if(!u)return 0; 50 if(L<=l&&r<=R){ 51 return t[u].v; 52 } 53 int mid=(l+r)/2,ret=0; 54 pd(u,l,r); 55 if(L<=mid)ret=query(l,mid,t[u].ls,L,R); 56 if(mid<R)ret+=query(mid+1,r,t[u].rs,L,R); 57 return ret; 58 } 59 void ins(int l,int r,int s){ 60 for(;s<=N;s+=lb(s)){ 61 updata(1,n,rts[s],l,r,1); 62 } 63 } 64 ll ask(int l,int r,int s){ 65 ll ret=0; 66 for(;s;s-=lb(s)){ 67 ret+=query(1,n,rts[s],l,r); 68 } 69 return ret; 70 } 71 void work(int l,int r,ll s){ 72 ll ret=ask(l,r,N),tmp=0; 73 if(ret<=s){ 74 printf("1 %lld\n",ret); 75 ans=1; 76 return; 77 } 78 s=ret-s; 79 for(int bt=bit;bt;bt/=2){ 80 if(tmp+bt<=N){ 81 int q=query(1,n,rts[tmp+bt],l,r); 82 //printf("%lld %lld %lld\n",tmp,bt,q); 83 if(s>q)s-=q,tmp+=bt; 84 } 85 } 86 ll rt=ask(l,r,tmp+1); 87 rt=ret-rt; 88 printf("%lld %lld\n",tmp+2,rt); 89 ans=tmp+2; 90 } 91 int main(){ 92 scanf("%d%d",&n,&m); 93 for(int i=1;i<=n;i++){ 94 scanf("%d",&TheseNumbersDoNotMakeAnySense); 95 } 96 for(int i=1;i<=m;i++){ 97 scanf("%d%d%d%lld",&op,&l,&r,&s); 98 if(op==0){ 99 s^=ans; 100 ins(l,r,s); 101 }else{ 102 work(l,r,s); 103 } 104 } 105 return 0; 106 }?
轉(zhuǎn)載于:https://www.cnblogs.com/dcdcbigbig/p/9696111.html
總結
以上是生活随笔為你收集整理的(2016北京集训十二)【xsy1542】疯狂求导的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: Vuex 环境配置
- 下一篇: express添加权限拦截