- 博客(84)
- 资源 (2)
- 收藏
- 关注
RSS订阅
原创 ACM中的数学问题合集
数学在ACM竞赛中通常会占很大的比重,例如这次暑假的多校比赛中,几乎所有题目都和数学有着很大的联系所以统计一下一些常用的知识点数论欧拉定理威尔逊定理中国剩余定理费马小定理...
2019-08-20 22:30:35
11843
3
原创 Factories（树上背包转移 树形dp）
original link - https://nanti.jisuanke.com/t/41291题意：给出一棵树，你需要选择KKK个叶子，使得两两距离和最小。解析：考虑从儿子到父亲的转移，用dp[i][j]dp[i][j]dp[i][j]记录iii的子树中，选取jjj个叶子的预期最小距离。转移的时候用背包进行转移，将dp[son][h]dp[son][h]dp[son][h]当作一个...
2019-08-31 21:58:24
300
原创 Relief grain（树链剖分 权值线段树 离线）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=5029题意：给出一棵树，kkk个操作：u,vu,vu,v路径上的点的多重集合内放入一个ppp。求所有操作后，每个点的集合内最多的元素是哪个。解析：这题是所有操作结束后再查询，所以应该不需要实时更新，考虑离线。先树链剖分，将所有的链转化为区间。考虑区间更新[L,R](add...
2019-08-31 20:33:05
258
原创 Vertex Covers（折半 高维前缀和 顶点覆盖）
original link - https://nanti.jisuanke.com/t/A1771题意：给出n≤36n\leq36n≤36个点，mmm条边，点权，一个集合的权值为点权和。求所有完全覆盖所有边的集合权值和。解析：因为nnn卡在363636这个特殊的位置，所以考虑折半做法。枚举两边的点的所有状态，对于左边的状态XXX和右边的YYY，首先要求XXX覆盖左边内部的边，YYY覆盖...
2019-08-31 20:21:14
416
4
原创 4836: [Lydsy1704月赛]二元运算（cdq分治 FFT）
original link - https://www.lydsy.com/JudgeOnline/problem.php?id=4836题意：给出一个nnn的序列aaa，mmm的序列bbb，每次查询给出ccc，求有多少对(i,j)(i,j)(i,j)使得ai  op  bj=ca_i\;op\;b_j=caiopbj=c。a&am...
2019-08-30 20:59:38
263
原创 4555: [Tjoi2016&Heoi2016]求和（NTT 第二类斯特林数）
original link - https://www.lydsy.com/JudgeOnline/problem.php?id=4555题意：解析：第二类斯特林数先拆开，得到：f(n)=∑i=0n∑j=0i1j!∑k=0j(−1)kCjk(j−k)i∗2j∗(j!)f(n)=∑i=0n∑j=0i∗2j∗∑k=0j(−1)kCjk(j−k)if(n)=∑i=0n∑j=0i∗2j∗j!∗∑...
2019-08-30 17:27:56
208
原创 组合数学总结
by 段ls1. 排列与组合1.1 排列（1）在没有其他条件的情况下，从 nnn 个不同元素中选取 rrr 个不同的元素的排列数为Anr=n!(n−r)!A_{n}^{r} = \frac{n!}{\left( n - r \right)!}Anr=(n−r)!n!，当 r＞nr＞nr＞n 时，AnrA_{n}^{r}Anr=0（2）在 nnn 个不同元素中选取 rrr 个元素的...
2019-08-30 09:41:25
1100
原创 NTT模板
#include<cstdio>#define swap(x,y) x ^= y, y ^= x, x ^= y#define LL long long const int MAXN = 3e6 + 10, P = 998244353, G = 3, Gi = 332748118; char buf[1<<21], *p1 = buf, *p2 = buf;inl...
2019-08-29 14:23:31
200
原创 K-Inversions（FFT）
original link - https://vjudge.net/problem/Gym-101002E题意：给出一个长度为nnn的ABABAB序列，对于每个k∈[1,n−1]k\in[1,n-1]k∈[1,n−1]，求出有多少个BABABA对，满足posA−posB=kpos_A-pos_B=kposA−posB=k。解析：FFTFFTFFT的裸题吧，两种数列，计数目标为距离差固...
2019-08-29 13:41:26
394
原创 Programming Team（01分数规划 树形依赖背包）
original link - https://vjudge.net/problem/Gym-101002E题意：给出一棵树，到根节点的路径上的点都被选择时才能选择一个点。根节点必须选择。每个点有p,sp,sp,s，要求选择kkk个点使得∑p∑s\dfrac{\sum p}{\sum s}∑s∑p最大。解析：先按照01分数规划，二分v=∑p∑sv=\dfrac{\sum p}{\sum ...
2019-08-29 13:32:12
229
原创 huntian oy（杜教筛 欧拉函数）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6706题意：求f(n,a,b)=∑i=1n∑j=1igcd(ia−ja,ib−jb)[gcd(i,j)=1]%(109+7)f(n,a,b)=\sum_{i=1}^n \sum_{j=1}^i gcd(i^a-j^a,i^b-j^b)[gcd(i,j)=1]\%(10^9+7)...
2019-08-25 16:50:30
556
3
原创 2200+专项：D. Rectangle Painting 1（棋盘dp）
original link - http://codeforces.com/contest/1198/problem/D题意：n∗nn*nn∗n矩阵有些点为#\##，你可以选择一块区域，染成∗*∗，花费为区域的长和宽的最大值。求染完所有#\##点的最少花费。解析：显然，我只会染正方形。原来的花费至多只有nnn。如果有更少的，那么一定会有一行、一列没有被染色。所以我们可以枚举每行或者每列（取...
2019-08-25 11:03:51
272
原创 2019中国大学生程序设计竞赛（CCPC） - 网络选拔赛（主席树，后缀数组）
第7题想出来了但是没过，挺烦的吧~Path题意：任意路径中的第KKK短路径的长度。首先是将所有边都放到setsetset，然后每次选择最短的往外扩。解析：由于我们只要求第KKK条就行，所以队列内的最大数量应该是K−K-K−已经遍历的数量，实时维护即可。（当当前路径劣于setsetset内最长路径时，且已经到达上限，那么就无视这条。否则加入，若加入后超出上限必须立刻删除尾部原始并更新最...
2019-08-24 20:44:26
391
原创 ACM知识清单
桶排序排序数列分布较密，值域较小，例如1e61e61e6个数属于∈[1,10]\in[1,10]∈[1,10]。斜率优化dpdpdp的增量具有单调性时，可以确定转移的最优态，实现O(n)O(n)O(n)转移。矩阵快速幂O(logn)O(logn)O(logn)求解线性递推式第nnn项值。容斥欲正则反，虚虚实实。位操作将一个数看成二进制位，考虑每个位的贡献。逆元可由费马小定理或...
2019-08-22 22:09:52
535
原创 B. Russo's Russian（字符串处理）
original link - http://codeforces.com/gym/102299/problem/B题意：M = H '|' P 或者 '|' M 或者 PH = M 或者 '$'P = P ':' T 或者 TT = '{' M '}' 或者 数字串可以在任意连接处塞空格，问给出的串是否可以视为MMM串。解析：各种情况判一判就行了。注意MMM和PPP的分隔符要在...
2019-08-22 19:13:28
226
原创 严格次小生成树（LCA法）
original link - https://www.luogu.org/problem/P4180题意：给出一个图，求严格次小生成树，即边权和严格小于MSTMSTMST。解析：定理： 若存在次小生成树，则必然存在一个次小生成树与MSTMSTMST只有一条边的差异。所以我们先做一遍MSTMSTMST，然后去判断每条边加进去后的情况。显然加进去后形成一个环，要在环上原来的边中删除一条边（...
2019-08-22 11:24:11
462
原创 Walk Through Squares（AC自动机 dp）
original link - https://cn.vjudge.net/contest/320257#problem/K题意：求包含串A,BA,BA,B的串的个数，要求满足num(′R′)=n,num(′D′)=mnum('R')=n,num('D')=mnum(′R′)=n,num(′D′)=m解析：先把两个串做...
2019-08-21 21:17:32
272
原创 一次同余方程（ax%p=b）
最小非负解求 ax≡b(mod  p)ax\equiv b(mod\; p)ax≡b(modp)的最小非负解。公式化为：ax+py=bax+py=bax+py=b此时用扩展欧几里得算法求解即可。/* * Author : Jk_Chen * Date : 2019-08-21-10.26.44 */#include<bits/stdc++....
2019-08-21 11:01:49
967
原创 多项式插值法
Lagrange多项式插值给出样例插值点，求出XXX的函数值F(X)F(X)F(X)。其中X∈RX\in RX∈R时间复杂度为O(k2)O(k^2)O(k2)，kkk为插值点个数代码：#include <bits/stdc++.h>#define rep(i,a,b) for(ll i=a;i<=b;i++)using namespace std;typedef l...
2019-08-20 22:28:09
2258
原创 博客格式化
介绍 - 全英文段：加粗，加字体Courier - 代码块：无变化 - 保留空行 - 全*号行为分界线，无变化 - 出现汉字段加字体楷体 - 纯加粗段，加size=4 - 字体标签加在标题符号#后面代码：#include <bits/stdc++.h>#define rep(i,a,b) for(int i=a;i<=b;i++)using namespac...
2019-08-20 20:55:29
183
原创 BSGS，扩展BSGS（求离散对数 b^x%p=n）
应用大步小步法（Baby−Step−Giant−StepBaby-Step-Giant-StepBaby−Step−Giant−Step，简称BSGSBSGSBSGS），可以在O(logP)O(logP)O(logP)的时间内求出basex≡n(mod  P)base^x\equiv n(mod\;P)basex≡n(modP)的解xxx。过程其实就是一个暴力而已...
2019-08-20 20:37:26
607
原创 Yaoge’s maximum profit（树链剖分 维护两个方向的最值）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=5052题意：给出一棵树，点权。每次查询要从xxx走到yyy，可以在路径之前的地方取一个值v1v_1v1，之后的地方v2v_2v2，你可以得到v2−v1v_2-v_1v2−v1，当然可以不取为0。求每次得到的最大值。走完后，路径上的点权值加上vvv。解析：考虑用树链剖...
2019-08-20 15:59:12
249
原创 Z Algorithm（后缀匹配）
介绍一个和马拉车算法很像的东西，可以O(n)O(n)O(n)求出一个串的所有后缀与前缀的匹配长度。因为往后扩不会返回，右端点一直加，所以时间复杂度为O(N)O(N)O(N)两个串的情况如果要匹配两个串，串xxx的所有后缀匹配串yyy的前缀，只需要使用特殊字符将xxx连到yyy后面即可。orginal link - http://acm.hdu.edu.cn/showproblem....
2019-08-20 13:16:39
369
原创 Rikka with Travels（树形dp）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6686题意：给出一棵树，两条路径不相交定义为不能有重复点，一个有向关系<p1,p2><p1,p2><p1,p2>，表示两天不相交的路径的点数为p1,p2p1,p2p1,p2，求有多少对。解析：要维护以下变量...
2019-08-20 09:34:16
212
原创 codeblocks快捷键——高效编程从编译器开始
基础：Ctrl+Z,X,C,V,S,A，Ctrl+Shfit+ZF9：Build+RunF8：Debug不广为人知的比较牛的快捷键：Ctrl+O：打开项目
2019-08-19 19:25:48
1071
原创 Rikka with Coin（想法）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6685题意：有4种硬币：10,20,50,10010,20,50,10010,20,50,100，每种无限个，现在你要挑选最少的个数，使得可以组成给出的n个价格。解析：开始想到超过100的应该都是直接100比较少。所以模了10到90，发现只需要10,20,20,5010,2...
2019-08-19 18:50:35
631
6
原创 java BigInteger 大数版矩阵快速幂模板
import java.math.*;import java.util.*; public class Main { public static BigInteger __ng=BigInteger.valueOf(-1); public static BigInteger __0=BigInteger.valueOf(0); public static BigInt...
2019-08-19 11:55:16
382
原创 Mindis 2019百度之星初赛第一轮（建图 最短路）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6670题意：给出n个矩阵，起点到终点，可以往四个方向移动，被kkk个矩形覆盖时，速度为k+1k+1k+1，边也算覆盖。求最短时间。解析：将所有坐标离散化，建网格图（400*400）。两个点之间的距离已经可以计算了，关键是被多少个矩形覆盖。首先我们二维差分一下，处理出每个...
2019-08-17 22:40:28
328
原创 2019牛客暑期多校训练营（第十场）——Wood Processing（斜率优化dp）
original link - https://ac.nowcoder.com/acm/contest/890/J题意：给出n个矩形的高和宽，你可以自己排列这些矩形。现在要将其合并为j堆，每一堆的所有矩形都要砍到这堆的最小值为止，求最少砍掉的面积。解析：先从小到大排序，显然最优解一定要先排序。设dp[i][j]dp[i][j]dp[i][j]为前iii堆合并成jjj堆的最少砍掉面积，h...
2019-08-17 22:01:39
281
原创 数据生成器（带根树，无根树，序列，有向图）
#include<bits/stdc++.h>using namespace std;#define LL long long#define rep(i,a,b) for(int i=a;i<=b;i++)const int maxn=1e5+5;namespace RootTree{ /* output : the ...
2019-08-16 22:24:30
817
原创 Tree（可持久化字典树 树上区间 LCA）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=4757题意：给出一棵点权树，每个查询要求出两个点之间的路径上的一个点，与给出的值异或最大。解析：区间单点异或===字典树维护一个可持久化字典树，那么(x,y)(x,y)(x,y)分解为(x,lca)+(lca,y)(x,lca)+(lca,y)(x,lca)+(lca,y...
2019-08-16 22:01:37
265
原创 Query on A Tree（可持久化字典树）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6191题意：给出一棵点权树，每次查询求出一棵子树RtRtRt中的一个点，与valvalval异或后数最大，求这个最大的数。解析：首先，子树是一个典型的转化，按照dfsdfsdfs序排到数组中，那么一棵子树所在位置为一段连续区间。问题变为区间内，单点值问题。而单点与给出...
2019-08-16 21:31:44
265
原创 2019牛客暑期多校训练营（第九场）——Cutting Bamboos（主席树区间更新 标记持久化）
original link - https://ac.nowcoder.com/acm/contest/889/H题意：有很多个竹子，给出高度，每次查询一个区间，求一个高度HHH使得：区间内的竹子，低于HHH的部分长度比上高于的部分为y−x:xy-x:xy−x:x。解析：考虑从左到右每次加入一个竹子，高度hhh，我们更新[1,h][1,h][1,h]区间加1。那么对于查询[L,R][L,R...
2019-08-16 19:08:42
219
原创 2019牛客暑期多校训练营（第九场）——The power of Fibonacci（中国剩余定理 循环节拆分）
orginal link - https://ac.nowcoder.com/acm/contest/889/A题意：给出两个数n,mn,mn,m，求∑i=1nFim\sum_{i=1}^n F_i^m∑i=1nFim，模数为1e91e91e9解析：我们将1e91e91e9拆分为29∗592^9*5^929∗59，分别求出斐波那契数列在模数为292^929和595^959下的循环节：7...
2019-08-15 22:35:12
254
5
原创 2019牛客暑期多校训练营（第九场）——All men are brothers（逻辑推理）
original link - https://ac.nowcoder.com/acm/contest/889/E题意：给出n个人，开始每个人独立。然后每次操作选择两个人，连接两个人所在的集合。问每次操作后，我选择4个人不在一个集合的方案数。解析：最开始的时候，A4A4A4表示选择4个人的方案数，为：Cn4C_n^4Cn4，A3,A2,A1A3,A2,A1A3,A2,A1同理。然后维护...
2019-08-15 19:22:08
517
原创 2019牛客暑期多校训练营（第九场）——Quadratic equation（二次同余方程求解）
original link - https://ac.nowcoder.com/acm/contest/889/B题意：给出b,cb,cb,c，求出满足x+y≡b,x∗y≡c,mod  1e9+7x+y\equiv b,x*y\equiv c,mod\;1e9+7x+y≡b,x∗y≡c,mod1e9+7的一对x,yx,yx,y。解析：(x+y)2≡x2+y2+2...
2019-08-15 19:03:46
245
原创 阶、原根、指标（离散对数）
阶定义： 最小的ttt使得：at≡1(mod  p)a^t\equiv1(mod \;p)at≡1(modp)，则称ttt为aaa在模ppp下的阶（(a,p)=1(a,p)=1(a,p)=1），记作：δp(a)=t\delta_p(a)=tδp(a)=t理解： a%pa\%pa%p下的阶ttt相当于a%pa\%pa%p的循环节，每乘上ata^tat，值不变。拓展...
2019-08-14 22:30:11
1414
原创 2019 Multi-University Training Contest 8——Acesrc and Good Numbers（数学 想法）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6659题意：定义F(d,n)F(d,n)F(d,n)为1,2...n1,2...n1,2...n的十进制下位中，数字d出现的个数，找出小于等于nnn的最大的数xxx使得F(d,x)=xF(d,x)=xF(d,x)=x。解析：神奇的做法，可惜比赛的时候在敲其他题目没有好好想过...
2019-08-14 18:42:52
610
原创 2019 Multi-University Training Contest 8——Acesrc and Travel（树形dp）
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6662题意：一棵树，点权为a−ba-ba−b，你先选择一个点作为起点，然后从对面开始依次走，每个点只能走一次。你要值最大，对面要值最小。 最后可以得到的最大值。解析：显然固定好位置后，对面先走。有两种走法，向下和向上。向下很好处理，使用down[i][0]down[i][...
2019-08-14 17:46:58
377
原创 python的奇淫巧计
python做题的一些细节： 1. 当负数取模时，模数为正数时与其他语言的取模不同，需要将模数改为对应负数再取模2. 二维数组可以append多个list得到，但是效率极低3. 输入一行多个数字时，先使用split分成多个字符串，再一个一个赋值示例：http://codeforces.com/gym/101933/problem/K# 会TLE~def Mod(a): if a...
2019-08-14 11:16:52
500
mnist手写数字数字集
2019-12-11
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人