捣鼓了好长时间才终于把zkw线段树的区间最值修改给A掉了
= =zkw神犇的课件里竟然有错!
(╯▽╰)下面的伪代码里被注释的是我修改的><
Func Add_x(s,t,x)for (s=s+M-1,t=t+M+1;s^t^1;s>>=1,t>>=1) {if (~s&1) T[s^1]+=x;if ( t&1) T[t^1]+=x;A=min(T[s],T[s^1]),T[s]-=A,T[s^1]-=A, T[s>>1]+=A;A=min(T[t],T[t^1]),T[t]-=A,T[t^1]-=A,T[t>>1]+=A;}while (s!=1){A=min(T[s],num[s^1]);T[s]-=A;T[s^1]-=A;T[s>>=1]+=A;}Func Max(s,t)/*for (s=s+M-1,t=t+M+1;s^t^1;s>>=1,t>>=1) {LAns+=T[s],Rans+=T[t];if (~s&1) LAns=max(LAns,T[s^1]);if ( t&1) RAns=max(RAns,T[t^1]);}*/if(s!=t)for (s=s+M,t=t+M;s^t^1;s>>=1,t>>=1) {LAns+=T[s],Rans+=T[t];if (~s&1) LAns=max(LAns,T[s^1]);if ( t&1) RAns=max(RAns,T[t^1]);}*/Ans = max(LAns/*+T[s]*/,RAns/*+T[t]*/);while (s>1) Ans += T[s>>=1];
第一个是在Add_x函数里,zkw神犇没有进行上传操作,这样单次修改没有任何问题,多次操作就会出现错误。
第二个是在Max函数里,zkw神犇没有上传到根节点,那么就会在上浮一层后出现混乱。
按照我的代码A掉了NOIP2012年提高组的借教室,感觉挺高兴= =
Description
在大学期间,经常需要租借教室。大到院系举办活动,小到学习小组自习讨论,都需要向学校申请借教室。教室的大小功能不同,借教室人的身份不同,借教室的手续也不一样。面对海量租借教室的信息,我们自然希望编程解决这个问题。
我们需要处理接下来n天的借教室信息,其中第i天学校有ri个教室可供租借。共有m份订单,每份订单用三个正整数描述,分别为dj,sj,tj,表示某租借者需要从第sj天到第tj天租借教室(包括第sj天和第tj天),每天需要租借dj个教室。
我们假定,租借者对教室的大小、地点没有要求。即对于每份订单,我们只需要每天提供dj个教室,而它们具体是哪些教室,,每天是否是相同的教室则不用考虑。
借教室的原则是先到先得,也就是说我们要按照订单的先后顺序依次为每份订单分配教室。如果在分配的过程中遇到一份订单无法完全满足,则需要停止教室的分配,通知当前申请人修改订单。这里的无法满足指从第sj天到第tj天中有至少一天剩余的教室数量不足dj个。现在我们需要知道,是否会有订单无法完全满足。如果有,需要通知哪一个申请人修改订单。
Input Description
第一行包含两个正整数n,m,表示天数和订单的数量。
第二行包含n个正整数,其中第i个数为ri,表示第i天可用于租借的教室数量。
接下来有m行,每行包含三个正整数dj,sj,tj,表示租借的数量,租借开始、结束分别在第几天。每行相邻的两个数之间均用一个空格隔开。天数与订单均用从1开始的整数编号。
Output Description
如果所有订单均可满足,则输出只有一行,包含一个整数0。否则(订单无法完全满足)输出两行,第一行输出一个负整数-1,第二行输出需要修改订单的申请人编号。
Sample Input
43
2543
213
324
424
Sample Output
-1
2
Data Size & Hint
对于10%的数据,有1≤n,m≤10;
对于30%的数据,有1≤n,m≤1000;
对于70%的数据,有1≤n,m≤105;
对于100%的数据,有1≤n,m≤10^6,0≤ri,dj≤10^9,1≤sj≤tj≤n。
/*Name:wxjlzbcdMemory:17MBTime:1574ms*/#include<cstdio>#include<cstdlib>#include<cctype>#include<cstring>#include<algorithm>using namespace std;int x,a,b,c,n,m,pre=1;int num[1048576];inline int read(){int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}inline void add(int a,int b,int c){int tmp=0;a+=pre-1;b+=pre+1;while(a^b^1){if(~a&1)num[a^1]+=c;if(b&1)num[b^1]+=c;tmp=min(num[a],num[a^1]);num[a]-=tmp;num[a^1]-=tmp;num[a>>1]+=tmp;tmp=min(num[b],num[b^1]);num[b]-=tmp;num[b^1]-=tmp;num[b>>1]+=tmp;a>>=1;b>>=1;}while(a!=1){tmp=min(num[a],num[a^1]);num[a]-=tmp;num[a^1]-=tmp;num[a>>=1]+=tmp;}}inline int query(int a,int b){int lans=0,rans=0,ans=0;a+=pre;b+=pre;if(a!=b)while(a^b^1){lans+=num[a];rans+=num[b];if(~a&1)lans=min(lans,num[a^1]);if(b&1)rans=min(rans,num[b^1]);a>>=1;b>>=1;}ans=min(lans+num[a],rans+num[b]);while(a>1)ans+=num[a>>=1];return ans;}int main(){int i;n=read(),m=read();while(pre<=n)pre<<=1;for(i=pre+1;i<=pre+n;++i)num[i]=read();for(i=pre-1;i>=1;–i){num[i]=min(num[i<<1],num[(i<<1)+1]);num[i<<1]-=num[i];num[(i<<1)+1]-=num[i];}for(i=1;i<=m;++i){c=read(),a=read(),b=read();if(query(a,b)>=c)add(a,b,-c);else{printf("-1\n%d",i);return 0;}}putchar('0');return 0;}
君子无故,玉不去身。