woshidalaocxy的博客

原创 扩展欧几里得 学习笔记（扩欧 || exgcd 入门）

什么是欧几里得算法？简单点说，(a,b)=(b,a(a,b)=(b,a(a,b)=(b,a modmodmod b)b)b)什么是扩展欧几里得算法？扩展欧几里得算法是欧几里得算法（又叫辗转相除法）的扩展。除了计算aaa和bbb两个整数的最大公约数，此算法还能找到整数xxx yyy（其中一个很可能是负数）。通常谈到最大公因子时, 我们都会提到一个非常基本的事实: 给予二整数 aaa与 bb...

原创 洛谷 P3743 kotori的设备 题解

题目地址题目大意：kotori有n个设备，每个设备使用速度为aia_iai​每秒,一开始的电量为bib_ibi​，一旦有一个设备电量为0即停止使用，有一个充电宝可以连续给一个设备每秒冲ppp的电量，求最多这套设备能用多久，如果用之不竭输出-1数据范围：1≤n,p,ai,bi≤1e51\le n,p,a_i,b_i\le1e51≤n,p,ai​,bi​≤1e5读完题之后，我们会发现这道题珂以...

原创 [转载]CSP-S2019 树的重心

题面分析 原题面 小简单正在学习离散数学，今天的内容是图论基础，在课上他做了如下两条笔记： 一个大小为 nn 的树由 n ...

原创 $\max$卷积优化

