浅谈CDQ分治

前置知识：分治

步骤：

           1、把原问题分为若干个形式相同，规模更小的子问题，直到子问题足够小

           2、递归求解子问题

           3、合并子问题的解，得到原问题的解

__________________________________________________________________

 

CDQ分治：一种特殊的分治，可以代替各种高级数据结构，常数较小，但是必须离线处理

它可以把时间看做一层维度，把插入询问的操作变成先插入后删除的操作

思路如下：          

         1、把区间分成

         2、递归求解

         3、处理对的影响

例1：逆序对

所谓的归并排序模板题，也就是一种CDQ分治

#include<cstdio>
#define ll long long
using namespace std;

const int N=1e6+5;
int n,a[N],t[N];
ll ans;

void work(int l,int r)
{
	if(l==r) return;   //足够小的问题 
	int mid=l+r>>1;
	work(l,mid),work(mid+1,r);   //分解 
	int i=l,j=mid+1,k=l;
	while(i<=mid&&j<=r)
	{
		if(a[i]<=a[j]) t[k++]=a[i++];
			else t[k++]=a[j++],ans+=mid-i+1;
	}
	while(i<=mid) t[k++]=a[i++];
	while(j<=r) t[k++]=a[j++];
	for(int i=l;i<=r;i++) a[i]=t[i];    //处理前面对于后面的影响 
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	work(1,n);
	printf("%lld\n",ans);
	return 0;
} 

  如果把时间轴提取出来，就要求，满足且有几对

这类问题就叫二维偏序

时间默认有序，所以是有序的

在处理这类问题时，把排序，然后再求顺序对或者逆序对就可以了

 

练1：树状数组模板题

题解：想清楚O_O

 

例2：陌上花开

多了一维，也就是三维偏序

注意离散化+去重

按二维的套路：对排序，对进行cdq分治，就可以抛去对答案的影响（？一定要想清楚）

问题就变成了给定两个序列，求且的数量有多少（i为前一个序列的下标，j为后一个序列的下标）

这不就是二维CDQ？

注意标记，因为把b重排后a的顺序将发生变化

不得不说，这题离散化真的恶心

两个CDQ，所以时间复杂度为

#include<cstdio>
#include<algorithm>
using namespace std;

const int N=1e6+5;
int n,m,num[N],ans[N],f[N],n_n;
struct A{ int id,a,b,c; }q[N],q1[N],q2[N];
bool cmpa(A i,A j){ return i.a<j.a||i.a==j.a&&i.b<j.b||i.a==j.a&&i.b==j.b&&i.c<j.c; }

void work2(int l,int r)
{
    if(l==r) return;
    int mid=l+r>>1;
    work2(l,mid),work2(mid+1,r);
    int i=l,j=mid+1,k=l,y=0;
    while(i<=mid&&j<=r)
    {
        if(q1[i].c<=q1[j].c) q2[k++]=q1[i],y+=q1[i].a?num[q1[i].id]:0,i++; 
            else 
            {
                q2[k++]=q1[j];
                if(!q1[j].a) f[q1[j].id]+=y;
                j++;
            }
    }
    while(i<=mid) q2[k++]=q1[i++];
    while(j<=r) 
    {
        q2[k++]=q1[j];
        if(!q1[j].a) f[q1[j].id]+=y;
        j++;
    }
    for(int i=l;i<=r;i++) q1[i]=q2[i];
}

void work1(int l,int r)
{
    if(l==r) return;
    int mid=l+r>>1;
    work1(l,mid),work1(mid+1,r);
    int i=l,j=mid+1,k=l;
    while(i<=mid&&j<=r)
    {
        if(q[i].b<=q[j].b) q1[k]=q[i++],q1[k].a=1,k++; //标记是前面的还是后面的 
            else q1[k]=q[j++],q1[k].a=0,k++;
    }
    while(i<=mid) q1[k]=q[i++],q1[k].a=1,k++;
    while(j<=r) q1[k]=q[j++],q1[k].a=0,k++;
    for(int i=l;i<=r;i++) q[i]=q1[i];
    work2(l,r);
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d%d%d",&q[i].a,&q[i].b,&q[i].c);
    sort(q+1,q+n+1,cmpa);
    n_n=1; num[1]=1; q[1].id=1;
    for(int i=2;i<=n;i++) 
        if(q[i].a==q[i-1].a&q[i].b==q[i-1].b&q[i].c==q[i-1].c) num[n_n]++;
            else q[++n_n]=q[i],q[n_n].id=n_n,num[n_n]++;
    work1(1,n_n);
    for(int i=1;i<=n_n;i++) ans[f[i]+num[i]-1]+=num[i];
    for(int i=0;i<n;i++) printf("%d\n",ans[i]);
    return 0;
} 

 

那是不是可以一直叠下去？

理论上来说是的，但是到后来都没有暴力来的优秀

 

练2：动态逆序对

题解：想清楚O

 

四维偏序   留着坑等以后再填

 

dp优化CDQ分治 

练3：货币兑换

题解：想清楚O