wannafly day4 G 置置置换

题目链接：https://www.zhixincode.com/contest/17/problem/G?problem_id=254

题目大意：求合法的1-n的排列数满足奇数为峰，偶数为谷（n>=2）

 

法1 DP

状态表示有点不太好理解，解法参考https://blog.csdn.net/jayye1994/article/details/12361481 转自JayYe

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#define ll long long
using namespace std;
const int maxn=1005;
const int mod=1e9+7;
ll dp[maxn][maxn];

int main()
{
    memset(dp,0,sizeof(dp));
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        dp[1][i]=1;
    }
    for(int i=2;i<=n;i++)//前i个数
    {
        if(i%2==1)
        {
            dp[i][1]=0;
            for(int j=2;j<=i;j++)
            {
                dp[i][j]=(dp[i][j-1]+dp[i-1][j-1])%mod;
            }
        }
        else
        {
            dp[i][i]=0;
            for(int j=i-1;j>=1;j--)
            {
                dp[i][j]=(dp[i][j+1]+dp[i-1][j])%mod;
            }
        }
    }
    ll ans=0;
    for(int i=1;i<=n;i++)
    {
        ans+=dp[n][i];
        ans%=mod;
    }
    cout<<ans<<endl;

    return 0;
}

法2：组合数学法

参考当晚上台讲题的同学 (流弊！)

题解思想：ans[i]表示前i个数合法排列的结果

ans[x+1]就是将x+1固定，在其左边和右边各插入x个数

1.要满足左右独自合法，这里代入ans[i]，数字大小改变了，但是相对大小关系不变，选i个数的合法排列答案就是ans[i]

2.要满足x+1这个数所在的位置合法，这个数一定最大，所以必为峰，只能在偶数位，*inv[2]就可以了

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#define ll long long
using namespace std;
const int N =1005;
const int mod =(int)1e9+7;
int comb[N][N];
ll ans[N];
ll f[N],finv[N],inv[N];//f中存的是阶乘，finv中存的是逆元的阶乘（分母）
void init()
{
    inv[1]=1;
    for(int i=2;i<N;i++)
    {
        inv[i]=(mod-mod/i)*1ll*inv[mod%i]%mod;
    }
    f[0]=finv[0]=1;
    for(int i=1;i<N;i++)
    {
        f[i]=f[i-1]*1ll*i%mod;
        finv[i]=finv[i-1]*1ll*inv[i]%mod;
    }
}
int cal(int n,int m)
{
    if(m<0||m>n)return 0;
    return f[n]*1ll*finv[n-m]%mod*finv[m]%mod;
}



int main()
{
    init();
    int n;
    cin>>n;
    ans[1]=ans[0]=1;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<=i;j++)
        {
            ans[i+1]=(ans[i+1]%mod+(cal(i,j)*ans[j]%mod*ans[i-j])%mod)%mod;
        }
        ans[i+1]=(ans[i+1]*inv[2])%mod;
    }
    cout<<ans[n]<<endl;


    return 0;
}