- 博客(43)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 LOJ #6360. 复燃「恋之埋火」（最小高维圆覆盖）
题目前置知识1：最小圆覆盖前置知识2：nnn维圆上n+1n+1n+1点求圆心然后就可以解决此题了。首先最小圆覆盖的时间复杂度分析套用至高维是没有问题的。其次可以同样证明按照每次都从小到大的顺序枚举点是不会出现让你求共线的三个点，共面的四个点之类的外接高维圆的。最后求nnn维kkk个点都在圆上的最小圆覆盖，把这kkk个点互减得到k−1k-1k−1个向量基底，把圆心用这几个向量做基底表示出...
2019-12-27 18:03:22
401
原创 LOJ #6034. 「雅礼集训 2017 Day2」线段游戏(李超线段树)
题目主席树套李超树。一开始以为空间是O(nlog2n)O(n\log^2n)O(nlog2n)的。然后发现这个神奇的数据结构（李超树）是O(n)O(n)O(n)空间的。。。。注意不要往李超树里面加斜率不存在的线段，好像处理不了的样子，在外面简单维护即可。时间复杂度O(nlog2n)O(n\log ^2n)O(nlog2n)AC Code\mathrm {AC \ Code...
2019-12-26 20:33:47
251
原创 最小圆覆盖
讲得很好的一篇博客具体来讲就是把一个引理运用的淋漓尽致。引理：如果一个点集的最小圆覆盖没有覆盖这个点，那么这个点在，这个点集加上这个点的新点集的最小覆盖圆（边界）上。证明：反证，如果不在，那么显然新点集的最小覆盖圆就应该是原点集的最小覆盖圆。（这个。。。卡圆来调整证明吧，我显然一下。）所以我们可以求出前i−1i-1i−1个点的最小覆盖圆，如果第iii个点不在其中，那么第iii个点一定在前i...
2019-12-25 22:09:34
470
原创 LOJ #2159. 「POI2011 R1」收缩点 Plot（最小圆覆盖，倍增代替二分）
题目首先可以二分答案ansansans。然后再从左往右在最小圆半径<=ans<=ans<=ans的情况下尽量多分点在同一组。因为最小圆覆盖需要随机，所以不能直接从左往右加。这个过程也不能二分，因为checkcheckcheck的复杂度是O（二分答案）O（二分答案）O（二分答案）的。但是发现我们有总答案=n=n=n于是可以用倍增代替二分，那么我们checkcheckch...
2019-12-25 22:08:54
211
原创 LOJ #2158. 「POI2011 R1」移方块 Shift（构造）
题目我太难了。发现把第三个放到开头有点鬼畜。我们来把这个操作规约到我们熟悉的操作。做一次AABAABAAB我们得到了什么？把[2,n][2,n][2,n]每个数向后移两个位置。所以我们可以从i=[2,n]i=[2,n]i=[2,n]依次把iii移到111，然后把1.2.3..i−11.2.3..i-11.2.3..i−1这个之前已经排好序的序列用AABAABAAB移到整个序列的最后，如...
2019-12-25 20:28:36
185
原创 LOJ #2156. 「POI2011 R1」棒棒糖 Lollipop（构造）
题目我们把价格为222看做是长度为222,构成一个新的有很多个分界点为原来棒棒糖中不同口味的分界点的棒棒糖。考虑暴力，我们是从左往右扫，如果相距为kkk的两个点都是分界点，那么这就是一组答案。考虑如果1.当前两个点都不是分界点，那么右移长度1就一定两个都是。2.如果有一个不是，、那么如果是分界点那一个点后面是长度为1的口味，、那么右移长度1也可以合法，、否则：、、右移长度1，重复...
2019-12-25 20:15:08
422
原创 LOJ #2155. 「POI2011 R1」同谋者 Conspiracy（深入性质分析）
题目性质一：两个合法方案中，在一个方案AAA是后勤组且在另一个方案BBB中不是的人最多有一个，同谋者同理。证明：反证，如果有两个，那么在AAA中同为后勤组保证他们间有边，与他们在BBB中同为同谋者矛盾。那么拿到这个性质就可以按照网上的大多数题解直接2−sat2-sat2−sat一波找出特解然后一个人一个人的尝试切换组之类的就可以了。但是我们初步算一下，这个题没给边数的范围，所以我们是可能...
2019-12-25 20:06:46
246
原创 LOJ #2021. 「AHOI / HNOI2017」大佬（DP）
题目感觉这题和劈配差不多恶心人的脑筋。能造成的伤害只和你有几天来打伤害有关，至于是那几天，这不重要。所以我们就可以通过dpdpdp求出你最多能苟多少天来打伤害（注意即使你苟最多天的方案会导致你最后丧失斗争能力你也要接受）。然后BFS（搜）出苟iii天用一次怼大佬能打出的伤害vvv的所有方案(v,i)(v,i)(v,i)注意那个iii要取最小的iii。然后对于每个大佬，我们对于需要找出两...
2019-12-24 15:39:27
342
原创 LOJ #2023. 「AHOI / HNOI2017」抛硬币（组合计数，拓展Lucas）
题目有人讲得比我好嗯，总的来说就是把情况通过反转每一位的值来建立一一映射。建立映射之后就可以简单+不能建立映射的方案再/2了，这个题对于2的逆元的处理十分巧妙，也可以当做一个不错的exLucasexLucasexLucas板题。AC Code\mathrm {AC \ Code}AC Code#include<bits/stdc++.h>#define...
2019-12-24 15:27:36
308
原创 LOJ #2022. 「AHOI / HNOI2017」队长快跑 （单调队列，转化思想，计算几何）
题目题意：平面上有nnn条射线，求从(0,0)(0,0)(0,0)到(Tx,Ty)(Tx,Ty)(Tx,Ty)最短合法路径长度（保证存在合法路径），合法路径是指不经过射线（可以经过不被其他射线覆盖的端点）的，任意时刻走的方向在xxx轴正方向上的投影为正的一条路径。题解：首先，可以发现对于每条射线，我们都是从下或从上绕过去，那么可以把每条射线调整为竖直向上或向下。然后我们想到可以dpdpdp...
2019-12-24 15:16:00
400
原创 多项式的各种算法\cdot极
极\cdot多项式假之前的一篇感觉太不系统了。这次将用链套链的方式写好这些博文造好这些轮子。0.前置知识：0.1 极限0.2 生成函数1.FFT\mathrm {FFT}FFT2.NTT\mathrm {NTT}NTT3.MTT暂且跳过。4.多项式求逆5.多项式对数函数6.多项式指数函数7.多项式开根8.多项式快速幂9.多项式除法10....
2019-12-22 17:04:08
2918
原创 多项式多点求值
给定一个nnn次多项式f(x)f(x)f(x),现在请你对于 i∈[1,m]i\in [1,m]i∈[1,m] ，求出 f(ai)(mod998244353)f(a_i)\pmod {998244353}f(ai)(mod998244353)。好像有一个啥定理：f(ai)=f(x)mod (x−ai)f(a_i) = f(x) \mod (x-a_i)f(ai)=f(x)mod(x−ai...
2019-12-22 17:03:27
340
原创 多项式除法
已知A(x),B(x)A(x),B(x)A(x),B(x)，求A(x)=B(x)C(x)+R(x)A(x) = B(x)C(x)+R(x)A(x)=B(x)C(x)+R(x)且degR<degB\deg R < \deg BdegR<degBdegf(x)\deg f(x)degf(x)为f(x)f(x)f(x)的最高次数。...
2019-12-21 17:18:10
435
原创 生成函数
1.普通生成函数OGF\texttt{OGF}OGF对于无顺序的计数f(x)=∑i=0aixif(x) = \sum_{i=0} a_ix^if(x)=i=0∑aixi这里的xxx是形式幂级数，我们不关心xxx的具体取值，我们只关心系数，可以通过复分析证明可以无视收敛条件化简：∑i=0xi=11−x\sum_{i=0} x^i = \frac 1{1-x}∑i=0xi=1−x1Ex...
2019-12-21 15:15:50
315
原创 极限
limx→∞f(x)=A\lim_{x\rightarrow \infty} f(x) = Alimx→∞f(x)=A即为limx→+∞f(x)=limx→−∞f(x)=A\lim_{x\rightarrow +\infty} f(x) =\lim_{x\rightarrow -\infty} f(x) = Alimx→+∞f(x)=limx→−∞f(x)=A费马定理：设f(x)...
2019-12-21 15:12:59
272
原创 多项式快速幂
已知f(x)f(x)f(x)，求g(x)=f(x)k(modxn)g(x) = f(x)^k\pmod{x^n}g(x)=f(x)k(modxn)版本1:1:1:qk≤1e6qk\leq1e6qk≤1e6，其中qqq是f(x)f(x)f(x)的最高次数。将f(x)f(x)f(x)FFTFFTFFT后直接对每个点值求kkk次方后再IDFTIDFTIDFT回去即可。非常快，不要小看了这种做法。...
2019-12-21 09:50:09
1525
5
原创 多项式开根
已知f(x)f(x)f(x)，求g(x)=f(x)12(modxn)g(x) = f(x)^{\frac 12} \pmod{x^n}g(x)=f(x)21(modxn)保证常数项为111算法111：g(x)=exp(lnf(x)2)g(x) = \exp(\frac {\ln f(x)}2)g(x)=exp(2lnf(x))O(nlogn)O(n\log n)O(nlogn)算法...
2019-12-20 22:41:49
405
原创 多项式指数函数
已知f(x)f(x)f(x)，求g(x)=exp(f(x))(modxn)g(x) = \exp(f(x)) \pmod{x^n}g(x)=exp(f(x))(modxn)lng(x)−f(x)=0\ln g(x) - f(x) = 0lng(x)−f(x)=0直接牛顿迭代g1(x)=g0(x)−lng0(x)−f(x)1g0(x)g_1(x) = g_0(x) - \frac{\ln...
2019-12-20 22:22:28
547
原创 多项式对数函数
已知f(x)f(x)f(x),求g(x)=lnf(x)(modxn)g(x) = \ln f(x) \pmod {x^n}g(x)=lnf(x)(modxn)保证常数项为111g(x)′=f(x)′f(x)g(x)' = \frac {f(x)'}{f(x)}g(x)′=f(x)f(x)′g(x)=∫f(x)′f(x)g(x)=\int \frac{f(x)'}{f(x)}g(x)=∫f...
2019-12-20 21:58:07
292
原创 多项式求逆
求B(x)B(x)B(x)使A(x)B(x)=1(modxn)A(x)B(x) = 1 \pmod {x^n}A(x)B(x)=1(modxn)即A(x)−1B(x)=0(modxn)A(x)-\frac 1{B(x)}=0 \pmod{x^n}A(x)−B(x)1=0(modxn)直接牛顿迭代即可。B1(x)=B0(x)−A(x)−1B0(x)(A(x)−1B(x))′∣B(x)=B0(...
2019-12-20 21:37:32
298
原创 快速数论变换（NTT）
其实对于看了FFT的人来说，也没有什么好讲的。对于没有看过FFT的人来说，讲的内容和FFT也没什么两样。所以就直接上板子吧。#include<bits/stdc++.h>#define maxn 3000005#define mod 998244353#define rep(i,j,k) for(int i=(j);i<=(k);i++)#define per(i,...
2019-12-20 21:07:17
508
转载 【转载】STL中的Complex
原文链接:https://blog.csdn.net/jingza/article/details/77334417这expexpexp和sqrtsqrtsqrt感觉好强的样子。
2019-12-20 20:05:20
181
原创 老魔杖（博弈，不能用SG）
你有aaa个1，bbb个222，ccc个333,ddd个444操作1：把一个数字分成两个正数满足总和不变。操作2：把nnn个nnn删除。不能操作者输。求一个(a,b,c,d)(a,b,c,d)(a,b,c,d)是否先手必胜。a,b,c,d<=1e10000a,b,c,d<=1e10000a,b,c,d<=1e10000这个题不能用SGSGSG函数，因为你把大的数字分了，...
2019-12-17 21:29:35
224
原创 和箭头一起游走（主席树+倍增）
n,m,q<=1e5,T<=1e15n,m,q<=1e5,T<=1e15n,m,q<=1e5,T<=1e15题目保证箭头不会相交。一道不卡精度的计算几何发现一个箭头上的人只会走到另一个确定的箭头上或走出去，那么我们可以把这个走的关系建成图。求出每一步的距离后即可分三类简单讨论得到终点坐标。对于T<=1e15T<=1e15T<=1e...
2019-12-17 20:49:08
127
原创 ARC 062 F AT2143 AtCoDeerくんとグラフ色塗り / Painting Graphs with AtCoDeer（burnside引理+点双）
题目分单边，单环，复环来讨论。发现复环中可以交换任意两条边的颜色，直接组合数即可。单环用burnside引理。单边直接乘。AC Code\rm AC \ CodeAC Code#define maxn 55#define maxm 205#define mod 1000000007using namespace std;int n,m,k;int i...
2019-12-15 23:00:32
324
原创 【UOJ#422】【集训队作业2018】小Z的礼物（min-max反演，简单插头DP）
题目一看数据范围我以为要插头DP。然后发现不维护连通性的插头DP不是随便打吗。。。。。。直接min−maxmin-maxmin−max反演。求E(∑T⊂S(−1)∣T∣+1min(T))=∑T⊂S(−1)∣T∣+1E(min(T))E(\sum_{T\sub S} (-1)^{|T|+1}\min(T)) = \sum_{T\sub S} (-1)^{|T|+1}E(\min(T))E...
2019-12-15 22:21:14
191
原创 #6544. 联通块（套路--线段树）
题意： 给你一颗nnn个节点的树，问有多少个区间 [L,R](L<=R)[L,R] (L<=R)[L,R](L<=R) 满足 区间内的点在树上恰好是一个联通块。n<=3e5n<=3e5n<=3e5题解：好套路啊。树上一个点集为联通块−>->−>树上一个点集的导出子图中点数-边数=1=1=1−>->−>按序加入动态维护区...
2019-12-15 11:50:54
414
原创 #6543. 大葱的神谕（杜教筛）
It’s Obvious that the Chinese Green Onion is the STRONGest MAN in the universe.It’s Such A Universal truth that Even woman and children all know it !!求∑i=1n(<=1e9)φ(n)(−1)n+1\sum_{i=1}^{n(<=...
2019-12-15 11:15:09
132
原创 杨氏矩阵（杨表）
杨氏矩阵更像一个维护最长上升子序列的数据结构，实际上最长上升子序列的O(nlogn)O(n\log n )O(nlogn)的单调栈解法就是维护杨表的第一行。。。所以杨表就是最长上升子序列的究极工具。推荐阅读IOI2019IOI2019IOI2019国家集训队论文 袁方舟的那篇。例题：CTSC2017最长上升子序列CTSC2017最长上升子序列CTSC2017最长上升子序列我们可以只维护...
2019-12-14 21:36:09
1400
原创 CF280C Game on Tree（期望的线性性）
题目给出一棵树，每次随机等概率选择一未染黑的点，将它及其子树染黑。问期望多少次操作可以将树全部染黑。n<=1e5n<=1e5n<=1e5期望操作次数 -> 每个点的期望操作次数之和 -> 每个点被操作的概率之和 -> ∑i1depi\sum_{i} \frac 1{dep_i}∑idepi1我真的没有在水博客。AC Code\rm AC...
2019-12-13 17:30:02
105
原创 20191212 C 简单点分治
给定一棵有n个节点的无根树，树上的每个点有一个非负整数点权。定义一条路径的价值为路径上的点权和-路径上的点权最大值。 给定参数PPP，我们想知道，有多少不同的树上简单路径，满足它的价值恰好是PPP的倍数。 注意：单点算作一条路径；u!=vu!=vu!=v时，(u,v)(u,v)(u,v)和(v,u)(v,u)(v,u)只算一次。n<=1e5,P<=1e7n<=1e5 , P&l...
2019-12-12 21:43:47
181
原创 20191212 A 状压DP
有一个数组aaa，其中有nnn个元素，现在请构造一个数组bbb，其中也有nnn个元素，使得数组bbb中所有元素的异或值为000，且aia_iai到bib_ibi的距离之和最小。T<=6,n<=15,ai<=1e9T<=6 , n <=15 ,a_i<=1e9T<=6,n<=15,ai<=1e9这个。从高位到低位DPDPDP。当所有...
2019-12-12 21:41:10
103
原创 BZOJ 3681 Arietta(可持久化线段树合并优化建边|网络流)
题目尽管是一道板题但是还是很难的。中途输出了一下结果，4e54e54e5的点1e61e61e6的边。。。SAPSAPSAP在一条链的时候TTT了，换成DinicDinicDinic就过了。。。应该是DinicDinicDinic比SAPSAPSAP更能创造奇迹。AC Code:\rm AC \ Code:AC Code:#include<bits/stdc++...
2019-12-10 20:57:23
172
原创 [SHOI2007]善意的投票 （最小割|niubility的贪心）
题目…最小割建图显然.AC Code\rm AC\ CodeAC Code#include<bits/stdc++.h>#define maxn 305#define maxm maxn * maxn * 2#define inf 0x3f3f3f3f#define Copy(a,b) memcpy(a,b,sizeof a)using names...
2019-12-10 15:44:21
88
原创 BZOJ 1305: [CQOI2009]dance跳舞 (二分答案网络流|错误的贪心)
题目循规蹈矩的做法是比较显然的。有着看似严谨的贪心最后的评论一语道出贪心的错误。有些的度数看似是有很多，但是因为对方太忙，都是废的度数。。。所有人度数一致才可以套结论。。。。。。AC Code\rm AC \ CodeAC Code#include<bits/stdc++.h>#define maxn 250#define maxm 50005#...
2019-12-10 15:22:53
79
原创 Comet OJ 序列（分解质因数卡常+平衡复杂度（分块））
题面首先考虑把分解质因数的复杂度从O(x+σ0(x))O(\sqrt x+\sigma_0(x))O(x+σ0(x))降低。可以写PollardRhoPollard_RhoPollardRho但是没必要。把x\sqrt xx以下的质因数预处理即可O(xlnx+σ0(x))O(\frac {\sqrt x}{\ln x} + \sigma_0(x))O(lnxx+σ0(x))...
2019-12-05 17:35:30
264
原创 Stoer_Wagner算法
浅谈无向图最小割问题的一些算法及应用 王文涛定理证明详见论文。HDU3691Nubulsa Expo代码：#include<bits/stdc++.h>#define maxn 305#define LL long longusing namespace std;int n,m,S,c[maxn][maxn],d[maxn],vis[maxn],id[maxn],an...
2019-12-04 22:49:00
849
原创 LOJ #2146. 「SHOI2017」寿司餐厅（最大权闭合子图）
题目注意是CCC种。大概就这样一个图。AC Code\rm AC\ CodeAC Code#include<bits/stdc++.h>#define maxn 10005#define maxm 40005#define inf 0x3f3f3f3f#define Clear(a,b) memset(a,b,sizeof a)#define Co...
2019-12-03 21:16:23
186
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人