BZOJ4517

https://www.lydsy.com/JudgeOnline/problem.php?id=4517

这题竟然60s

我日

不过我的代码带上读入挂跑了14s就出来了

题目很简单

预处理加错排问题公式

其实完全可以跑的更快

把inv提前预处理就行了

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1e6,M=1e9+7;
ll f[N+10],JC[N+10];
ll quickmod(ll a,ll b,ll m){ll ans=1;while(b)
{if(b&1)ans=(ans*a)%m;b>>=1;a=a*a%m;}return ans;}

ll inv(ll a,ll p){return quickmod(a, p-2, p);}
inline void read(ll &x){char ch;bool ok;
for(ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if(ch=='-') ok=1;
for(x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());if(ok) x=-x;}
ll C(ll a,ll b,ll m){
    return JC[a]*inv(JC[b],m)%m
            *inv(JC[a-b],m)%m;
}
int main(){
    ll T;
    read(T);
    ll ans,n,m;
    f[0]=0,f[1]=0,f[2]=1;
    for(int i=3;i<=N;++i){
        f[i]=(f[i-1]+f[i-2])*(i-1)%M;
    }
    JC[0]=1;
    for(ll i=1;i<=N;++i){
        JC[i]=JC[i-1]*i%M;
    }
    while(T--){
        read(n),read(m);
        if(n==m)
            ans=1;
        else
            ans=f[n-m]*C(n,m,M)%M;
        printf("%lld\n", ans);
    }
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值