- 博客(415)
- 收藏
- 关注
RSS订阅原创 【GDOI2019模拟2019.4.13】数据结构
Description:1<=n,m<=1e5时限：-O2,1s空限：128MB题解：把询问区间的左端点视为x，右端点视为y。转换模型后不难发现就是对二维平面上一些东西做一个东西。设s表示一段时间内这个区间和的增量，v表示和的增量的历史版本最小值(v>=0)。这个玩意显然是可以合并的，即有两个相邻的时间段的s和v，就可以合并出新的。那么用二维线段树去维护这个吗...
2019-04-13 21:38:00
364
原创 2019十二省联考（SH） 游记
考完一段考终于可以出去玩了，吼鸡冻!这篇比赛主要是游记，题解待我研究透题目后会写（D2T3不知道会咕咕到什么时候）。Day0：早上到了学校后发现我居然没有带身—份—证于是冲回家再直奔机场，一个小时好算是赶上了。飞机上唯一的感受是空调好冷，飞机餐要了两份居然吃饱了。到酒店时还挺早的，瞄了眼房价牌感慨魔都物价真的高，三星酒店居然一晚……700+，帝都同地段四星也就这个价。G2的跑去外滩浪...
2019-04-11 22:27:59
421
原创 单位根反演学习小记
裸题：bzoj3328在学习这个之前建议对FFT中的求和引理有一定的了解。单位根反演就是对：∑i∣kCni(mod p)\sum_{i|k}C_{n}^i(mod~p)∑i∣kCni(mod p)进行变换，p是质数，k∣p−1k|p-1k∣p−1，n比较大设mod pmod~pmod p下的原根为g,w=g(p−1)/kw=g^{(p-1)/...
2019-04-11 21:49:17
641
原创 半平面学习小记
从前的我看zzy那篇论文，真的是麻烦，后面发现其实半平面交其实是只比凸包难一点的东西。在学习半平面交之前，你需要对向量计算几何有初步的认识，还有就是atan2(y,x)函数的运用。atan2(y,x)函数见这篇博客：https://blog.csdn.net/chen825919148/article/details/7582510网上的这篇博客应该是比较易懂的：https://blog...
2019-04-01 20:53:44
491
原创 agc031 C Differ by 1 Bit
传送门：https://atcoder.jp/contests/agc031/tasks/agc031_c题解：考虑把x xor y，问题变成从0->x xor y因为2n−12^n-12n−1是奇数，所有x xor yx~xor~yx xor y有偶数个1无解。由于我们任意交换一下位是没有关系的，所以考虑就是要构造出：000…−&...
2019-03-18 21:17:29
413
原创 CodeForces 1056 G（非旋可持久化treap）
传送门。题解：设Ans（x,t）Ans（x,t）Ans（x,t）表示在x点，时间为t的答案假设我们能预处理出所有Ans(x,n)Ans(x,n)Ans(x,n)，那么这题就很简单了。考虑倒着想，假设现在倒着走了若干步，f[i]表示i到了哪里。每次相当于从f中copy四段再接起来，那么用可持久化treap维护。然后你发现treap的深度炸了：Samjia的解释是分出的四段有相同的部分，...
2019-03-17 22:50:53
524
原创 多次剩余原根法
多次剩余：求xa=b(mod p)x^{a}=b(mod~p)xa=b(mod p)，p是质数的所有解。找到p的原根g，设gB=bg^{B}=bgB=b，这个用bsgs求。xa=b(mod p)x^a=b(mod~p)xa=b(mod p)等价于ga∗x=gB(mod p)g^{a*x}=g^{B}(mod~p)ga∗x=gB(mod&nb...
2019-03-10 16:37:17
297
原创 二次剩余Cipolla法
翻着翻着居然翻到了两个退役选手zy和mw的博客，想到自己可能也会像他们一样省选失利……由于我比较水，目前只会模数是质数的情况，之后会了可能会补。很多具体的证明我也不会，限于感性认识，所以请看zy的博客：https://blog.csdn.net/a_crazy_czy/article/details/51959546二次剩余要做什么：给出n,pn,pn,p，p为一质数。要求xxx,...
2019-03-09 09:02:06
326
原创 【NOI2018】你的名字
题目大意：一开始给出字符串S。Q次询问，每次给出l,r和字符串T。问T有多少不同子串不是S[l…r]的子串。∑∣T∣&lt;=1e6,∣S∣,∣T∣&lt;=5e5\sum|T|&lt;=1e6,|S|,|T|&lt;=5e5∑∣T∣<=1e6,∣S∣,∣T∣<=5e5题解：不同子串一类问题优先考虑SAM。感觉这题考的更多是对SAM的理解。...
2019-03-04 16:43:47
971
原创 2019.3.2 B题 排行（交互）
题目大意：今年是 3345 年，你打算还原千年前这场比赛的名次，但是主办方称由于技术原因排名遗失了，连用于排名的 a 也丢失了。所幸，在 Byteland 中生活的都是机器人，所以你可以询问这些千年前的参赛选手。为了谨慎起见，你可以每次询问一个参赛者，某一个参赛者的比赛成绩比它好还是比它差。但是，机械心理学家告诉你，这些选手不一定愿意回答你的提问。具体地：名次小于 1 的选手由于耿耿于...
2019-03-02 19:41:57
517
原创 自适应辛普森积分法
推荐博客：浅谈自适应Simpson法今天做比赛时遇到一道题，就是求两个球的体积并，在XHM的指导下，我临时学习了这个方法并AC掉了这题（我知道可以百度到公式）。既然这玩意儿有公式，那么我们不妨换个题面：求平面上若干个圆的并的面积，这总没有公式了吧。我们想一下积分中最正常的思想，即把整个x轴划分成若干（->∞）段，然后每一段用左边的点的值去替代整一段的值。划分为无限段的时候，那就是答...
2019-03-01 22:50:41
1417
原创 原根相关
原根，定义为，对一给定模数p（p是奇质数），若有一g，满足g1−p−1g^{1-p-1}g1−p−1两两不同，则称g为p的原根。不会证明的结论QAQ：奇质数p一定存在原根，而且通过打表我们可以发现最小的原根都很小。task1：如何判断一个数x是否为p的一个原根？费马小定理：xp−1=1(p为质数)x^{p-1}=1(p为质数)xp−1=1(p为质数)假设x1−p−1x^{1-p-1}x1...
2019-02-28 22:51:30
918
原创 JZOJ 幽幽子与森林
题目大意：迷途竹林可以看成是一个n个点的森林，幽幽子定义dis(u,v)为u到v路径上的边的数量，若u和v不连通则为m。她定义整个森林的危险度为为了去拜访永琳师匠，幽幽子需要提前知道迷途竹林的危险度。但迷途竹林的形态是时刻变化着的，所以幽幽子希望知道危险度的期望是多少。为了避免浮点运算，答案对998244353取模。1<=n<=2e5题解：这个计数写得真的是累。很容易想到...
2019-02-28 16:48:15
13324
原创 【GDOI2019模拟2019.2.23】字符串
Description：题解：首先思考一个经典的问题怎么做？就是求一个串的子串的不同子串数，这玩意儿冬令营上《字符串前沿算法》讲过。考虑在把整个串的SAM建出来，注意建的时候记录每个结尾点。然后一个一个把结尾点激活，之所以要先建出来的原因是树的形态就稳固了。假设要查询[l…r]的不同子串数，先把r以前的结尾点激活。如果知道每个点的rmax，这个点的深度代表的是[x…y]，那么就给...
2019-02-24 22:35:54
307
原创 【GDOI2019模拟2019.2.23】签到
题目：题解：设是a[i]的倍数的限制为a，是b的倍数限制为b，f=1的限制为c。相当于求：a|b|c=a+b+c-a&b-a&c-b&c+a&b&c发现我们求a并不好求，考虑把所有的a改成!a=(n-!a)+b+c-(b-!a&b)-(c-!a&c)-(b&c)+(b&c-!a&b&c)=n...
2019-02-24 18:33:59
425
原创 CodeForces 1110 F Nearest Leaf
传送门.翻译：给出一棵带权树，树的dfs序是确定的。每次询问给出x,l,r，问dfs序属于[l…r]的叶子点到x最近的距离是多少。题解：不难想到一个暴力的做法：每次x往父亲跳，查询其它子树里的叶子节点。复杂度O(n2)O(n^2)O(n2)可以强行套上树链剖分搞这个，时间复杂度为(nlog2n)(n log^2n)(nlog2n)，空间复杂度通过离线询问可以达到O(nlogn)O(...
2019-02-22 21:13:09
204
原创 CodeForces 1110 D. Jongmah
传送门.看错题交了11发才过。给出n个数，分最多组，使得每组标号相邻递增或全部相同。1<=n<=10^6首先自然地想到dp设前两个有多少个，当然这个dp会超时。不难想到如果一个数出现的次数比较多，那么应该先把他自己取了。那么这个数具体是多少呢？如果我们取了（a,b,c）三次，是可以转成取a,b,c分别三次。那么一个数最多作为前中后分别2次，最多的次数就是2*3=6次也...
2019-02-20 22:14:13
485
原创 CodeForces1109D Sasha and Interesting Fact from Graph Theory（扩展Cayley公式）
传送门.题目翻译：给出n,m,a,b要求生成一棵有标号有边权的n个点的树，边权∈[1…m]是a->b的距离=m问方案数，对1e9+7取模。n,m<=10^6题解：很容易想到枚举a,b之间的边数设为x，然后开始统计：首先要选x-1个点出来，P(n−2,x−1)P(n-2,x-1)P(n−2,x−1)边的和为m，挡板问题（stars and bars method），C...
2019-02-17 22:39:14
510
原创 WC2019打铁记
真想不到就这样打铁了，快乐至极。讲课到现在已经忘的差不多了，印象里只有一场杂题选讲听懂了个大概，其它的断线重连失败，有些还是有价值的。考试日：刚T1刚了三个小时，无果，只有暴力48分T2搞了半小时的编译，然后捞了前3个点就走了。T3一开始以为grader是赋个权值然后比较，结果发现下发的grader是随机的，然后很自然的以为正式的grader也是这样的，于是写了个随机下的做法，顺利爆零。...
2019-02-03 22:31:10
650
原创 Codeforces 1103B Game with modulo
传送门.翻译：交互有一固定a，每次你可以询问，给出参数x,y，返回[x%a>=y%a]求a，询问次数<=60次。题解：B题这么难可以退役了。60次那肯定是什么二分乱搞。假设现在确定a∈[x…y]设mid=(x+y)/2我们模拟一下会发现做不了。其难点在于2*x不一定大于mid于是不如倍增（那些10分钟切的人怎么想得这么快？），先查询[2i,2i+1][2^i,2...
2019-01-23 10:59:43
348
原创 PKUWC2019游记
因为在本校举办，THU那边又不让去，所以作为“体验最菜营员”来划水。day0：两点钟坐从XL坐轻轨到NL接人，发现我居然和要接的人在同一班车上。于是下了车马上屁颠屁颠跑到出口去等他们。送走这班人，发现下一班人3H后来，于是和JY商量去看电影。适合看的只有《大黄蜂》，评分有9+，然并卵。还是很快乐地度过了3H。晚上CF...
2019-01-20 19:47:03
727
4
原创 2019.1.17模拟
T1:求∑i=1nti∗Cn−ik−1\sum_{i=1}^nt^i*C_{n-i}^{k-1}∑i=1nti∗Cn−ik−11&lt;=k&lt;=1e7,1&lt;=n&lt;=1e91&lt;=k&lt;=1e7,1&lt;=n&lt;=1e91<=k<=1e7,1<=n<=1e9也不知道怎么搞...
2019-01-17 18:46:18
305
原创 常系数齐次线性递推学习小记
此篇文章可能写的很不标准，请读者见谅。常系数齐次线性递推，可以理解为dp，大概是给出f[0..k−1]f[0..k-1]f[0..k−1]的初值，对于i&amp;amp;amp;amp;amp;gt;=k,f[i]=∑j=1k−1f[i−j]∗c[j]i&amp;amp;amp;amp;amp;gt;=k,f[i]=\sum_{j=1}^{k-1}f[i-j]*c[j]i&amp;amp;amp;amp;gt;=k,f[i]=∑j=1k−1f[i
2019-01-16 16:42:54
354
原创 【PKUWC2019模拟2019.1.15】Square（结论题）
Description:题解：怎么说呢，一道联赛题没有做出来。以后还是要在一开始就想清楚，后面真的是越想越迷。首先需要判断答案是否存在，这非常简单，每次找到一行或一列全是一样颜色的，然后把他们变成通配颜色，如果矩阵能全部变成通配颜色，就可行。不然发现要不每一行都染了，要不每一列都染了，不然会有点染不到。假设每一行都染了，我们不用关心每一行染了什么，只用看最多有多少列完全相同，行的染色...
2019-01-16 16:11:17
336
原创 【WC2019模拟2019.1.12】二分的代价（状态互换dp）
Description:题解：非常容易想到一个O(n3)O(n^3)O(n3)的dp，于是做不动了。考虑一个数变大，这个序列的答案只会更大，所以答案不会超过9∗log n9*log~n9∗log n因此状态互换。设fi,jf_{i,j}fi,j表示从i开始，在花费j的代价最远能走到哪里。枚举一个点p：fi,j=max(fi,j−1,[fi,j−a[p]&...
2019-01-15 22:05:39
454
原创 JZOJ 5990. Bear（斜线dp）
Description：n<=12,m<=30题解：很容易想到一种一行一行的dp，我们觉得n太小了，想一列一列做，发现不行。然后脑洞打开，直接斜着dp，发现就没有了。Code：#include<cstdio>#include<cstring>#define pp printf#define ll long long#define fo(i,...
2019-01-15 21:55:06
261
原创 【WC2019模拟2019.1.14】选数
Description：题解：有gcd那先反演掉。∑∑gcd(a,b)\sum\sum gcd(a, b)∑∑gcd(a,b)=∑dd∗[d∣a]∗[d∣b]∗[(a/d,b/d)=1]=\sum_{d}d*[d|a]*[d|b]*[(a/d,b/d)=1]=∑dd∗[d∣a]∗[d∣b]∗[(a/d,b/d)=1]=∑dd∗∑d′μ(d′)∗∑dd′∣a∑dd′∣bI=\sum_...
2019-01-15 21:48:59
607
原创 FWT模板
我觉得那几个证明真的有点猥琐，以后有时间再搞懂吧（或许搞懂了也不会补的）推荐博客：yyb我觉得除了异或，结论都可以感性出来，异或也很好记。IFWT的话就是直接反解，不需要记。or：FWT(A)=(FWT(A0),FWT(A0+A1))FWT(A)=(FWT(A0),FWT(A0+A1))FWT(A)=(FWT(A0),FWT(A0+A1))IFWT(A)=(IFWT(A0),IFWT(...
2019-01-15 21:31:31
488
1
原创 2019.1.15模拟
T1:给你一个长度为n的序列a，求一个最小n的排列p，使maxi=1n−1(a[p[i]] xor a[p[i+1]])\max_{i=1}^{n-1}(a[p[i]]~xor~a[p[i+1]])maxi=1n−1(a[p[i]] xor a[p[i+1]])最小。题解：不知道为什么想了这么久就是不觉得这题是个贪心。按照位运算贪心的基本思想，肯...
2019-01-15 20:42:09
290
原创 【北大2019冬令营模拟2019.1.1】多边形
Description:1&lt;=T&lt;=105,1&lt;=k,m&lt;=n&lt;=1061&lt;=T&lt;=10^5,1&lt;=k,m&lt;=n&lt;=10^61<=T<=105,1<=k,m<=n<=106题解：这题很好的考察了观察能力、组合数乱搞
2019-01-03 22:50:29
404
原创 浅谈我对（扩展）中国剩余定理Crt的理解
我感觉中国剩余定理我并不能顺推出来，这是一个很巧妙的构造，这个构造思想值得学习，但是扩展的可以顺推。问题：有一xxx，满足：x≡c1(mod&amp;nbsp;m1)x≡c1(mod~m1)x≡c1(mod&amp;nbsp;m1)x≡c2(mod&amp;nbsp;m2)x≡c2(mod~m2)x≡c2(mod&amp;nbsp;m2)…x≡ck(mod&amp;nbsp;mk)x≡ck(mod~m
2019-01-03 12:23:39
746
原创 浅谈我对（扩展）Lucas定理的理解
可能是太太太久没有用这个定理了，考场上都不记得有这个东西，然后自己推了个组合数取模（模数是一个质数p，且p比较小），考后发现就是Lucas定理，可以说是顺推出了Lucas定理，而不是有了结论再去证明。目标：求Cnm(mod&amp;amp;nbsp;p)C_{n}^m(mod~p)Cnm(mod&amp;amp;nbsp;p)，p是质数，且比较小。step1：考虑如何求n!(mod&amp;amp;nbsp;p)n!(mod~p)n...
2019-01-03 10:29:59
535
原创 有上下界的网络流
不知道为什么一直不会这个东西，可能是太懒了无源汇上下界可行流：我们知道最基本的网络流就是有上界的网络流，那么加了一个下界有什么不一样呢？一个大胆的想法：默认每条边流到了下界，然后把边的流量限制设为上限-下限，直接跑。这显然是不对的，不然可行流不就一定存在了吗？我们发现这样不能满足网路流的流量守恒。如果默认每条边流到了下界，那么对于一个点来说，它就可能不满足流量守恒了。设q[i]q...
2018-12-28 16:32:43
297
原创 Codeforces 1087 E. Vasya and Templates
传送门.真的是被细节恶心到了，是不是我的方法太麻烦了？要直接s’>=a且s’<=b显然是不好做的，因为要同时考虑两个限制。于是就要睿智，求>=a的最小的s’。那么这就很水了。高位要尽可能的相等，如果不行的话，就回溯，找到最右的一个位置，使它有大于a[x]的选择，然后再把后面填最小的了，注意没有用过的字符给它们匹配一下。Code：#include<cstdio...
2018-12-26 21:56:05
547
原创 Codeforces 1087F Rock-Paper-Scissors Champion
传送门。题解：比赛要回宿舍去睡觉没时间做了。GD选手表示只会强行数据结构。首先先来猜一下结论。对一个人，左右分开考虑，只有都能赢，才能赢。一边能赢的条件是，要么不存在克制它的，要么存在克制它的同时存在它克制的。这个是很显然的。既然带修，强行套上一个线段树，对于线段树的每个区间维护f[l][r]表示左边传过来的状态是l，右边传过来的状态是r，这个区间有多少个可行的，状态是二进制状态，...
2018-12-25 21:46:37
314
原创 Samjia Contest 12.22
传送门：https://floj.tech/contest/449题解：注意到每次只会把一段往前copy，那么一个点每次都会只会向前走，假设现在x∈[yi,yi+1)，hi&amp;gt;0,那么x一次向左会跳yi-hi的距离，直到跳到yi左边。这个可以直接除一下来得到跳出这段后的位置，倒着扫，x最终会定在一个点上。如果两个点最终汇聚在一个点上，它们就是一样的。时间复杂度O((a+q)b)O(...
2018-12-24 21:57:26
392
原创 JZOJ 4583【NOI2016模拟7.5】求和
1<=n,m<=500000题解：为什么这么多题目的名字叫求和？这题其实有很多种做法，但是用一个非常简单的性质就可以做出来。。。s(pq)=∑i∣p∑j∣q[(i,j)=1]∗i∗q/js(pq)=\sum_{i|p}\sum_{j|q}[(i,j)=1]*i*{q/j}s(pq)=∑i∣p∑j∣q[(i,j)=1]∗i∗q/j证明之前写过：https://blog....
2018-12-21 21:09:18
199
原创 【清华2019冬令营模拟12.15】排列
Description:n<=5000题解：考虑对一个值i它的区间的长度是多少。区间可以只考虑右半部分的，左边的倒过来做一遍就行了。若i是固定的，找到右边第一个固定的j>i那么区间肯定不能超过j所在的位置。可以直接枚举区间的右端点，这个右端点选的大于i，它左边的小于i，剩余的乱排，用排列数算即可。这一部分O(n2)O(n^2)O(n2)若i不是固定的，找到整个序列固...
2018-12-15 15:27:45
1465
4
原创 SnackDown 2019 - Online Elimination Round Adi and the Matrix(二维Polya计数)
原题链接.刚了这么久就是没有往burnside靠。题目大意：求有多少个不同的n*m的01矩阵。两个矩阵相同，定义为：通过每次交换两行或两列若干次之后完全相同。n*m&lt;=550,答案模1e9+7题解：burnsideburnsideburnside引理告诉我们：不同数=∑每种置换不动点个数置换数不同数={\sum_{每种置换}不动点个数\over 置换数}不同数=置换数∑每种置...
2018-12-10 22:33:03
379
原创 SnackDown 2019 - Online Elimination Round Suffix Palindromes
原题链接.没上OEIS真的是亏了。题目大意：问有多少个字符集大小为|S|的长度为n的字符串，不存在一个长度大于1的回文后缀。1<=n<=10001<=|S|<mo=1e9+7题解：这种题一看就要容斥原理首先总方案数是S^n减去后缀是回文串的，但是我们发现直接减肯定减重。不妨设f[i]f[i]f[i]表示长度为i的回文串，且它也不存在长度大于1的回文后缀的...
2018-12-10 22:04:59
354
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人