日记

日记

Mon May 26 2025
2 分钟

这 gbc 看的我有点 emo 怎么个事?

GCD - Extreme (II) Link to

我们发现他让我们求这个东西:i=1nj=i+1ngcd(i,j)\sum\limits_{i = 1}^{n}\sum\limits_{j = i + 1}^{n}\gcd(i,j),显然可以转化一下:

i=1nj=i+1ngcd(i,j)=x=1nxj=1ni=1j1[gcd(i,j)=x]=x=1nxj=1nxi=1j=1[gcd(i,j)=1]=x=1nxj=1nxφ(j)\begin{aligned} & \sum_{i = 1}^{n}\sum_{j = i + 1}^{n}\gcd(i,j) \\ & = \sum_{x = 1}^{n}x\sum_{j = 1}^{n}\sum_{i = 1}^{j - 1}[\gcd(i,j) = x] \\ & = \sum_{x = 1}^{n}x\sum_{j = 1}^{\lfloor \frac{n}{x} \rfloor}\sum_{i = 1}^{j = 1}[\gcd(i,j) = 1] \\ & = \sum_{x = 1}^{n}x\sum_{j = 1}^{\lfloor \frac{n}{x} \rfloor} \varphi(j) \end{aligned}

到这个时候已经可以 O(nn)O(n \sqrt{n}) 做了,但是我们发现 n106n \le 10^6O(nn)O(n \sqrt{n}) 的复杂度显然过不了。我们发现:一车 nx\lfloor \frac{n}{x} \rfloor 都是一样的。于是我们直接数论分块。复杂度直接降到 O(n)O(\sqrt{n})

GCD Link to

也像上面一样,枚举每个质数,然后变成 φ(x)\varphi(x) 的形式。

报数 Link to

我们发现:i=1np(i)4×106\sum\limits_{i = 1}^{n}p(i) \approx 4 \times 10^6,套个 log\log 也能过。于是直接 O(nlogn)O(n \log n) 筛就行。

Trailing Loves (or L’oeufs?) Link to

我都能做出来的题肯定是不用说了。

后日谈 Link to 后日谈

今天做数论做疯了。

溜冰溜爽了。给我看的差点哭出来怎么个事?

照理来说这里是日记,应该写点感想,但是我写日记的时候忘了当时的感觉了,于是不写了。再听一遍《熙熙攘攘,我们的城市》就下。