【CDQ分治】Non-Decreasing Subsequences

CDQ分治好题，至少对于我这个CDQ分治初学者

loj3247

题目：

求区间不下降子序列个数
且$a_i\le k,k\le 20,N\le 5\times 10^4,Q\le 2\times 10^5$

题解：

很难想到，要用CDQ分治
要分了之后怎么进行维护复杂度尽量小，这是CDQ分治很难受的地方。

• 定义数组
  在每个$solve(l,r)$中都有
  $$f_{i,j}表示$$
  $$(l\le j\le mid）j为左端点，k\in [j,mid]中所有a_k=i的k为右端点，不下降子序列个数$$
  $$(mid<j\le r）j为右端点，k\in (mid,j]中所有a_k=i的k为左端点，不下降子序列个数$$
• 维护数组
  尽可能的找数组之间的联系
  首先因为$mid$将$f_{i,j}$分为了两个部分，意义不同（其实泛泛来说差不多）,所以我们也分两部分讨论
  这里我们先只分析右边的，左边的同理理解
  $$f_{i,j}=[a_j==i]+\sum _{mid\le k<j,a_k\le a_j}{f}_{i,k}$$
  然后我开始看到这式子，是懵逼的
  $[a_j==i]$指如果$a_j==i$的话，$j->j$就是一个合法区间，$f_{i,j}+1$
  ${\sum }_{mid\le k<j,a_k\le a_j}{f}_{i,k}$指与前面的已合法区间拼接，因为还要满足拼接后不下降，所以$a_k\le a_j$
  于是我们可以想到树状数组维护（略）
  意右边循环顺序是$mid+1->r$，但左边循环顺序是$mid->l$,因为是逆着来的，所以注意树状数组也要逆一下

注意：

因为空集也算合法区间，所以要把$f_{1,mid},f_{k,mid+1}$都要+1

• 询问求值
  yy一下我们想到，我们需要将我们的$f_{i,j}$前缀和一下（左边是后缀和）
  我们先只考虑跨过$mid$的询问区间（ql,qr分别为询问的左右端点）
  $$ans=\sum _{1\le i\le k}{F}_{i,ql}\sum _{i\le j\le k}{F}_{j,qr}$$
  理解下，就是左右两端拼接，乘法原理
  （我曾迷惑：$ql->s,s->qr$，$mid<s$的情况没有考虑啊，但其实将$s$移到 $\overline{s}<mid$，$ql->\overline{s},\overline{s}<qr$这样是等效的）
• 复杂度 $O(nklognlogk+qk)$

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll mod=1e9+7;
const int N=5e4+10,K=25,Q=2e5+10;
int n,k,q;
int a[N];
struct qus{int l,r;}qs[Q];
ll ans[Q];
ll f[K][N];
//---
ll t[K];
int lowbit(int x){return x&(-x);}
void clear(){memset(t,0,sizeof(t));}
void add(int i,int x){for(;i<=k;i+=lowbit(i))t[i]=(t[i]+x)%mod;}
ll qur(int i){ll res=0;for(;i;i-=lowbit(i))res=(res+t[i])%mod;return res;}
//----
ll ls[N];
void solve(int l,int r,vector<ll>v)
{
	if(l==r)
	{
		for(int i=0;i<v.size();i++)ans[v[i]]=2;
		return ;
	}
	int mid=(l+r)>>1;
	for(int i=1;i<=k;i++)
	{
		clear();
		for(int j=mid;j>=l;j--)
		{
			f[i][j]=qur(k+1-a[j]);
			if(a[j]==i)f[i][j]=(f[i][j]+1)%mod;//自己->自己 
			add(k+1-a[j],f[i][j]);//逆转树状数组 
		}
		for(int j=mid-1;j>=l;j--)f[i][j]=(f[i][j]+f[i][j+1])%mod;//前缀和 
	}
	for(int j=l;j<=mid;j++)f[1][j]=(f[1][j]+1)%mod;//空集
	for(int i=1;i<=k;i++)
	{
		clear();
		for(int j=mid+1;j<=r;j++)
		{
			f[i][j]=qur(a[j]);
			if(a[j]==i)f[i][j]=(f[i][j]+1)%mod;
			add(a[j],f[i][j]);
		}
		for(int j=mid+2;j<=r;j++)f[i][j]=(f[i][j]+f[i][j-1])%mod;
	}
	for(int j=mid+1;j<=r;j++)f[k][j]=(f[k][j]+1)%mod;
	vector<ll>vl,vr;
	vl.resize(0);vr.resize(0); 
	for(int i=0;i<v.size();i++)
	{
		int l=qs[v[i]].l,r=qs[v[i]].r;
		if(r<=mid)	  vl.push_back(v[i]);
		else if(mid<l)vr.push_back(v[i]);
		else
		{
			for(int j=k;j>=1;j--)ls[j]=(ls[j+1]+f[j][r])%mod;//前缀又后缀。。。
			for(int j=1;j<=k;j++)ans[v[i]]=(ans[v[i]]+ls[j]*f[j][l]%mod)%mod;
		}
	}
	if(vl.size())solve(l,mid,vl);
	if(vr.size())solve(mid+1,r,vr);
}
vector<ll>v;
int main()
{
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	scanf("%d",&q);
	for(int i=1;i<=q;i++)scanf("%d%d",&qs[i].l,&qs[i].r);
	for(int i=1;i<=q;i++)v.push_back(i);
	solve(1,n,v);
	for(int i=1;i<=q;i++)printf("%lld\n",ans[i]);
}

从这位巨佬博客看懂的