自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

KsCla

なかない君と嘆きの世界 輪廻の刻は暮れゆけど

  • 博客(160)
  • 收藏
  • 关注

原创 手动友链

没有摘要

2019-12-27 17:50:31 499

原创 咕咕咕

没有摘要

2019-07-09 23:24:08 821 3

原创 BZOJ4538:[Hnoi2016]网络 (整体二分+Lca+树状数组/线段树+路径交/树链剖分+Heap)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4538题目分析:这题网上好多人写树剖啊,都是把一条路径的区间搞出来之后取反更新,删除的话就套个Heap或者写线段树CDQ分治blablabla,时间复杂度O(nlog3(n))O(nlog^3(n)),好像因为树剖和Heap的常数特别小所以根本不虚。网上某大神用这种方法8200ms就过了,而

2019-12-28 01:25:53 627

原创 CSP2019Day2T3 洛谷P5666:树的重心 (树上倍增)

题目传送门:https://www.luogu.org/problem/P5666很早以前就觉得,凡是考树的重心相关的题,到最后都变成一道模拟题。树的重心有许多优秀的性质,比如:结论一:记f(node)表示以node为根的最大子树的大小。从无根树上的任意一个点x出发,向相邻的点走一步。假设某个点y的f值比x的小,那么x向y走就相当于向树的重心移动了一步。假设不存在这样的y,那么x就是重心。...

2019-11-24 00:59:31 823

原创 BZOJ5495:[2019省队联测]异或粽子 (Trie+贪心)

题目传送门BZOJ(没题面)洛谷(有题面)题目大意给出一个长为n的序列,选择k个不完全重合的区间使得每个区间的异或值的总和最大。题解作为高三退役狗康复训练的第一题,我发现我的做法有点特立独行……4月份Ghastlcon跟我说这题的时候,我在想为什么要用可持久化Trie。而现在我YY出了一个既不用可持久化也不用堆的方法。首先把所有的异或前缀和扔进Trie,接下来目标就是找最大的k...

2019-06-17 18:03:24 347

原创 UOJ#348:【WC2018】州区划分 (FMT优化DP)

题目传送门:http://uoj.ac/problem/348题目分析:题面就是要求将n个点划分为若干个集合,使得刚好包含某个集合的点以及它们之间的边的子图不存在欧拉回路。然后题面给出了一种方法计算某种方案的贡献。由于n很小,不妨先用h[s]表示s中的点能不能刚好划分为一个集合。算出h[]的时间是O(n22n)O(n22n)O(n^22^n)的。对于题面给出的式子,我们发现每个集合的贡...

2018-04-17 22:06:42 808 1

原创 CodeChef Counting D-sets (容斥原理+组合数学)

vjudge题面传送门:https://cn.vjudge.net/problem/CodeChef-CNTDSETS(PS:vjudge上中文版的题面有误,一个点集的直径应该定义为其中点对的切比雪夫距离的最大值。切比雪夫距离是两个点各个维度之差的绝对值取max。这一点看回英文版题面就能知道)题目分析:一道思维难度较大,代码量极少的题。直径=d的点集数 = 直径<=d的点集数...

2018-03-29 21:40:41 618

原创 BZOJ4767:两双手 (组合数学+DP+容斥原理)

题目传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=4767题目分析:一开始看题目名还以为是两只手,后来感觉有些不对劲……题面保证了给出的两个向量叉积为0,就是说它们不平行。不平行的两个向量可以作为一组基底,这样原先平面上的所有点就获得了一个新坐标。于是问题变成了:从(0,0)走到(n,m),中间不能经过指定的k个点,求方案数。...

2018-03-29 16:28:10 849

原创 circle (容斥原理+数据结构)

题目大意:在一条直线上有2*n个点,点与点之间两两配对成n组。现在要你选出三组点对,使得这三组点对满足112233,122331,123123的其中一种形式,问方案数。n≤105n≤105n\leq 10^5。题目分析:多年前的老坑,昨天晚上想填一下,发现还是不会做,而且我还是看不懂题解。懵逼了一整晚,最后翻出标程来看,终于看懂了做法。从n个点对中选取3个点对,有C3nCn3C_n^3...

