MyException - 我的异常网
当前位置:我的异常网» C++ » 洛谷 P2023 [AHOI2009]维护序列 || 线段树加法跟乘

洛谷 P2023 [AHOI2009]维护序列 || 线段树加法跟乘法运算

www.MyException.Cn  网友分享于:2013-11-16  浏览:0次
洛谷 P2023 [AHOI2009]维护序列 || 线段树加法和乘法运算
"); $("#cnblogs_post_body").prepend("
"); $("#cnblogs_post_body").prepend("
如果行号影响了复制,请点击代码框左上角的按钮。
"); //$("#deng4").after("
"); $("#backTop").click(function () { var speed=500;//滑动的速度 $('body,html').animate({ scrollTop: 0 }, speed); return false; });}); //$("body").append("

洛谷 P2023 [AHOI2009]维护序列 || 线段树加法和乘法运算

原理倒是非常简单。设原数为x,加法的lazytag为b,乘法的lazytag为a,操作数为c,那么原式为ax+b,乘上c后(ax+b)c=(ac)*x+b*c,加上c后(ax+b)+c=ax+(b+c),因此加法时只需要更新加法的lazytag,乘法的时候就需要同时乘乘法和加法的lazytag。(乘法时的操作曾经搞错)

有些细节却很难受啊...比如很容易忘记加取模(一般都是程序打完了之后再补上,但是这个程序要补的地方太多,容易漏,要小心),把乘法的lazytag初始化为1出错,很容易少几句话...

  1 //以下代码中:加important标记的是曾经缺少的,注释掉的是曾经多的
  2 #include<cstdio>
  3 #include<cstring>
  4 #define lc (num<<1)
  5 #define rc (num<<1|1)
  6 #define mid ((l+r)>>1)
  7 typedef long long LL;
  8 LL a[300100];
  9 LL tree[1200100],laz[1200100],laz2[1200100];
 10 LL L,R,x,n,m,md;
 11 void build(LL l,LL r,LL num)
 12 {
 13     laz2[num]=1;//important
 14     if(l==r)
 15     {
 16         tree[num]=a[l]%md;
 17         return;
 18     }
 19     build(l,mid,lc);
 20     build(mid+1,r,rc);
 21     tree[num]=(tree[lc]+tree[rc])%md;
 22 }
 23 void pushdown(LL l,LL r,LL num)
 24 {
 25     if(laz2[num]!=1)
 26     {
 27         laz2[lc]=(laz2[lc]*laz2[num])%md;
 28         laz2[rc]=(laz2[rc]*laz2[num])%md;
 29         /*important*/
 30         laz[lc]=(laz[lc]*laz2[num])%md;
 31         laz[rc]=(laz[rc]*laz2[num])%md;
 32         tree[lc]=(tree[lc]*laz2[num])%md;
 33         tree[rc]=(tree[rc]*laz2[num])%md;
 34         /*-important*/
 35         /*
 36         tree[lc]=(tree[lc]+(mid-l+1)*laz2[num])%md;
 37         tree[rc]=(tree[rc]+(r-mid)*laz2[num])%md;
 38         */
 39         laz2[num]=1;
 40     }
 41     if(laz[num])
 42     {
 43         laz[lc]=(laz[lc]+laz[num])%md;
 44         laz[rc]=(laz[rc]+laz[num])%md;
 45         tree[lc]=(tree[lc]+(mid-l+1)*laz[num]%md)%md;
 46         tree[rc]=(tree[rc]+(r-mid)*laz[num]%md)%md;
 47         laz[num]=0;
 48     }
 49 }
 50 void update(LL l,LL r,LL num)
 51 {
 52     if(L<=l&&r<=R)
 53     {
 54         tree[num]=(tree[num]+(r-l+1)*x)%md;
 55         laz[num]=(laz[num]+x)%md;
 56         return;
 57     }
 58     pushdown(l,r,num);
 59     if(L<=mid)    update(l,mid,lc);
 60     if(mid<R)    update(mid+1,r,rc);//if(mid+1<=R)//等价
 61     /*important*/tree[num]=(tree[lc]+tree[rc])%md;
 62 }
 63 void update2(LL l,LL r,LL num)
 64 {
 65     if(L<=l&&r<=R)
 66     {
 67         tree[num]=(tree[num]*x)%md;
 68         /*importaant*/laz[num]=(laz[num]*x)%md;/*important*/
 69         laz2[num]=(laz2[num]*x)%md;
 70         return;
 71     }
 72     pushdown(l,r,num);
 73     if(L<=mid)    update2(l,mid,lc);
 74     if(mid<R)    update2(mid+1,r,rc);//if(mid+1<=R)//等价
 75     /*important*/tree[num]=(tree[lc]+tree[rc])%md;
 76 }
 77 LL query(LL l,LL r,LL num)
 78 {
 79     if(L<=l&&r<=R)    return tree[num];
 80     pushdown(l,r,num);
 81     LL ans=0;
 82     if(L<=mid)    ans=(ans+query(l,mid,lc))%md;
 83     if(mid<R)    ans=(ans+query(mid+1,r,rc))%md;
 84     return ans;
 85 }
 86 int main()
 87 {
 88     LL i,t;
 89     scanf("%lld%lld%lld",&n,&m,&md);
 90     for(i=1;i<=n;i++)
 91         scanf("%lld",&a[i]);//不应该将laz2[]=1写在此处important
 92     build(1,n,1);
 93     for(i=1;i<=m;i++)
 94     {
 95         scanf("%lld",&t);
 96         if(t==2)
 97         {
 98             scanf("%lld%lld%lld",&L,&R,&x);
 99             update(1,n,1);
100         }
101         else if(t==3)
102         {
103             scanf("%lld%lld",&L,&R);
104             printf("%lld\n",query(1,n,1)%md);
105         }
106         else
107         {
108             scanf("%lld%lld%lld",&L,&R,&x);
109             update2(1,n,1);
110         }
111     }
112     return 0;
113 }

