首页 > 其他分享 >闲话 23.01.02

闲话 23.01.02

时间:2023-01-02 22:13:11浏览次数:59  
标签:02 right frac 闲话 sum 23.01 times sj binom

闲话

有人说不如让我把学术部分去掉,我就把闲话部分去掉了。

image

多项式杂题

凑数?凑数就凑数吧。

[HEOI2016/TJOI2016]求和

给定 \(n\),计算

\[\sum_{i=0}^n \sum_{j=0}^i {{i}\brace{j}}\times j!\times 2^j \]

\(n\le 10^5\)。

我们发现斯特林数不是很好做,于是先把它放在内层循环。

\[\sum_{j=0}^n j!\times 2^j \times \sum_{i=0}^n {{i}\brace{j}} \]

然后自然地展开斯特林数:

\[\sum_{j=0}^n j!\times 2^j \times \sum_{i=0}^n \sum_{k=0}^j (-1)^k \frac{(j - k)^i}{(j - k)! k!} \]

把 \(i\) 放在里面求和:

\[\sum_{j=0}^n j!\times 2^j \times \sum_{k=0}^j (-1)^k \frac{(j - k)^{n + 1} - 1}{(j - k - 1)(j - k)! k!} \]

然后就可以卷积了。复杂度 \(O(n\log n)\)。

当然这题存在 \(O(n)\) 做法。具体数学的公式 \((6.19)\) 提供了一种新的展开式,套用得到

\[\sum_{j=0}^n 2^j \sum_{i=0}^n \sum_{k=0}^j \binom{j}{k} k^i(-1)^{j-k} \]

然后用 \(\text{q-analog}\) 化简一下

\[\sum_{k=0}^n [n + 1]_k\sum_{j=k}^n 2^j \binom{j}{k} (-1)^{j-k} \]

考虑一个

\[\sum_{m = i}^n \binom{m}{i} a_m q^{n - m} \]

的形式。

我们发现,若 \(F(x)\) 生成 \(\langle a_m \rangle\),则这个东西的生成函数就是 \(F(x + q)\)。现在 \(q = -1, F(x) = \sum_{i=0}^n (2x)^i = [n + 1]_{2x}\),因此我们得到了 \(F(x - 1) = [n + 1]_{2x - 2}\)。这个东西微分有限,可以递推。

然后 \([n]_k\) 容易在 \(O(n)\) 复杂度内求得。因此总复杂度为 \(O(n)\)。代码没实现。



[HAOI2018]染色

为了报答小 C 的苹果,小 G 打算送给热爱美术的小 C 一块画布,这块画布可以抽象为一个长度为 \(N\) 的序列,每个位置都可以被染成 \(M\) 种颜色中的某一种。然而小 C 只关心序列的 \(N\) 个位置中出现次数恰好为 \(S\) 的颜色种数,如果恰好出现了 \(S\) 次的颜色有 \(K\) 种,则小 C 会产生 \(W_k\) 的愉悦度。小 C 希望知道对于所有可能的染色方案,他能获得的愉悦度的和对 \(1004535809\) 取模的结果是多少。

\(n\le 10^7, m\le 10^5, S\le 150\)。

我们考虑使用 egf 来刻画这个性质。假设 \(F, G\) 分别生成确定颜色的位置和不确定颜色的位置,取第 \(n\) 项表示用特定的一种颜色填 \(n\) 个位置的种类数。由于 \(F\) 对应着已经确定了位置数的颜色,因此 \(F(x) = x^s / s!\),\(G\) 则正好相反,是 \(\sum_{i = 0}^{\infty} x^i / i! - x^s / s!\),即 \(e^x - x^s / s!\)。取颜色数为幂次,并分配位置得到这部分的总方案数

\[[x^n / n!] \left(\frac {x^s}{s!}\right)^K\left(e^x - \frac {x^s}{s!}\right)^{m - K} \]

考虑化简。

\[\begin{aligned} &\sum_{i=0}^m \binom mi\times w_i\times [x^n / n!] \left(\frac {x^s}{s!}\right)^i\left(e^x - \frac {x^s}{s!}\right)^{m - i} \\ = \ & n!\times \sum_{i=0}^m \binom mi\times w_i\times [x^n] \left(\frac {x^s}{s!}\right)^i\left(e^x - \frac {x^s}{s!}\right)^{m - i} \\ = \ & n!\times \sum_{i=0}^m \binom mi\times w_i\times [x^{n - si}] \frac{1}{(s!)^i}\left(e^x - \frac {x^s}{s!}\right)^{m - i} \end{aligned}\]

先看提取系数。

\[\begin{aligned} &[x^{n - si}] \frac{1}{(s!)^i}\left(e^x - \frac {x^s}{s!}\right)^{m - i} \\ = \ & [x^{n - si}] \frac{1}{(s!)^i}\sum_{j = 0}^{m - i} (-1)^{j} \binom{m - i}{j} e^{x(m - i - j)}\left(\frac {x^s}{s!}\right)^{j} \\ = \ & \sum_{j = 0}^{m - i} (-1)^{j} \binom{m - i}{j} \frac{1}{(s!)^{i + j}} [x^{n - s(i + j)}]e^{x(m - i - j)} \\ = \ & \sum_{j = i}^{m} (-1)^{j-i} \binom{m - i}{j - i} \frac{1}{(s!)^{j}} [x^{n - sj}]e^{x(m - j)} \\ = \ & \sum_{j = i}^{m} (-1)^{j-i} \binom{m - i}{j - i} \frac{1}{(s!)^{j}} \frac{(m - j)^{n-sj}}{(n-sj)!} \end{aligned}\]

