P4071 [SDOI2016]排列计数

错排问题
在这里插入图片描述
最开始以为是组合数取余结果样例兜不上去,才发现是错排问题

#include<bits/stdc++.h>
using namespace std;
#define ll long long
long long mod = 1e9+7;
ll up[1000010],d[1000010];
ll Pow(ll a,ll b)
{
	ll ans=1;
	while(b){
		if(b&1) ans=(ans*a)%mod;
		a=(a*a)%mod;
		b>>=1;
	}
	return ans;
}
ll C()
{
	//if(m==0) return 1;
	//if(m>n-m) m=n-m;
	//ll up=1,down=1;
	up[0]=1;
	for(ll i=1;i<=1e6;i++){
		up[i]=(up[i-1]*i)%mod;
		//down=(down*i)%p;
	}
	d[2]=1,d[3]=2;
	for(ll i=4;i<=1e6;i++)
	{
		d[i]=(i-1)*(d[i-1]+d[i-2])%mod;
	}
//	cout<<up<<" "<<down<<" ";
	//return up*Pow(up,p-2)%p;
}
ll c(ll n,ll m)
{
	if( m==0 )	return 1;
//	cout<<up[n]<<up[m]<<up[n-m]<<endl;
	ll bb=up[n]*Pow(up[m]*up[n-m]%mod,mod-2)%mod;
	return bb;
}
signed main()
{
	int t;
	cin>>t;
	C();
	while(t--)
	{
		ll n,m;
		cin>>n>>m;
		if(n==m)cout<<"1"<<endl;
		else cout<<c(n,m)*d[n-m]%mod<<endl;
	}
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值