Codeforces Round #FF (Div. 1) C. DZY Loves Fibonacci Numbers

 用线段树维护

记录三个值，一个是区间和，一个是该区间 1 , 2 位置的斐波那契数列前两项大小（因为两个斐波那契数列相加得到的还是一个斐波那契数列，且该数列只与两个首项和项数有关，另外斐波那契数列前N 项和  与  该斐波那契数列第N+2减去第二项  的值 相同）

该线段树的pushdown有一些改变，直接利用已有的update函数将懒标记下传，下传时改变的区间是传入pushdown的区间大小（虽然不太明白这个姿势有何好处，但是觉得特别方便）

#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#define CLR(x,y) memset(x,y,sizeof(x))
#define mp(x,y) make_pair(x,y)
#define eps 1e-9
#define INF 0x3f3f3f3f
#define LLINF 1LL<<62

using namespace std;

typedef long long ll;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef complex<ld> point;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;

template<class T>
inline bool read(T &n)
{
    T x = 0, tmp = 1; char c = getchar();
    while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
    if(c == EOF) return false;
    if(c == '-') c = getchar(), tmp = -1;
    while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
    n = x*tmp;
    return true;
}
template <class T>
inline void write(T n)
{
    if(n < 0)
    {
        putchar('-');
        n = -n;
    }
    int len = 0,data[20];
    while(n)
    {
        data[len++] = n%10;
        n /= 10;
    }
    if(!len) data[len++] = 0;
    while(len--) putchar(data[len]+48);
}
//-----------------------------------
#define MID(x,y) (x+((y-x)>>1))

const int MOD=1000000009;
const int MAXN=300010;

struct Node
{
	int l,r;
	ll v,x,y;
}t[MAXN<<2];

int n,m,q,l,r;
ll f[MAXN];

inline ll fib(ll a,ll b,int n)
{
    if(n==0)
        return a;
    if(n==1)
        return b;
    return (a*f[n-1]%MOD+b*f[n]%MOD)%MOD;
}

void update(int id,int l,int r,ll a,ll b);
void build(int id,int l,int r)
{
    t[id].l=l,t[id].r=r;
    t[id].v=t[id].x=t[id].y=0;
    if(l==r)
    {
        read(t[id].v);
        return;
    }
    build(id<<1,l,MID(l,r));
    build(id<<1|1,MID(l,r)+1,r);
    t[id].v=(t[id<<1].v+t[id<<1|1].v)%MOD;
}

void pushdown(int id)
{
    if (t[id].x||t[id].y)
    {
        update(id<<1,t[id].l,t[id].r,t[id].x,t[id].y);
        update(id<<1|1,t[id].l,t[id].r,t[id].x,t[id].y);
        t[id].x=0;
        t[id].y=0;
    }
}

void update(int id,int l,int r,ll a,ll b)
{
//    cout<<id<<"->("<<t[id].l<<","<<t[id].r<<") -> "<<MID(t[id].l,t[id].r)<<"-> (" <<l<<","<<r<<")"<<endl;
    if(l<=t[id].l && t[id].r<=r)
    {
        ll na,nb;
        na=fib(a,b,t[id].l-l);
        nb=fib(a,b,t[id].l-l+1);
        a=na;
        b=nb;
        t[id].x=(t[id].x+a)%MOD;
        t[id].y=(t[id].y+b)%MOD;
        t[id].v=(t[id].v+fib(a,b,t[id].r-t[id].l+2)-b)%MOD;
        return;
    }
    pushdown(id);
    int mid=MID(t[id].l,t[id].r);
    if(l<=mid)
        update(id<<1,l,r,a,b);
    if(r>mid)
        update(id<<1|1,l,r,a,b);
    t[id].v=(t[id<<1].v+t[id<<1|1].v)%MOD;
}

ll query(int id,int l,int r)
{
//    cout<<id<<"->("<<t[id].l<<","<<t[id].r<<") -> "<<MID(t[id].l,t[id].r)<<"-> (" <<l<<","<<r<<")"<<endl;
    if(l<=t[id].l && t[id].r<=r)
    {
        return t[id].v;
    }
    pushdown(id);
    ll ret=0;
    int mid=MID(t[id].l,t[id].r);
    if(l<=mid)
        ret+=query(id<<1,l,r);
    if(r>mid)
        ret+=query(id<<1|1,l,r);
    return ret%MOD;
}

int main()
{
	read(n),read(m);
	f[1]=f[2]=1;
	for(int i=3;i<MAXN;i++)
        f[i]=(f[i-1]+f[i-2])%MOD;
	build(1,1,n);
    while(m--)
    {
        read(q),read(l),read(r);
        if(q==1)
        {
            update(1,l,r,1LL,1LL);
        }
        else
        {
            ll ans=(query(1,l,r)+MOD)%MOD;
            write(ans),putchar('\n');
        }
    }
    return 0;
}