- 博客(138)
- 收藏
- 关注
RSS订阅原创 HDU2065："红色病毒"问题（快速幂+指数生成函数）
题面 题意：有A，B，C，D 4种字母 A，B只能用偶数个 问组成长为n的串的方案数可重组合问题用指数生成函数 无限制的指数生成函数为exexe^x 只能选偶数的指数生成函数为(ex+e−x)/2(ex+e−x)/2(e^x+e^{-x})/2乘起来为e4x+2e2x+14e4x+2e2x+14\frac{e^{4x}+2e^{2x}+1}{4}泰勒展开，第n项系数大概是4n...
2018-03-13 22:11:38
416
原创 BZOJ3027：Sweet/BZOJ1272：Gate Of Babylon（生成函数+广义二项式定理+卢卡斯定理）
题面 首先问题可以差分对于每一罐都搞一个生成函数 大概是1−xm+11−x1−xm+11−x\frac{1-x^{m+1}}{1-x} 总方案数的生成函数就是 1(1−x)n∏i=1n(1−mi)1(1−x)n∏i=1n(1−mi)\frac{1}{(1-x)^n}\prod_{i=1}^{n}(1-m_i) =(1−x)−n∏i=1n(1−mi)=(1−x)−n∏i=1n(1−mi...
2018-03-13 21:36:26
311
原创 BZOJ3510：首都（LCT+启发式Link）
题面 题意：一个森林，两个操作，①连边。②问一个连通块的重心。由于它只有连边，所以我们可以搞个类似启发式的东西考虑一个小树的u和一棵大树的v连边 重心一定打大树里面 设大树的重心为ro 则新的重心一定在ro到v的路径上 并且移动的步数不会超过小树的siz用LCT维护子树大小 Evert(ro)，Access(v)后 在Splay上拎出一条长为小树的siz的链，在上面一步一...
2018-03-13 16:58:50
183
原创 51nod1355：斐波那契的最小公倍数（数论）
题面 题意给出n个a，问LCM{ f(a) }，f为斐波那契数。知乎靠谱的题解记住这两个路人性质就好 ①容斥求LCM lcm{S}=∏T⊆S,T≠∅gcd{T}(−1)|T|+1lcm{S}=∏T⊆S,T≠∅gcd{T}(−1)|T|+1\text{lcm}\{S\}=\displaystyle\prod_{T\subseteq S,T\ne \,\emptyset} \gcd\{T...
2018-03-13 09:56:56
250
原创 洛谷3401：洛谷树（树链剖分+机巧的位运算）
题面 题意：一棵树，支持边修改，问路径的子路径异或和。异或运算每位独立 gdkoi魔卡烧酒Sakura的套路 区间的子区间的异或和为区间异或前缀和 中0的个数乘1的个数 因为在前缀和中的0和1可以任意配对，也只有01配对能产生贡献 用线段树维护这个就行了根据Kscla的扫盲教育 长为len的区间有len+1个异或前缀和 若区间的异或和为0，则前缀和0/1的个数和后缀和0/1的...
2018-03-13 09:11:31
282
原创 BZOJ4652：循环之美（数学+杜教筛）
题面 这题天秀 题意：求有多少个互质的x，y，x≤n，y≤m 使得xyxy\frac{x}{y}的小数在k进制下从第一位开始无限循环我看到这题，突然发现和我小学5年级奥数学的内容很像（我小学好厉害）当时讲了怎么把一个循环小数化成分数 从第一位开始循环的小数 循环节长度就是分母中9的个数 分子就是循环节 (证明就是等比数列求和，然后求极限，好像当时也讲了） 比如0.12345...
2018-03-12 20:34:22
283
原创 BZOJ2127：happiness（最小割）
题面 题意：一群人，每人选文选里都有收益，相邻两个人同时选文和选理也有收益 问最大收益二选一，应该是网络流 网络流忘得差不多了 我刚看到的时候什么都想不出二选一能处理的问题 ①两个人选不同的有代价 ②一群人选相同某个的收益（一群人中的某人不选有代价）属于S集则选文，否则选理只考虑两个人A，B 套①，看似无法处理共同的问题 可以把共同选文的收益平分加在与S相连的边上 ...
2018-03-12 10:05:02
241
1
原创 BZOJ3209：花神的数论题（数位dp）
题面 题意：s(n)为n的二进制1的个数。 求s的前缀积，n≤1e15即对于每个i 求小于等于n的数中 二进制有i个1的有几个大概就一个数位dp了#include <iostream>#include <fstream>#include <algorithm>#include <cmath>#include <cti...
2018-03-11 20:19:59
156
原创 BZOJ4180：字符串计数（后缀自动机+Floyd+倍增+二分）
这题是打错题号找到的 题面 题意：给出一个串，用它的子串前后拼接来生成串 问所以长为n的串最少拼接次数的最大值拼接次数最少就要求每个子串都要是极长的 即它不能和后一个子串的任何前缀形成新的子串由于子串的某种性质，所以只和首字母有关 即它不能和后一个子串的首字母形成新的子串 在后缀自动机表现为该状态没有对应儿子由于要拼接次数最大，故每个状态要在关注首字母的情况下取最短 设f...
2018-03-09 21:41:58
410
原创 洛谷4260：博弈论与概率统计（组合数学+莫队/分块）
题面 题意：小L在玩游戏，赢了n场，输了m场 赢一场得1分，输一场扣1分 若当前为0分，则不会扣 问期望得分前置技能 有一个n个1和m个-1的序列，求前缀和最小值≥0的方案数 考虑不合法的 找到第一个和为-1的前缀 将其1与-1翻转 得到一个有n+1n+1n+1个1和m−1m−1m-1个-1的序列 恰好与不合法的方案一一对应类比得前缀和最小值恰好为−i−i-i的方案数为...
2018-03-09 19:28:21
530
原创 BZOJ3771：Triple（生成函数+FFT+容斥）
题面 题意：一些数，对于每个x，选3个以下的数，不计顺序，和为x的方案数。顺序就是阶乘，先不管，最后乘上就行了设生成函数为A 考虑选两个数，直接卷积后会算重复的，减去两个相同的方案数 三个的话，会算重选了两个一样的 减去后还要加上三个一样的 可以强行对于每种都选了两个一样的，再做一个生成函数B 选三个一样的，生成函数C容斥一下发现是大概这样的 A3−3A∗B+2∗CA3−...
2018-03-08 21:13:14
245
原创 BZOJ4259：残缺的字符串（FFT+卡常）
题面 题意：两个串A和B，有通配符，问A在B中出现了几次。好题（除了卡常以外） 将通配符设为0 枚举开头x 设f[x]=∑m−1i=0A[i]∗B[x+i]∗(A[i]−B[x+i])2f[x]=∑i=0m−1A[i]∗B[x+i]∗(A[i]−B[x+i])2f[x]=\sum_{i=0}^{m-1}A[i]*B[x+i]*(A[i]-B[x+i])^2若f[x]为0，则A在B[...
2018-03-08 19:54:26
244
原创 POJ3691：DNA repair（AC自动机+DP）
题面 题意：有很多个只包含4个字母的短串和一个长串 问长串至少改几个字符就不包含任何短串用短串建AC自动机，标记好所有结束状态 设f[i][j]f[i][j]f[i][j]为长度为i，后缀为状态j的最小更改数 每次枚举4个字母转移就可以了#include <iostream>#include <fstream>#include <algorithm...
2018-03-08 15:19:35
188
原创 BZOJ1408：Robot（数论）
题面 题意：以分解质因数形式给出m 考虑m的所有无平方因子的约数d 对有偶数个质因子的d和有奇数个质因子的d 分别求它们φφ\varphi的和 （不会说了，反正很简单）φ(p)=p−1φ(p)=p−1\varphi(p)=p-1 一算下去就好了#include <iostream>#include <fstream>#include <algo...
2018-03-07 21:36:52
161
原创 BZOJ3211：花神游历各国（线段树）
题面 题意：序列，区间开根，区间求和 题解：见题意为什么做这题？ 炉石玩多了，只会做一眼题了#include <iostream>#include <fstream>#include <algorithm>#include <cmath>#include <ctime>#include <cstdio>...
2018-03-07 20:42:57
197
原创 BZOJ2466：树（树形dp）
题面 题意：一棵树，选最少的点使得所有点都被覆盖奇数次 选一个点会覆盖它自身和相邻的点为点x设4个状态，分别表示x选或不选，仅考虑x的子树，x是否被覆盖 若某个点不选，则每个儿子都要被覆盖，反之亦然就可以转移了，虽然有一丢丢烦#include <iostream>#include <fstream>#include <algorithm>...
2018-03-06 22:02:41
345
原创 BZOJ2342：双倍回文（回文自动机+树上差分）
题面 题意：给你一个串，求出它的最长子串 满足该子串有两个长度相等且为偶数的回文串拼接而成该串显然是回文串，就对应回文自动机上的一个状态若某个状态的长度为4的倍数 且存在某个祖先的长度为其一半 则该状态可以贡献答案因为祖先即回文后缀和前缀 长度为其一半的前缀和后缀都是回文串 显然符合题意我的做法是树上差分#include <iostream>#in...
2018-03-05 19:20:16
195
原创 BZOJ2565：最长双回文串（回文自动机）
题面 题意：给出一个串，求出其最长子串T，使得T是两个回文串拼接而成。我的做法是在建出回文自动机时，求出每个前缀的最长回文后缀。 然后用反串在自动机上跑，求出每个后缀的最长回文前缀 就可以算答案了而串在回文自动机上跑 由于每次转移是增加了两个字符 不仅要判有么有这个儿子 还要判是不是回文#include <iostream>#include <fstre...
2018-03-05 15:37:16
239
原创 BZOJ4480：快乐的jyy（回文自动机）
题面 题意：给出两个串，问所有回文串在两个串出现次数之积之和。应该是回文自动机果题 插完一个串后重置last，再插入另一个串 每个状态对于两个串分别统计Right集的大小 就可以统计答案了#include <iostream>#include <fstream>#include <algorithm>#include <cmath&g...
2018-03-04 21:27:53
3642
原创 BZOJ3676：[Apio2014]回文串（回文自动机）
很水的我刚刚打出了听说半年的回文自动机，就开始吹了长为n的串不同的回文子串只有n个 回文自动机有两个初始状态分别表示长度为0和-1的空串 其余状态表示一个回文子串有两个属性 Right集：表示其出现位置 len：子串的长度三个指针 son[x]：在串的两端加一个字母x所形成的状态 fa：son的逆指针 pre：后缀中的最长回文串的状态与sam有些差不多的性质 Rig...
2018-03-02 22:06:36
212
原创 BZOJ4033：[HAOI2015]树上染色（树形dp）
题面 题意：n个点有边权的树，让你选择k个点为白色，其余为黑色。 问同色点两两之间距离的和的最大值。设f[x][i]f[x][i]f[x][i]为子树x里i个白点的最大贡献 包括x的父边边权的贡献然后一个一个儿子合并，字面意思转移看似n3n3n^3，其实转移时相当于在lca处枚举点对，故是n2n2n^2的#include <iostream>#include &...
2018-03-02 16:23:08
215
原创 BZOJ3295：[Zjoi2015]地震后的幻想乡（状压dp）
题面我理解的题意：n个点，m条边的图。将m条边随机排序，对于每种序列， 找到一个最小的k使得编号为1~k的边能时图联通，问k的期望除以m+1。根据大佬的理解，期望就是个积分。故可用算积分的方法算期望 设p[x]p[x]p[x]为选了x条边还不联通的概率 则k的期望为p的反常积分（大概就这个意思，乱吹的，应该是假的）而求概率就是求方案数，设cnt[S]cnt[S]cnt[S]为点集...
2018-03-02 09:22:07
230
原创 洛谷4233：射命丸文的笔记 （dp+多项式求逆）
题面 题意：求n个点的强联通竞赛图的哈密顿回路条数的期望。和城市规划那题很像，应该打表扔oeis里也会有结果吧。f[n]f[n]f[n]为n个点的强联通竞赛图的个数，根据题解，有 f[n]=2C2n−∑i=1n−1f[i]⋅C(n,i)⋅2C2n−if[n]=2Cn2−∑i=1n−1f[i]⋅C(n,i)⋅2Cn−i2f[n]=2^{C_n^2}-\sum_{i=1}^{n-1}{f[...
2018-02-28 10:02:21
412
原创 BZOJ3309 ：DZY Loves Math（莫比乌斯反演）
题面f(i)f(i)f(i)表示将iii分解质因数后的最高次幂。 求∑i=1a∑j=1bf(gcd(i,j))∑i=1a∑j=1bf(gcd(i,j))\sum_{i=1}^a \sum_{j=1}^bf(gcd(i,j)) 老年人数论题，好像没什么好说的一轮推导，有ans=∑i=1a⌊ai⌋⌊bi⌋∑d|if(d)μ(id)ans=∑i=1a⌊ai⌋⌊bi⌋∑d|if(d)μ(id)...
2018-02-27 19:57:42
205
原创 BZOJ3561：DZY Loves Math VI（莫比乌斯反演+线性筛）
题面题意：∑i=1n∑j=1mlcm(i,j)gcd(i,j)∑i=1n∑j=1mlcm(i,j)gcd(i,j)\sum_{i=1}^n \sum_{j=1}^m lcm(i,j)^{gcd(i,j)} 以下内容十分丑陋 就像背包连着胆假就是平底锅一样 这个O(n)O(n)O(n)的外层枚举连着三个幂和，它这就是个nlog(n)nlog(n)nlog(n) 这就是我用线性筛打出的...
2018-02-27 15:02:53
251
原创 BZOJ3601：一个人的数论（莫比乌斯反演+伯努利数）
题面题意：给出d和n（n以分解质因数给出），问所有与n互质的数的d次幂之和，即∑xxd[gcd(x,n)==1]∑xxd[gcd(x,n)==1]\sum_xx^d[gcd(x,n)==1]套路推♂倒=∑i|nμ(i)∗id∗∑x=1nixd=∑i|nμ(i)∗id∗∑x=1nixd=\sum_{i|n}\mu(i)*i^d*\sum_{x=1}^{\frac{n}{i}}x^d fa♂...
2018-02-26 21:03:20
546
1
原创 BZOJ5093：图的价值（数学+ntt）
题面 题意：n个有标号点的图，问所有方案所有点度数的k次方的和。前置技能：用ntt求一行第二类斯特林数 S(k,p)=∑q=0p(−1)qq!∗(p−q)k(p−q)!S(k,p)=∑q=0p(−1)qq!∗(p−q)k(p−q)!S(k,p)=\sum_{q=0}^p\frac{(-1)^q}{q!}*\frac{(p-q)^k}{(p-q)!}设点的度数为d，则它的系数为n∗2(n...
2018-02-26 15:34:27
380
原创 BZOJ4596: [Shoi2016]黑暗前的幻想乡（矩阵树定理+容斥）
题面 题意：有n个点，n-1组边，问每组中选一条边，构成一棵树的方案数。 n≤17n≤17相当于有n-1个条件，每个条件要求要在这一组中选边，求满足n-1个要求的方案数。我们可以一本正经的上容斥 设f[S]f[S]为集合S中的条件不满足（其他任意）的方案数 g[S]g[S]为只有集合S中的条件不满足的方案数f[S]f[S]能用矩阵树定理n3n^3求出来 根据套路 g[S]=
2018-01-24 19:07:40
240
原创 洛谷3790：文艺数学题（莫比乌斯反演+矩阵树定理）
题面 题意：一个nn个点，mm条带权边的图。问所有(生成树边权GCD)的和。 n≤60，m≤3000n≤60，m≤3000某个大佬曾说过，看到GCD就知道是个反演。 所以这题就是个反演。设f(x)f(x)为gcd为xx的生成树个数 g(x)g(x)为gcd为xx倍数的生成树个数 我们可以通过枚举边权是xx倍数的边，通过矩阵树定理n3n^3地算出g(x)g(x)套路反演 an
2018-01-24 15:34:13
342
原创 BZOJ3684：大朋友和多叉树（拉格朗日反演+多项式逆元+ntt）
题面 题意：求出满足两个性质的有根多叉树的数量（结点无标号，孩子有顺序） ①共有 n 个叶子结点（n ≤ 1e5） ②每个非叶结点的儿子数量∈ S（1∉S）设答案为fi，f_i，生成函数为FF 它要么是叶子，f1=1f_1=1 要么有s∈Ss∈S个儿子，则F(x)=x+∑s∈SFsF(x)=x+\sum_{s∈S}F^s设集合SS的一般生成函数为TT 则F=x+T(F)F=x+
2018-01-23 10:44:05
398
原创 51nod 1258：序列求和 V4（伯努利数+多项式求逆+NTT+CRT）
题面 题意：求Sk(n)S_k(n) n≤1e18，k≤5e4。用多项式逆元求出伯努利数，代公式就可以了。 上维基查了下伯努利数，B+B^+与B−B^-的差别只是B1B_1的正负而已.但我依然不知道什么时候用哪个，两个都试一下，发现这题用B+B^+。所以我做这题是为了学CRT的。假设题意要求模P，三模数为p1,p2,p3。我起初还天真地以为可以算出分别模p1,p2,p3的
2018-01-22 16:59:15
435
原创 洛谷3711：仓鼠的数学题（NTT+伯努利数）
题面 题意：给出a数组，求 ∑nk=0Sk(x)ak\sum_{k=0}^nS_k(x)a_k 所表示多项式的每一项系数。额，直接将伯努利数带进S里，得 =∑k=0nakk+1∑g=0kCgk+1Bgxk+1−g=\sum_{k=0}^n\frac{a_k}{k+1}\sum_{g=0}^k C_{k+1}^g B_g x^{k+1-g}=∑k=0nakk!∑g=0nBgg!xk+1−
2018-01-20 10:36:32
480
原创 51nod1228：序列求和（伯努利数）
题面题意：求幂和∑ni=0ik\sum_{i=0}^ni^k 模1e9+7 n≤1e18，k≤2e3，5000组数据曾记得L指导koi后教我们杜教筛，普及了下幂和怎么求 一本正经地推试子 原来那就是伯努利数 还说实在不会就写高斯消元由于我还很菜，还在学，就随便吹一些。当某个数列Bi{B_i}指数生成函数为xex−1\frac{x}{e^x-1}， 我们定义它叫伯努
2018-01-19 21:43:14
399
原创 HDU5659：CA Loves Substring（后缀自动机）
题面 题意：一个字符集0~9的串，将第i个字符和i+1个字符分开，得到两个串。 问这两个串中有多少个不同的子串。对于每个i都要知道答案。看到切开成两个，很自然地去考虑少了哪些串 设f[i]为切开位置为i时所少的串的个数 （大佬们说得好，正难则反）先用SAM求出总数 少的串显然是每次出现都跨过了i的。对于每次出现，当然就是Right集所反映的。对于长度为len的子串x R
2018-01-19 15:58:56
379
1
原创 HDU6063：RXD and math（数论）
题面 题意： ∑i=1nkμ2(i)∗nki−−−√\sum_{i=1}^{n^k}\mu^2(i)*\sqrt{\frac{n^k}{i}}1~n的每个数都可以写成a∗b2a*b^2，都对应1个a。a无平方因子。 循环就相当于枚举a，剩下的是算出有多少个b。 没问题，妥妥的。#include #include #include #include #include #
2018-01-19 11:11:13
241
原创 BZOJ3328：PYXFIB（数论+矩阵快速幂）
我觉得这是一道非常细的数学题。题面 题意：模P下，求∑k|i,inCin∗Fi\sum_{k|i,i F为斐波那契数，n≤1e18，p≤1e9，k≤2e4，k|(p−1)n≤1e18，p≤1e9，k≤2e4，k|(p-1)。其实数学题也是有套路的（根据本题现学） 看到菲波那契就应该是矩阵 看到整除P-1就应该是原根有关 看到组合数就应该是二项式先搞个算斐波那契的矩阵AA={1
2018-01-19 10:37:58
273
原创 [AGC005F]：Many Easy Problems（NTT+容斥）
题面 题意：给出一棵树，若S为某个点集，f[S]为最小的包含S的联通块大小，对于所有1≤k≤n，求出在n个中选k个的CknCnkC_n^k种方案中，f之和。膜924844033 。额，根据套路，应该考虑每个点对每个k算贡献。 对于点x，它不属于联通块仅当所选的k个点没有路径跨过x。 设ai为以x为根的某棵子树的大小，故贡献为Ckn−∑ai∗CikCnk−∑ai∗CkiC_n^k-\sum...
2018-01-19 08:24:16
276
原创 BZOJ2298：[HAOI2011]problem a（DP）
题面每句话等价于： ①排名属于[a+1,n-b]的人分数相等。 ②那个人排名属于[a+1,n-b]。原题变成了序列上选若干条不相交的线段，使得权值最大。 且每条线段的权值不能超过其长度。#include #include #include #include #include #include #include #include using namespace s
2018-01-18 17:01:13
254
原创 BZOJ1042：[HAOI2008]硬币购物（DP+容斥原理）
题面 题意：4种硬币。面值分别为c1,c2,c3,c4。某人去商店买东西，去了tot次。 每次带di枚ci硬币，买价值s，问多少种付款方法。若没有限制，就是一个完全背包的计数。额，考虑容斥，直观说来就是 ans=没限制的方案数-1种硬币超限的方案数+2种硬币超限的方案数-3种硬币超限的方案数+4种硬币超限的方案数写成柿子有f[S]=∑T⊇Sg[T]⇒g[S]=∑T⊇S(−1)|T
2018-01-17 22:09:45
230
原创 BZOJ3560：DZY Loves Math V（数论，线性筛）
也许我应该下放一下我的懒惰标记了。 题面题意：n≤1e5，a≤1e7，求 由于积性，对每个质数p求贡献。φ(pe)=pe∗p−1p\varphi(p^e)=p^e*\frac{p-1}{p}，对于某个a，我们掏空它里面p的幂，得到贡献 ∑ki=0φ(pi)=p−1p∗(−1+∑ki=0pi)+1\sum_{i=0}^k\varphi(p^i)=\frac{p-1}{p}*(-1+
2018-01-17 20:04:07
246
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人