博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
3398: [Usaco2009 Feb]Bullcow 牡牛和牝牛
阅读量:5766 次
发布时间:2019-06-18

本文共 1478 字,大约阅读时间需要 4 分钟。

3398: [Usaco2009 Feb]Bullcow 牡牛和牝牛

Time Limit: 1 Sec  Memory Limit: 128 MB
Submit: 243  Solved: 167
[][][]

Description

    约翰要带N(1≤N≤100000)只牛去参加集会里的展示活动,这些牛可以是牡牛,也可以是牝牛.牛们要站成一排.但是牡牛是好斗的,为了避免牡牛闹出乱子,约翰决定任意两只牡牛之间至少要有K(O≤K<N)只牝牛.
    请计算一共有多少种排队的方法.所有牡牛可以看成是相同的,所有牝牛也一样.答案对5000011取模

Input

    一行,输入两个整数N和K.

Output

 
    一个整数,表示排队的方法数.

Sample Input

4 2

Sample Output

6
样例说明
6种方法分别是:牝牝牝牝,牡牝牝牝,牝牡牝牝,牝牝牡牝,牝牝牝牡,牡牝牝牡

HINT

 

Source

/*以下牡牛为a,牝牛为b。学完排列计数后试着来写这题,“至少”一词可以给我们提示,我们可以枚举a为x头(x>1),然后算出对应的排列累计起来。对于x头a,首先我们先缩掉必要的k头牛(x-1)*k,然后这时可以特判可以先结束(因为单调的),然后在缩好后的x个点和n-x-(x-1)*k个点进行多重排列就行了。逆元:a^(phi(n)-1) mod n */#include
#include
using namespace std;const int mod=5000011;typedef long long ll;int n,k;ll ans;ll fpow(ll a,ll p){ ll res=1; for(;p;p>>=1,a=a*a%mod) if(p&1) res=res*a%mod; return res;}ll C(ll n,ll m){ m=min(m,n-m);ll r1=1,r2=1; for(ll i=n-m+1;i<=n;i++) r1=r1*i%mod; for(ll i=1;i<=m;i++) r2=r2*i%mod; return r1*fpow(r2,mod-2);}int main(){ scanf("%d%d",&n,&k); for(ll i=0;i<=n;i++){ ll t=n-(i-1)*k; if(t

 

/*DP:设 f[i]表示取的最后一个数是i的方案数则 f[i]=siama(f[j]) i-j>k*/#include
#define mod 5000011using namespace std;int n,k,f[(int)1e5+5];int main(){ scanf("%d%d",&n,&k); f[0]=1; int sum=1,ans=1; for(int i=1;i<=n;i++){ if(i>k+1) sum=(sum+f[i-k-1])%mod; f[i]=sum; ans=(ans+f[i])%mod; } printf("%d\n",ans); return 0;}

 

转载于:https://www.cnblogs.com/shenben/p/6286551.html

你可能感兴趣的文章
Mongodb启动命令mongod参数说明
查看>>
TCP&UDP压力测试工具
查看>>
Android 最简单的自定义Dialog之一
查看>>
磨刀不误砍柴 - 配置适合工作学习的桌面环境
查看>>
redux v3.7.2源码解读与学习之 applyMiddleware
查看>>
【React】为什么我不再使用setState?
查看>>
Git原理与高级使用(3)
查看>>
从JDK源码看Writer
查看>>
Express 结合 Webpack 实现HMRwi
查看>>
基于protobuf的RPC实现
查看>>
我的友情链接
查看>>
HAProxy负载均衡原理及企业级实例部署haproxy集群
查看>>
开源中国动弹客户端实践(三)
查看>>
Win 8创造颠覆性体验:预览版关键更新
查看>>
vim在多文件中复制粘贴内容
查看>>
Android ContentObserver
查看>>
疯狂java学习笔记1002---非静态内部类
查看>>
ISA2006实战系列之一:实战ISA三种客户端部署方案(上)
查看>>
TCP服务器
查看>>
AC旁挂三层交换机管理ap,二层接入ap心得
查看>>