bzoj 2839 集合计数 - 二项式反演

35 篇文章 0 订阅
18 篇文章 0 订阅

傻逼题

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define lint long long
#define mod 1000000007
#define N 1000010
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
int fac[N],facinv[N];
inline int mul(int x) { return x<0?x+mod:x; }
inline int sol(lint x,int s) { return x%=mod,((s&1)?(x?mod-x:0):x); }
inline int C(int n,int m) { return (lint)fac[n]*facinv[m]%mod*facinv[n-m]%mod; }
template<int p>inline int fast_pow(int x,int k,int ans=1)
{   for(;k;k>>=1,x=(lint)x*x%p) (k&1)?ans=(lint)ans*x%p:0;return ans; }
inline int prelude(int n)
{
    for(int i=fac[0]=1;i<=n;i++) fac[i]=(lint)fac[i-1]*i%mod;
    facinv[n]=fast_pow<mod>(fac[n],mod-2);
    for(int i=n-1;i>=0;i--) facinv[i]=facinv[i+1]*(i+1ll)%mod;
    return 0;
}
int main()
{
    int n,k,ans=0;scanf("%d%d",&n,&k);prelude(n);
    for(int i=k,fi;i<=n;i++)
        fi=(lint)C(n,i)*mul(fast_pow<mod>(2,fast_pow<mod-1>(2,n-i))-1)%mod,
        (ans+=sol((lint)C(i,k)*fi,i+k))%=mod;
    return !printf("%d\n",ans);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值