摘要：
IV.LOJ#572. 「LibreOJ Round #11」Misaka Network 与求和 首先先考虑莫反推一波式子。 \(\begin{aligned}&\sum\limits_{i=1}^n\sum\limits_{j=1}^nf^k\Big(\gcd(i,j)\Big)\\=&\sum 阅读全文
posted @ 2021-04-05 23:05
Troverld
阅读(119)
评论(0)
推荐(0)
摘要：
III.UOJ#188. 【UR #13】Sanrd 题意：求 \(\sum\limits_{i=l}^rf(i)\)，其中 \(f(i)\) 为 \(i\) 的次大质因子。 显然其可以被转为两个前缀和相减的形式。 明显 \(f(i)\) 并非积性函数，所以常规min25筛处理不了。但是我们可以用非 阅读全文
posted @ 2021-04-05 23:01
Troverld
阅读(125)
评论(0)
推荐(0)
摘要：
II.LOJ#6053. 简单的函数 重申一下min25筛应用的条件： 是积性函数。 质数处取值是低阶多项式。 质数次幂处取值可以快速求出。 满足以上三点的任意函数均可以min25筛。 现在看到这题。乍一看 \(p\operatorname{xor}a\) 这种东西看上去一脸非多项式的样子；但是因为 阅读全文
posted @ 2021-04-05 22:59
Troverld
阅读(145)
评论(0)
推荐(0)
摘要：
4.min25筛 听说这玩意能干杜教筛干不了的事？ 同杜教筛一样，这也是用来求积性函数前缀和的东西。其复杂度为 \(O(\dfrac{n^{0.75}}{\log n})\)，大部分时候要略优于杜教筛。 min25筛作用的积性函数，应保证对于一切质数 \(p\)，\(f(p)\) 均是有关 \(p\ 阅读全文
posted @ 2021-04-05 22:36
Troverld
阅读(76)
评论(0)
推荐(0)
摘要：
VII.[NOI2016] 循环之美 依据小学数论知识，我们要求 \(\sum\limits_{i=1}^n\sum\limits_{j=1}^m[\gcd(i,j)=1][\gcd(j,k)=1]\) 因为后面的 \(k\) 是个常数，所以我们就想把它搞出来。 \(\begin{aligned}& 阅读全文
posted @ 2021-04-05 21:55
Troverld
阅读(106)
评论(0)
推荐(0)
摘要：
VI.LJJ爱数数 题目给出要求这样的东西 \(\dfrac{1}{a}+\dfrac{1}{b}=\dfrac{1}{c}\) 开始胡搞 \(\begin{aligned}\dfrac{1}{a}+\dfrac{1}{b}&=\dfrac{1}{c}\\\dfrac{a+b}{ab}&=\dfra 阅读全文
posted @ 2021-04-05 21:53
Troverld
阅读(83)
评论(0)
推荐(0)
摘要：
V.Product 要求这个东西： \(\prod\limits_{i=1}^n\prod\limits_{j=1}^n\dfrac{\operatorname{lcm}(i,j)}{\gcd(i,j)}\) 开始推式子。 \(\begin{aligned}\\&\prod_{i=1}^n\prod 阅读全文
posted @ 2021-04-05 21:49
Troverld
阅读(307)
评论(0)
推荐(0)
摘要：
IV.Lcm 既然上一道题中的DZY都能自定义函数，那我们为什么不能呢？ 定义$f(x)$为$x$中是否含有平方项。没有则为$1$，有则为$0$。显然，它是积性函数。而我们要求的，就是 \(\sum_{i=1}^n\sum_{j=1}^m\dfrac{ij}{\gcd(i,j)}f(\gcd(i,j 阅读全文
posted @ 2021-04-05 21:46
Troverld
阅读(221)
评论(0)
推荐(0)
摘要：
III.DZY Loves Math 题意：求 \(\sum\limits_{i=1}^n\sum\limits_{j=1}^mf(\gcd(i,j))\)，其中$f(x)$表示$x$的所有质因数中次数最高的一个的次数。 近乎套路的一堆操作后，我们得到了 \(\sum\limits_{i=1}^{\ 阅读全文
posted @ 2021-04-05 21:44
Troverld
阅读(59)
评论(0)
推荐(0)
摘要：
II.[CQOI2015]选数 我们要求这个东西： \(\sum\limits_{a_1=L}^R\sum\limits_{a_2=L}^R\dots\sum\limits_{a_n=L}^R[\gcd\limits_{i=1}^n(a_i)=k]\) 老套路，除一下，得到 \(\sum\limit 阅读全文
posted @ 2021-04-05 21:42
Troverld
阅读(51)
评论(0)
推荐(0)
摘要：
I.简单的数学题 在做这题之前，我们先来见一位老朋友： \(\sum\limits_{i=1}^n\sum\limits_{j=1}^n\gcd(i,j)\) 我们在1.v.[NOI2010]能量采集中就已经接触到了这道题。当时我们运用了$\sum\limits_{d|n}\mu(d)=[n=1]\ 阅读全文
posted @ 2021-04-05 21:38
Troverld
阅读(144)
评论(0)
推荐(0)
摘要：
3.杜教筛 之前在做莫反的题时，有很多题都需要用到杜教筛，因而我非常不爽。因此便来研究杜教筛了。 杜教筛可以干什么？ 在非线性时间内（准确说，\(O(n^{\frac{2}{3}})\)）求出某些积性函数的前缀和。例如，\(\sum_{i=1}^n\mu(i)\)。 怎么办呢？ 假设我们要求$S(n 阅读全文
posted @ 2021-04-05 21:36
Troverld
阅读(56)
评论(0)
推荐(0)
摘要：
2.狄利克雷卷积与数论函数 在1.v.[NOI2010]能量采集中，我们第一次认识到了狄利克雷卷积这个概念。下面我们将介绍它的更多性质。 我们之前得到了如下性质： \(\boxed{h(n)=(f*g)(n)\Leftrightarrow h(n)=\sum_{d|T}f(d)*g(\dfrac{T 阅读全文
posted @ 2021-04-05 21:34
Troverld
阅读(279)
评论(0)
推荐(0)
摘要：
ix.[51Nod1222]最小公倍数计数 求 \(\sum\limits_{i}\sum\limits_{j}\Big[\operatorname{lcm}(i,j)\in[a,b]\Big]\)。 考虑差分，问题转换为 \(\sum\limits_{i}\sum\limits_{j}\Big[\ 阅读全文
posted @ 2021-04-05 21:32
Troverld
阅读(65)
评论(0)
推荐(0)
摘要：
viii.[SDOI2017]数字表格 题意：求出 \(\prod\limits_{i=1}^n\prod\limits_{j=1}^mf_{\gcd(i,j)}\)，其中$f$是斐波那契数列。 就算是积，我们也一样能反演，只是反演到了指数头上。 \(\begin{aligned}\prod_{i= 阅读全文
posted @ 2021-04-05 21:28
Troverld
阅读(75)
评论(0)
推荐(0)
摘要：
vii.[SDOI2014]数表 仍然是线性筛筛各种东西。我们引出一个东西$\sigma(x)=\sum\limits_{d|n}d$，也就是$x$的约数和。这个东西明显是积性函数。设$x=\prod\limits_n(P_i)\(，则\)\sigma(x)=\prod\limits_n(\sum\ 阅读全文
posted @ 2021-04-05 21:23
Troverld
阅读(59)
评论(0)
推荐(0)
摘要：
vi.于神之怒加强版 在这之前，我们引出一个数论函数$idk(x)=xk$。这个函数就是整数域上的$k$次函数。很明显，它是积性函数，准确地说，是完全积性函数。 它的两个特例，一是$k=1$，就是我们之前提到的$id$函数。二是$k=0$，即$id0$函数。因为$\forall x\in\mathb 阅读全文
posted @ 2021-04-05 21:18
Troverld
阅读(81)
评论(0)
推荐(0)
摘要：
v.[NOI2010]能量采集 真正自己做出来的第一道莫反题祭~~~~ 题意： 求$\sum_^n\sum_^m(2\gcd(i,j)-1)$。 开始推式子： \(\begin{aligned}\sum_{i=1}^n\sum_{j=1}^m(2\gcd(i,j)-1) & =2\sum_{i=1} 阅读全文
posted @ 2021-04-05 21:14
Troverld
阅读(52)
评论(0)
推荐(0)
摘要：
iv.[SDOI2015]约数个数和 完蛋了，我们前几题里面都有$\gcd(\dots)$，但是这道题没有，怎么办呢？ 引理： \(\boxed{d(ij)=\sum\limits_{x|i}\sum\limits_{y|j}[\gcd(x,y)==1]}\) 换句话说，两个数$(i,j)$积的因数 阅读全文
posted @ 2021-04-05 21:12
Troverld
阅读(67)
评论(0)
推荐(0)
摘要：
iii.[HAOI2011]Problem b 第一道自己做出来的莫比乌斯反演题祭~~~ 实际上就是对上一道题套上一个类似于二维前缀和的东西。 把上一道题的东西的答案设为$calc(a,b,d)$， 则依据容斥原理，本题答案即为$calc(b,d,k)-calc(a-1,d,k)-calc(b,c- 阅读全文
posted @ 2021-04-05 21:10
Troverld
阅读(43)
评论(0)
推荐(0)
摘要：
ii.[POI2007]ZAP-Queries 如果前一道题没有听懂的话，是我的锅。毕竟这道题应该放在第一道题，上一道题明显是这道题的升级版。 首先，观察一下题目，发现这道题让我们求的就是上一道题中的$f(d)$。 我们再来推一下$f(d)$： 设$f(x)\(为\)\gcd(i,j)=x$的个数， 阅读全文
posted @ 2021-04-05 21:07
Troverld
阅读(56)
评论(0)
推荐(0)
摘要：
I.YY的GCD 这就是莫比乌斯反演？咋长得不像呢？ 我们看一下式子： \(ans=\sum\limits_{i=1}^n\sum\limits_{j=1}^m[\gcd(i,j)\ is\ prime]\)。其中方括号相当于强制把方括号内的东西转成$bool$形。 完蛋了，这个里面看不到任何函数， 阅读全文
posted @ 2021-04-05 21:05
Troverld
阅读(100)
评论(0)
推荐(0)
摘要：
1.莫比乌斯函数与莫比乌斯反演 O.约定 \(\color{white}\colorbox{red}{本blog中所有的分数，无论有无下取整符号，均默认下取整。}\) 主要是因为我太懒了，下取整符号的$\LaTeX$表达式太长了 I.作用 设有一函数$g(n)=\sum\limits_{d|n}f( 阅读全文
posted @ 2021-04-05 21:02
Troverld
阅读(136)
评论(0)
推荐(0)
摘要：
XXX.[HDU3151]Cave Crisis 一眼看上去跟XII.[ABC181F]Silver Woods完全一致，因此考虑一样的思路。 于是我们现在问题变为求出两个多边形间的距离。 首先先考虑如何判断它们是否有交。有交只有一种可能，就是边有交。于是我们枚举两个多边形所有的边，然后判断它们是否 阅读全文
posted @ 2021-04-05 21:00
Troverld
阅读(56)
评论(0)
推荐(0)
摘要：
XXIX.CF1195F Geometers Anonymous Club 闵可夫斯基和是关于两个凸包的运算，其几何意义是所有来自两个凸包内部的向量之和所构成的集合。 可以被证明的是，两个凸包的闵可夫斯基和，可以通过对两个凸包上的边按照极角大小排序后依次首尾相接得到。 回到本题。依照我们上述理论，我 阅读全文
posted @ 2021-04-05 20:57
Troverld
阅读(74)
评论(0)
推荐(0)
摘要：
XXVIII.[TopCoder12693]EnclosingTriangle 经典套路是固定一个点，求出所有合法的剩余两个点。 为了方便，我们将环状的图形拆开，拆成 \(4n\) 个点。然后，我们枚举一个点 \(i\) 明显发现，剩下两个点必定位于 \(i\) 两侧的一端区间内，不妨设一半是 \( 阅读全文
posted @ 2021-04-05 20:55
Troverld
阅读(87)
评论(0)
推荐(0)
摘要：
XXVII.CF611G New Year and Cake 做题时居然忘记了叉积满足分配律/jk 我们先将图形翻转成为逆时针排布。 首先，我们发现，若总图形的面积是 \(area\)，切完后，较小一半的面积是 \(nowarea\)，则贡献是 \(area-2nowarea\)。 我们记点 \(p 阅读全文
posted @ 2021-04-05 20:53
Troverld
阅读(56)
评论(0)
推荐(0)
摘要：
XXVI.[SHOI2012]信用卡凸包 一种异端想法是因为只需保证两位精度所以直接在每个半圆上取 \(100\) 个点代表该半圆，没有试过，但说不定也能过…… 书归正传。 我们考虑画出最终所得到的图形，发现其就是一堆小扇形，再加上中间的“裁去边角后的信用卡”的部分。大眼观察可得那堆小扇形拼一起就得 阅读全文
posted @ 2021-04-05 20:50
Troverld
阅读(63)
评论(0)
推荐(0)
摘要：
XXV.CF598F Cut Length 题解 阅读全文
posted @ 2021-04-05 20:48
Troverld
阅读(52)
评论(0)
推荐(0)
摘要：
XXIV.「SWTR-04」Taking a Walk 题解 阅读全文
posted @ 2021-04-05 20:45
Troverld
阅读(47)
评论(0)
推荐(0)
摘要：
XXIII.[HNOI2011]数矩形 题解 阅读全文
posted @ 2021-04-05 20:41
Troverld
阅读(45)
评论(0)
推荐(0)
摘要：
XXII.[USACO10OPEN]Triangle Counting G 题解 阅读全文
posted @ 2021-04-05 20:39
Troverld
阅读(35)
评论(0)
推荐(0)
摘要：
XXI.最小圆覆盖 随机增量法。 引理1.对于任意一组点集$\mathbb$和某点$P$，则$P$要么在$\mathbb\(的外接圆内，要么在\){\mathbb\cup P}$的外接圆上。 于是我们可以设计出如下的解法： 我们枚举一个$1\sim n$的变量i，并判断当前点是否在当前外接圆内。如果 阅读全文
posted @ 2021-04-05 20:37
Troverld
阅读(84)
评论(0)
推荐(0)
摘要：
XX.[COCI2009-2010#6] XOR 题解 阅读全文
posted @ 2021-04-05 20:35
Troverld
阅读(25)
评论(0)
推荐(0)
摘要：
XIX.[POI2007]OSI-Axes of Symmetry 题解 阅读全文
posted @ 2021-04-05 20:33
Troverld
阅读(47)
评论(0)
推荐(0)
摘要：
XVIII.[POI2010]OWC-Sheep 题解 阅读全文
posted @ 2021-04-05 20:31
Troverld
阅读(49)
评论(0)
推荐(0)
摘要：
XVII.[HNOI2012]射箭 强烈谴责本道卡精度屑题。 首先，乍一看，二次函数$ax^2+b$在$x=x_0$处的值要在$[y_0,y_1]$之内？带进去不就是一个关于$a,b$的半平面交吗？ 然后再一看，要找到半平面交非空的最大位置， 难不成要用动态凸包？ 后来想想动态凸包什么的完全没有必要 阅读全文
posted @ 2021-04-05 20:30
Troverld
阅读(64)
评论(0)
推荐(0)
摘要：
XVI.[JOISC2014]二人の星座 这题乍一看，和之前X.[POI2008]TRO-Triangles好像思想差不多；但是实际操作一番并进行了很多失败的尝试后，发现并不能简单应用。 后来知道了一种判两个（三点不共线的）三角形相离的做法：它们一定存在且只存在两条相同的外割线。（一个三角形的外割线 阅读全文
posted @ 2021-04-05 20:28
Troverld
阅读(71)
评论(0)
推荐(0)
摘要：
XV.CF1045E Ancient civilizations 神题。 我们先考虑如果凸包上只有三个点时的情形。 假如该三个点是同色的，我们考虑能否在该三角形内部找到一个异色点。假如能找到，我们便可以将这个大三角形拆分成三个小三角形，每个小三角形以该异色节点和凸包上两个点为顶点，这就使得小三角形的 阅读全文
posted @ 2021-04-05 20:26
Troverld
阅读(86)
评论(0)
推荐(0)
摘要：
XIV.CF70D Professor's task 没错，我们提到之前的Andrew算法的目的，就是为了运用在这题上——动态凸包。 明显这题不能用Graham，因为我们并不能找到一个固定的、一直在凸包上的点。 我们使用两个std::set来分别维护上下凸壳，当插入一个节点时，就向前向后不断删节点即 阅读全文
posted @ 2021-04-05 20:24
Troverld
阅读(66)
评论(0)
推荐(0)
摘要：
XIII.[USACO5.1]圈奶牛Fencing the Cows /【模板】二维凸包 相信大家都知道Graham算法（gift wrapping）。但是这个算法不好推广。这里使用便于推广的Andrew算法。 具体而言，本算法将所有点按照$x$为第一关键字，$y$为第二关键字排序，然后从前往后扫一 阅读全文
posted @ 2021-04-05 20:22
Troverld
阅读(61)
评论(0)
推荐(0)
摘要：
XII.[ABC181F]Silver Woods 一种很蠢的思路是把平面三角剖分然后建图，然后二分，但这很明显是个很蠢的主意 我们考虑反向思考——对于所有的有序点集$p_1,\dots,p_k$，球的直径一定$\leq\min\Big{dis(y=100,p_1),dis(p_1,p_2),\do 阅读全文
posted @ 2021-04-05 20:20
Troverld
阅读(116)
评论(0)
推荐(0)
摘要：
XI.[JSOI2016]炸弹攻击2 题解 阅读全文
posted @ 2021-04-05 20:18
Troverld
阅读(36)
评论(0)
推荐(0)
摘要：
X.[POI2008]TRO-Triangles 本题介绍两种做法。 一种是我自己的做法： 考虑某$\triangle ABC$，\(2S_{\triangle ABC}=\Big|\vec{AB}\times\vec{AC}\Big|=\Big|(\vec{B}-\vec{A})\times(\v 阅读全文
posted @ 2021-04-05 20:17
Troverld
阅读(93)
评论(0)
推荐(0)
摘要：
IX.[ZJOI2008]瞭望塔 ZJOI2008咋总喜欢考这种计算几何的题啊…… 同我们的IV.监视摄像机[CTSC1998]一样，本题相当于求多边形的核，可以使用半平面交。但因为和VI.[HNOI2008]水平可见直线一样，所有半平面都是向上的，故直接求凸包即可。 发现求完凸包后，答案的横坐标一 阅读全文
posted @ 2021-04-05 20:15
Troverld
阅读(68)
评论(0)
推荐(0)
摘要：
VIII. 题解 阅读全文
posted @ 2021-04-05 20:14
Troverld
阅读(53)
评论(0)
推荐(0)
摘要：
VII.[BalticOI 2005] Polygon 题解 阅读全文
posted @ 2021-04-05 20:12
Troverld
阅读(32)
评论(0)
推荐(0)
摘要：
VI.[HNOI2008]水平可见直线 一开始以为这是半平面交模板；后来一想，直接求出凸包来就行了。 我们仍然将所有直线按照斜率从小到大排序。不过这时，我们只需要使用单调栈维护即可。 具体而言，设栈顶次位、首位直线分别为$y=k_1x+b_1,y=k_2x+b_2$ 则其交点位于 \(\Bigg(\ 阅读全文
posted @ 2021-04-05 20:10
Troverld
阅读(59)
评论(0)
推荐(0)
摘要：
V.[NOI2004]降雨量 本题思路就和I.[HNOI2012]三角形覆盖问题差不多了。 首先，我们特判掉有长度为$W$的伞的情况——此时答案即为$0$。 否则，对于每一时刻，我们计算出下面三种情况中，最先来到的一个： 有一把伞撞到了边缘 有两把伞，它们的某两个边缘相遇了（不管是哪两个边缘） 时刻 阅读全文
posted @ 2021-04-05 20:08
Troverld
阅读(151)
评论(0)
推荐(1)
摘要：
IV.[CTSC1998]监视摄像机 这题就是半平面交模板。因为能看到一条边某侧的所有位置的一个点必定处于此边所在直线的内侧，故直接求半平面交即可。 另外，这题#5的第39个测试点似乎出了问题，得特判掉。 代码： #include<bits/stdc++.h> using namespace std 阅读全文
posted @ 2021-04-05 20:05
Troverld
阅读(80)
评论(0)
推荐(0)
摘要：
III.[CQOI2006]凸多边形 /【模板】半平面交 半平面交开始~~ 这里介绍一种做法：随机增量法，其可以在$O(n\log n)$的时间内完成半平面交的求解。 首先，我们可以用一个向量来表示直线；之所以使用向量来表示，是因为我们将强制该向量的左方表示半平面。 接着，我们考虑将所有直线按照向量 阅读全文
posted @ 2021-04-05 20:02
Troverld
阅读(77)
评论(0)
推荐(0)
浙公网安备 33010602011771号