bzoj4402: Claris的剑（组合数学）

最新推荐文章于 2020-01-13 20:23:43 发布

SC.ldxcaicai

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量145

点赞数 1

分类专栏： # 组合数学

本文链接：https://blog.csdn.net/dreaming__ldx/article/details/96039332

版权

组合数学专栏收录该内容

35 篇文章 0 订阅

订阅专栏

传送门
思路：
考虑怎么才不会算重。
显然你的序列可以变化成： $1\ 2\ 1\ 2....x\ x+1\ x\ x+1...max$ 和
$1\ 2\ 1\ 2....x\ x+1\ x\ x+1...max\ max-1$
中的任意一种。
证明显然。
于是我们枚举 $m a x$ 是多少就成了简单组合计数问题。
代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int N=2e6+5;
typedef long long ll;
const int mod=1e9+7;
inline int add(const int&a,const int&b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(const int&a,const int&b){return a>=b?a-b:a-b+mod;}
inline int mul(const int&a,const int&b){return (ll)a*b%mod;}
inline void Add(int&a,const int&b){a=a+b>=mod?a+b-mod:a+b;}
inline void Dec(int&a,const int&b){a=a>=b?a-b:a-b+mod;}
inline void Mul(int&a,const int&b){a=(ll)a*b%mod;}
inline int ksm(int a,int p){int ret=1;for(;p;p>>=1,a=mul(a,a))if(p&1)Mul(ret,a);return ret;}
int fac[N<<1],ifac[N<<1],n,m;
inline void init(){
    fac[0]=fac[1]=ifac[0]=ifac[1]=1;
    for(ri i=2;i<=n+m;++i)fac[i]=mul(i,fac[i-1]),ifac[i]=mul(ifac[mod-mod/i*i],mod-mod/i);
    for(ri i=2;i<=n+m;++i)Mul(ifac[i],ifac[i-1]);
}
inline int C(int n,int m){return (ll)fac[n]*ifac[m]%mod*ifac[n-m]%mod;}
int main(){
    cin>>n>>m;
    init();
    int ans=1;
    for(ri a,b,i=2;i<=m;++i){
        if(n<i)break;
        a=(n-i)>>1,b=i-1;
        Add(ans,C(a+b,b));
        if(n<i+1)break;
        a=(n-i-1)>>1,b=i-1;
        Add(ans,C(a+b,b));
    }
    cout<<ans;
    return 0;
}