交换题解

最新推荐文章于 2023-11-11 15:01:36 发布

LV24twx

最新推荐文章于 2023-11-11 15:01:36 发布

阅读量490

点赞数

分类专栏：日常刷题记录

本文链接：https://blog.csdn.net/weixin_44351585/article/details/103532004

版权

日常刷题记录专栏收录该内容

8 篇文章 0 订阅

订阅专栏

先考虑暴力，暴力枚举 $q$ 序列，然后代入统计答案。

时间复杂度： $O(2^n\times n)$ 即 $O\ \left(不可过\right)$

考虑倒推法，求能把 $p$ 变成 $\{0,1,\cdots,n-1\}$ 的 $q$ 序列。

如果第一次交换了 $p [i]$ 和 $p [i + 1]$ ，那么容易得到 $p [i]$ 和 $p [i + 1]$ 不会再次交换。

所以 $p[0,\cdots,i]$ 和 $p[i+1,\cdots,n-1]$ 会被 $(i, i + 1)$ 这次交换独立出来。

$p[0,\cdots,i]$ 和 $p[i+1,\cdots,n-1]$ 是两个规模更小的子问题。

设 $a s d [n] [s t]$ 表示 $n$ 个数，由 $\{st,st+1,\cdots,st+n-1\}$ 组成的排列 $H$ 为 $p$ 的一个子序列，将 $H$ 排成从左到右递增的方案数。

枚举 $i$ ，表示第一次交换 $H [i - 1]$ 和 $H [i]$ 。

1. $H$ 会被 $(i - 1, i)$ 这次交换独立成两个区间 $H[0,\cdots,i-1]$ 和 $H[i,\cdots,n-1]$

2.注意交换 $H [i - 1]$ 和 $H [i]$ 之后只有 $H[0,\cdots,i-1]$ 包含了 $[s t, s t + i - 1]$ 内的数且 $H[i,\cdots,n-1]$ 包含了 $[s t + i, s t + n - 1]$ 内的数，第一次交换 $(i - 1, i)$ 才合法。即注意判断合法的 $i$ 。

3.左区间 $H[0,\cdots,i-1]$ 是一个 $s t$ 到 $s t + i - 1$ 的序列，方案数为 $a s d [i] [s t]$ 。

4.右区间 $H[i,\cdots,n-1]$ 是一个 $s t + i$ 到 $s t + n - 1$ 的序列，方案数为 $a s d [n - i] [s t + i]$ .

5.在除第一次交换后剩下的 $n - 2$ 次交换中选出 $i - 1$ 次交换分配给 $H[0,\cdots,i-1]$ ，剩下的 $n - i - 1$ 次交换分配给 $H[i,\cdots,n-1]$ ，方案数就是 $C_{n-2}^{i-1}$ 。

6.转移方程为 $asd[n][st]+=asd[i][st]\times asd[n-i][st+i]\times C(n-2,i-1)$ (交换 $H [i - 1]$ 和 $H [i]$ 也合法)

由于需要 $O (n)$ 的时间来判断交换是否合法，所以时间复杂度为 $O(n^4)$

#include<bits/stdc++.h>
using namespace std;
#define f1(a,b,c) for(ll c=a;c<=b;c++)
#define f2(a,b,c) for(ll c=a;c>=b;c--)
#define f3(a,b,c) for(ll c=a;c;c=b)
#define so1(a,n) sort(a+1,a+n+1,mycmp);
#define so2(a,n) sort(a+1,a+n+1);
#define ll long long
#define itn ll
#define ubt ll 
#define mp make_pair
#define pii pair<ll,ll>
#define pll pair<ll,ll>
const ll twx=55;
const ll inf=0x3f3f3f3f;
const ll MOD=1e9+7;
ll read()
{
    ll sum=0;
    ll flag=1;
    char c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-')
        {
            flag=-1;
        }
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        sum=((sum*10)+c-'0');
        c=getchar();
    }
    return sum*flag;
}
ll n;
ll p[twx];
ll qu[twx][twx][twx];
ll asd[twx][twx];
ll fac[twx];
ll inv[twx];
ll add(ll x,ll y)
{
    return x+y>=MOD?x+y-MOD:x+y;
}
ll sub(ll x,ll y)
{
    return x-y<0?x-y+MOD:x-y;
}
ll mul(ll x,ll y)
{
    return 1LL*x*y%MOD;
}
ll C(ll a, ll b)
{
    return 1ll*fac[a]*inv[b]%MOD*inv[a-b]%MOD;
}
void pre()
{
    fac[0]=1;
    f1(1,n,i)
    {
        fac[i]=mul(fac[i-1],i);
    }
    inv[0]=inv[1]=1;
    f1(2,n,i)
    {
        inv[i]=MOD-(ll)MOD/i*inv[MOD%i]%MOD;
    }
    f1(2,n,i)
    {
        inv[i]=(ll)inv[i]*inv[i-1]%MOD;
    }
}
void init()
{
	n=read();
    f1(0,n-1,i)
    {
        p[i]=read();
    }
    memset(asd,-1,sizeof asd);
}
void alize()
{
    f1(0,n-1,i)
    {
        f1(i,n-1,j)
        {
            ll tot=0;
            f1(0,n-1,k)
            {
                if(i<=p[k]&&p[k]<=j)
                {
                    tot++;
                    qu[i][j][i+tot-1]=p[k];
                }
            }
        }
    }
}
bool check(ll l,ll r,ll mid)
{
    f1(l,mid-1,i)
    {
        if(l>qu[l][r][i]||qu[l][r][i]>mid)
        {
            return 0;
        }
    }
    if(l>qu[l][r][mid+1]||qu[l][r][mid+1]>mid)
    {
        return 0;
    }
    f1(mid+2,r,i)
    {
        if(mid+1>qu[l][r][i]||qu[l][r][i]>r)
        {
            return 0;
        }
    }
    if(mid+1>qu[l][r][mid]||qu[l][r][mid]>r)
    {
        return 0;
    }
    return 1;
}
ll dp(ll nn,ll st)
{
    if(asd[nn][st]!=-1)
    {
        return asd[nn][st];
    }
    if(nn==1)
    {
        return asd[nn][st]=1;
    }
    ll res=0;
    f1(1,nn-1,i)
    {
        if(check(st,st+nn-1,st+i-1))
        {
            res=add(res,mul(mul(dp(i,st),dp(nn-i,st+i)),C(nn-2,i-1)));
		}
    }
    return asd[nn][st]=res;
}
void print()
{
    printf("%lld\n",dp(n,0));
}
int main()
{
    init();
    alize();
    pre();
    print();
	return 0;
}