文章评论

Google伦敦新总部 犹如星级庄园
Google伦敦新总部 犹如星级庄园
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
那些性感的让人尖叫的程序员
那些性感的让人尖叫的程序员
如何成为一名黑客
如何成为一名黑客
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
Web开发者需具备的8个好习惯
Web开发者需具备的8个好习惯
Java程序员必看电影
Java程序员必看电影
程序员都该阅读的书
程序员都该阅读的书
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
漫画:程序员的工作
漫画:程序员的工作
旅行,写作,编程
旅行,写作,编程
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
程序员必看的十大电影
程序员必看的十大电影
老程序员的下场
老程序员的下场
 程序员的样子
程序员的样子
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
我是如何打败拖延症的
我是如何打败拖延症的
初级 vs 高级开发者 哪个性价比更高?
初级 vs 高级开发者 哪个性价比更高?
要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
总结2014中国互联网十大段子
总结2014中国互联网十大段子
亲爱的项目经理,我恨你
亲爱的项目经理,我恨你
鲜为人知的编程真相
鲜为人知的编程真相
2013年中国软件开发者薪资调查报告
2013年中国软件开发者薪资调查报告
为啥Android手机总会越用越慢?
为啥Android手机总会越用越慢?
Web开发人员为什么越来越懒了?
Web开发人员为什么越来越懒了?
10个调试和排错的小建议
10个调试和排错的小建议
程序员应该关注的一些事儿
程序员应该关注的一些事儿
我的丈夫是个程序员
我的丈夫是个程序员
老美怎么看待阿里赴美上市
老美怎么看待阿里赴美上市
程序员周末都喜欢做什么?
程序员周末都喜欢做什么?
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
程序员的鄙视链
程序员的鄙视链
5款最佳正则表达式编辑调试器
5款最佳正则表达式编辑调试器
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
写给自己也写给你 自己到底该何去何从
写给自己也写给你 自己到底该何去何从
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
60个开发者不容错过的免费资源库
60个开发者不容错过的免费资源库
每天工作4小时的程序员
每天工作4小时的程序员
我跳槽是因为他们的显示器更大
我跳槽是因为他们的显示器更大
不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
编程语言是女人
编程语言是女人
“肮脏的”IT工作排行榜
“肮脏的”IT工作排行榜
为什么程序员都是夜猫子
为什么程序员都是夜猫子
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
程序员和编码员之间的区别
程序员和编码员之间的区别
10个帮程序员减压放松的网站
10个帮程序员减压放松的网站
软件开发程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有