【51nod 1251】 Fox序列的数量(以及带限制插板法讲解)
为什么网上没有篇详细的题解【雾
可能各位聚聚觉得这道题太简单了吧 /kk
题意
首先题目是求满足条件的序列个数,条件为:出现次数最多的数仅有一个
分析
感谢 刚睡醒的 JZ姐姐在咱写题解忽然陷入自闭的时候为咱知道迷津 QwQ
那么套路来说第一想到的其实就是容斥辣 然后懒得想去网上黈力【雾
在然后,发现网上要么没推导公式(虽说是直接列不用化简的)要么就是推导公式是错的...
于是坐下来冷静分析:首先,我们考虑最高的出现次数为 x 的情况下,有多少个序列是满足条件的
容斥一发发现式子长这个样子(用到了隔板法,然后容斥顺便也能把有数字出现次数高于 x 的情况给容掉)
\]
大概是说,我们确定最高次数为 x ,那么我们先钦定一个数为恰好出现 x 次的数,然后剩下 m-1 个数插板容斥求方案 QwQ ,于是我们就可以递推求解辣:
\]
然后我们发现复杂度爆炸了呢 Σ(⊙▽⊙"a 管他呐,交上去就好辣
然后咱发现后面那个式子如果 \(n-i*(j+1)+m-2\) 大于 \(m-2\) 的话贡献为 0 ,一顿计算发现 j 要小于 n/i 才有贡献,那么第二个枚举的上界的 m-1 改一改就好咯...
\]
然后调和级数分析一波,发现复杂度是 \(n·H(n)\) 的,能过...
顺便稍微讲讲带限制的板怎么插好了 QAQ
注意这里说的限制是对于任意两个隔板之间的元素数量而言,如果隔板间元素数量的限制是不同,可能需要什么黑科技(反正咱只会指数级别容斥【逃)
那么假设现在有 n 个元素,要分成 m 堆元素,也就是有 m-1 个插板
0.不带限制的插板
显然是 \(\binom{n+m-1}{m-1}\)
意思就是我们把 n 个元素和 m-1 个插板排成一列,计算方案数
那么我们只要计算 m-1 个插板有多少种摆法就好了,隔板确定完了,元素的放法自然就确定了
说白了就是求了个 m-1 个插板在所有 n+m-1 个物品中的位置方案数
1.带下界的插板
假设下界为 x ,即隔出来的每堆元素数量不小于 x
答案是 \(\binom{n+m-m·x-1}{m-1}\)
我们发现这玩意儿比较简单,只要把原本要拿来插板的元素先减去 \(m*x\) (也就是预先给每堆元素分掉这个下界 ) 然后再去插板就好辣
2.(JZ 姐姐教的QwQ)带上界的插板
假设上界为 x ,即隔出来的每堆元素数量不大于 x
答案是 \(\sum_{i=0}^{m} (-1)^i \binom{m}{j} \binom{n+m-i*(x+1)-1}{m-1}\)
有容斥呢 QwQ 解释一下,就是说我们枚举超出上限的元素堆个数,然后强制他们选了 (x+1) 个,这样不管它们接下来怎么选都已经超限了,然后我们加个 \((-1)\) 的幂次容斥一下就好辣
如果康不懂的话请重修容斥 QvQ
这样咱发现上面的那个式子一下子就很明了了呢 QvQ
Code
注意阶乘及逆元的 1e5 数组是假的,要两倍的!花了咱 5 点头盾 get 到的细节 【枯
//by Judge
#define HGS_AK_IOI true
#include<bits/stdc++.h>
#define mod 1000000007
#define Rg register
#define fp(i,a,b) for(Rg int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(Rg int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(Rg int i=head[u],v=e[i].to;i;v=e[i=e[i].nxt].to)
#define open(S) freopen(S".in","r",stdin),freopen(S".out","w",stdout)
#define ll long long
using namespace std;
const int M=2e5+3;
typedef int arr[M];
#ifndef Judge
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline T Min(T x,T y){return x<y?x:y;}
template<class T>inline T Max(T x,T y){return x>y?x:y;}
template<class T>inline bool cmax(T& a,T b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T& a,T b){return a>b?a=b,1:0;}
inline int read(){ int x=0,f=1; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;
} int n,m,ans; arr fac,inv;
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline void Pls(int& x,int y){if((x+=y)>=mod)x-=mod;}
inline int qpow(int x,int p=mod-2){ Rg int s=1;
for(;p;p>>=1,x=mul(x,x)) if(p&1) s=mul(s,x); return s;
}
inline void prep(int n=2e5){
fac[0]=1; fp(i,1,n) fac[i]=mul(fac[i-1],i);
inv[n]=qpow(fac[n]); fd(i,n,1) inv[i-1]=mul(inv[i],i);
}
inline int Ceil(int x,int y){ return x%y>0?x/y+1:x/y; }
inline int C(int n,int m){ if(n<m) return 0;
return mul(fac[n],mul(inv[m],inv[n-m]));
}
inline void Solv(){ n=read(),m=read(),ans=0;
if(m==1) return printf("1\n"),void();
if(n==1) return printf("%d\n",m),void();
fp(i,Ceil(n,m),n) fp(j,0,Min(n/i-1,m-1))
Pls(ans,mul(j&1?(mod-1):1, mul(C(m-1,j),C(n+m-i*(j+1)-2,m-2)) ));
printf("%d\n",mul(ans,m));
}
int main(){
#ifdef Judge
freopen("1.in","r",stdin);
#endif
int T=read(); prep();
while(T--) Solv(); return 0;
}
【51nod 1251】 Fox序列的数量(以及带限制插板法讲解)的更多相关文章
- 51nod 1251 Fox序列的数量 (容斥)
枚举最多数字的出现次数$k$, 考虑其他数字的分配情况. 对至少$x$种数出现$\ge k$次的方案容斥, 有 $\sum (-1)^x\binom{m-1}{x}\binom{n-(x+1)k+m- ...
- 【51nod】1251 Fox序列的数量
题解 容斥题 我们枚举出现次数最多的数出现了K次 然后我们需要计算的序列是所有数字出现个数都不超过K - 1次 我们枚举不合法的数字的数目j,说明这个排列里除了我们固定出现K次的数至少有j个数是不合法 ...
- (转)AS3正则:元子符,元序列,标志,数量表达符
(转)AS3正则:元子符,元序列,标志,数量表达符: AS3正则:元子符,元序列,标志,数量表达符 七月 4th, 2010 归类于 AS3前端技术 作者Linkjun 进行评论 as3正则:元子符, ...
- Java实现蓝桥杯模拟正整数序列的数量
问题描述 小明想知道,满足以下条件的正整数序列的数量: 1. 第一项为 n: 2. 第二项不超过 n: 3. 从第三项开始,每一项小于前两项的差的绝对值. 请计算,对于给定的 n,有多少种满足条件的序 ...
- Comparing Your Heros拓扑序列的数量
给出N行英雄的比较,每一行包含两个英雄的名字,代表第一个英雄比第二个英雄更受欢迎. 英雄的数目不超过16个.问有多少种可能的受欢迎程度的序列满足N行英雄的比较. 由于只有英雄数目不超过16个,可以用二 ...
- 51nod 1042数字0-9的数量
1042 数字0-9的数量 基准时间限制:1 秒 空间限制:131072 KB 分值: 10 难度:2级算法题 收藏 关注 给出一段区间a-b,统计这个区间内0-9出现的次数. 比如 10-19 ...
- 51nod 1042 数字0-9的数量 数位dp
1042 数字0-9的数量 基准时间限制:1 秒 空间限制:131072 KB 分值: 10 难度:2级算法题 收藏 关注 给出一段区间a-b,统计这个区间内0-9出现的次数. 比如 10-1 ...
- 【BZOJ】【1251】序列终结者
Splay 还是splay序列维护,这题我WA了的原因是:在Push_up的时候,当前子树的max我是直接取的L.R和v[x]的最大值,但是如果没有左/右儿子,默认是会访问0号结点的mx值,而这个值没 ...
- 51nod 1522 上下序列
题目描述 现在有1到n的整数,每一种有两个.要求把他们排在一排,排成一个2*n长度的序列,排列的要求是从左到右看,先是不降,然后是不升. 特别的,也可以只由不降序列,或者不升序列构成. 例如,下面这些 ...
随机推荐
- spring-boot-devtools 热部署
一.简单介绍 spring-boot-devtools是一个为开发者服务的一个模块,其中最重要的功能就是修改代码后自动启动springboot服务,速度比手动停止后再启动要快,节省出来的并不是手工操作 ...
- 快学Scala 第十三课 (类型层级,对象相等性)
Scala 类型层级: 对象相等性: 和Java一样要重写equals方法和hashcode方法 class Student(val id: Int, val name: String) { over ...
- 快学Scala 第十二课 (抽象类, 抽象字段, 提前定义)
抽象类: Scala 抽象类中,抽象方法不需要使用abstract. 在子类中重写超类抽象方法时,不需要使用override. abstract class Person { def say(s: S ...
- MongoDB 学习笔记之 DBRef
DBRef: MongoDB建模有两种方式,一种是内嵌(Embed),另一种是连接(Link).内嵌比较好理解,就是字段内容是个数组,数组内再包含文档,而我们今天介绍的是另一种,称为链接DBRef.由 ...
- 《构建之法》MSF&需求分析
第七章 MSF MSF基本原则 推动信息共享与沟通 为共同的远景而工作 充分授权和信任 各司其职,对项目共同负责 交付增量的价值 保持敏捷,预期和适应变化 投资质量 学习所有的经验 与顾客合作 MSF ...
- 第二篇 python进阶
目录 第二篇 python进阶 一 数字类型内置方法 二 字符串类型内置方法 三 列表类型内置方法(list) 四 元组类型内置方法(tuple) 五 字典内置方法 六 集合类型内置方法(self) ...
- 从键盘录入输入3 个数num1,num2,num3,按从大到小进行输出
本题分别使用嵌套分支语句以及三目运算符来实现,两种方法,可以对比看看. import java.util.Scanner; /** * 从键盘录入输入3 个数a,b,c,按从大到小进行输出 * @au ...
- 微信小程序路由跳转
微信小程序路由跳转 1.wx.switchTab(Object object) 这里的tabBar是底下的导航栏指定的页面, 跳转到 tabBar 页面,并关闭其他所有非 tabBar 页面 参数 O ...
- unittest介绍
unittest框架是python中一个标准的库中的一个模块,该模块包括许多的类如 test case类.test suit类.texttest runner类.texttest resuite类.t ...
- 渗透测试-基于白名单执行payload--Cmstp
0x01 Cmstp简介 Cmstp安装或删除“连接管理器”服务配置文件.如果不含可选参数的情况下使用,则 cmstp 会使用对应于操作系统和用户的权限的默认设置来安装服务配置文件. 微软官方文档: ...