题目描述：http://www.lydsy.com/JudgeOnline/problem.php?id=4621
题解：
观察到原排列进行k次操作后同种数字必成为连续的一段，我们可以定义f[i,j]表示前i个数字中进行了j次操作。
DP时先暴力找出可以把当前a[i]扩展到的区间找到，记为[l, r]。那么可以得出一个非常显然的DP方程：
dp[i][j]+=Σdp[k][j-1], l-1<=k<r （具体情况画一下图应该比较就能得出来了）
特别的，这样计算时会将本次操作区间长度为1的算进去，即有：f[i][j]+=f[i-1][j-1]
但这样的操作是没有意义的，对原数列没有一点改变，所以要将这一步的影响消除：
f[i][j-1]=f[i][j-1]+f[i-1][j-1];f[i][j]=f[i][j]-f[i-1][j-1];（不计算这一步操作）
时间复杂度为O（n^3）
参考代码：
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
#define rep(i,j,k) for(int i=(int)j;i<=(int)k;i++)
#define per(i,j,k) for(int i=(int)j;i>=(int)k;i--)
#define mod 1000000007
int n,m,a[510],f[510][510];
int main(){
scanf("%d%d",&n,&m);
rep(i,1,n)scanf("%d",&a[i]);
f[0][0]=1;
int l,r,tmp;
rep(i,1,n){
l=r=i;
while(l>1 && a[l-1]<a[i])l--;
while(r<n && a[r+1]<a[i])r++;
f[i][m]=(f[i][m]+f[i-1][m])%mod;
per(j,m,1){
tmp=f[l-1][j-1];
rep(k,l,r){
f[k][j]=(f[k][j]+tmp)%mod;
tmp=(tmp+f[k][j-1])%mod;
}
f[i][j-1]=(f[i][j-1]+f[i-1][j-1])%mod;
f[i][j]=(f[i][j]-f[i-1][j-1]+mod)%mod;
}
}
int ans;
rep(i,0,m)ans=(ans+f[n][i])%mod;
printf("%d\n",ans);
return 0;
}