给定一个长度为 \(n\) 的置换 \(p\)，以及一个正整数 \(k\)．
对于一个置换 \(q\)，要求对于所有满足 \(1\leq i<j\leq n\) 的 \(i,j\)，有以下不等式成立：
\[p_{q_i}\leq p_{q_j}+k \]现在请求出满足条件的置换 \(q\) 中，逆序对数最小的 \(q\)，它逆序对数是多少。
\(1\leq n\leq 5000,1\leq k\leq 8\)．
相当于要对 \(p\) 进行重排，如果将某个前缀出现过的 \(p\) 看作一个二进制数，那么即为要求最高的 \(1\) 和最低的 \(0\) 之间相差不能超过 \(k\)．所以我们知道如果填了 \(i\)，那么 \(<i-k\) 的那些数一定已经填过了。所以先 \(\mathcal{O}(n^2)\) 把 \(<i-k\) 的那些逆序对给统计上。
那么就好 dp 了，\(f_{i,S}\) 表示从低到高填到了 \(i\)，然后 \(i\) 前面 \(k\) 个数填的状态是 \(S\)，只记前 \(k\) 个是因为由于条件的限制 \(<i-k\) 的都一定已经填过了。
这样时间复杂度就是 \(\mathcal{O}(n^2+nk^22^k)\)，因为转移的时候枚举下一个填啥是 \(\mathcal{O}(k)\) 的，还有一个暴力 \(\mathcal{O}(k)\) 计算逆序对贡献。
#include<cstdio>
#include<vector>
#include<queue>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<ctime>
#include<random>
#include<assert.h>
#define pb emplace_back
#define mp make_pair
#define fi first
#define se second
#define dbg(x) cerr<<"In Line "<< __LINE__<<" the "<<#x<<" = "<<x<<'\n';
#define dpi(x,y) cerr<<"In Line "<<__LINE__<<" the "<<#x<<" = "<<x<<" ; "<<"the "<<#y<<" = "<<y<<'\n';
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int>pii;
typedef pair<ll,int>pli;
typedef pair<ll,ll>pll;
typedef pair<int,ll>pil;
typedef vector<int>vi;
typedef vector<ll>vll;
typedef vector<pii>vpii;
typedef vector<pil>vpil;
template<typename T>T cmax(T &x, T y){return x=x>y?x:y;}
template<typename T>T cmin(T &x, T y){return x=x<y?x:y;}
template<typename T>
T &read(T &r){
r=0;bool w=0;char ch=getchar();
while(ch<'0'||ch>'9')w=ch=='-'?1:0,ch=getchar();
while(ch>='0'&&ch<='9')r=r*10+(ch^48),ch=getchar();
return r=w?-r:r;
}
template<typename T1,typename... T2>
void read(T1 &x,T2& ...y){read(x);read(y...);}
const int mod=998244353;
inline void cadd(int &x,int y){x=(x+y>=mod)?(x+y-mod):(x+y);}
inline void cdel(int &x,int y){x=(x-y<0)?(x-y+mod):(x-y);}
inline int add(int x,int y){return (x+y>=mod)?(x+y-mod):(x+y);}
inline int del(int x,int y){return (x-y<0)?(x-y+mod):(x-y);}
const int inf=0x3f3f3f3f;
inline int bit(int x){
return 1<<(x-1);
}
inline int all(int x){
return (1<<x)-1;
}
const int N=5010;
int n,k,p[N],vis[N],c[N];
int f[N][260],ans;
int calc(int x,int s,int i){
int sum=vis[x]>vis[x-i];
for(int j=1;j<=k&&x-j>=1;j++)
if(bit(j)&s)
if(vis[x-j]>vis[x-i])
++sum;
for(int j=k+1;x-j>=1&&j<=k+i;j++)
if(vis[x-j]>vis[x-i])
++sum;
return sum;
}
signed main(){
#ifdef do_while_true
// assert(freopen("data.in","r",stdin));
// assert(freopen("data.out","w",stdout));
#endif
read(n,k);
for(int i=1;i<=n;i++)read(p[i]),vis[p[i]]=i;
for(int i=1;i<=n;i++)
for(int j=1;j<i-k;j++)
if(vis[j]>vis[i])
++ans;
memset(f,0x3f,sizeof(f));
f[0][(1<<k)-1]=0;
for(int i=0;i<=n;i++)
for(int s=0;s<(1<<k);s++)if(f[i][s]!=inf){
for(int j=1;j<=k;j++)
if(!(bit(j)&s)){
cmin(f[i][s|(1<<(j-1))],f[i][s]+calc(i,s,j));
}
for(int j=1;j<=k+1&&i+j<=n;j++){
if(j<=k && !(bit(k-j+1)&s))break;
int to=((s<<j)|bit(j))&all(k);
cmin(f[i+j][to],f[i][s]+calc(i+j,to,0));
}
}
cout << f[n][(1<<k)-1]+ans << '\n';
#ifdef do_while_true
// cerr<<'\n'<<"Time:"<<1.0*clock()/CLOCKS_PER_SEC*1000<<" ms"<<'\n';
#endif
return 0;
}