[2018.11.03 T2] 执行任务

暂无链接

执行任务

有一支部队,由 nn 个成员组成,编号分别为 1,2,,n1,2,…,n。现在要求选出若干人去执行一项 任务,由于任务很艰巨,现在要求选出的人数不少于 mm

现在要求你求出有多少种选人的方案。

输入:

一行两个正整数 nmn,m

输出:

一行一个正整数 xx,表示答案。由于答案可能很大,输出答案对 10000000071000000007 取模的结果。

样例输入:

6 3

样例输出:

42

数据范围:

对于 20%20\%的数据,1mn201≤m≤n≤20

对于 40%40\%的数据, 1mn10001≤m≤n≤1000

对于 60%60\%的数据,1mn1000001≤m≤n≤100000

对于额外 10%10\%的数据,m=1,1n1000000000m=1,1≤n≤1000000000

对于 100%100\%的数据,mn,1m100000,1n1000000000m≤n,1≤m≤100000,1≤n≤1000000000

题解

数学题:
i=mn(ni)=i=1n(ni)i=1m1(ni)=2n1i=1m1(ni) \begin{aligned} &\sum_{i=m}^n\binom{n}{i}\\ =&\sum_{i=1}^n\binom{n}{i}-\sum_{i=1}^{m-1}\binom{n}{i}\\ =&2^n-1-\sum_{i=1}^{m-1}\binom{n}{i} \end{aligned}
因为(ni)=n×(n1)××(ni+1)1×2××i\binom{n}{i}=\frac{n\times(n-1)\times\cdots\times(n-i+1)}{1\times 2\times\cdots\times i},意思就是说分号上面和下面的项数不会超过mm项,是可以预处理的,此题完结。。。

代码
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=1e5+5,mod=1e9+7;
ll fac[M],inv[M],sum;
int n,m;
ll power(ll x,ll p){ll r=1;for(;p;p>>=1,x=x*x%mod)if(p&1)r=r*x%mod;return r;}
ll C(int n,int m){return fac[m]*inv[m]%mod;}
void in(){scanf("%d%d",&n,&m);}
void ac()
{
	fac[0]=1;for(int i=1;i<=m;++i)fac[i]=fac[i-1]*i%mod;
	inv[m]=power(fac[m],mod-2);
	for(int i=m-1;i>=0;--i)inv[i]=inv[i+1]*(i+1)%mod;
	fac[1]=n;for(int i=2;i<=m;++i)fac[i]=fac[i-1]*(n-i+1)%mod;
	for(int i=1;i<m;++i)(sum+=C(n,i))%=mod;
	printf("%lld",(power(2,n)-1-sum+mod)%mod);
}
int main(){in(),ac();}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章