51nod2717 小Biu的序列

2717 小Biu的序列

用N个不同的字符(编号1 - N),每种字符最多使用一次,组成一个字符串,有如下要求:
(1) 对于编号为i的字符,如果2 * i > n,则该字符可以作为结尾字符。如果不作为结尾字符而是中间的字符,则该字符后面可以接任意字符。
(2) 对于编号为i的字符,如果2 * i <= n,则该字符不可以作为结尾字符。作为中间字符,那么后面接的字符编号i一定要满足 2*i>n。
问有多少长度为M且符合条件的字符串,由于数据很大,只需要输出该数Mod 10^9 + 7的结果。

输入

输入2个数,N, M中间用空格分割,N为不同字符的数量,M为字符串的长度。(2 <= N, M <= 10^6)

输出

输出符合条件的字符串的数量。由于数据很大,只需要输出该数Mod 10^9 + 7的结果。

输入样例

3 3

输出样例

4

 解析:

      如果用k表示“前半段字符的个数”,那么n-k可以表示“后半段字符的个数”,下面我们考虑如何计数。

      设表示字符串时用到了t个“前半段字符”,先选出这t个字符,即C(k,t)。

      然后,考虑到 前半段字符后面必须是后半段字符 ,我们选出k个"后半段字符"任意放在刚才选出来的那t个”前半段字符“的后面,组成一组,即A(n-k,t)。

      这样,我们就有了t个组,然后再从剩下的”后半段字符“里挑出m-2*t个字符自成一组,即C(n-k-t,m-2*t),然后我们就有了m-t个组,总共m个字符。
      最后,我们再把这m-t个组任意排序,即A(m-t,m-t).
      最终式子为  ΣC(k,t)*A(n-k,t)*C(n-k-t,m-2*t)*A(m-t,m-t)。

放代码:

#include <bits/stdc++.h>
#define ll long long
const ll M=1e9+7;
using namespace std;
ll fac[1000006]={1},n,m, ans;
ll pw(ll x,int ti=M-2){
    ll ret=1;
    while(ti){
        if(ti&1)ret=ret*x%M;
        ti>>=1;
        x=x*x%M;
    }
    return ret;
}
ll A(int a, int b){
	if(a<0||b<0)return 0;
    return fac[a]*pw(fac[a-b])%M;
}
ll C(int a, int b){
    return A(a,b)*pw(fac[b])%M;
}
int main() {
	for(int i=1; i<=1e6;i++)fac[i]=fac[i-1]*i%M;
	cin>>n>>m;
	ll k=n>>1, lt=max(m-n+k,0ll),rt=min(k,m/2);
	for(int t=lt;t<=rt;t++)
	    ans+=C(k,t)*A(n-k,t)%M*C(n-k-t,m-2*t)%M*A(m-t,m-t)%M,ans%=M;
	cout<<ans;
	return 0;
}

上一篇:带你学习DSP系列一:走进TMS320F28377


下一篇:java课堂测试2第一阶段:方法运用