E  DZY Loves Fibonacci Numbers

In mathematical terms, the sequence Fn of Fibonacci numbers is defined by the recurrence relation

F1 = 1; F2 = 1; Fn = Fn - 1 + Fn - 2 (n > 2).

DZY loves Fibonacci numbers very much. Today DZY gives you an array consisting of nintegers: a1, a2, ..., an. Moreover, there are m queries, each query has one of the two types:

  1. Format of the query "1 l r". In reply to the query, you need to add Fi - l + 1 to each element ai, where l ≤ i ≤ r.
  2. Format of the query "2 l r". In reply to the query you should output the value of  modulo 1000000009 (109 + 9).

Help DZY reply to all the queries.

Input

The first line of the input contains two integers n and m (1 ≤ n, m ≤ 300000). The second line contains n integers a1, a2, ..., an (1 ≤ ai ≤ 109) — initial array a.

Then, m lines follow. A single line describes a single query in the format given in the statement. It is guaranteed that for each query inequality 1 ≤ l ≤ r ≤ n holds.

Output

For each query of the second type, print the value of the sum on a single line.

Example

Input
4 41 2 3 41 1 42 1 41 2 42 1 3
Output
1712

Note

After the first query, a = [2, 3, 5, 7].

For the second query, sum = 2 + 3 + 5 + 7 = 17.

After the third query, a = [2, 4, 6, 9].

For the fourth query, sum = 2 + 4 + 6 = 12.

题目的大意就是,定义了fib数列前两项F[1]=1,F[2]=1。给出整数n,m,再给出n个数和m项操作。

每个操作包含三个数op,l,r,若op=1,就在[l,r]的区间相应加上F[1,r-l+1](每个数对应着加),若op=2,就将[l,r]中所有数统计和并输出。

首先,我们要知道,fib数列有个特殊的性质_两个fib数列相加,还是fib数列,只是前两项的值变动了而已(从而导致后面的数变动).

那么,根据这个性质,就可以通过累计的方法实现_用线段树进行维护.

首先,我们需要知道一些性质:

设fib[1]=a,fib[2]=b,则fib[n]=fib[n-1]*b+fib[n-2]*a.(可推)

fib[1]+fib[2]+fib[3]+......+fib[n]=fib[n+2]-fib[2].(可推)

通过这两个性质,我们可以巧妙的记录线段树上每一个节点的前两个(fib[1],fib[2])的值,然后计算出这个节点上[L,R],所有数的和.

怎么来?

我们对于每一个线段树上的节点,都记录两个值,ta,tb,分别表示这个区间前两项的值,在更新时,

 void upt(int now,int fir,int sec,int len){
     T[now].ta=(T[now].ta+fir)%TT,T[now].tb=(T[now].tb+sec)%TT;
     T[now].key=((]-(long long)sec)%TT;
 }

其中前两句是将标记累计,最后以句就是顺带计算出这个节点的值.

同时,在pushdown函数中,也要对此进行更新.

 void push_down(int now){
     if (!T[now].ta&&!T[now].tb) return;
     int L=T[now].L,R=T[now].R;
     upt(now<<,T[now].ta,T[now].tb,mid-L+);
     ,a,b;
     a=((]+(long long)T[now].tb*fib[len])%TT;
     b=((])%TT;
     upt(now<<|,a,b,R-mid);
     T[now].ta=T[now].tb=;
 }

其中,对于某个节点的子节点,有两段.一段[L,mid],一段[mid+1,R].

那么更新[L,mid]时,它的前缀和当前节点是一样的,所以不需改动信息,直接下传.

更新[mid+1,R]时,需要重新计算出这段区间前两项的值(这可以用上面说的性质算),然后才能下传.

最后下传完了记得把标记清空.

不过不知道为什么,我的线段树开4,5倍空间会炸,开10倍才不炸,很玄学......

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define mid ((L+R)>>1)
 using namespace std;
 ,maxn=;
 int n,m,a[maxn],fib[maxn];
 ];
 int read(){
     ; char ch=getchar();
     ') ch=getchar();
     +ch-',ch=getchar();
     return x;
 }
 void make(int now,int L,int R){
     T[now].L=L,T[now].R=R;
     if (L==R){T[now].key=a[L]; return;}
     make(now<<,L,mid),make(now<<|,mid+,R);
     T[now].key=(T[now<<].key+T[now<<|].key)%TT;
 }
 void upt(int now,int fir,int sec,int len){
     T[now].ta=(T[now].ta+fir)%TT,T[now].tb=(T[now].tb+sec)%TT;
     T[now].key=((]-(long long)sec)%TT;
 }
 void push_down(int now){
     if (!T[now].ta&&!T[now].tb) return;
     int L=T[now].L,R=T[now].R;
     upt(now<<,T[now].ta,T[now].tb,mid-L+);
     ,a,b;
     a=((]+(long long)T[now].tb*fib[len])%TT;
     b=((])%TT;
     upt(now<<|,a,b,R-mid);
     T[now].ta=T[now].tb=;
 }
 void alter(int now,int aimL,int aimR){
     int L=T[now].L,R=T[now].R;
     if (L>aimR||R<aimL) return;
     ],fib[L-aimL+],R-L+); push_down(now); return;}
     push_down(now);
     alter(now<<,aimL,aimR),alter(now<<|,aimL,aimR);
     T[now].key=(T[now<<].key+T[now<<|].key)%TT;
 }
 int seek(int now,int aimL,int aimR){
     int L=T[now].L,R=T[now].R;
     ;
     if (L>=aimL&&R<=aimR) return T[now].key;
     push_down(now);
     ,aimL,aimR)+seek(now<<|,aimL,aimR))%TT;
 }
 int main(){
     n=read(),m=read(),memset(T,,sizeof T);
     ; i<=n; i++) a[i]=read();
     make(,,n);
     fib[]=,fib[]=fib[]=;
     ; i<=n+; i++) fib[i]=(fib[i-]+fib[i-])%TT;
     ; i<=m; i++){
         int tp=read(),L=read(),R=read();
         ) alter(,L,R);
         ,L,R));
     }
     ;
 }

