以1,1为起点
设dp[i][j] 为 以i个数全排列中j上升段的个数。
特殊的dp[1][1] 为 1
转移方程 看上一次某些状态
假设上一次的一个状态为 dp[i-1][j]
可以把当前的i这个数每次分别插入到每一段的末尾, 即dp[i][j] = j*dp[i-1][j];
假设上一次的一个状态为 dp[i-1][j-1]
可以把当前的i这个数每次分别插入到除每一段的末尾的地方, 即dp[i][j] = (i-j+1)*dp[i-1][j-1];
对于其它状态dp[i-1][j-2]不可能转移到dp[i][j].
故有dp[i][j] = j * dp[i-1][j] + (i-j+1)*dp[i-1][j-1];
与欧拉数三角形有关系,这个看着有数字三角形那个感觉。
代码:
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
#define int long long
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const double eps = 1e-4;
const int mod = 999911659;
const int N = 110;
int dp[N][N];
/*
以1,1为起点
设dp[i][j] 为 以i个数全排列中j上升段的个数。
特殊的dp[1][1] 为 1
转移方程 看上一次某些状态
假设上一次的一个状态为 dp[i-1][j]
可以把当前的i这个数每次分别插入到每一段的末尾, 即dp[i][j] = j*dp[i-1][j];
假设上一次的一个状态为 dp[i-1][j-1]
可以把当前的i这个数每次分别插入到除每一段的末尾的地方, 即dp[i][j] = (i-j+1)*dp[i-1][j-1];
对于其它状态dp[i-1][j-2]不可能转移到dp[i][j].
故有dp[i][j] = j*dp[i-1][j] + (i-j+1)*dp[i-1][j-1];
与欧拉数三角形有关系
*/
signed main(){
// IOS;
#ifdef ddgo
freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
#endif
int n,m; cin>>n>>m;
dp[1][1] = 1;
for(int i=2;i<=n;i++)
for(int k=1;k<=i;k++)
dp[i][k] = k*dp[i-1][k]+(i-k+1)*dp[i-1][k-1];
cout<<dp[n][m]<<endl;
return 0;
}