牡牛和牝牛

牡牛和牝牛

有n个0或者1，进行全排列，要求任意两个0间至少有k个1,询问其方案数%5000011。
对于全部数据，对于全部数据，\(1≤N≤10^5,0≤K<N\)。

解：

显然为排列组合问题，考虑方向自然为通项与递推方程。

法一（通项公式）：

首先0决定了1的摆放，其次数据范围支持对0的枚举，于是枚举0的个数，设其x，于是每个间隔至少要有k个1,不妨先构造让其满足条件，需要(x-1)k个1,显然有\((x-1)k+x\leq n\Rightarrow [\frac{n+k}{k+1}]\),于是问题即变成对于剩下的1的插入0之间的间隔，
共有x+1个间隔,剩下的1的个数为\(n-(x-1)k-x\)。

思路一：

思路二：

参考代码：

#include <iostream>
#include <cstdio>
#define il inline
#define ri register
#define ll long long
#define yyb 5000011
using namespace std;
ll jc[100001],jv[100001],lsy(1);
il ll pow(ll,ll),c(int,int);
int main(){
    int n,i,j,k,li;scanf("%d%d",&n,&k);
    for(i=jc[0]=1;i<=n;++i)jc[i]=jc[i-1]*i%yyb;
    jv[n]=pow(jc[n],yyb-2),li=(n+k)/(k+1);
    for(i=n,jv[0]=1;i>1;--i)jv[i-1]=jv[i]*i%yyb;
    for(i=1;i<=li;++i)(lsy+=c(n-(i-1)*k,i))%=yyb;
    printf("%lld",lsy);
    return 0;
}
il ll c(int n,int r){
    if(n<r)return 0;
    return jc[n]*jv[r]%yyb*jv[n-r]%yyb;
}
il ll pow(ll x,ll y){
    ll ans(1);
    while(y){
        if(y&1)ans=ans*x%yyb;
        x=x*x%yyb,y>>=1;
    }return ans;
}

法二（递推方程）

经验告诉我们以序列长度为状态，于是设\(f[i]\)表示填到第i个位置的方案数，显然策略为填0或者1,填1恒满足累加\(f[i-1]\)，填0导致前k个都不能填0,故累加\(f[i-k-1]\),于是有:
\[f[i]=f[i-1]+f[i-k-1]\]

参考代码：

#include <iostream>
#include <cstdio>
#define il inline
#define ri register
#define yyb 5000011
using namespace std;
int dp[100001];
int main(){
    int n,k,i;
    scanf("%d%d",&n,&k);
    for(i=1;i<=k+1;++i)dp[i]=i+1;
    for(i=k+2;i<=n;++i)
        dp[i]=(dp[i-k-1]+dp[i-1])%yyb;
    printf("%d",dp[n]);
    return 0;
}

原文地址：https://www.cnblogs.com/a1b3c7d9/p/10780306.html