Statement

带标号仙人掌计数问题.

\(n< 131072\).

Solution

设\(x\)个点的仙人掌个数的生成函数为\(C(x)\)

  • 对于与根相邻的块, 还是仙人掌, 生成函数为\(C(x)\)

  • 包含根的环, 生成函数为\(\sum_{i>1}\frac{C(x)^i}{2}\)

组合起来:
\[
C(x) = x \exp{\frac{2C(x)-C(x)^2}{2-2C(x)}}
\]
设\(G(C(x)) = x\exp{\frac{2C(x)-C(x)^2}{2-2C(x)}}-C(x)\), 那么:
\[
\small{
\begin{aligned}
G'(C(x)) &= x\left(\exp{\frac{2C(x)-C(x)^2}{2-2C(x)}}\right)'-1 \\
&= x \exp\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)' - 1 \\
&= x \exp\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)
\left(\frac{\left(2-2C(x)\right)^2-\left(2C(x) - C(x)^2\right)(-2)}{(2-2C(x))^2}\right)
- 1\\
&= x \exp\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)
\left(1+\frac{4C(x) - 2C(x)^2}{(2-2C(x))^2}\right)
- 1
\end{aligned}
}
\]
牛顿迭代:
\[
\begin{aligned}
C_1(x) &= C(x) - \frac{G(C(x))}{G'(C(x))} \\
&= C(x) - \frac{2x\exp\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)-2C(x)}
{x \exp\left(\frac{2C(x)-C(x)^2}{2-2C(x)}\right)
\left(1+\frac{1}{(C(x)-1)^2}\right)
- 2}
\end{aligned}
\]

[LOJ161] 仙人掌计数的更多相关文章

  1. [LOJ6569] 仙人掌计数

    Statement 带标号仙人掌计数问题. \(n< 131072\). Solution 设\(x\)个点的仙人掌个数的生成函数为\(C(x)\) 对于与根相邻的块, 还是仙人掌, 生成函数为 ...

  2. WinterCamp2017 游记

    Winter is coming! Day0 Day0前一天打了一轮CF,做完了ABCD,Div2 Rank59.然后就去开开心心的睡觉,准备第二天的行程. 快到一点的时候躺在了床上,睡不着,翻来覆去 ...

  3. [日常] NOIWC2019 冬眠记

    NOIWC 2019 冬眠记 辣鸡rvalue天天写意识流流水账 Day 0 早上没有跑操(极度舒服.png) 和春哥在博客颓图的时候突然被来送笔电的老爹查水表(捂脸) 母上大人骗我说这功能机不能放存 ...

  4. WC2017游记

    Day0 到杭州之后出了点锅换了辆车,等了好久才开= =到宿舍发现路由器就在房门口,稳啊,过了一会儿就连不上了= =而且只有门口那个连不上,可以连上楼下的= =之后干了啥也忘了…… Day1 上午直接 ...

  5. REHの收藏列表

    搬运自本人的AcWing,所以那里的文章会挺多. 友链(同类文章) :bztMinamoto 世外明月 mlystdcall 新人手册:AcWing入门使用指南 前言 有看到好文欢迎推荐(毛遂自荐也可 ...

  6. $dy$讲课总结

    字符串: 1.广义后缀自动机(大小为\(m\))上跑一个长度为\(n\)的串,所有匹配位置及在\(parent\)树上其祖先的数量的和为\(min(n^2,m)\),单次最劣是\(O(m)\). 但是 ...

  7. UOJ#290. 【ZJOI2017】仙人掌 仙人掌,Tarjan,计数,动态规划,树形dp,递推

    原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ290.html 题解 真是一道好题! 首先,如果不是仙人掌直接输出 0 . 否则,显然先把环上的边删光. ...

  8. 【做题】ZJOI2017仙人掌——组合计数

    原文链接 https://www.cnblogs.com/cly-none/p/ZJOI2017cactus.html 给出一个\(n\)个点\(m\)条边的无向连通图,求有多少种加边方案,使得加完后 ...

  9. [BZOJ]1016 JSOI2008 最小生成树计数

    最小生成树计数 题目描述 现在给出了一个简单无向加权图.你不满足于求出这个图的最小生成树,而希望知道这个图中有多少个不同的最小生成树.(如果两颗最小生成树中至少有一条边不同,则这两个最小生成树就是不同 ...

随机推荐

  1. 从Android手机中取出已安装的app包,导出apk

    从Android手机中取出已安装的app包,导出apk TAG:Android,提取,apk,adb,pm,root,导出apk 假设有这样一个场景,A君看到你手机上一个实用APP,想要安装到自己手机 ...

  2. node和数据库建立连接

    var express = require('express') , app = express(); var querystring = require('querystring'); var ut ...

  3. virutalenv一次行安装多个requirements里的文件

  4. 【ABAP系列】SAP ABAP ALV中的TOP_OF_PAGE添加任意图标

    公众号:SAP Technical 本文作者:matinal 原文出处:http://www.cnblogs.com/SAPmatinal/ 原文链接:[MM系列]SAP ABAP ALV中的TOP_ ...

  5. 工具 - VNC

    安装 ubuntu下vnc客户端的安装命令sudo apt-get install xtightvncviewer 重置密码 cd /root/.vnc/ rm -rf passwd vncserve ...

  6. spring boot 将对象转换为json返回

    Spring Boot默认使用Jackson将对象转换为json,需要配置以下依赖: compile group: 'com.fasterxml.jackson.core', name: 'jacks ...

  7. 模块内高内聚?模块间低耦合?MVC+EF演示给你看!

    前言 在软件项目开发过程中,我们总能听见“高内聚,低耦合”,即使这种思想在我们学习编程的过程中就已经耳濡目染.可一旦当我们上项目,赶进度的时候我们就会“偷懒”,能省时间就省.管他什么设计模式,什么软件 ...

  8. 2019寒假作业二:PTA7-1币值转换

    7-1 币值转换 (20 分) 输入一个整数(位数不超过9位)代表一个人民币值(单位为元),请转换成财务要求的大写中文格式.如23108元,转换后变成“贰万叁仟壹百零捌”元.为了简化输出,用小写英文字 ...

  9. [BNDSOJ] #1106代码

    #include<bits/stdc++.h> using namespace std; ]; ][]; int n; bool check(int i,int j) { ]==]==]= ...

  10. CSU-1110 RMQ with Shifts (单点更新+区间最小值 zkw线段树)

    In the traditional RMQ (Range Minimum Query) problem, we have a static array A. Then for each query ...