女神(goddess)——组合数学
出自某模拟赛。
题目大意:

对1e9+7取模。
数据范围
20 % : n<=300
40 % : n<=2,000
50 % : n<=10,000
70 % : n<=1,000,000
100 % : n<=1,000,000,000
题解
方法众多。
然而我太蒟了,2h43min愣是写了一个40pts暴力走人。
开始指着想正解,但是发现想不出来,然后写n^2暴力,结果总是要么漏算要么算重要么取模出错。。。最后才调出来。
1.如果比较菜,请尝试打表:
$n^3$暴力显然。然后打表。
$n=1,ans_1=1$
$n=2,ans_2=6$
$n=3,ans_3=24$
$n=4,ans_4=80$
然后,套路的先观察相邻两项的关系。
然后尝试和项数的下标套上关系。
$ans_2=ans_1\times \frac{6}{1}$
$ans_3=ans_2\times \frac{8}{2}$
$ans_4=ans_3\times \frac{10}{3}$
诶,然后发现了规律!!
可以递推。
$f_{n+1}=f_n\times \frac{4+2\times n}{n}$
$f_{n+1}=f_n\times \frac{2\times(n+2)}{n}$
迭代下去得:
$f_{n+1}=f_1 \frac{2^n\times (n+2)!}{2n!}$
$f_{n+1}=\frac{2^{n}\times (n+1)\times(n+2)}{2}$
$f_{n}=\frac{2^{n-1}\times n\times(n+1)}{2}$
即可出结果
2.正解:
有点意思的是,i+n-i-1=n-1,k-j+j-1=k-1对于任意的i,j恒成立。
这就是突破口
考虑组合数的意义。
$\sum_{k=1}^n(k\times\sum_{j=1}^k\sum_{i=0}^{n-1}(C_i^{k-j}\times C_{n-i-1}^{j-1})$
$=\sum_{k=1}^n(k\times\sum_{i=0}^{n-1}\sum_{j=1}^k(C_i^{k-j}\times C_{n-i-1}^{j-1})$
j和i换了位置之后,
发现,其实i就是枚举的一个分割点,
然后对于选择的k-1个数,在1~i个数中选择k-j个,
在i+1~n-i-1个数中选择j-1个。
好像和$C_{n-1}^{k-1}$有些关系。
发现,对于$C_{n-1}^{k-1}$中的每个方案。
其实都可以找出0~n-1这n个分界点,然后统计一次。
每个方案被统计了n次。
所以,
原式
$=\sum_{k=1}^nk\times n\times C_{n-1}^{k-1}$
已经可以O(n)递推了。
我们可以用刚才的打表中方法,搞出递推式,然后迭代出来通项公式。
即可O(logn)求解。
3.但是这个还不够漂亮!!!
这个可是组合数啊!!不是一般的数。
组合数毕竟有实际的意义。
观察这个式子的组合意义。
$\sum_{k=1}^nk\times n\times C_{n-1}^{k-1}$
这个k-1和k,n有点麻烦。
提出来:
$=n\times (\sum_{k=1}^n(k-1)\times C_{n-1}^{k-1})+n\times 2^{n-1}$
第一个括号里面是什么意义?
对于n-1个数的集合中,所有子集的大小的和。
套路地,我们转化研究对象。
考虑每个元素被统计了几次。
就是:$2^{n-1-1}=2^{n-2}$
因为每个数自己必须出现一次,然后其他的数爱出现不出现。
所以,
$=n\times( n\times2^{n-2})+n\times 2^{n-1}$
然后就可以O(logn)计算了。
总结:
0.这个式子,我们尝试用数学知识、组合数公式化简,发现不容易化简。然后就要考虑组合数的意义。
1.组合数是一个有意义的数。这样的数学式子的推导,可以通过寻找式子的意义来进行化简。
往往起到立竿见影的效果。
2.打表找规律,要考虑把结果,递推关系和项的编号放在一起。
女神(goddess)——组合数学的更多相关文章
- JDK1.8新特性之Optional
概念 Optional 是JDK1.8中出现的一个容器类,代表一个值存在或者不存在.原来使用null表示一个值不存在,现在Optional可以更好的表达这个概念.并且可以避免空指针异常. 场景分析 需 ...
- FJNU 1154 Fat Brother And His Love(胖哥与女神)
FJNU 1154 Fat Brother And His Love(胖哥与女神) Time Limit: 2000MS Memory Limit: 257792K [Description] [ ...
- [ACdream] 女神教你字符串——三个气球
Problem Description 女神邀请众ACdream开联欢会,显然作为ACM的佼佼者,气球是不能少的~.女神准备了三种颜色的气球,红色,黄色,绿色(交通信号灯?) 有气球还不能满足女神,女 ...
- kali linux之窥看女神上网隐私(ettercap+wireshark+zenmap +dsniff)
作者:小波 http://www.cnblogs.com/xiaobo-Linux/ 无聊就玩了玩,不要干一些坏事哟~~网上也有一些文章关于kali linux的,就实战了一番.kali是用的debi ...
- 海洋女神建新installshield交流群了,原来的老群都满了,请加新群哦,记得认真填写验证信息
群号511751143 海洋女神installshield群
- 南邮CTF隐写之丘比龙的女神
刚开始下载下图片来 习惯性的binwalk一下 没发现东西 formost一下也没分离出来 扔进c32asm中发现有nvshen.jpg 于是改后缀名字为.zip 解压nvshen.jpg发现无法解压 ...
- 找女神要QQ号码
引言 我们组来了个美女程序员,我心里窃喜,哈哈这下机会来了.我在想怎么下手呢?好吧,还是从QQ号码开始,找到女神要到QQ号,哈哈,我真是个天才~~~ 是这样子滴 想法是美好的,现实是残酷的,找女神要Q ...
- poj 3734 Blocks 快速幂+费马小定理+组合数学
题目链接 题意:有一排砖,可以染红蓝绿黄四种不同的颜色,要求红和绿两种颜色砖的个数都是偶数,问一共有多少种方案,结果对10007取余. 题解:刚看这道题第一感觉是组合数学,正向推了一会还没等推出来队友 ...
- 美女程序员是如何将QQ转换成题目中那串数字的--读博文《找女神要QQ号码》
我只能说好好的端午节你们不约么?,还在这里写代码?我也是够无聊的,下班了不走也在这跟风写着玩!<找女生要QQ号码原文>原文链接http://www.cnblogs.com/iforever ...
随机推荐
- 搭建docker 私有镜像仓库
前期准备 服务器:centos 7.3 docker-ce: 18.06.1-ce docker-compose: 1.22.0 docker 安装 首先,更新系统 yum update yum up ...
- kNN--近邻算法
kNN--近邻算法 kNN算法的核心思想是如果一个样本在特征空间中的k个最相邻的样本中的大多数属于某一个类别,则该样本也属于这个类别,并具有这个类别上样本的特性. 在机器学习中常用于分类. 数学内容: ...
- [python]序列的重复操作符
当你需要需要一个序列的多份拷贝时,重复操作符非常有用,它的语法如下: sequence * copies_int In [1]: a = [1,2,3,4] In [2]: a * 5 Out[2]: ...
- JUnit initializationError错误
一.JUnit Test 测试 initializationError错误 MyMaincom.test.sunc.MyMaininitializationError(com.test.sunc.My ...
- WeakHashMap介绍
WeakHashMap简介 WeakHashMap 继承于AbstractMap,实现了Map接口. 和HashMap一样,WeakHashMap 也是一个散列表,它存储的内容也是键值对(key ...
- 用 Python 3 的 async / await 做异步编程
前年我曾写过一篇<初探 Python 3 的异步 IO 编程>,当时只是初步接触了一下 yield from 语法和 asyncio 标准库.前些日子我在 V2EX 看到一篇<为什么 ...
- 404 Note Found -选题报告
目录 NABCD分析引用 N(Need,需求): A(Approach,做法): B(Benefit,好处): C(Competitors,竞争): D(Delivery,交付): 初期 中期 个人贡 ...
- POJ题目分类推荐 (很好很有层次感)
著名题单,最初来源不详.直接来源:http://blog.csdn.net/a1dark/article/details/11714009 OJ上的一些水题(可用来练手和增加自信) (POJ 3299 ...
- dubbo源码分析1——负载均衡
dubbo中涉及到的负载均衡算法只要有四种:Random LoadBalance(随机均衡算法).RoundRobin LoadBalance(权重轮循均衡算法).LeastAction LoadBa ...
- DS01--抽象数据类型
一.作业内容 二.数据结构.函数说明 1.头文件 common.h 2.数据结构 Rational.h 三.代码实现说明 1.构造有理数T 2.销毁有理数T 3.e返回有理数的分子或分母 4.用e改变 ...