bzoj 4517: [Sdoi2016]排列计数【容斥原理+组合数学】

第一个一眼就A的容斥题!

这个显然是容斥的经典问题------错排,首先考虑没有固定的情况,设\( D_n \)为\( n \)个数字的错排方案数。

\[
D_n=n!-\sum_{t=1}^{n}(-1)^{t-1}\sum_{i_1<i_2<...<i_t}(n-t)!
\]

\[
D_n=n!+\sum_{t=1}^{n}(-1)^tC_{n}^{t}(n-t)!
\]

\[
D_n=n!+\sum(-1)^t\frac{n!}{t!}
\]

推到这一步就可以了,然后观察数据范围显然是要线性预处理,于是计算递推式:

\[
D_{(n+1)}=(n+1)!+\sum_{t=1}^{n+1}(-1)^t\frac{(n+1)!}{t!}
\]

\[
D_{(n+1)}=(n+1)!+(n+1)\sum_{t=1}^{n+1}(-1)^t\frac{n!}{t!}
\]

\[
D_{(n+1)}=(n+1)!+(n+1)(\sum_{t=1}^{n}(-1)^t\frac{n!}{t!}+(-1)^{n+1}\frac{n!}{(n+1)!})
\]

\[
D_{(n+1)}=(n+1)!+(n+1)\sum_{t=1}^{n}(-1)^t\frac{n!}{t!}+(-1)^{n+1}
\]

\[
D_{(n+1)}=(n+1)(n!+(n+1)\sum_{t=1}^{n}(-1)^t\frac{n!}{t!})+(-1)^{n+1}
\]

\[
D_{(n+1)}=(n+1)D_n+(-1)^{n+1}
\]

\[
D_i=i*D_{i-1}+(-1)^i
\]

然后考虑有\( m \)的限制,就相当于\( m \)个数字固定,剩下\( n-m \)个数字错排,直接从预处理的\( D \)里面查即可,最后乘上选出\( m \)个固定位的方案数,对组合数预处理阶乘、逆元。由此可得答案为:

\[
ans=D_{(n-m)}*C_{n}^{m}
\]

这东西推起来真刺激

#include<iostream>
#include<cstdio>
using namespace std;
const long long N=1000005,mod=1e9+7;
long long T,n,m,inv[N],fac[N],cp[N];
int read()
{
    int r=0;
    char p=getchar();
    while(p>‘9‘||p<‘0‘)
        p=getchar();
    while(p>=‘0‘&&p<=‘9‘)
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r;
}
long long ksm(long long a,long long b)
{
    long long r=1ll;
    while(b)
    {
        if(b&1)
            r=r*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return r;
}
long long C(long long n,long long m)
{
    return fac[n]*inv[n-m]%mod*inv[m]%mod;
}
int main()
{
    fac[0]=1;
    for(long long i=1;i<=N-5;i++)
        fac[i]=fac[i-1]*i%mod;
    inv[N-5]=ksm(fac[N-5],mod-2);
    for(long long i=N-6;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod;
    cp[0]=1;//这里的cp数组即是上文提到的D(cuopai 23333)
    for(long long i=1;i<=N-5;i++)
        cp[i]=(i*cp[i-1]+((i&1)?-1:1))%mod;
    T=read();
    while(T--)
    {
        n=read(),m=read();
        printf("%lld\n",(cp[n-m]*C(n,m)%mod+mod)%mod);
    }
    return 0;
}

原文地址:https://www.cnblogs.com/lokiii/p/8215044.html

时间: 2024-11-05 17:30:47

bzoj 4517: [Sdoi2016]排列计数【容斥原理+组合数学】的相关文章

BZOJ 4517: [Sdoi2016]排列计数 错排+逆元

4517: [Sdoi2016]排列计数 Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. Input 第一行一个数 T,表示有 T 组数据. 接下来 T 行,每行两个整数 n.m. T=500000,n≤1000000,m≤1000000 Output 输出 T 行,每行一个数,表示

数学(错排):BZOJ 4517: [Sdoi2016]排列计数

4517: [Sdoi2016]排列计数 Time Limit: 60 Sec  Memory Limit: 128 MBSubmit: 693  Solved: 434[Submit][Status][Discuss] Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. Input 第

BZOJ 4517: [Sdoi2016]排列计数

4517: [Sdoi2016]排列计数 Time Limit: 60 Sec  Memory Limit: 128 MBSubmit: 911  Solved: 566[Submit][Status][Discuss] Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. Input 第

bzoj-4517 4517: [Sdoi2016]排列计数(组合数学)

题目链接: 4517: [Sdoi2016]排列计数 Time Limit: 60 Sec  Memory Limit: 128 MBSubmit: 846  Solved: 530[Submit][Status][Discuss] Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. I

4517: [Sdoi2016]排列计数

Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. Solution 答案是:\(C(n,m)*D(n-m)\) \(D(n)\) 是长度为\(n\)的错排的方案数 \(D(n)=n!*(1-\frac{1}{1!}+\frac{1}{2!}-\frac{1}{3!}+(-1)^n\fr

[SDOI2016] 排列计数 (组合数学)

[SDOI2016]排列计数 题目描述 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7109+7 取模. 输入输出格式 输入格式: 第一行一个数 T,表示有 T 组数据. 接下来 T 行,每行两个整数 n.m. 输出格式: 输出 T 行,每行一个数,表示求出的序列数 输入输出样例 输入样例#1: 5 1 0 1 1

BZOJ 2111 Perm 排列计数(满二叉树)

题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=2111 题意:求1到n有多少种排列满足:A[i]>A[i/2](2<=i<=n). 思路:形式类似二叉树.建模之后其实就是n个节点的不同的满二叉树有多少种?用f[i]表示i个节点的满二叉树个数,则f[n]=f[L]*f[R]*C(n-1,L).其中L和R对于确定的n来说是确定的.比如n=10时,左右子树分别有6.3个点. i64 a[N],n,p,f[N]; void init(

【BZOJ4517】[Sdoi2016]排列计数 组合数+错排

[BZOJ4517][Sdoi2016]排列计数 Description 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. Input 第一行一个数 T,表示有 T 组数据. 接下来 T 行,每行两个整数 n.m. T=500000,n≤1000000,m≤1000000 Output 输出 T 行,每行一个

[Sdoi2016]排列计数

问题 A: [Sdoi2016]排列计数 时间限制: 3 Sec  内存限制: 512 MB 题目描述 求有多少种长度为 n 的序列 A,满足以下条件: 1 ~ n 这 n 个数在序列中各出现了一次 若第 i 个数 A[i] 的值为 i,则称 i 是稳定的.序列恰好有 m 个数是稳定的 满足条件的序列可能很多,序列数对 10^9+7 取模. 输入 第一行一个数 T,表示有 T 组数据. 接下来 T 行,每行两个整数 n.m. T=500000,n≤1000000,m≤1000000 输出 输出