CF633H Fibonacci-ish II

VIII.CF633H Fibonacci-ish II

强烈谴责卡常屑题(

莫队。于是要支持插入数、查询与斐波那契数列的点积和。

离散化后就直接在线段树上处理。于是变成单点插入/删除(这个很 easy)以及后缀斐波那契数列前移/后移。

直接上矩阵维护。但是会TLE(因为常数是 \(2^3\))。预处理出斐波那契数列转移矩阵的正负 \(n\) 次幂,就能将转移优化到 \(2^2\),实际效果还要更好(时间是原来的 \(1/4\))。

代码:

#include<bits/stdc++.h>
using namespace std;
int n,mod,a[30100],buc[30100],res[30100],m,lim;
const int BBB=240;
struct Matrix{
	int a[2][2];
	int*operator[](const int &x){return a[x];}
	Matrix(){memset(a,0,sizeof(a));}
	friend Matrix operator*(register Matrix&u,register Matrix&v){
		Matrix w;
		for(register int i=0;i<2;i++)for(register int j=0;j<2;w[i][j++]%=mod)for(register int k=0;k<2;k++)w[i][j]+=u[i][k]*v[k][j];
		return w;
	}
	void print()const{for(register int i=0;i<2;i++){for(register int j=0;j<2;j++)printf("%d ",a[i][j]);puts("");}}
}F,G,I,f[30100],g[30100];
Matrix ksm(register Matrix x,register int y){register Matrix z=I;for(;y;y>>=1,x=x*x)if(y&1)z=z*x;return z;}
vector<int>v;
struct SegTree{int s[2],sz,tag;}seg[120100];
inline void MUL(register int *x,register Matrix M){
	register int z[2];z[0]=z[1]=0;
	for(int i=0;i<2;i++)for(int j=0;j<2;j++)z[j]+=x[i]*M[i][j];
	x[0]=z[0]%mod,x[1]=z[1]%mod;
}
#define lson x<<1
#define rson x<<1|1
void SUM(int&x,int&y,int&z){x=y+z;if(x>=mod)x-=mod;}
inline void MUL(register int x,register int y){MUL(seg[x].s,y>=0?f[y]:g[-y]),seg[x].tag+=y;}
inline void pushdown(register int x){if(seg[x].tag)MUL(lson,seg[x].tag),MUL(rson,seg[x].tag),seg[x].tag=0;}
inline void pushup(register int x){for(register int i=0;i<2;i++)SUM(seg[x].s[i],seg[lson].s[i],seg[rson].s[i]);seg[x].sz=seg[lson].sz+seg[rson].sz;}
inline void turnon(register int x,register int l,register int r,register int P,register int val){
	if(l==r){seg[x].s[0]=0,seg[x].s[1]=v[P-1],MUL(x,val),seg[x].sz=1;return;}
	register int mid=(l+r)>>1;pushdown(x);if(P<=mid)turnon(lson,l,mid,P,val);else turnon(rson,mid+1,r,P,val);pushup(x);
}
inline void turnoff(register int x,register int l,register int r,register int P){
	if(l==r){seg[x].s[0]=seg[x].s[1]=seg[x].sz=0;return;}
	register int mid=(l+r)>>1;pushdown(x);if(P<=mid)turnoff(lson,l,mid,P);else turnoff(rson,mid+1,r,P);pushup(x);
}
inline int querysz(register int x,register int l,register int r,register int P){
	if(l>P)return 0;if(r<=P)return seg[x].sz;
	register int mid=(l+r)>>1;return querysz(lson,l,mid,P)+querysz(rson,mid+1,r,P);
}
inline void modifymatrix(register int x,register int l,register int r,register int P,register int tp){
	if(r<P)return;if(l>=P){MUL(x,tp);return;}
	register int mid=(l+r)>>1;pushdown(x),modifymatrix(lson,l,mid,P,tp),modifymatrix(rson,mid+1,r,P,tp),pushup(x);
}
struct QUERY{int l,r,id;friend bool operator<(const QUERY&u,const QUERY&v){if(u.l/BBB!=v.l/BBB)return u.l<v.l;return(u.l/BBB)&1?u.r<v.r:u.r>v.r;}}q[30100];
inline void Push(register int x){if(!buc[a[x]]++)turnon(1,1,lim,a[x],querysz(1,1,lim,a[x])),modifymatrix(1,1,lim,a[x],1);}
inline void Pop(register int x){if(!--buc[a[x]])turnoff(1,1,lim,a[x]),modifymatrix(1,1,lim,a[x],-1);}
char buf[1<<20],*ss=buf,*tt=buf;
#define gc() (tt==ss&&(tt=(ss=buf)+fread(buf,1,1<<20,stdin),ss==tt)?EOF:*ss++)
inline void read(register int &x){x=0;register char c=gc();while(c>'9'||c<'0')c=gc();while(c>='0'&&c<='9')x=(x<<3)+(x<<1)+(c^48),c=gc();}
inline void print(register int x){if(x<=9)putchar('0'+x);else print(x/10),putchar('0'+x%10);}
int main(){
	read(n),read(mod);
	F[0][0]=0,F[0][1]=1,F[1][0]=1,F[1][1]=1;
	G[0][1]=1,G[0][0]=mod-1,G[1][0]=1,G[1][1]=0;
	I[0][0]=1,I[1][1]=1,I[1][0]=0,I[0][1]=0;
	f[0]=g[0]=I;for(int i=1;i<=n;i++)f[i]=f[i-1]*F,g[i]=g[i-1]*G;
	for(register int i=1;i<=n;i++)read(a[i]),v.push_back(a[i]);
	sort(v.begin(),v.end()),v.resize(lim=unique(v.begin(),v.end())-v.begin());
	for(register int i=1;i<=n;i++)a[i]=lower_bound(v.begin(),v.end(),a[i])-v.begin()+1;
	for(int i=0;i<v.size();i++)v[i]%=mod;
	read(m);
	for(register int i=1;i<=m;i++)read(q[i].l),read(q[i].r),q[i].id=i;
	sort(q+1,q+m+1);
	register int L=1,R=0;
	for(register int i=1;i<=m;i++){
		while(L>q[i].l)Push(--L);
		while(R<q[i].r)Push(++R);
		while(L<q[i].l)Pop(L++);
		while(R>q[i].r)Pop(R--);
//		iterate(1,1,lim);puts("");
		res[q[i].id]=seg[1].s[0];
	}
	for(int i=1;i<=m;i++)print(res[i]),putchar('\n');
	return 0;
}

上一篇:seata1.4.1简单记录


下一篇:Django框架之forms组件