兔子 （dp）

兔子(rabbit)

10.26

第1,2个测试点:
O(n!). STL中的next_permutation()就可以得到这20分
第3,4个测试点:
O(n!)打表,O(1)回答. STL中的next_permutation()跑个几分钟就出来了
第5,6,7个测试点:
O(2n*n2)状压DP,f[S][i][0/1]表示已经使用了数字集合S中的数字,最后一个数字是i,最后一个数字比倒数第二个数字大/小的方案数.
第8,9个测试点:
定义f[i][j]表示长度为i,第一个元素为j,且第一个元素大于第二个元素的所有排列中,”抖动排列”的数目.g[i][j]表示长度为i,第一个元素为j,且第一个元素小于第二个元素的所有排列中,”抖动排列”的数目.考虑长度为i,第一个元素为j,且第一个元素大于第二个元素的一个排列,我们删掉它的第一个元素,然后把其他所有大于i的元素都减1,可以得到一个1到i-1的排列,这个排列是第一个元素小于第二个元素的.我们枚举删掉第一个元素之后的新排列的第一个元素的大小,发现f[i][j]可以从g[i-1][1,2,3…(j-1)]转移过来.g[i][j]的求法类似.
DP方程:
f[i][j]=g[i-1][1]+g[i-1][2]+g[i-1][3]+…+g[i-1][j-1]
g[i][j]=f[i-1][j]+f[i-1][j+1]+…+f[i-1][i-1]
直接转移,复杂度为O(n3)
第10个测试点:
在第8,9个测试点的DP上加入前缀和优化即可.
注意模数的范围是int范围内,两个int范围内的数字相加可能会炸int,在第10个测试点中设置了这个坑点,增加AC此题的难度.

思路：
f[i][0/1] 表示前i个数1~i，偶数位置比奇数位置大（奇数位置比偶数位置大）的方案数，然后就考虑i，可以怎么放进前i-1个数的排列中，枚举放进的位置k，然后在i-1个数中选k-1个放在前面（组合数）再按照题意转移（前后两段）。（边乘边模呀~~~）

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define LL long long
#define N 2010
using namespace std;

int n;
LL mod;
LL f[N][2], comb[N][N];

void init(){
    comb[0][0] = 1;
    for(int i=1; i<=n; i++){
        comb[i][0] = 1;
        for(int j=1; j<=i; j++){
            comb[i][j] = (comb[i-1][j-1] + comb[i-1][j]) % mod;
        }
    }
}

int main(){
    freopen ("rabbit.in", "r", stdin);
    freopen ("rabbit.out", "w", stdout);
    cin >> n >> mod;
    init();
    f[0][1] = f[0][0] = 1;
    f[1][1] = f[1][0] = 1;
    for(int i=2; i<=n; i++){
        for(int k=1; k<=i; k++){
            if((k & 1) == 0) f[i][0] = (f[i][0] + f[k-1][0] * f[i-k][0] % mod * comb[i-1][k-1]) % mod;
            else f[i][1] = (f[i][1] + f[k-1][1] * f[i-k][0] % mod * comb[i-1][k-1]) % mod;
        }
    }
    LL ans = (f[n][0] + f[n][1]) % mod;
    cout << ans << endl;
}