Description

DJL為了避免成為一只鹹魚，來找czgj學習Fibonacci數列。
通過czgj的諄諄教導，DJL明白了Fibonacci數列是這樣定義的：
F(1)=1;F(2)=1;F(n)=F(n-1)+F(n-2)(n>2)
Czgj深諳熟能生巧的道理，於是他給了DJL一個數列，並安排了如下的訓練計劃：
1、“1 L r”，表示給ai 加上F(i-L+1) ，其中L<=i<=r ；
2、“2 L r”，表示詢問 的值。
DJL經過長時間的學習，感覺身體被掏空，他希望你能幫他解決這個問題。

4 4
1 2 3 4
1 1 4
2 1 4
1 2 4
2 1 3

17
12
樣例解釋
經過第一次操作，數列變為a=[2,3,5,7] ；
 
第二次詢問，sum=2+3+5+7=17 ；
經過第三次操作，數列變為a=[2,4,6,9] ；
第四次詢問，sum=2+4+6=12 。

平衡規劃

對於區間修改操作先不直接修改，每當儲存到了√n個後，再進行重構。每次詢問的話，先查詢重構了的數組，再遍歷儲存下來的操作，然後再單獨做，可以優化時間復雜度

題解

• 那麽有了上面的思想，套用到這題目裏
• 每次修改就是O(1)的，然後維護一個數組d
• 當一次區間[l,r]修改時
• d[l]=d[l]+1,d[r+1]=d[r+1]-f[r-l+2],d[r+2]=d[r+2]-f[r-l+1],d[i]=d[i]+d[i-1]+d[i-2]
• 然後重構a數組，並求出前綴和
• 最後時間復雜度就是O(m√n)

代碼

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <cmath>
 4 #include <cstring>
 5 using namespace std;
 6 struct edge {int l,r;}e[100010];
 7 int n,m,a[100010],cnt,feng;
 8 long long sum[100010],k[100010],tot[100010],f[100010],ans,mo=1e9+9;
 9 int main()
10 {
11     scanf("%d%d",&n,&m); feng=sqrt(m);
12     for (int i=1;i<=n;i++) scanf("%d",&a[i]),sum[i]=(sum[i-1]+a[i])%mo;
13     f[1]=tot[1]=1;
14     for (int i=2;i<=n;i++) f[i]=(f[i-1]+f[i-2])%mo,tot[i]=(tot[i-1]+f[i])%mo;
15     for (int j=1;j<=m;j++)
16     {
17         int d,l,r;
18         scanf("%d%d%d",&d,&l,&r);
19         if (d==1) e[++cnt].l=l,e[cnt].r=r; 
20         else 
21         {
22             ans=0;
23             for (int i=1;i<=cnt;i++)
24             {
25                 int a=e[i].l,b=e[i].r;
26                 if (r<a||l>b) continue;
27                 ans=(ans+tot[min(b,r)-(a-1)]-tot[max(a,l)-a]+mo)%mo;
28             }
29             printf("%lld\n",(ans+sum[r]-sum[l-1]+mo)%mo);
30         }
31         if (cnt==feng)
32         {
33             memset(k,0,sizeof(k));
34             for (int i=1;i<=cnt;i++)
35             {
36                 int l=e[i].l,r=e[i].r;
37                 k[l]=(k[l]+1)%mo;
38                 k[r+1]=(k[r+1]-f[r-l+2]+mo)%mo,k[r+2]=(k[r+2]-f[r-l+1]+mo)%mo;
39             }
40             cnt=0;
41             for (int i=1;i<=n;i++)
42             {
43                 k[i]=(k[i]+k[i-1]+k[i-2])%mo;
44                 a[i]=(a[i]+k[i])%mo;
45                 sum[i]=(sum[i-1]+a[i])%mo;
46             }
47         }
48     }
49 }

[平衡規劃][模擬][前綴和] Jzoj P4724 斐波那契