bzoj 5125: [Lydsy12月赛]小Q的书架

题意

给你一个序列，要你把它分成k段，使得每段的逆序对和最小
k≤10

题解

这题还是很棒的题啊
考虑到k很小，于是我们可以考虑一个十分暴力的DP
$f[i][j]$表示前i位，分了k段的最优值
转移和很简单，枚举这一段开头在哪里，加上逆序对个数就可以了
显然地，这个决策点的移动是单调的
就是随着i的右移，决策点肯定会右移
但是会有一个问题，就是他的决策值不是单调的，所以我们不可以使用单调队列
因此我们考虑使用对决策分治
我们对于分治，每一次传入代计算区域以及决策区域
取$mid$值，对他决策
就找到决策点了
至于怎么求所要区间的逆序对
有一个十分巧妙的写法，我们可以向莫队一样移端点
可以发现，如果我们每一次，先递归左儿子，再递归右儿子（反过来也行，只要不是随机访问就好）
那么你会发现，左端点每一次是从1移到n，一共移动O(n)次
右端点在不断往返，但是从一层来看，移动了n次
然后一共有log层，每一层访问n个元素，因此这里是$nlogn$的
但是由于要维护逆序对，如果使用树状数组的话，就会多一个log
Claris题解上说可以进行可持久化分块，但是我不会。。

CODE:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int N=40005;
int n,k;
int a[N];
int tr[N];//树状数组
int lb (int x){return x&(-x);}
void change (int x,int y){while (x<=n) {tr[x]+=y;x+=lb(x);}}
int find (int x)
{
	int lalal=0;
	while (x>=1){lalal=lalal+tr[x];x-=lb(x);}
	return lalal;
}
int g[N],f[N];
int L,R,now;
void go (int l,int r)
{
	while (R<r) now=now+(R-L+1-find(a[R+1])),change(a[++R],1);
	while (L<l) now=now-find(a[L]-1),change(a[L++],-1);
	while (L>l) now=now+find(a[L-1]-1),change(a[--L],1);
	while (R>r) now=now-(R-L+1-find(a[R])),change(a[R--],-1);
}
void solve (int l,int r,int dl,int dr)//要处理啊的区间    可以转移的区间
{
	int mid=(l+r)>>1,dm=dl;
	for (int u=dl;u<=min(mid-1,dr);u++)
	{
		go(u+1,mid);
		int t=g[u]+now;
		if (t<f[mid]) f[mid]=t,dm=u;
	}
	if (l<mid) solve(l,mid-1,dl,dm);
	if (mid<r) solve(mid+1,r,dm,dr);
}
int main()
{
	scanf("%d%d",&n,&k);
	for (int u=1;u<=n;u++) scanf("%d",&a[u]);
	for (int u=1;u<=n;u++)
	{
		change(a[u],1);
		f[u]=f[u-1]+u-find(a[u]);
	}
	for (int u=2;u<=k;u++)
	{
		memcpy(g,f,sizeof(g));
		memset(f,63,sizeof(f));
		memset(tr,0,sizeof(tr));
		L=1;R=0;now=0;
		solve(u,n,u-1,n);
	}
	printf("%d",f[n]);
	return 0;
}