[Luogu4609][FJOI2016]建筑师

Posted zhoushuyu

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[Luogu4609][FJOI2016]建筑师相关的知识,希望对你有一定的参考价值。

luogu

description

一个(1...n)的排列,其前缀最大值有(A)个,后缀最大值有(B)个,求满足要求的排列数。
一个位置(i)满足前缀最大当且仅当不存在(j<i)使得(a_j>a_i)。后缀最大亦然。
(Tle2 imes10^5,nle5 imes10^4,A,Ble100)

sol

考虑(n)这个数。它一定既是前缀最大又是后缀最大,所以它的前面还有(A-1)个前缀最大,(B-1)个后缀最大。
考虑每个前缀/后缀最大值,它一定挡住了若干个数(可能是零个)使得它们不能成为前缀最大。不妨将其视作一个整体,其中最大的那个数是前缀/后缀最大值,且位置在最左/右边,其余的数可以随意排列。
这等价于把除了(n)以外的(n-1)个数分成(A+B-2)个环排列。而这就是第一类斯特林数的定义。
所以答案就是(s(n-1,A+B-2) imesinom{A+B-2}{A-1})
(O(n(A+B)))预处理,然后(O(1))回答。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int gi(){
    int x=0,w=1;char ch=getchar();
    while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
    if (ch=='-') w=0,ch=getchar();
    while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return w?x:-x;
}
const int N = 50005;
const int M = 205;
const int mod = 1e9+7;
int S[N][M],C[M][M];
int main(){
    S[0][0]=1;
    for (int i=1;i<N;++i)
        for (int j=1;j<=i&&j<M;++j)
            S[i][j]=(1ll*S[i-1][j]*(i-1)+S[i-1][j-1])%mod;
    C[0][0]=1;
    for (int i=1;i<M;++i){
        C[i][0]=1;
        for (int j=1;j<=i;++j)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
    }
    int Case=gi();while (Case--){
        int n=gi(),a=gi(),b=gi();
        printf("%lld
",1ll*S[n-1][a+b-2]*C[a+b-2][a-1]%mod);
    }
    return 0;
}

以上是关于[Luogu4609][FJOI2016]建筑师的主要内容,如果未能解决你的问题,请参考以下文章

P4609 [FJOI2016]建筑师

[FJOI2016]建筑师

[FJOI2016]建筑师(斯特林数)

[FJOI2016] 建筑师

loj2173 「FJOI2016」建筑师

FJOI2016 建筑师 和 CF960G Bandit Blues