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经过长时间的学习,感觉身体被掏空,他希望你能帮他解决这个问题。
Input
第一行两个整数n和m,表示原始数列的长度,和总的训练次数。
第二行n个整数a1,a2,…,an(1<=ai<=10^9) ,表示czgj给DJL的原始数列。
接下来m行,每一行给出三个整数,表示问题描述中的两种训练的一种。保证1<=L<=r<=n 。
Output
对于每一种形如“2 L r”的训练,输出一行一个整数值。
Sample Input
4 4
1 2 3 4
1 1 4
2 1 4
1 2 4
2 1 3
Sample Output
17
12
样例解释
经过第一次操作,数列变为a=[2,3,5,7] ;
第二次询问,sum=2+3+5+7=17 ;
经过第三次操作,数列变为a=[2,4,6,9] ;
第四次询问,sum=2+4+6=12 。
Data Constraint
对于20%的数据,1≤n, m≤100;
对于40%的数据,1≤n, m≤1000;
对于100%的数据,1≤n, m≤100000。
思路
一开始的想法是线段树,但是。。。好麻烦啊!!!
题解说,分块就可以过啦!
由于是斐波那契数列,所以我们可以弄一个类似于差分数组的东西
每次修改就是O(1)的
每次查询就扫一遍之前的询问,把所有书加上,再加上数组a即可。
但如果单纯暴力肯定会TLE
所以考虑把操作分成k段,每一段就对a数组进行重构
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn=1e5+77,mod=1e9+9;
int n,m,a[maxn],f[maxn],sf[maxn],s[maxn],c1[maxn],c2[maxn];
struct B
{
int x,l,r;
}b[maxn];
int main()
{
scanf("%d%d",&n,&m);
for(int i=1; i<=n; i++) scanf("%d",&a[i]),s[i]=(s[i-1]+a[i])%mod;
int k=sqrt(m)+1,t=0;
f[1]=1; f[2]=1;
for(int i=3; i<=n; i++) f[i]=(f[i-1]+f[i-2])%mod;
for(int i=1; i<=n; i++) sf[i]=(sf[i-1]+f[i])%mod;
while(m--)
{
t++;
scanf("%d%d%d",&b[t].x,&b[t].l,&b[t].r);
if(b[t].x==1)
{
int l=b[t].l,r=b[t].r;
c1[l]=(c1[l]+1)%mod;
c2[r+1]=(c2[r+1]-f[r-l+1])%mod;
c1[r+1]=(c1[r+1]-f[r-l+2])%mod;
}
if (b[t].x==2)
{
int l=b[t].l,r=b[t].r;
long long ans=(s[r]-s[l-1]+mod)%mod;
for(int i=1; i<t; i++)
{
if (b[i].x==1)
{
if (b[i].r<l||b[i].l>r) continue;
int ll=max(l,b[i].l),rr=min(r,b[i].r);
ans=(ans+sf[rr-b[i].l+1]-sf[ll-b[i].l]+mod)%mod;
}
}
printf("%lld\n",ans);
}
if(t==k)
{
for (int i=1;i<=n;i++)
{
c2[i]=((c2[i]+c1[i-1])%mod+mod)%mod;
c1[i]=((c1[i]+c2[i-1]+c1[i-1])%mod+mod)%mod;
a[i]=(a[i]+c1[i])%mod;
s[i]=(s[i-1]+a[i])%mod;
}
memset(c1,0,sizeof(c1)); memset(c2,0,sizeof(c2));
t=0;
}
}
}