题目的意思就是一个序列本来是0,1,2...n-1.然后要你变成p[0],p[1],p[2]...p[n-1],通过交换相邻的两个数每两个相邻的数都可以交换一次,你可以选择交换的顺序.问你有多少种交换顺序?
很有意思的一道记忆化搜索题目,emm,大佬教我做的...
具体怎么做呢?就是一个典型的分治,加记忆化搜索..分治就是枚举要交换的顺序,具体在哪个点交换,很显然的一个结论,假如那个点交换了,就相当于断开了,然后就可以分治了...具体再记忆化搜索下就可以了~dp[l][r][st][ed].表示区间l-r中开头为st,结尾是ed,然后交换的好序列的方案数是多少.然后复杂度是O(n^5)
代码如下:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int N=52;
ll dp[N][N][N][N];
bool mp[N][N][N][N];
ll p[N];
ll fact[N];
inline ll qp(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1) res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res;
}
inline ll C(ll n,ll m)
{
ll res=fact[n];
return res*qp(fact[m],mod-2)%mod*qp(fact[n-m],mod-2)%mod;
}
void init()
{
fact[1]=1;fact[0]=1;
for(int i=2;i<=N-1;i++)
{
fact[i]=i*fact[i-1]%mod;
}
}
inline ll dfs(ll l,ll r,ll st,ll ed)
{
ll res=0,x1,y1,x2,y2;
if(l==r)//递归出口.
{
if(st==p[l]) return 1;
return 0;
}
if(mp[l][r][st][ed]) return dp[l][r][st][ed];
mp[l][r][st][ed]=true;
for(ll i=l;i<r;i++)//枚举交换位子. st i+1 i+2 ... ed
{
if(l==r-1)
{
x1=ed,y1=ed;
x2=st,y2=st;
}
else if(i==l)
{
x1=i+1,y1=i+1;
x2=st,y2=ed;
}
else if(i+1==r)
{
x1=st,y1=ed;
x2=i,y2=i;
}
else
{
x1=st,y1=i+1;
x2=i,y2=ed;
}
//cout<<x1<<' '<<y1<<' '<<x2<<' '<<y2<<endl;
res=(res+dfs(l,i,x1,y1)*dfs(i+1,r,x2,y2)%mod*C(r-l-1,i-l))%mod;
}
return dp[l][r][st][ed]=res;
}
int main()
{
ll n;
scanf("%lld",&n);
init();
for(int i=1;i<=n;i++)
{
scanf("%lld",&p[i]);
p[i]++;
}
printf("%lld\n",dfs(1,n,1,n));
return 0;
}

京公网安备 11010502036488号