现在有\[h_i = \sum _ {\max(j,k) = i} f_j\times g_k\]求\(h\)。设\[F_i = \sum_{j = 1} ^ i f_j\\G_i = \sum_{j = 1} ^ i g_j\]则\[\sum_{i = 1} ^ n h_i = F_n \times G_n\\h_n = F_n \times G_n - F_{n - 1} \...

原创 [ABC077D] Small Multiple

同余最短路，把增加的数位和作为代价，\(\times 10\)显然代价为零，然后加上\(i(0 \le i < K)\)的代价是\(i\)，这里的运算都是在\(\bmod K\)意义下。跑最短路即可。...

原创 P4630 [APIO2018] 铁人两项

题目大意给一个\(n\)个点\(m\)条边的无向图，求符合条件的三元组\((s,c,f)\)个数：存在一条\(s \to f\)的路径，经过点\(c\)，并且每个点最多经过一次。\(n \le 10^5,m \le 2 \times 10^5\)。Solution圆方树板题。建完圆方树后，枚举每个点为\(c\)，计算\(s,f\)的个数。考虑每个点的权值，对于圆点权值为\(-1\)，方...

原创 Atcoder ABC 266 EF

E题目大意有一个游戏，你可以玩\(n\)次，每次投一个骰子，若数字为\(X\)，则：若这把是第\(n\)把，那么你的分数为\(X\)，游戏结束否则，你可以选择继续游戏，或者立刻停止游戏，分数为\(X\)，游戏结束求最大的得分期望。\(n \le 100\)(???)Solution设\(f(d,x)\)为第\(d\)次游戏，骰子数为\(x\)的最大期望得分。\[f(d,x) ...

原创 BZOJ 4771 七彩树 | scz 8.5 逆向思维（强制在线）

题目大意给一颗\(n\)个节点的有根树，每个节点有一个颜色\(c_i(1 \le c_i \le n)\)，设节点\(x\)的深度为\(deep_x\)，回答\(m\)次询问，每次询问形如：x d 询问\(x\)子树里深度不超过\(deep_x + d\)的所有节点的颜色种数\(T\)组测试，每组测试相互独立，强制在线。\(T \le 500,n,m \le 1 \times 10^5...

原创 Atcoder ABC 265DEF

D题目大意给定一个序列\(A = (A_0,\cdot,A_{N - 1})\)，判断是否能找到一个四元组\((x,y,z,w)\)满足：\(0 \le x < y < z < w \le N\)\(\sum_{i = x} ^ {y - 1} A_i = P\)\(\sum_{i = y} ^ {z - 1} A_i = Q\)\(\sum_{i = z} ^ ...

原创 COCI 2021-2022 Round 2 选做

A.Kaučuk根据题意模拟即可。B.Kutije对于每个\(p_i\)，连一条\(i \to p_i\)的有向边，这样对于每次询问\(a,b\)，相当于询问\(a\to b\)是否联通。传递闭包bitset优化即可。时间复杂度\(\mathcal{O}\left(\dfrac{n^3}{64}\right)\)...

原创 Atcoder ABC 263E 期望，数学

题意有\(n\)个地方，编号为\(1\sim n\)，每个地方有一个骰子，骰子上标有整数\(0,1,\cdots , A_i\)，一个人在\(i\)掷骰子到\(j\)，那么他会走到编号为\(i+j\)的地方。若一个人不在编号为\(n\)的地方，那么他会一直投骰子。求投骰子的期望次数。\(n \le 2 \times 10^5,A_i \le n - i\).Solution根据套路，设\(...

原创 scz 7.30 B

A.regia线段树随意维护。全场唯一挂分B.papy我们沿着这个词从左到右迭代，并保留我们还没有配对的字母数组。假设我们遇到一个字母A.我们将检查最后一个不成对的字母是否为A：如果是，我们可以将其与新的字母配对，并将其从不成对的字母数组中删除。如果最后一个不成对的字母是B，我们必须将新的字母A添加到不成对的字母数组(我们不能与以前的一个A配对)，因为如果我们尝试连接B后，它要形成的弧线将...

原创 多组询问不同指数询问$a^x$的值的trick

一般做法：对于每个询问\(x\)，用快速幂求出\(a^x\)。复杂度\(\mathcal{O}(T\log{x})\)优化：考虑到\(a\)为定值，每次只改变\(x\)的值，设\(x\)的最大值为\(N\)，设\(A = \sqrt{N}\)，预处理出\(a^0,a^1,\cdots,a^{A - 1}\)和\((a^A)^1,(a^A)^2,\cdots,(a^A)^A\)。对于每次询问，易...

原创 scz 7.28 B

A.number给一个 \(n\) 位数，将每一位重新排序，使得新的数尽量大，并且为 \(30\) 的倍数，如果没有 \(30\) 的倍数，输出 \(-1\)。\(n \le 10^5\)A.solution若一个数是 \(3\) 的倍数，那么它的各位数和也是 \(3\) 的倍数。末位是 \(0\) 的数是 \(10\) 的倍数。先判断各位数和是否为 \(3\) 的倍数；再确定末位...

原创 Codeforces Round #805 (Div. 3)（A~F，G待更）

ABCDEFG+++2+++-2A签到题不讲。B模拟题，每次要记录到的字母超过3个时就加一天。时间复杂度\(\mathcal{O}(\sum|S|)\)C不难发现要把\(u\)先离散化一波，设车站\(x\)所处最靠前的位置为\(Min_x\)，最靠后的为\(Max_x\)，每次询问\(a_i,b_i\)时，若\(Min_{a_i} \l...

原创 CSP-S2 2021 场外

说在前面由于我初赛太逊了，导致我今年没去S2。不过我在复赛过后还是看了卷子，以下内容大概都是在独立思考下想出来的。T1感觉这个题比排水系统/julian良心。不难发现飞机肯定是优先排在前面，那么也就是说我们先求出来有无限多廊桥的情况，然后作前缀和即可，国内和国际单独做，最后枚举0~n作为国内廊桥个数取max。那么无限多廊桥的情况怎么弄呢？一开始属实把我整不会了。后来想想，主要是解决几个...

原创 P1965 [NOIP2013 提高组] 转圈游戏

Problem给定\(n,m,k,x\)，\(x\)每次会变成\((x + m) \bmod n\)，称为1次变换，求经过\(10^k\)次变换后\(x\)的值。\(n \le 10^6,m < n,k,x \le 10^9\)。Solution看见\(n\)数据范围显然可以想到整循环节，但是我们不会推，咋办，发现求循环节至多\(\mathcal{O}(n)\)，求完之后相当于求\(...

原创 2021 暑期总结

说在前面这个暑假大概是我过的最充实的一个暑假。整个暑期时间可以分三段：强基前，强基后 ~ 8.27,8.27 ~ 结束。Part 1强基前基本上每天就是刷数物的阶段。分类精粹和小蓝书几本有用的基本上都做了一遍(但事实证明小蓝书好像做了没考)，文科基本可以说是没有复习(因为我也不知道怎么搞)。考前非常自信地认为会考组合数学(flag)和数论(flag)，所以考前甚至还做了不少竞赛类型的数...

原创 牛客IOI周赛28-普及组 题解

说在前面被毒瘤题虐了，看到昨天牛客有普及模拟赛，来vp。大概花了1h AK。比赛链接A发现做完\(n\)次操作等于没做，所以做\(x\)次操作等同于做\(x \bmod n\)次操作，然后先输出[x + 1,n]，再输出[1,x]即可。# include <bits/stdc++.h>using namespace std;const int N = 1e5 + 5;...

原创 P4198 楼房重建

Problem\(1 \le X_i \le N,1 \le Y_i \le 10^9,1 \le N,M \le 10^5\)Solution定义\(\text{slope}(i)\)表示\((0,0)\)与\((i,H_i)\)的连线的斜率。不难发现，如果第\(i\)个楼要被看见，当且仅当\(\text{slope}(i)\)是前\(i\)个最大的。我们把\(\text{slope}...

原创 P5278 算术天才⑨与等差数列

Problem给一个序列，要求支持2种操作：单点修改查询\(a_l \cdots a_r\)从小到大排序是否是一个公差为\(k\)的等差数列\(n,q,k \le 3 \times 10^5\)Solution看见题目，感觉正解不太会，想通过多种等差数列的限制乱搞做题。第一个想到的是区间\(\min\)和区间\(\max\)。\(\max - \min = (r - l)k\)才...

原创 P6327 区间加区间sin和

Problem维护一个序列\(a\)，每次操作有两种类型：1 l r v 将\(a_l \sim a_r\)的所有数加\(v\)。2 l r 求\(\sum_{i = l} ^ r \sin (a_i)\)。输入数据均小于\(2 \times 10^5\)。Solution上午刚做完一个Ynoi，以为是最简单的Ynoi题了，结果这题更是重量级。不难发现通过两个诱导公式可以做：...

原创 分层图总结

分层图分层图主要就是通过将图分为几层，每层内部都有边，然后各层之间也有边，用于解决对边权进行操作且操作数少的题目。Example给一个无向图，求\(1 \to n\)的最短路径，其中可以选择不超过\(k\)条边让边的边权变为一半。 \(n \le 50,m \le 1000,k \le 50\)。考虑建立\(0,1,2,\cdots,k\)层共\(k + 1\)层图。定义\(P(x,...

原创 P7515 [省选联考 2021 A 卷] 矩阵游戏

原创 P4745 [CERC2017]Gambling Guide

Problem\(1 \le n,m \le 300000\)Solution考虑设\(dp_x\)为\(x\)到\(n\)的最小硬币期望值。\(dp_n = 0\)\(dp_x = \left(\sum \min\{dp_v,dp_x\}\right) \cdot \dfrac{1}{deg_x} + 1\)然后搞搞。...

原创 P3573 [POI2014]RAJ-Rally

Problem给一个\(n\)个点，\(m\)条边的DAG，找到一个点，使得删去这个点后的最长路径最短。\(2 \le n \le 500000,1 \le m \le 1000000\)。Solution看到DAG就想到拓扑。求出以\(x\)为起点的最长路径长度\(ds_x\)和以\(x\)为终点的最长路径长度\(dt_x\)。可以通过跑正反图拓扑实现。设第\(x\)个点的拓扑序为\...

原创 P5304 [GXOI/GZOI2019]旅行者

Problem给一个\(n\)个点，\(m\)条边的有向图。给\(k\)个特殊点\(K_1,K_2,\cdots,K_k\)求\(k\)个特殊点中两两最短路的最小值。数据范围：SolutionThinking 1Floyd...好像暴力都没得打。/kkThinking 2DAG!可以拓扑搞。然而并不会！Thinking 3思路真的很妙。建超级起点\(S\)和超级终点\(T...

原创 最短路径树（Shortest Path Tree,SPT）

什么是SPT？SPT是原图的一颗生成树。具有生成树的所有性质。从根节点到任意节点的距离与其在原图的最短路距离相等。如何求出SPT结合单源最短路径算法。这里使用堆优化Dijkstra.记录一个\(pre_i\)表示\(i\)的上一条边。这样只要有根就可以遍历出整棵树。void dij(void){ priority_queue <pair<int,int&gt...

原创 P4159 [SCOI2009] 迷路

Problem给一个\(n\)个节点\(m\)条边的带权有向图。求从\(1\to n\)的长度为\(t\)的路径条数。对\(2009\)取模。\(2 \le n \le 10,1 \le t \le 10^9\)值得一提的是，此题的输入格式：说明边权\(\le 9\)。SolutionThinking 1考虑边权都为\(1\)怎么做。设\(f[t][i][j]\)为长度为\(t\)...

原创 P6617 Search

Problem给一个长度为\(n\)的序列\(a\)和一个正整数\(w\)，\(m\)次操作，每次操作为：1 x y 修改\(a_x\)为\(y\)2 l r 查询是否存在\(i,j\)满足\(l \le i < j \le r\)，且\(a_i + a_j = w\)。\(1 \le x \le n,0 \le y \le w,1 \le l \le r \le n \le ...

原创 P1850 [NOIP2016 提高组] 换教室

ProblemSolution设\(dp[i][j][0/1]\)表示在前\(i\)节课中，申请\(j\)节，同时第\(i\)节课申不申请的期望值。首先需要Floyd预处理出各个教室之间最短路。不难发现\(dp[i][..][..]\)一定由\(dp[i - 1][..][..]\)转移来。分类讨论：\(i\)不申请\(i - 1\)不申请\(dp[i][j][0] = dp...

原创 P2831 [NOIP2016 提高组] 愤怒的小鸟

Problem平面上有\(n\)个在第一象限的点，求至少要用几个类似于\(y = ax^2 + bx(a < 0,a,b \in \mathbb{R})\)的抛物线能将其全部覆盖。\(n \le 18\)SolutionThinking 1\(1 \le n \le 18\)，状压dp.设\(dp[S]\)为S状态下至少要用几个抛物线。Thinking 2发现抛物线必经\((...

原创 P4211 [LNOI2014]LCA

Problem给一棵\(n\)个节点的有根树，\(m\)次询问，每次询问给\(l,r,z\)，求\(\sum_{i = l}^r dep[LCA(i,z)]\)\(1 \le n,m \le 500000\)SolutionThinking 1我会暴力！树剖LCA，\(\mathcal{O}(n\log{n} + nq)\)属于是。你要树剖LCA不会只会倍增的话那还得加个\(\log{...

原创 P4137 Rmq Problem / mex 强制在线做法

Problem给定一个数组\(a\)，每次询问给定一个区间\([l,r]\)，求区间\(\operatorname{mex}\)。\(n \le 2 \times 10 ^ 5,a_i \le 10 ^ 9\)。Solution考虑用主席树做。每个节点记录它代表的区间权值在当前最早出现的位置。查询的时候直接在\(Root_r\)上查询$ < l$的即可。# include <...

原创 P2272 [ZJOI2007]最大半连通子图

Problem当一个有向图\(G= (V,E)\)满足：\(\forall u,v \in V,u \to v\)或\(v \to u\)。求这个图的最大半联通子图(节点数最多)节点数和最大半联通子图个数模\(X\)的值。\(n \le 10^5,m \le 10^6,X \le 10^8\)Solution显然SCC也是强连通子图。所以先用Tarjan缩点，然后将图变成一个DAG(套路...

原创 P5664 [CSP-S2019] Emiya 家今天的饭

很久之前的题今天才做。Solution发现每种主要食材不能超过\(\left\lfloor \frac{k}{2} \right\rfloor\)，那考虑不合法的情况，有也仅有一种主要食材会出现\(> \left\lfloor \frac{k}{2} \right\rfloor\)的情况。于是，我们考虑如何求出总情况数和不合法情况数，然后再相减就是答案。考虑枚举不合法的主要食材\(...

原创 全程NOIP计划 模拟赛1

说在前面这场比赛打炸了，听说是套路题，但是除了会骗点t1的分啥都不会。预期：100+75+10+10实际：90+36+0+10(捆绑测试\fn)yzhang和dead_X差点ak，orz。A 平面几何考察：写代码能力，小思维(？)我考场上尝试\(\mathcal{O}(n^2)\)(没跑满)+极限卡常，可惜还是失败了。考虑作一条\(Y = a\)的直线，其中\(a\)是一个很大...

原创 AtCoder Beginner Contest 211

ScoreABCDEFACACAC(1)AC//A送分题。B送分题。C设\(dp[i][j]\)为\(S\)的前\(i\)个字符组成目标串的前\(j\)个字符的总数，则：\[dp[i][j] = \begin{cases}\sum_{k < i,S[k] = St[j - 1]} dp[k][j - 1] & (S[i]...

原创 P2569 [SCOI2010]股票交易

Problem初始有\(\infty\)元钱，有\(T\)天，每天可以买卖股票，告诉你第\(i\)天股票买入价\(AP_i\)，股票卖出价\(BP_i\)(保证\(AP_i \ge BP_i\))，至多买\(AS_i\)股，至多卖出\(BS_i\)股。股票交易至少要间隔\(W\)天，每天一个人手里的股票数不得超过\(MaxP\)。求\(T\)天后最大赚钱数。\(1 \le BP_i \le A...

原创 全程NOIP计划 题目选做

Day 1P1314二分答案P2367差分板子P2882找性质，然后模拟，异或差分维护。P1719转换成最大子段和问题。P3017二分，先考虑水平方向，定两个点\(l,r\)，在\([l,r]\)行内，再挂两个点\(l_0,r_0\)，表示当前竖直方向的区间，如果矩阵和\(\ge mid\)，那么继续，否则继续延伸。最后若分的块数\(< B\)，则不行。水平同理。P1...