ixRic-CSDN博客

原创 [CQOI 2015] 选数（莫比乌斯反演 + 杜教筛） | 错题本

文章目录题目分析题目[CQOI 2015] 选数分析∑a1=LH∑a2=LH⋯∑aN=LH[gcd⁡i=1Nai=K]=∑a1=⌈LK⌉⌊HK⌋∑a2=⌈LK⌉⌊HK⌋⋯∑aN=⌈LK⌉⌊HK⌋[gcd⁡i=1Nai=1]=∑d=1⌊HK⌋μ(d)f(d)\begin{aligned}&\sum_{a_1= L}^{H}\sum_{a_2 = L}^{H} \cdots \sum_{a_N = L}^{H}[\gcd_{i = 1}^N a_i = K] \\ =&\sum_{a_

2020-09-07 12:02:39 174

原创杜教筛及其时间复杂度分析

文章目录杜教筛方法举例莫比乌斯函数欧拉函数时间复杂度杜教筛杜教筛用于求一类积性函数的前缀和，时间复杂度可以做到 O(n23)O(n^{\frac{2}{3}})O(n32)。方法设我们要求的是积性函数函数 f(x)f(x)f(x) 的前缀和 S(n)=∑i=1nf(i)S(n) = \sum_{i = 1}^n f(i)S(n)=∑i=1nf(i)。我们取一个积性函数 ggg（ggg 具体是什么由 fff 决定，但这不影响推式子）则 ∑i=1n(f∗g)(i)=∑i=1n∑d∣ig(d)f(i

2020-09-07 10:44:15 981 3

原创 [HAOI 2011] Problem b（莫比乌斯反演） | 错题本

文章目录题目分析代码题目[HAOI 2011] Problem b分析转换成求四块的答案，然后和 YY的GCD 一样。代码#include <bits/stdc++.h>typedef long long LL;const int MAXN = 50000;int Mu[MAXN + 5];bool Vis[MAXN + 5];std::vector<int> Primes;void Init(int n) { Mu[1] = 1; for (in

2020-09-05 11:34:33 216

原创 [洛谷 P2257] YY的GCD（莫比乌斯反演） | 错题本

文章目录题目分析代码题目[洛谷 P2257] YY的GCD分析不妨设 N≥MN \geq MN≥M，PPP 是素数集，枚举 d=gcd⁡(x,y),i=xd,j=ydd = \gcd(x, y), i = \frac{x}{d}, j = \frac{y}{d}d=gcd(x,y),i=dx,j=dy，要求的即 ∑d=1M[d∈P]∑i=1⌊Nd⌋∑j=1⌊Md⌋[gcd⁡(i,j)=1]=∑d=1M[d∈P]∑i=1⌊Nd⌋∑j=1⌊Md⌋∑k∣gcd⁡(i,j)μ(k)=∑d=1M[d∈P]

2020-09-05 11:17:00 198

原创 [CodeForces 839D] Winter is here（莫比乌斯反演） | 错题本

文章目录题目分析代码题目[CodeForces 839D] Winter is here分析设 mmm 是最大的数，枚举 gcd⁡=d\gcd = dgcd=d，子序列长度为 lll，要求的就是 ∑d=2m∑l=1nd⋅l∑1≤k1<k2<⋯<kl≤n[gcd⁡i=1laki=d]=∑d=2m∑l=1nd⋅l∑1≤k1<k2<⋯<kl≤n(∏i=1l[d∣aki])[gcd⁡i=1lakid=1]\begin{aligned} &\sum_{d = 2}^

2020-09-05 09:47:33 205

原创 [CodeForces 547C] Mike and Foam（莫比乌斯反演） | 错题本

文章目录题目分析代码题目[CodeForces 547C] Mike and Foam分析设 cic_ici 表示 iii 的个数，要求 i<ji < ji<j 比较麻烦，考虑求 t=∑i=1n∑j=1n[gcd⁡(ai,aj)=1]t = \sum_{i = 1}^n\sum_{j = 1}^{n}[\gcd(a_i, a_j) = 1]t=i=1∑nj=1∑n[gcd(ai,aj)=1] 则答案为 t−c12\frac{t - c_1}{2}2t−c1。设 mmm

2020-09-04 22:14:21 235

原创 [CodeForces 439E] Devu and Birthday Celebration（莫比乌斯反演） | 错题本

文章目录题目分析代码题目[CodeForces 439E] Devu and Birthday Celebration分析莫比乌斯函数比较重要的性质：μ∗1=ε\mu * 1 = \varepsilonμ∗1=ε 即 ∑d∣nμ(d)=[n=1]\sum_{d | n}\mu(d) = [n = 1]d∣n∑μ(d)=[n=1]证明：设 n=∏i=1mpiαin = \prod_{i = 1}^{m} p_i^{\alpha_i}n=∏i=1mpiαi，n′=∏i=1mpin' = \

2020-09-04 12:33:29 164

原创 [洛谷 P3935] Calculating（整数分块） | 错题本

文章目录题目分析代码题目[洛谷 P3935] Calculating分析f(i)f(i)f(i) 是约数个数函数。令 g(n)=∑i=1nf(i)g(n) = \sum_{i = 1}^{n} f(i)g(n)=∑i=1nf(i)，则 ∑i=lrf(i)=g(r)−g(l−1)\sum_{i = l}^{r}f(i) = g(r) - g(l - 1)∑i=lrf(i)=g(r)−g(l−1)。g(n)=∑i=1n∑j=1i[i mod j=0]=∑j=1n⌊nj⌋\be

2020-09-04 11:01:13 221

原创 [LOCAL] 过河（数学） | 错题本

文章目录题目分析代码题目一只青蛙想从河的一侧跳到另一侧。河中有一些石子，青蛙必须要经过这些石子才能够顺利过河。由于青蛙一次跳过的距离是总是正整数，我们可以把河中青蛙可能到达的点看成数轴上的一串整点：0,1,⋯ ,L0,1,\cdots,L0,1,⋯,L（其中 LLL 是河的宽度）。当青蛙的坐标大于 LLL 或者小于 000 时，我们就认为青蛙跳出了这条河。一旦青蛙跳出这条河，青蛙就不会选择回到河中。注意青蛙不一定最后要跳到河的另一侧。现在，青蛙将会从 000 坐标开始，每次可以随意选择向正方向跳跃

2020-09-01 10:07:02 230

原创 [LOCAL] 画画图（结论 + DP + 常数分析） | 错题本

文章目录题目分析代码题目题目描述小胡同学是一个热爱画画的好（xiao）孩（shen）子（tong）。小胡同学画了一棵树，小胡同学的树是这么画的：在纸上画 nnn 个点，从 111 到 nnn 标号。随机选两个点 x,yx, yx,y，若 xxx 和 yyy 已经连通则继续随机，直到 x,yx, yx,y 不连通为止。在 x,yx, yx,y 之间连一条边。若 nnn 个点已经连通则停止，否则返回 222。小胡同学画完这棵树之后觉得似乎少了点什么，于是他给每条边都加了一个权值，这样他就得

2020-08-21 11:56:30 343

原创 [CodeForces 494C] Helping People（树形 DP） | 错题本

文章目录题目分析错因代码题目[CodeForces494C] Helping People分析区间不相交意味着各个区间之间形成了一个树形结构。设原数组最大值为 mmm，我们只需要保留 [m,m+q][m, m + q][m,m+q] 之间的数，又 q∈[1,5000]q \in [1, 5000]q∈[1,5000] 意味着我们可以将“最大值”最为下标。于是定义状态 dp[i][j]dp[i][j]dp[i][j] 表示 iii 子树中最大值与 mmm 的差小于等于 jjj 的概率，因为最大值与

2020-08-21 11:35:12 162

原创 [CodeFoeces 1025D] Recovering BST（DP 状态优化） | 错题本

文章目录题目分析代码题目[CodeFoeces 1025D] Recovering BST分析一个子树必然对应一个连续的区间，记为 [l,r][l, r][l,r]，发现一个子树根的父亲必定是 l−1l - 1l−1 对应的节点或 r+1r + 1r+1 对应的节点，仅两种情况。因此我们将状态 dp[i][j][k]dp[i][j][k]dp[i][j][k]（左端点、右端点、根节点）优化成 dp[i][j][0/1]dp[i][j][0/1]dp[i][j][0/1]（左端点、右端点、根节点的父亲

2020-08-21 10:56:33 153

原创 [CodeFoeces 1051E] Vasya and Big Integers（单调性 DP + Hash） | 错题本

文章目录题目分析代码题目[CodeFoeces 1051E] Vasya and Big Integers分析用双指针维护 dp[i]dp[i]dp[i] 的转移范围，问题的难点就是快速比较两个字符串的大小。比较两个字符串的大小转化为比较 LCP 的后一个字符的大小，因此 Hash + 二分即可。要注意处理前导零的问题。找 LCP 的时候先暴力比前面的几位，否则会被卡 Hash。真就调了整整一天。代码调得很暴躁，代码不好看。#include <algorithm>#incl

2020-08-21 10:49:23 201

原创 [洛谷 P1341] 无序字母对（欧拉路径） | 错题本

文章目录题目分析错因代码题目[洛谷 P1341] 无序字母对分析找一条字典序最小的欧拉路径即可。判无解：奇点个数只能为 000 或 222 ，否则都是无解的。并且如果有奇点的话，起点是奇点。具体方法就是走一条边删一条边，退栈的时候储存答案。错因欧拉路径要退栈时存！原因见欧拉路径为什么要退栈时储存并逆序输出。代码#include <bits/stdc++.h>const int MAXN = 60;int N;int Deg[MAXN + 5];bool G[

2020-08-21 10:24:54 187

原创欧拉路径为什么要退栈时储存并逆序输出

考虑一个有两个奇点的图：找字典序最小的欧拉路时从 111 开始，如果在进栈时存储会出现 1 2 3 4 1，显然是不对的，因为 DFS 会连续回溯导致路径中间“断层”。但是如果在退栈时储存，答案是 1 4 3 1 2，因为存退栈路径不会有“断层”，退栈肯定是一步一步走的。...

2020-08-21 10:16:52 358 1

原创 [CodeForces 1119F] Niyaz and Small Degrees（堆优化多次 DP） | 错题本

文章目录题目分析代码题目[CodeForces 1119F] Niyaz and Small Degrees分析先考虑对于一个给定的 xxx DP：dp[u][0/1]dp[u][0 / 1]dp[u][0/1] 表示 uuu 结点的度数为 (x−1)/x(x - 1) / x(x−1)/x，且其子树合法的删边代价。dp[u][0]dp[u][0]dp[u][0] 就意味着 uuu 与其父亲的边保留，dp[u][1]dp[u][1]dp[u][1] 意味着 uuu 与其父亲的边删除。设 vvv 是

2020-08-04 17:07:53 265 1

原创 [AGC 004B] Colorful Slimes（DP） | 错题本

文章目录题目分析代码题目[AGC 004B] Colorful Slimes分析枚举位移的次数 kkk，显然 0≤k<N0 \leq k < N0≤k<N，接下来只需要知道每个史莱姆在什么时候被取走，因此答案为 min⁡0≤k<N{kx+∑1≤i≤Nmin⁡{ai,ai+1,⋯ ,ai+k}}\min_{0 \leq k < N} \left\{kx + \sum_{1 \leq i \leq N} \min\{a_i, a_{i + 1}, \cdots, a_{i

2020-08-03 11:28:10 173

原创 [洛谷 4317] 花神的数论题（数位 DP） | 错题本

文章目录题目分析错因代码题目[洛谷 4317] 花神的数论题分析考虑统计二进制含 jjj 个 111 的数有多少个，于是设计 dp[i][j][0/1]dp[i][j][0/1]dp[i][j][0/1] 表示前 iii 位中出现了 jjj 个 111，且前 iii 位比 nnn 小 /// 和 nnn 相等的数个数，转移刷表，分填 000 和填 111 即可。本题可以算个板子。错因DP 的时候不要取模。md 指数当然不能取模代码#include <algorithm>#i

2020-08-01 09:49:53 175

原创 [SPOJ 1026] FAVDICE - Favorite Dice（期望入门题） | 错题本

文章目录题目分析代码题目[SPOJ 1026] FAVDICE - Favorite Dice分析定义一回合表示经过若干次投掷，掷出一个之前没出现过的数字。考虑现在进行第 iii 回合，那么掷出一个之前没有出现过的数字的概率是 pi=n−i+1np_i=\frac{n - i + 1}{n}pi=nn−i+1，于是掷出一个之前没出现过的数字的期望次数为 Ei=1pi=nn−i+1\mathbb E_i=\frac{1}{p_i} = \frac{n}{n - i + 1}Ei=pi1=n−

2020-08-01 08:40:46 237

原创 [CodeForces 743E] Vladik and cards（二分答案 + 状压 DP） | 错题本

文章目录题目分析题目[CodeForces 743E] Vladik and cards分析

2020-07-31 08:20:15 173

原创 [CodeForces 165E] Compatible Numbers（状压 DP） | 错题本

文章目录题目分析代码题目[CodeForces 165E] Compatible Numbers分析简单状压 DP，将 S−{i}S-\{i\}S−{i} 的答案给 SSS 即可。代码#include <algorithm>#include <cstdio>#include <cstring>const int MAXN = 1000000;const int MAXA = 4000000;const int LOG = 22;int N;i

2020-07-30 14:57:50 212 1

原创 [AtCoder Petrozavodsk Contest 001F] XOR Tree（巧妙的转化 + 状压 DP） | 错题本

文章目录题目分析代码题目[AtCoder Petrozavodsk Contest 001F] XOR Tree分析一条路径上的边权全部异或一个值比较恶心，于是有一个神仙转化：考虑到路径上所有非端点的度都为 222，也就是说进入一个点和出去一个点都异或了一个值，所以我们将点权设为与它相连的边的边权异或和。容易证明所有点权均为 000 是所有边权均为 000 的充分必要条件。必要性显然。证明充分性只需要考虑不断找到度为 111 的点，其对应的边边权为 000，于是可以删掉它们，不断删掉后即可证明

2020-07-30 12:13:31 192

原创 [AtCoder ARC058C] 和風いろはちゃん / Iroha and Haiku（Hash + 状压 DP） | 错题本

文章目录题目分析错因代码题目[AtCoder ARC058C] 和風いろはちゃん / Iroha and Haiku分析显然只用存末尾的一段和小于等于 171717 的状态，先深搜一波。#include <algorithm>#include <cstdio>#include <cstring>int tot;void Dfs(int sum) { if (sum > 17) return; tot++; for (int i = 1

2020-07-30 11:10:11 635

原创 [HAOI 2015] 树上染色（树上背包） | 错题本

文章目录题目分析错因题目[HAOI 2015] 树上染色（树上背包）分析考虑 dp[u][i]dp[u][i]dp[u][i] 表示 uuu 子树中有 iii 个黑点时，对答案的贡献，转移过程中计算边 (u,v)(u, v)(u,v) 的贡献即可。注意贡献有白色贡献和黑色贡献。错因...

2020-07-29 19:40:51 167

原创 [CodeForces 1101D] GCD Counting（树形 DP） | 错题本

文章目录题目分析错因代码题目[CodeForces 1101D] GCD Counting分析将 aua_uau 分解质因数得到不同的质因子 p[u]1,p[u]2,⋯ ,p[u]tp[u]_1, p[u]_2, \cdots, p[u]_tp[u]1,p[u]2,⋯,p[u]t。dp[u][i]dp[u][i]dp[u][i] 表示从 uuu 向下只经过含因子 pip_ipi 的数能延伸的最长距离，转移的时候暴力枚举 x,y (p[u]x=p[v]y)x, y\ (p[u]_x

2020-07-29 15:51:21 195

原创 [CTSC1997] 选课（树上背包） | 错题本

文章目录题目分析代码题目[CTSC1997] 选课（树形背包）分析发现自己好像之前没做过树形背包其实是树形 DP 的基础上套一个对每个结点的背包。设 dp[u][i][j]dp[u][i][j]dp[u][i][j] 表示从结点 uuu 的前 iii 个儿子中选 jjj 门课的最大学分，设 cvc_vcv为 vvv 的儿子个数，可以得到转移 dp[u][i][j]=max⁡1≤k<jdp[u][i−1][j−k]+dp[v][cv][k]dp[u][i][j] = \max_{1 \leq

2020-07-29 14:55:01 379 1

原创 [BZOJ 4712] 洪水（动态 DP + 树链剖分） | 错题本

文章目录题目分析代码题目[洛谷 P6021] 洪水分析记 fuf_ufu 为 uuu 子树的答案，vvv 是 uuu 的儿子，则 fu=max⁡{∑ufv,au}f_u = \max\left\{\sum_{u}f_v, a_u\right\}fu=max{u∑fv,au} 接下来的过程与【模板】“动态 DP“&动态树分治几乎一致，找到转移方程的矩阵形式：(guau−∞0)×(fh0)=(fu0)\begin{pmatrix}g_u & a_u \\ -\infty &

2020-07-28 10:27:15 193

原创 [洛谷 P4719] 【模板】“动态 DP“&动态树分治（树链剖分 + 矩阵） | 错题本

文章目录题目分析代码题目【模板】“动态 DP”&动态树分治分析首先写出朴素的 DP 式：f[u][0/1]f[u][0/1]f[u][0/1] 表示不选 /// 选 uuu 点，uuu 的子树内的最大权独立集的权值和。设 vvv是 uuu 的儿子（下同），则 f[u][0]=∑vmax⁡{f[v][0],f[v][1]}f[u][1]=∑vf[v][0]+au\begin{aligned}f[u][0] &= \sum_{v} \max\{f[v][0], f[v][1]\} \\

2020-07-28 08:57:52 279 2

原创 [CodeForces 750E] New Year and Old Subsequence（动态 DP） | 错题本

文章目录题目分析错因代码题目CF750E New Year and Old Subsequence分析设 dp[i][0/1/2/3/4]dp[i][0/1/2/3/4]dp[i][0/1/2/3/4] 表示前 iii 个数字只能依次匹配到 ϕ/{2}/{2,0}/{2,0,1}/{2,0,1,7}\phi/\{2\}/\{2, 0\}/\{2, 0, 1\}/\{2, 0, 1, 7\}ϕ/{2}/{2,0}/{2,0,1}/{2,0,1,7} 的最小删除代价，则 dp[i][0]={dp[i−1

2020-07-27 13:03:31 217 1

原创 [SPOJ 1716 GSS3] Can you answer these queries III（动态 DP） | 错题本

文章目录题目分析错因代码题目Can you answer these queries III分析本题当然可以用传统线段树做记 fif_ifi 表示以 aia_iai 结尾的最大子段和，gig_igi 表示前 iii 个数的最大子段和。则 fi=max⁡{fi−1+ai,ai}gi=max⁡{fi,gi−1}=max⁡{fi−1+ai,ai,gi−1}\begin{aligned} f_i &= \max\{f_{i - 1} + a_i, a_i\} \\ g_i &= \m

2020-07-27 10:08:36 300 1

原创 [NOIP 2018 Plus] 奇怪的题（贪心 + DP） | 错题本

文章目录题目分析题目题目描述我们称一个长度为 nnn 的排列 {pn}\{p_n\}{pn} 为“match”，当且仅当对于任意 iii，pi≠i∧ppi=ip_i\not=i\land p_{p_i}=ipi=i∧ppi=i；称两个“match“不同，当且仅当对于任意 iii，pi≠qip_i\not=q_ipi=qi。设 p,qp,qp,q 是两个不同的”match“，给定序列 {an}\{a_n\}{an}，求：min⁡{∑i=1n∣ai−api∣+∣ai−aqi∣

2020-07-23 08:53:57 241

原创 [APIO 2018] Duathlon 铁人两项（圆方树 + 算贡献） | 错题本

文章目录题目分析错因代码题目[APIO 2018] Duathlon 铁人两项分析考虑已知 s,fs, fs,f，则 ccc 大致等于圆方树上 s→fs \to fs→f 的简单路径经过的点双的大小（即方点的度）之和。注意如果直接全部加起来会算多，因为圆方树上的非叶子圆点都是割点，即属于多个点双，而我路径上会将这些割点重复统计。解决办法很巧妙：方点的权值不变（所在点双的大小即度数），圆点的权值设置为 −1-1−1。这样一来统计 s,fs, fs,f 在圆方树上的简单路径的权值和即为已知 s,fs,

2020-07-22 20:51:43 222

TA关注的人

C++编写打飞机游戏

几何画板中文破解版+工具包

《父与子的编程之旅》Python安装包+源代码

C++编写扫雷

Markdown在小于号后面输字母后面的文本就隐藏了

C++怎么在运行时等待输入

GUIDE C++怎么输出中文？