BZOJ 4517 [Sdoi2016]排列计数
2017-06-08 11:35
281 查看
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 行,每行一个数,表示求出的序列数Sample Input
51 0
1 1
5 2
100 50
10000 5000
Sample Output
01
20
578028887
60695423
HINT
Source
鸣谢Menci上传~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
组合数问题~
从中选出m个, c(n,m),剩下的要用到神奇的公式,学自SilverN,并不知道原理:
f[0]=1,f[1]=1,f[2]=1,f[i]=(i-1)*(f[i-1]+f[i-2]),序列的错序种类数。
#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long
#define mod 1000000007
int t,n,m,sheng[1000001],jiang[1000001],f[1000001];
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
while(ch>='0' && ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
int mi(int u,int v)
{
int now=1;
for(;v;v>>=1,u=(ll)u*u%mod) if(v&1) now=(ll)now*u%mod;
return now;
}
int c(int n,int m)
{
return (ll)sheng
*jiang[m]%mod*jiang[n-m]%mod;
}
int main()
{
sheng[0]=1;
for(int i=1;i<=1000000;i++) sheng[i]=(ll)sheng[i-1]*i%mod;
jiang[1000000]=mi(sheng[1000000],mod-2);
for(int i=1000000-1;~i;i--) jiang[i]=(ll)jiang[i+1]*(i+1)%mod;
f[0]=f[2]=1;
for(int i=3;i<=1000000;i++) f[i]=(ll)(i-1)*(f[i-1]+f[i-2])%mod;
t=read();
while(t--)
{
n=read();m=read();
if(n<m)
{
puts("0");continue;
}
printf("%d\n",(ll)c(n,m)*f[n-m]%mod);
}
return 0;
}
相关文章推荐
- bzoj 4517: [Sdoi2016]排列计数
- BZOJ 4517: [Sdoi2016]排列计数
- BZOJ 4517: [Sdoi2016]排列计数 组合数 错排公式
- bzoj 4517: [Sdoi2016]排列计数【容斥原理+组合数学】
- BZOJ_4517_[Sdoi2016]排列计数_组合数学
- BZOJ 4517 [Sdoi2016]排列计数
- BZOJ 4517|SDOI 2016|排列计数|错排公式|组合数学
- BZOJ 4517 [Sdoi2016]排列计数
- bzoj-4517 4517: [Sdoi2016]排列计数(组合数学)
- 【BZOJ】4517 [Sdoi2016]排列计数
- bzoj 4517: [Sdoi2016]排列计数 排列组合+动态规划
- bzoj 4517: [Sdoi2016]排列计数
- bzoj 4517: [Sdoi2016]排列计数
- bzoj 4517: [Sdoi2016]排列计数 递推
- BZOJ 4517: [Sdoi2016]排列计数 错排公式
- BZOJ.4517.[SDOI2016]排列计数(错位排列 逆元)
- Bzoj 4517: [Sdoi2016]排列计数(排列组合)
- [组合 错排] BZOJ 4517 [Sdoi2016]排列计数
- BZOJ 4517: [Sdoi2016]排列计数 错排+逆元
- bzoj 4517 [Sdoi2016]排列计数