闲话 23.1.26
闲话
下午补。
今日推歌:FREELY TOMORROW - Mitchie M feat. 初音ミク
有兴趣的可以自己了解。
今日模拟赛 T3
小恶心。原题似乎是 pe444,但是好像不让写题解的样子?
反正不是原题啦,随便!
首先观察 \(E(111) = 5.2912\) 这个数值就很不对劲。扔进计算器发现 \(H(111) = 5.291243607252986\)。这里的 \(H\) 就是调和级数,常这么记。
然后猜了这个结论,也就是 \(S_0(n) = H(n)\)。
具体是咋证的呢?
我们假设现在翻开的牌上数的集合是 $ A$,它的补是 $ B$。有一个结论是如果 \(\max A > \min B\) 就和前面的最大值换,反之就不换。
比较显然?因为你可以选择的一定是 \(\max A\),选了一定是不劣的。这点也是需要自己看不见自己的牌才能决定的。
然后发现只有自己是后缀最大值时才会不换。假设当前的排列是 \(p\),我们设 \(suff(i)\) 是 \(i\sim n\) 的最大值,我们要的就是 \(\sum[p_i = suff(i)]\) 的期望。
这个东西就是 \(H(n)\)。
看到前缀和就想到用生成函数来刻画。假设 \(H(n)\) 的生成函数是 \(F(x)\),我们考虑表出 \(F(x)\)。
\(H(n)\) 的差分是 \(\dfrac{1}{n}\),而这个东西是好表示的,就是 \(-\ln(1-x)\),证明考虑对 \((1-x)^{-1}\) 做积分。场上以为是 x((1-x)^{-1})' 所以推了老半天错解
于是 \(F(x) = \dfrac{-\ln(1-x)}{(1-x)}\)。\(S_k(n)\) 的生成函数 \(F_k(x)\) 也好写出了,就是
我们想要得到递推,考虑分析 ode。
我们容易提取上面的式子两边的系数,也就是
我们希望对 \(k\) 作递推,化简得到
到这里我们就可以快速地递归到 \(S_0(n)\) 后计算了,组合数可以 \(O(k)\) 地求得。
那我们如何得到 \(S_0(n)\) 呢?
\(n\) 较小时可以预处理 \(H(n)\)。\(n\) 较大时呢?可以发现答案对精度的要求并不高。因此我们容易想到
其中 \(\epsilon_n = O\left(\frac 1n\right)\),\(\gamma\) 是 Euler-Mascheroni 常数,约为 0.57721566490153286060651209
。背不过也可以直接算一个出来,几秒钟就能得到小数点后 ⑨ 位精度!
code
#include <bits/stdc++.h>
using namespace std;
using pii = pair<int,int>; using vi = vector<int>; using vp = vector<pii>; using ll = long long;
using ull = unsigned long long; using db = double; using ld = long double; using lll = __int128_t;
template<typename T1, typename T2> T1 max(T1 a, T2 b) { return a > b ? a : b; }
template<typename T1, typename T2> T1 min(T1 a, T2 b) { return a < b ? a : b; }
#define multi int T; cin >> T; while ( T -- )
#define timer cerr << 1. * clock() / CLOCKS_PER_SEC << '\n';
#define iot ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
#define file(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define rep(i,s,t) for (register int i = (s), i##_ = (t) + 1; i < i##_; ++ i)
#define pre(i,s,t) for (register int i = (s), i##_ = (t) - 1; i > i##_; -- i)
#define eb emplace_back
#define pb pop_back
const int N = 1e6 + 10;
const int inf = 0x3f3f3f3f;
const ll infll = 0x3f3f3f3f3f3f3f3fll;
#define gamma 0.57721566490153286060651209
#define H(n) (__builtin_logl(n) + gamma)
ll k, n;
ld h[1000003];
ld S(int k, ll n) {
if (k == 0 and n <= 1e6) return h[n];
if (k == 0 and n > 1e6) return H(n);
ld ret = 1. / k;
rep(i,1,k) ret = ret / i * (n + i);
return S(k - 1, n + 1) / k * (n + 1) - ret;
}
signed main() {
iot; cout.unsetf(ios::fixed); cout.setf(ios::scientific); cout.precision(9);
cin >> k >> n;
rep(i,1,min(1e6, n + k + 1)) h[i] = h[i - 1] + (ld)1. / i;
cout << S(k, n) << '\n';
}
一个小级数
闲的没事看 b 站,然后看到一个视频。
没看 up 主的讲解,自己搞了一会儿出来了,发现挺好玩的,写写。
我们想要的就是证明
考虑从左往右证。首先这个 \(i\) 的范围挺难受的,不如从 \(0\) 开始,也就是
提出来 \(x^2\),再合并一下同类项得到
用原丁真,鉴定为可以直接套广义二项级数。
根据 定理 附 \(3.1.3\),我们可以写出后面的封闭形式,也就是
众所周知,\(\mathcal B _2(z)\) 是卡特兰数的生成函数,也就是 \(\dfrac{1 - \sqrt{1 - 4z}}{2z}\)。带入 \(z = x - x^2\) 可以得到化简的答案即为 \(\dfrac{2 - 2p}{2p(1-p)} = \dfrac{1}{p}\)。
你说根号咋确定正负?还是考察是否合理。
带回原式,就是
结束啦!
以下是博客签名,与正文无关。
请按如下方式引用此页:
本文作者 joke3579,原文链接:https://www.cnblogs.com/joke3579/p/chitchat230126.html。
遵循 CC BY-NC-SA 4.0 协议。
请读者尽量不要在评论区发布与博客内文完全无关的评论,视情况可能删除。