回代得到

\[\begin{aligned} & n!\times \sum_{i=0}^m \binom mi\times w_i\times \sum_{j = i}^{m} (-1)^{j-i} \binom{m - i}{j - i} \frac{1}{(s!)^{j}} \frac{(m - j)^{n-sj}}{(n-sj)!} \\ = \ & m! n!\times \sum_{i=0}^m \frac{w_i }{i!(m - i)!}\times \sum_{j = i}^{m} (-1)^{j-i} \frac{(m - i)!}{(j - i)!(m-j)!}\frac{(m - j)^{n-sj}}{(n-sj)!(s!)^{j}} \\ = \ & m! n!\times \sum_{i=0}^m \frac{w_i}{i!}\times \sum_{j = i}^{m} \frac{(-1)^{j-i}}{(j - i)!}\frac{(m - j)^{n-sj}}{(m-j)!(n-sj)!(s!)^{j}} \end{aligned}\]

然后后面那个部分就可以卷积了。具体地,我们设后面那个是 \(g_{j-i} \times f_j\) 的形式。我们反转 \(f\),然后令 \(j = m - j\),新的形式就是 \(\sum_{j=0}^{m-i} g_{m-j-i}f_j\)。

直接做即可。总时间复杂度 \(O(n\log n)\)。

听说有个什么容斥还是反演的做法,我永远不会去看的!

Submission.



bzoj5093 图的价值

我们分别考虑每个点的贡献。假设这个节点的度数是 \(i\),它可以向外连边的方案数是 \(C_{n-1}^i\) 种,每种带贡献 \(2^{C_{n-1}^2}\times i^k\)。又因为有 \(n\) 个点,总贡献即为

\[n 2^{\binom{n-1}{2}} \sum_{i=0}^{n-1} i^k \binom {n-1}{i} \]

我们只看求和号。把幂用斯特林数化开,我们得到

\[\sum_{i=0}^{n-1} \sum_{j=0}^i\binom {n-1}{i} \binom{i}{j} {{k}\brace{j}} j! =\sum_{j=0}^{n-1} \sum_{i=j}^{n-1} \binom {n-1}{j} \binom{n-1-j}{i-j} {{k}\brace{j}} j! \]

于是把 \(i\) 消去,得到

\[\sum_{j=0}^{n-1} 2^{n-1-j} \binom {n-1}{j}{{k}\brace{j}} j! \]

然后我们发现 \(j \le \min(n-1, k)\)。这样直接求一行斯特林数,我们就能枚举 \(j\) 得到答案了。

没了。总时间复杂度为 \(O(n\log n)\)。

标签:02,right,frac,闲话,sum,23.01,times,sj,binom
From: https://www.cnblogs.com/joke3579/p/chitchat230102.html

相关文章

  • 2023.1.2周报
    2023.1.2周报本周总结:阳了,所以一直处于虚弱状态,现在基本好了,接下去可以恢复正常训练了。大致学了下矩阵的运用,但大部分时间在看文档,没什么精力写题目。矩阵构造有非常多......
  • 代码随想录算法训练营第五天LeetCode242,349,202,1
    代码随想录算法训练营第五天|LeetCode242,349,202,1LeetCode242有效的字母异位词题目链接:https://leetcode.cn/problems/valid-anagram/description///开了两个哈希数组,......
  • HTML培训课程-------Day02(表格和框架)
    表格在网页中表格是一种经常使用到得设计结构,就像表格的内容中可以包含任何的数据,如文字、图像、表单、超链接、表格等等,所有在HTML中可以使用的数据,都可以被设置在表格中,所......
  • Codeforces Good Bye 2022 CF 1770 F Koxia and Sequence 题解
    题目链接注意题目要求的是所有好序列的所有元素的XOR之和的XOR之和。我一开始以为是所有XOR之和的加法和导致不知道官方题解在讲什么。假设我们枚举一个位置\(1\lei\le......
  • the eleventh——2023.1.2
    scanf()函数一般只读取字符串中的一个单词,而不是一句话。例如:scanf("%s",name);printf("Hello,%s!",name) NingBabaHello,Ning!(后面的Baba在scanf这读取不到,在遇......
  • 2023.1.2周报
    本周总结:学习了《算法竞赛》第六章数论6.7-6.9、第七章组合数学7.1-7.6内容,牛客组合数学课程,做书上例题和习题。准备新手课堂文档和讲课,顺便出结训赛题目。大方向组合数......
  • 报告分享|2022年中国新冠医药研发趋势报告
    2020年爆发的新冠肺炎是百年一遇的全球性公共卫生事件,至今疫情虽有缓和趋势,但影响仍然存在。面对严峻疫情,中国全力投入新冠相关的研发创新,力图降低疫情影响、终结疫情蔓延......
  • 代码随想录day6 哈希表 LeetCode 242.有效的字母异位词 349. 两个数组的交集 202.
    当我们遇到了要快速判断一个元素是否出现集合里的时候,就要考虑哈希法242.有效的字母异位词 https://leetcode.cn/problems/valid-anagram/由于哈希表只有26个元素,故采......
  • 2023/01 LeetCode练习
    ......
  • 2023,整装待发!
    我是个怎么样的人呢.....高考过后,我常常断断续续的思索这些。但是对刚刚“解放”的我来说,思考这些还是过于疲惫和自找不快的一件事。所以就这样遗忘了,在混沌与半清醒......