celerity

摘要： 被卡常数的代码： #include<bits/stdc++.h> #define int long long #define N 300010 #define mo 1000000007 #pragma GCC optimize(3) int p[N],vi[N],ct,b[N],d[N],e[N] 阅读全文

该文被密码保护。 阅读全文

摘要： 先考虑怎么求lcm，考虑从序列后部插入一个数$y$的影响。 \(\frac{xy}{(x,y)}=x\frac{y}{(x,y)}\) $x$是容易计算的，边算边取模即可。 $\frac{(x,y)}$这个东西很难计算，因为$x$很大。 但是$(x,y)=(x\mod y,y)$，所以可以从左到右计 阅读全文

摘要： 根据题意，我们考虑从小到大枚举树上每个节点，并且让它的值最小。 一组已经确定的值的序列，序列上的每个点都会要求它的子树的值大于等于它。 已经确定的值已经满足条件，所以只需要考虑未被确定的值的限制，也就是要求若干子树（若干dfs序连续段）的值$\leq 某个数$。 考虑使用二分图完美匹配。 把每个数从 阅读全文

摘要： 给出一个时间复杂度为$O(d(k)\log_2n\log_2c(k))$（$c$是质因数个数函数）的算法。 原式就是$vk=\prod a_i$，$v,k$互质。 $k$已经给出，设唯一分解为$p_1^....p_n^$ 则$v$的$p_1,p_2...p_n$的指数必须是$0$，其他任意。 就是给 阅读全文

摘要： 1_998244353： 由于$998244353$是十分常见的模数，所以可以猜想到取模。 $1,19,361....$是$19^x$。 前两个点可以直接用快速幂计算。 第三个点由于数比较大，所以使用euler定理，把幂次模$998244352$后快速幂。 namespace s1{ int mo= 阅读全文

摘要： 对一个多项式$A$的快速幂可以在$O(n\log_2n)$的时间内计算多项式模$xn$的乘方。 如果多项式的项数$m$较少，则我们有一个$O(nm)$的算法。 \((A^k(x))'=kA^{k-1}(x)A'(x)\) \(A^k(x)'A(x)=kA^k(x)A'(x)\) \([x^n]A^k 阅读全文

摘要： 如果当$R-L+1$很小时怎么做？ 完全平方数的质因数分解的每个质数的质数肯定是$0$。 考虑每个数的唯一分解，把它的每个质因数模$2$后视为一个向量。 则判定是否合法，等于判定$[L,R]$中是否存在一个子集$S$使得$S$中的所有质因数分解向量在xor下$=0$。 答案就是$2^{自由基数量}$ 阅读全文

摘要： 研究两个图$G1$和$G2$的乘积的连通块的性质（agc011c） 根据定义，两个点$(a,b),(c,d)$连通的条件：存在$a,b$和$c,d$之间长度为$x$的路径 考虑$G1$的每个连通块和$G2$每个连通块合并后结果： 有任一是孤立点：显然不会连出任何边。 两个非二分图：会发现在$x$足够 阅读全文

摘要： 会发现，无论操作的顺序怎么样，每个节点形成的区域都是个矩形。 发现两个点连通要行/列同时连通。 假设行向左走$x$，右走$y$次，则$x+y\geq X$坐标差绝对值。 我们只需要关心$x+y$，所以把边权设为这个值。 $Y$坐标类似。 暴力连边边数达到$(nm)^2$，不可接受。 但是发现如果存在 阅读全文