[CodeForces - 447E] E - DZY Loves Fibonacci Numbers的更多相关文章

  1. codeforces 446C DZY Loves Fibonacci Numbers(数学 or 数论+线段树)(两种方法)

    In mathematical terms, the sequence Fn of Fibonacci numbers is defined by the recurrence relation F1 ...

  2. Codeforces 446-C DZY Loves Fibonacci Numbers 同余 线段树 斐波那契数列

    C. DZY Loves Fibonacci Numbers time limit per test 4 seconds memory limit per test 256 megabytes inp ...

  3. Codeforces Round #FF 446 C. DZY Loves Fibonacci Numbers

    參考:http://www.cnblogs.com/chanme/p/3843859.html 然后我看到在别人的AC的方法里还有这么一种神方法,他预先设定了一个阈值K,当当前的更新操作数j<K ...

  4. ACM学习历程—Codeforces 446C DZY Loves Fibonacci Numbers(线段树 && 数论)

    Description In mathematical terms, the sequence Fn of Fibonacci numbers is defined by the recurrence ...

  5. codeforces 446C DZY Loves Fibonacci Numbers 数论+线段树成段更新

    DZY Loves Fibonacci Numbers Time Limit:4000MS     Memory Limit:262144KB     64bit IO Format:%I64d &a ...

  6. Codeforces 446C —— DZY Loves Fibonacci Numbers(线段树)

    题目:DZY Loves Fibonacci Numbers 题意比較简单,不解释了. 尽管官方的题解也是用线段树,但还利用了二次剩余. 可是我没有想到二次剩余,然后写了个感觉非常复杂度的线段树,还是 ...

  7. 『题解』Codeforces446C DZY Loves Fibonacci Numbers

    更好的阅读体验 Portal Portal1: Codeforces Portal2: Luogu Description In mathematical terms, the sequence \( ...

  8. cf446C DZY Loves Fibonacci Numbers

    C. DZY Loves Fibonacci Numbers time limit per test 4 seconds memory limit per test 256 megabytes inp ...

  9. Codeforces446C - DZY Loves Fibonacci Numbers

    Portal Description 给出一个\(n(n\leq3\times10^5)\)个数的序列,进行\(m(m\leq3\times10^5)\)次操作,操作有两种: 给区间\([L,R]\) ...

随机推荐

  1. .NET下使用 Seq结构化日志系统

    前言 我们公司在日志管理方面一直没有统一,主要痛点有: 每个开发人员都是各用各的,存储日志的形式也是五花八门,如:本地文件,数据库,Redis,MongoDB 由于公司访问服务器要通过堡垒机,所以本机 ...

  2. 【Selenium2】【Shell】

    E:\test_object>Python all_test.py >> report/log.txt 2>&1

  3. unity shader base pass and additional pass

      [Unity Shaders]Shader中的光照,shadersshader 写在前面 自己写过Vertex & Fragment Shader的童鞋,大概都会对Unity的光照痛恨不已 ...

  4. python中常用的模块一

    一,常用的模块 模块就是我们将装有特定功能的代码进行归类,从代码编写的单位来看我们的程序,从小到大的顺序: 一条代码<语句块,<代码块(函数,类)<模块我们所写的所有py文件都是模块 ...

  5. python字符串用法

    一,数据类型的转换 1常用的数据类型有int,str,bool三种; int是整型:一般的操作是查找数据的二进制长度bit_length() 一般数据有: 1)二进制:逢二进一 2)十进制,(0-9) ...

  6. modal结合art-template

    内容div <div id="modal-cont"></div> 模板tpl <script id="modal-tpl" ty ...

  7. The content of element type "web-app" must match "(icon?,display-name?,description?,distributable?,context-param*,filter*,filter-mapping*,listener*,servlet*,servlet- mapping*,session-config?,mime-map

    修改了一下web.xml,加入了一个<filter>,然后就报这样的错??? The content of element type "web-app" must ma ...

  8. HDU 2569(简单的递推)

    彼岸 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submissi ...

  9. Codeforces 1077 F2 - Pictures with Kittens (hard version)

    F2 - Pictures with Kittens (hard version) 思路: 单调队列优化dp 代码: #pragma GCC optimize(2) #pragma GCC optim ...

  10. ASA与N6K对接

    ASA5545配置interface GigabitEthernet0/0 channel-group 10 mode active no nameif no security-level no ip ...