- 博客(388)
- 收藏
- 关注
RSS订阅原创 【BZOJ2733】【HNOI2012】永无乡（Splay启发式合并）
Descriptionclick meSolution每次合并时考虑将较小的Splay中的节点依次插入较大的Splay中，发现每个点最多会被插入log2nlog2nlog_2n次，时间复杂度为O(nlog22n)O(nlog22n)O(nlog^2_2n)Source/************************************************ * A...
2018-02-27 20:57:58
226
原创 【Luogu3898】【BZOJ3652】【湖南集训】大新闻（数位DP？）
Descriptionclick meSolution将题目看做两个子问题： - 加密 按位计算每一位的贡献即可 - 不加密 有个比较妙的方法，参考：https://www.cnblogs.com/TSHugh/p/8476934.html 我们先假设所有数都找到了他能找到的最好的匹配(就是异或后为二进制最高位与n-1相等的最大数)并且算出其异或后的总和,然后我们按...
2018-02-27 16:42:56
437
原创 【BZOJ4802】【Pollard-Rho模板题】欧拉函数
Description给定n(n≤1018)n(n≤1018)n(n\le 10^{18}),求φ(n)φ(n)\varphi(n)Solution直接用Pollard-Rho分解质因数，然后用φ(n)=n∏(1−1pi)φ(n)=n∏(1−1pi)\varphi(n)=n\prod(1-\frac{1}{p_i})（pipip_i表示nnn的不同质因数）Source...
2018-02-26 14:21:38
372
原创 【BZOJ3329】Xorequ（数位DP，矩阵乘法）
Descriptionclick meSolution发现x⊕3x=2xx⊕3x=2xx \oplus 3x = 2x即x⊕2x=3xx⊕2x=3xx \oplus 2x = 3x，考虑异或是不进位的加法，题目条件等价于xxx的二进制表示中不存在连续的111。 第一问数位dp，第二问直接矩乘优化dp即可。Source/***********************...
2018-02-25 17:10:19
281
原创 【AGC021 D】Reversed LCS（DP）
Description给定一个字符串，可以更改其中的kkk个字符，要求最大化该字符串与其翻转字符串的LCS长度，求该长度。Solution考虑DP，设fi,j,kfi,j,kf_{i,j,k}表示从iii到jjj的字符串更改kkk个字符可达到的最大长度，直接转移即可。Source/***********************************************...
2018-02-25 10:19:13
280
原创 快速数论变换（NTT）板子
原理讲解？反正我这么菜当然只能背板子了。。原理讲解？反正我这么菜当然只能背板子了。。\color{white}{原理讲解？反正我这么菜当然只能背板子了。。}/************************************************ * Au: Hany01 * Date: Feb 23th, 2018 * Prob: [UOJ34] NTT Template *...
2018-02-23 21:48:56
458
原创 【BZOJ1833】【ZJOI2010】count 数字计数 （数位DP）
Description给定a,ba,ba,b，求出aaa到bbb之间所有的整数中所有的数字出现个数和。Solution设fi,j,kfi,j,kf_{i,j,k}为以jjj开头、长度为iii的所有数（允许前导0）中数字kkk的出现次数。 预处理出fi,j,kfi,j,kf_{i,j,k}，然后数位DP即可。 感觉做复杂了点吧 ，应该有更简单的做法。Source/****...
2018-02-19 22:45:20
269
原创 【HDU2089】不要62（数位DP入门题）
Description给定mmm,nnn，求mmm到nnn的所有整数中及不含’4’也不含’62’的数的个数。Solution直接数位DP即可。。。Source/**************************** * Au: Hany01 * Prob: [HDU2089] 不要62 * Date: Feb 19th, 2018 * Email: hany01@...
2018-02-19 20:27:38
220
原创 【BZOJ2159】Crash的文明世界（树形DP，斯特林数的性质）
Description给定一棵nnn个节点的树和定值kkk，求每个点的S(u)S(u)S(u)，S(u)S(u)S(u)的定义如下： S(u)=∑v∈treedist(u,v)kS(u)=∑v∈treedist(u,v)kS(u)=\sum_{v\in tree} dist(u, v)^kSolution感觉很巧妙的一个题，好评（除了输入要解压外）先推一波式子： S(u)=∑i...
2018-02-19 17:54:47
457
原创 莫比乌斯函数、二项式、斯特林数以及它们的反演
莫比乌斯反演引入定理另一种莫比乌斯反演证明性质：求莫比乌斯函数值例题及题解二项式反演反演公式例题bzoj3622涂色问题DescriptionSolutionStirling反演第一类斯特林数递推方式第二类斯特林数求法递推容斥性质例题CF 932EDescriptionSolutionBZOJ2159...
2018-02-18 22:35:04
2893
原创 【BZOJ3622】已经没有什么好害怕的了（DP，二项式反演）
Descriptionclick meSolution先将{ai},{bi}{ai},{bi}\{a_i\},\{b_i\}排序，求出titit_i，表示有多少个bjbjb_j小于aiaia_i。 设dpi,jdpi,jdp_{i,j}为aaa中的前iii个有jjj个大于bbb，得到转移方程： dpi,j=dpi−1,j+dpi−1,j−1×max(0,ti−j+1)dpi,j...
2018-02-17 00:00:39
472
原创 【LA2218】【POJ1755】 Triathlon（半平面交）
Description铁人三项有333段：游泳、自行车和赛跑，给出每个选手每个项目的速度vi,uiwivi,uiwiv_i,u_iw_i，每个项目的长度不确定，求每个人是否有可能获得冠军。Solution非常巧妙的一道题。 设总长度为111，游泳长度为xxx，自行车长度为yyy，那么iii强于jjj的条件是： xvi+yui+1−x−ywi<xvj+yuj+1−x−ywjxv...
2018-02-14 14:06:31
300
原创 【LA3890】【POJ3525】 Most Distant Point from the Sea, Japan 2007（二分，半平面交）
Description给定一个多边形，要求在其内部找到一个点使得该点到边界的最短距离最长，求出该距离。Solution考虑二分距离ddd，剩下的问题只有怎么判断ddd可行： 将所有的边向内推进ddd，求半平面交即可。 记得空间要多开一点，不然会RE…Source/**************************** * Au: Hany01 * Prob: ...
2018-02-13 23:00:39
309
原创 【POJ 3130】How I Mathematician Wonder What You Are! ,Japan 2006（半平面交）
Description问给定的多边形中是否存在一点CCC，使得对于多边形内任意一点PPP都有线段CPCPCP在多边形内。Solution发现题目的实质是求该多边形的各边是否存在半平面交，上板子即可。Source/**************************** * Au: Hany01 * Prob: POJ3130 How I Mathematicia...
2018-02-13 20:59:35
329
原创 【LA 4728】Square, Seoul 2009 （凸包，旋转卡壳）
Description给定很多个正方形的定点，求这些点中两两之间距离的最大值。Solution直接求出凸包后旋转卡壳即可。。Source/************************************************ * Au: Hany01 * Date: Feb 13th, 2018 * Prob: [LA4728] Square, Seou...
2018-02-13 10:22:05
253
原创 【UVa 10652】Board Wrapping （凸包）
Description给定很多块木板，用面积尽量小的多边形将其包起来。Solution将木板的四个点全部加入点集，求凸包即可。Source/**************************** * Au: Hany01 * Prob: [UVa 10652] Board Wrapping * Date: Feb 12th, 2018 * Email: h...
2018-02-12 21:37:53
241
原创 【LA 2572】Viva Confetti, Kanazawa 2002 （计算几何，圆）
Description把nnn个圆盘依次摆到桌面上，给出顺序、圆心位置和半径，求最后有多少个圆盘可见？ 数据保证轻微扰动后不改变结果Code感觉这题挺巧妙的。 可以发现所有的可见部分的轮廓都是由圆弧组成的。 对于每个圆，求出边缘上的小圆弧，取每段圆弧的中点，判断其的内外部分的最上层分别属于哪个圆盘并计数即可。Source/***********************...
2018-02-12 15:40:00
253
原创 【BZOJ 2823】【AHOI2012】【最小圆覆盖模板题】信号塔
Descriptionclick meSolution最小圆覆盖问题模板： 1. 随机排列点 2. 依次将点加入 3. 若该点iii在当前圆内，则转2；否则该点一定在点1…i1…i1\dots i的最小圆边界上，以该点为圆心，半径为000，转4 4. 将1…i−11…i−11\dots i-1的点全部重新加入，每加入一个点jjj若不在当前圆内，那么i,ji,ji,j都在边界上，...
2018-02-12 11:47:59
237
原创 【UVa 12304】2D Geometry 110 in 1! （计算几何、圆）
Descriptionhttps://vjudge.net/problem/UVA-12304Solution都是一些最基本的操作，没什么技术含量，只是代码量有点大，，码了一天，摆上代码吧。Source/************************************************ * Happy Spring Festival! * Au: Hany...
2018-02-11 21:34:31
253
原创 【UVa 11796】Dog Distance （计算几何）
Description两条狗分别沿折线奔跑，同时出发、同时到达，且匀速奔跑，求奔跑过程中两条狗的最远距离、最近距离之差。Solution模拟整个过程，将时间分段使每段时间两条狗都是直线运动，将一条狗设为参照物，求另一条狗到该线段的距离即可。开始求点到线段距离时没有判断线段退化成一个点的情况，导致调了几个小时。。。Source/**********************...
2018-02-10 15:35:22
215
原创 【LA 3263】That Nice Euler Circuits (Shanghai 2004)（计算几何）
Description给定一个一笔画的图形，求将平面分成了多少个部分。Solution统计出所有的边、点（包括交点），由欧拉公式V+F=E+2V+F=E+2V+F=E+2求得答案即可。 震惊！调了一个小时的原因竟然是点积Dot(Vector A, Vector B写错了。。Source/***************************************...
2018-02-09 23:07:56
231
原创 【UVa 11178】Morley's Theorem （计算几何）
Description Morley定理：作三角形ABC每个内角的三等分线，相交成三角形DEF，则三角形DEF为等边三角形。给定A,B,CA,B,CA,B,C，求D,E,FD,E,FD,E,F坐标。Solution计算几何入门题。 先求出∠BAC,∠ABC∠BAC,∠ABC\angle BAC,\angle ABC，除以333，将AB旋转求出的角度，得到BF,AFBF,AF...
2018-02-09 20:38:59
234
原创 【BZOJ1415】【NOI2005】聪聪可可（期望dp、记忆化搜索）
Descriptionclick meSolution考虑设fi,jfi,jf_{i,j}为聪聪在iii节点、可可在jjj节点的期望剩余时间，发现记忆化搜索比较方便： 边界：fi,i=0fi,i=0f_{i,i}=0fi,j=1(disi,j=1或2)fi,j=1(disi,j=1或2)f_{i,j}=1 (dis_{i,j}=1或2) 转移：fi,j=fp,j+∑fp,k...
2018-02-08 15:40:16
243
原创 【BZOJ4318】OSU!（期望DP）
Descriptionosu 是一款群众喜闻乐见的休闲软件。 我们可以把osu的规则简化与改编成以下的样子: 一共有n次操作，每次操作只有成功与失败之分，成功对应1，失败对应0，n次操作对应为1个长度为n的01串。在这个串中连续的 X个1可以贡献X^3 的分数，这x个1不能被其他连续的1所包含（也就是极长的一串1，具体见样例解释） 现在给出n，以及每个操作的成功率，请你输出期望分数...
2018-02-08 09:52:31
280
原创 【BZOJ1419】Red is good（期望DP）
Descriptionclick meSolution用fi,jfi,jf_{i,j}表示还剩iii张红牌、jjj张黑牌的期望值，考虑倒推。 转移：fi,j=ii+j(fi−1,j+1)+ji+j(fi,j−1−1)fi,j=ii+j(fi−1,j+1)+ji+j(fi,j−1−1)f_{i,j}=\frac{i}{i+j}(f_{i-1,j}+1)+\frac{j}{i+j}(f_...
2018-02-08 09:23:42
173
原创 【BZOJ3566】【SHOI2014】概率充电器（树形期望DP）
Descriptionclick meSolution参考博客 直接求每个点充电的概率显然不太方便（其实还是很方便的，只是我太菜了。。），所以我们考虑求每个点不能充电的概率，这包括三个方面： 1. 节点的儿子没有传电（儿子没有电或者边不连通） 2. 自己没有电 3. 父亲没有传电（父亲没有电或者边不连通） 我们用fufuf_u表示该节点没有点且儿子没有传电的概率，那么有：...
2018-02-08 08:41:19
267
原创 【BZOJ4872】【六省联考2017】分手是祝愿（概率dp）
Descriptionclick meSolution设fifif_i为从还剩iii个位置没有操作到剩i−1i−1i-1个位置没有操作的期望操作次数。 那么有：fi=in+n−in(1+fi+fi+1)fi=in+n−in(1+fi+fi+1)f_i=\frac{i}{n}+\frac{n-i}{n}(1+f_i+f_{i+1}) 解方程得：fi=1+(n−i)(1+fi+1)...
2018-02-08 08:36:17
465
原创 【BZOJ2754】【SCOI2012】喵星球上的点名（后缀数组）
Descriptionclick meSolution据说AC自动机可做？反正我使用后缀数组做的。。把所有的串连在一起，用一个很大的数分割开来。 求出后缀数组，暴力地在height上找可以匹配的串并计数即可。 （应该是可以被卡的，但是没有出题人那么恶心的吧。。）Source/********************************** * Au: Hany
2018-02-06 21:16:49
309
原创 【UOJ131】【BZOJ4199】【NOI2015】品酒大会（后缀数组）
Descriptionclick meSolution先求出height数组XJB乱推式子算出第一问至于第二问，用st表维护每个区间的最大值、最小值，每次用可计入贡献的区间中的两个最大值的积、两个最小值的积更新即可。具体实现细节参看代码。Source/********************************** * Au: Hany01 * Prob:
2018-02-05 21:37:47
232
原创 【BZOJ1076】【SCOI2008】奖励关（DP、期望、状压）
Descriptionclick meSolution套路的状压期望DP题。。。考虑倒退期望：设fi,j" role="presentation" style="position: relative;">fi,jfi,jf_{i,j}为一直到第i−1" role="presentation" style="position: relative;">i−1i−1i-1轮
2018-02-04 13:27:51
314
原创 【BZOJ1101】【POI2007】ZAP-Queries（莫比乌斯反演）
Description给定n" role="presentation">nnn,m" role="presentation">mmm,k" role="presentation">kkk，求∑x≤n∑y≤m[gcd(x,y)=k]" role="presentation">∑x≤n∑y≤m[gcd(x,y)=k]∑x≤n∑y≤m[gcd(
2018-02-03 11:08:27
323
原创 【BZOJ3994】【SDOI2015】约数个数和（莫比乌斯反演）
Description给定nnn,mmm，求∑i=1n∑j=1md(i×j)∑i=1n∑j=1md(i×j)\sum_{i=1}^{n}\sum_{j=1}^{m}d(i\times j) 其中d(x)d(x)d(x)表示xxx的约数个数。Solution一个结论： Conclusion：d(nm)=∑i|n∑j|m[gcd(i,j)=1]d(nm)=∑i|n∑j|m...
2018-02-03 10:41:16
321
原创 【UVa 10601】Cubes （群论、Burnside引理、DP）
Description给定12" role="presentation" style="position: relative;">121212根木棍的颜色（用1…6" role="presentation" style="position: relative;">1…61…61\dots 6表示），求可以拼成多少个不同的正方体。Solution发现一共有6×
2018-02-02 14:12:41
261
原创 【BZOJ1004】【HNOI2008】Cards（群论、Burnside引理、背包dp）
Descriptionclick meSolution2.1由于题目有： 保证任意多次洗牌都可用这 m种洗牌法中的一种代替，且对每种洗牌法，都存在一种洗牌法使得能回到原状态。所以加上原状态，一共有m+1m+1m+1种置换，考虑用BurnsideBurnsideBurnside引理，剩下的问题只有计算当知道有多少个循环，每个循环的大小时，求不变元的个数。这个用将...
2018-02-02 11:07:55
236
原创 【LA 3641】Leonardo's Notebook （群论）
Descriptionclick meSolution发现规律： 长度为n" role="presentation" style="position: relative;">nnn的循环自乘时，如果n" role="presentation" style="position: relative;">nnn为奇数，那么得到的还是一个长度为n" role="presentati
2018-02-02 08:58:03
245
原创 【UVa 10294】Arif in Dhaka(First Love Part 2) （群论、Polya定理）
Descriptionclick me 默默复制了vjudge的链接（逃Solution可以XJB乱模拟算的题目为什么要用群论 题中有两种置换： 1. 旋转 对于跨度为iii的旋转，则有(i,n)(i,n)(i,n)个循环、每个循环有n(i,n)n(i,n)\frac{n}{(i,n)}个元素（对称性）。 则旋转的不动点总数为：A=∑i=0n−1t(i,n)A=∑i=...
2018-02-01 20:49:17
270
原创 【UOJ 219】【BZOJ 4650】【NOI 2016】优秀的拆分（后缀数组）
Descriptionclick meSolution2.1先考虑找到所有形如AA的串的位置： 枚举长度lenlenlen，找到关键点lenlenlen,2len2len2len,3len3len3len⋯⋯\cdots 发现所有长度为lenlenlen的A的AA串都经过了其中两个关键点。 令l=k×lenl=k×lenl=k\times len，r=(k+1)×le...
2018-02-01 15:01:16
248
原创 【BZOJ 1717】【USACO06DEC】牛奶模式Milk Patterns（后缀数组）
Descriptionclick meSolution后缀数组经典题，直接在height[]上二分即可。另外，ymy巨佬用后缀自动机也A掉了这题，真是太爷啦！！！Code/**************************** * Au: Hany01 * Prob: bzoj1717 * Date: Jan 31st, 2018 * Email:
2018-01-31 15:54:14
392
原创 【51nod 1244】莫比乌斯函数之和 （杜教筛）
Description给出a" role="presentation">aaa,b" role="presentation">bbb,求∑i=abμ(i)" role="presentation">∑bi=aμ(i)∑i=abμ(i)\sum_{i=a}^{b}\mu(i)Solution杜教筛板子。。 莫比乌斯函数前缀和的推导过程： 令M(n)" ro
2018-01-30 22:02:01
369
原创 【BZOJ3930】【CQOI2015】选数（递推 & 莫比乌斯反演）
Description我们知道，从区间[L,H][L,H]（LL和HH为整数）中选取NN个整数，总共有(H−L+1)N(H-L+1)^N种方案。小z很好奇这样选出的数的最大公约数的规律，他决定对每种方案选出的NN个整数都求一次最大公约数，以便进一步研究。然而他很快发现工作量太大了，于是向你寻求帮助。你的任务很简单，小z会告诉你一个整数KK，你需要回答他最大公约数刚好为KK的选取方案有多少个。由
2018-01-18 22:04:19
427
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人