2018-03-29 10:56:47 1255

原创 BZOJ4710:[Jsoi2011]分特产 (容斥原理+组合数学+DP)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4710题目分析:一开始看完全没有头绪,后来发现这不就是个和第二类stirling数很像的容斥吗?首先考虑没有“每个人至少要拿一个特产”这个条件怎么做。由于不同的特产之间是独立的,可以记h[i][j]表示前i个人拿了j件特产的方案数。转移方程为h[i][j]=∑jk=0h[i−1...

2018-03-28 10:54:59 407

原创 BZOJ1042:[HAOI2008]硬币购物 (容斥原理+DP)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1042题目分析:我看某个课件看到这题,一开始还以为每组询问都重新给出四个面值,导致我一直没有思路QAQ。由于四个面值是固定的,可以先做一次完全背包,将价值为1~maxs的答案记下来。每次询问的时候,记f(s)表示 只有 s集合中的硬币超过限制的方案数,记g(s)表示 至少有 s...

2018-03-28 09:23:41 511

原创 BZOJ4008:[HNOI2015]亚瑟王 (概率DP)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4008题目分析:一道很厉害的DP,做法和背包类似。记g[i][j]表示前i张卡牌,有j张卡牌发动了技能的概率,那么存在如下转移:g[i+1][j]=g[i][j]∗(1−p[i+1])r−jg[i+1][j]=g[i][j]∗(1−p[i+1])r−jg[i+1][j]=g[i...

2018-03-27 15:58:42 394

原创 BZOJ1444:[Jsoi2009]有趣的游戏 (AC自动机+概率DP+高斯消元)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1444题目分析:首先考虑静态的问题:如果已经生成一个字符串,如何让它跟所有模式串匹配?答案是建出所有模式串的AC自动机,然后让生成串在上面跑,如果跑到某个有endpos的节点就一直停在那里。然后考虑动态的问题:如果生成串无限长,如何求出它停在每个节点的概率?把AC自动机扩展成T...

2018-03-27 11:11:54 452

原创 BZOJ3143:[Hnoi2013]游走 (高斯消元+概率DP)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3143题目分析:很久之前就对这种高斯消元解DP值的题目有一种莫名的恐惧,因为它明明是DP却没有递推顺序因为我对概率论一窍不通。学了高消之后,我YY了一下这题的DP方程,发现一直过不了样例。最后居然被tututu在旁边随便口胡一句就过了?!本题要让我们自定每条边的权值,所以我们想...

2018-03-26 21:26:43 657

原创 BZOJ1923:[Sdoi2010]外星千足虫 (高斯消元+二进制压位)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1923题目分析:一开始还很好奇什么是异或高斯消元,后来发现就是用高消解异或方程组。普通的n3n3n^3高消理论上来说是过不了这题的,虽然实际上能过。观察一下,发现高斯消元的主要时间在于选定某一行作为主元后,将后面几行的对应位置变成0。而本题由于是异或操作,可以用bitset或者...

2018-03-26 16:52:50 374

原创 BZOJ1013:[JSOI2008]球形空间产生器sphere (高斯消元)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1013题目分析:设题面给出的第i个坐标为(ai,1,ai,2……ai,n)(ai,1,ai,2……ai,n)(a_{i,1},a_{i,2}……a_{i,n}),设答案为(b1,b2……bn)(b1,b2……bn)(b_1,b_2……b_n)。经过一番推导,令:pi=2∑j=1...

2018-03-26 15:12:36 326

原创 洛谷P3389:【模板】高斯消元法

题目传送门:https://www.luogu.org/problemnew/show/P3389题目分析:时隔多年(月),我终于入了高消这个坑。表示挂一发模板就跑,以后复习用。具体细节什么的还是自己YY吧,有益身心健康。CODE:#include<iostream>#include<string>#include<cstring&...

2018-03-26 11:07:46 299

原创 count (类插头DP+矩阵快速幂)

题目大意:有n个点,编号为1~n。第i个点和第j个点之间有一条无向边当且仅当|i-j|<=k。求这个图的生成树个数。k≤5,n≤1015k≤5,n≤1015k\leq 5,n\leq 10^{15}。题目分析:Coming在他初二时的资料里找到的一道题,是我校上古大神cdc给的。我不得不吐槽:难道前几届的dalao初二就能做这种题了吗?而且题面还很恶意地给出了怎么用矩阵树定理算无向图...

2018-03-22 20:19:00 328

原创 CodeChef Union on Tree (虚树+点分治)

vjudge题面传送门:https://cn.vjudge.net/problem/CodeChef-BTREE题目分析:sro wjmzbmr这是道码农神题。首先考虑简化版的问题:如果给出一个点x,再给出一个距离d,如何求出距离x不超过d的点的个数?这可以用点分治解决。先用点分治预处理出每个连通块的所有点到其分治中心mid的深度数组f。f[mid][son][dep]=num表示分治...

2018-03-22 14:38:41 871

原创 洛谷P4067:[SDOI2016]储能表 (数位DP)

题目传送门:https://www.luogu.org/problemnew/show/P4067题目分析:一道令我心态爆炸的数位DP。一调调一天,WA不花一分钱先说一下我理解的数位DP是什么。数位DP本质上还是个DP,它里面有很多重复的子问题。但现在题面给了DP的下标一个上界限制,而我们不能直接枚举下标,所以要贴着这个上界限制来DFS。形象地做个比喻就是走楼梯,DFS的时候要紧贴着楼...

2018-03-22 08:39:20 430

原创 洛谷P3345:[ZJOI2015]幻想乡战略游戏 (动态点分治)

题目传送门:https://www.luogu.org/problemnew/show/P3345题目分析:这题简直是噩梦,调了我一个下午。首先考虑暴力怎么做:在树上随机找一个点,然后计算以该点为补给点的答案,并算出它周围的点的答案。如果它周围存在一个点答案更优,就往该点走,直到当前点为最优点就停止。这样一次询问就是n2dn2dn^2d的,其中d=20d=20d=20,因为每个点的度数...

2018-03-20 19:21:28 757

原创 BZOJ2159:Crash 的文明世界 (第二类stirling数+组合数学+树形DP)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=2159题目分析:先说一下部分分怎么拿。20%20%20\%:直接以每个点为根DFS一遍,求出其它点的深度。预处理1k1k1^k~nknkn^k统计答案。时间复杂度O(n2+nk)O(n2+nk)O(n^2+nk)。50%50%50\%:令S[node][d]=∑u∈node...

2018-03-20 08:20:58 386

原创 POJ3696:The Luckiest number (欧拉定理)

题目传送门:http://poj.org/problem?id=3696题目大意:给出一个数LLL,请求出一个最小的数ans,使得ans全由8组成(即ans=88888……8),并且是LLL的倍数。要求输出ans的位数。多组数据。1≤L≤2∗1091≤L≤2∗1091\leq L\leq 2*10^9。题目分析:一道很神的题,要是思路错了就走远了。我一开始的想法是把8提出来,然...

2018-03-19 18:49:32 414

原创 洛谷P4253:[SCOI2015]小凸玩密室 (DP)

题目传送门:https://www.luogu.org/problemnew/show/P4253题目分析:模拟赛的时候因为T2的锅,就没仔细想这题,写了个20pts的暴力状压还错了。后来改的时候发现我的code在跳到当前节点i的第一个子树未满的祖先j时,忘记查看j自身是否被点亮了QAQ。这题的正解其实也不难,主要是一些细节比较烦。由于题面给出的本身就是棵二叉树,可以直接记f[node...

2018-03-19 15:37:38 496

原创 洛谷P4155:[SCOI2015]国旗计划 (贪心+倍增)

题目传送门:https://www.luogu.org/problemnew/show/P4155题目分析:考模拟赛的时候以为是个DP,结果一直没想出来,写了个O(n2)O(n2)O(n^2)的区间DP还错了。昨晚睡觉前想了想,发现贪心就可以了。因为区间之间没有包含关系,所以将其按左端点递增排序后,右端点必定单调右移,每一次选取尽可能靠右的和当前区间有交的区间作为下一个即可。很明显任何一个...

2018-03-19 08:52:56 897

原创 BZOJ3162:独钓寒江雪 (Hash判断树同构+树形DP+组合数学)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3162题目分析:一道很厉害的题,让我知道原来Hash还可以判断两棵树的形态是否相同。这题的具体做法还是看VFK的题解吧,我表示只能orz。用简单的话概括一下题解的内容就是:1.定义重心为树的直径的中点。如果树的直径长度为偶数,就在最中间的边上加一个虚点作为重心。 2.以重...

2018-03-16 11:16:22 440

原创 hdu5909:Tree Cutting (FWT/点分治优化树形DP)

题目传送门:http://acm.hdu.edu.cn/showproblem.php?pid=5909题目大意:给出一棵树,每个节点的权值都在[0,m)[0,m)[0,m)内。现在对于每个v∈[0,m)v∈[0,m)v\in[0,m),你都应给出权值为vvv的至少有一个点的子连通块个数。一个连通块的权值定义为其所有点的权值异或和。答案模109+7109+710^9+7,多组数据。...

2018-03-15 21:33:10 360

原创 UOJ#86:mx的组合数 (Lucas定理+原根+NTT+高精度)

题目传送门:http://uoj.ac/problem/86题目分析:高精度写死人系列,我写了一个晚上才写完QAQ。一开始拿到这题没什么头绪,然后从部分分开始想。上数学课的时候忽然间发现40分的部分分就是个暴力枚举+Lucas定理。根据:Cnm=C⌊np⌋⌊mp⌋∗CnmodpmmodpCmn=C⌊mp⌋⌊np⌋∗CmmodpnmodpC_{m}^{n}=C_{\lfloor\...

2018-03-15 19:00:02 518

原创 BZOJ1951:[Sdoi2010]古代猪文 (扩展欧拉定理+Lucas定理+CRT)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1951题目分析:又是一道把几个数论算法强行拼起来的题目。题面要求的其实是g∑d|nCdnmodpg∑d|nCndmodpg^{\sum_{d|n}C_n^d}\mod p。用扩展欧拉定理可以变成g(∑d|nCdn)modϕ(p)+ϕ(p)modpg(∑d|nCnd)modϕ(p...

2018-03-15 09:19:55 271

原创 洛谷P4245:【模板】MTT (拆系数+FFT)

题目传送门:https://www.luogu.org/problemnew/show/P4245题目分析:为什么同一道题的题解我写了两篇blog?当然是因为这是两种不同的方法骗访问量啊。今天终于写了传闻已久的拆系数+FFT。大概就是令sp=⌈p–√⌉sp=⌈p⌉sp=\lceil \sqrt{p}\rceil,然后将两个多项式按除以spspsp和模spspsp的结果分成两个多项式。最...

2018-03-14 22:07:39 1054

原创 多项式乘法优化 学习笔记

今早重新看了myy的论文,又掌握了一些多项式乘法的新姿势,于是写一篇blog巩固一下QAQ。①如何用一次DFT加一次IDFT求出两个实序列A和B的卷积?这里我们只要求卷积后的结果,不需要求DFT的值,所以有一种很简便的方法:令复数序列C的实部为A,虚部为B。将其自卷,所得结果虚部的值除以2就是要求的多项式。这个十分容易证明:C2[k]=∑j=0kC[j]C[k−j]C2[k]=...

2018-03-14 11:15:17 1362 1

原创 洛谷P4245:【模板】MTT (CRT+三模数NTT)

题目传送门:https://www.luogu.org/problemnew/show/P4245题目分析:一道任意模数多项式乘法的模板题。可以写拆项+FFT,或者三模数NTT。我暂时只写了后者。具体做法是这样:先选取三个乘积在1023102310^{23}以上的便于使用NTT的模数。在这里我选的是m1=998244353=223∗119+1m1=998244353=223∗119+1...

2018-03-13 22:10:19 1061 3

原创 POJ1006:Biorhythms (暴力/CRT)

题目传送门:http://poj.org/problem?id=1006题目大意:给出四个数a,b,c,d,要你求一个大于d的最小数字ans使得ans≡a(mod23),ans≡b(mod28),ans≡c(mod33)ans≡a(mod23),ans≡b(mod28),ans≡c(mod33)ans\equiv a(\mod 23),ans\equiv b(\mod 28),ans\eq...

2018-03-13 15:00:37 331

原创 CodeVS3990:中国余数定理 2 (CRT)

题目传送门:http://codevs.cn/problem/3990/题目分析:这是一道CRT入门题。关于CRT入门我是看MashiroSky大神的blog,个人认为写得很好。然后以下的东西是给自己复习用的,简记一下,大概写得很烂QAQ:现在要求kkk个形如x≡ci(modmi)x≡ci(modmi)x\equiv c_i(\mod m_i)的同余方程组的最小非负整数解xxx...

2018-03-13 11:13:29 350

原创 洛谷P3401:洛谷树 (位运算+树链剖分/LCT)

题目传送门:https://www.luogu.org/problemnew/show/P3401题目分析:今天是植树节。晚上闲暇之余,看见米娜桑都在种洛谷树,我也就种一棵休闲一下。结果居然种了2.5h,感觉自己写数据结构的能力都退化了。关于这题有一道弱化版的题目:GDKOI2016魔卡少女。那题是询问一个区间的子区间异或和。方法是拆开二进制的每一位,然后用线段树维护异或前缀和中0,1...

2018-03-12 22:11:48 498

原创 BZOJ4407:于神之怒加强版 (数论+线性筛)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4407题目分析:又是一道老年人数论题。不妨令n≤mn≤mn\leq m。经过一番推导倒,可以得到这个:ans=∑D=1n⌊nD⌋⌊mD⌋∑d|Ddkμ(Dd)ans=∑D=1n⌊nD⌋⌊mD⌋∑d|Ddkμ(Dd)ans=\sum_{D=1}^n\lfloor\frac{n...

2018-03-11 22:01:27 301

原创 BZOJ4180:字符串计数 (后缀自动机+二分答案+矩阵快速幂)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=4180题目分析:别人眼中的好题,然而我觉得就是把SAM强套上另一个算法的题。后缀自动机的特点使其能够很好地解决本题,因为它能够识别原串的所有子串。如果从Root开始匹配模式串,而走到的当前状态并不拥有模式串下一位的转移,那么原串便没有这个子串。接下来从静态的问题开始考虑:假...

2018-03-10 08:50:12 494

原创 洛谷P4260:[Code+#3]博弈论与概率统计 (组合数学+莫队/分块)

题目传送门:https://www.luogu.org/problemnew/show/P4260题目分析:一道很好的题,既不是无脑的算法套路题,也不是单纯的推式子题。因此我讲得详细一些。比赛的时候我因为时间问题没有看这题,后来补了题面,花了一节数学课自己推出了一些东西(O(Tn)O(Tn)O(Tn)的做法)。后来看了官方题解,发现了一种关于组合数前缀和的新姿势QAQ。首先,题面给你的...

2018-03-09 19:44:59 819

原创 BZOJ1095:Hide 捉迷藏 && SPOJ Qtree IV (树的三分治/括号序列)

题目传送门: BZOJ1095:http://www.lydsy.com/JudgeOnline/problem.php?id=1095 SPOJ Query on a tree IV vjudge题面:https://cn.vjudge.net/problem/SPOJ-QTREE4#题目分析:看完题就会发现后面那题是前面那题的加强版。这两题都可以用链剖,边分和点分实现。前两种解...

2018-03-08 21:20:12 324

原创 Codeforces 809E:Surprise me! (莫比乌斯反演+虚树)

题目传送门:http://codeforces.com/contest/809/problem/E题目分析:一道套路到极点的题目。公式推导直接见https://blog.sengxian.com/solutions/cf-809e,我懒得打了QAQ。最后推出:ans=∑T=1n∑d|Tdμ(Td)ϕ(d)∑d|ai∑d|ajϕ(ai)ϕ(aj)dis(i,j)ans=∑T=1n...

2018-03-07 21:11:35 557

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除