牡牛和牝牛
有n个0或者1,进行全排列,要求任意两个0间至少有k个1,询问其方案数%5000011。
对于全部数据,对于全部数据,\(1≤N≤10^5,0≤K
解:
显然为排列组合问题,考虑方向自然为通项与递推方程。
法一(通项公式):
首先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\)。
思路一:
间隔是有序的,而1是无序的,不好解决无序的元素放入有序的盒子的方案数的问题,反过来看则是有序的间隔放入无序的1,而间隔可以多次放入同一个1,即可重组合,于是方案数不难得知为\(C_{n-(x-1)k}^{x}\)。
思路二:
插入组合问题很困难,于是考虑组合转排列,转换模型,即等价于\(n-(x-1)k-x\)个1与\(x\)个0进行全排列,根据可重排列公式有:
\[\frac{(n-(x-1)k-x+x)!}{x!(n-(x-1)k-x)!}=\frac{(n-(x-1)k)!}{x!(n-(x-1)k-x)!}\]
\[=C_{n-(x-1)k}^{x}\]
所以只要枚举0的个数,代入公式累加即可。
参考代码:
#include
#include
#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>=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
#include
#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;
}