2018年12月_Cyhlnj

原创 LOJ6066：「2017 山东一轮集训 Day3」第二题

传送门二分答案 kkk，考虑如何 hashhashhash 使得做起来方便把每个点挂在 k+1k+1k+1 级祖先上，考虑在祖先上删除这道题巧妙在于其可以对于 dfsdfsdfs 序/括号序列 hashhashhash这样在 k+1k+1k+1 级祖先上暴力删除就好了# include <bits/stdc++.h>using namespace std;typedef ...

2018-12-31 10:50:47 513

原创 UOJ#316. 【NOI2017】泳池

传送门一道 DPDPDP 好题设 qqq 为一个块合法的概率套路一恰好为 kkk 的概率不好算，算小于等于 kkk 的减去小于等于 k−1k-1k−1 的那么设 fif_ifi 表示宽为 iii 的合法的泳池面积都小于等于 kkk 的概率设 gig_igi 表示宽为 iii 的合法的泳池面积都小于等于 kkk 且最下面一行都合法的概率那么考虑转移 fff套路二强制前面的满足一定的...

2018-12-30 17:30:24 265

原创 BZOJ4162：shlw loves matrix II

传送门利用Cayley-Hamilton定理，用插值法求出特征多项式 P(x)P(x)P(x)然后 Mn≡Mn(mod P(x))(mod P(x))M^n\equiv M^n(mod~P(x))(mod~P(x))Mn≡Mn(mod P(x))(mod P(x))然后就多项式快速幂+取模最后得到了一个关于 MMM 的多项式，代入 MiM^iMi 即...

2018-12-30 10:51:56 232

原创 BZOJ3168: [Heoi2013]钙铁锌硒维生素

传送门题意给定一个满秩的矩阵 AAA ，另一个矩阵 BBB对于 AAA 的每个行向量 AiA_iAi 找到一个匹配 BBB 的行向量 BpiB_{p_i}Bpi使得 AiA_iAi 替换成 BpiB_{p_i}Bpi 后的矩阵 A′A&#x27;A′ 仍然满秩Sol如果 BpiB_{p_i}Bpi 能替换 AiA_iAi，那么说明存在一个向量 λ\lambda...

2018-12-29 17:25:11 202

传送门如果大力推单位根反演就可以获得一个 k2lognk^2lognk2logn 的好方法但是可以发现这道题就是求 (1+x)n(1+x)^n(1+x)n 的循环卷积的系数而题目中 kkk 一定是 222 的幂，所以带入 wkiw_k^iwki 求出点值然后 IDFTIDFTIDFT 即可# include &amp;amp;lt;bits/stdc++.h&amp;amp;gt;using namespace st...

2018-12-29 15:42:11 189

原创牛客Wannafly挑战赛23F 计数(循环卷积+拉格朗日插值/单位根反演)

传送门直接的想法就是设 xkx^kxk 为边权，矩阵树定理一波后取出 xnkx^{nk}xnk 的系数即可也就是求出模 xkx^kxk 意义下的循环卷积的常数项考虑插值出最后多项式，类比 DFTDFTDFT 的方法假设我们要求Ci=∑j=0n∑k=0nAjBk[(j+k)&nbsp;mod&nbsp;n=i]C_i=\sum_{j=0}^{n}\sum_{k=0}^{n}A_jB_k[(...

2018-12-29 14:38:47 289

原创 N次剩余和二次剩余

N次剩余给定 N,a,PN,a,PN,a,P，且 PPP 最好为质数可以算出 xN≡a(mod p)x^N\equiv a(mod~p)xN≡a(mod p) 的解首先可以算出 PPP 的原根 ggg解方程 gy≡b(mod p)g^y\equiv b(mod~p)gy≡b(mod p)，这个直接 BSGSBSGSBSGS设 gz≡x(mod&nb...

2018-12-28 20:25:05 1042

原创 UOJ450. 【集训队作业2018】复读机

传送门d=1d=1d=1 输出 knk^nknd=2d=2d=2，构造生成函数，就是求(∑i=0∞[2∣i]eii!)k[xn]=(ex+e−x2)k(\sum_{i=0}^{\infty}[2|i]\frac{e^i}{i!})^k[x^n]=(\frac{e^x+e^{-x}}{2})^k(i=0∑∞[2∣i]i!ei)k[xn]=(2ex+e−x)k直接二项式定理展开求 nnn...

2018-12-27 22:32:49 443

原创 BZOJ3328: PYXFIB(单位根反演?)

传送门Sol设A=[1110]A=\begin{bmatrix}1 &amp;amp; 1 \\ 1 &amp;amp; 0\end{bmatrix}A=[1110]那么要求的相当于是∑i=0n[k∣i](ni)Ai\sum_{i=0}^{n}[k|i]\binom{n}{i}A^ii=0∑n[k∣i](in)Ai求出其中的 A0,0A_{0,0}A0,0 即可引入单位根(...

2018-12-27 21:23:42 400

原创 BZOJ4903: [Ctsc2017]吉夫特

传送门可以发现，(nm)≡1(mod 2)\binom{n}{m}\equiv 1(mod~2)(mn)≡1(mod 2) 当且仅当 m and n = mm~and~n~=~mm and n = m即 mmm 二进制下为 nnn 的子集那么可以直接写一个 3183^{18}318 的枚举子...

2018-12-27 09:13:50 323

原创 BZOJ4671：异或图

传送门直接求连通的不好做，考虑容斥设 gig_igi 表示至少有 iii 个连通块的方案数，fif_ifi 表示恰好有 iii 个的那么gx=∑i=xn{xi}fi&amp;ThickSpace;⟺&amp;ThickSpace;fx=∑i=xn(−1)i−x[xi]gig_x=\sum_{i=x}^{n}\begin{Bmatrix}x \\ i\end{Bmatrix}f_i\if...

2018-12-26 22:47:30 237

原创容斥与反演

反演设Fn∑i=0nAn,iGiF_n\sum_{i=0}^{n}A_{n,i}G_iFni=0∑nAn,iGiGn∑i=0nBn,iFiG_n\sum_{i=0}^{n}B_{n,i}F_iGni=0∑nBn,iFi下面的直接带入到上面Fn=∑i=0nAn,i∑j=0iBi,jFj=∑i=0nFi∑j=inAn,jBj,i=FnF_n=\sum_{i=0}^{n}A_{...

2018-12-26 19:35:51 282 1

原创 Luogu3307：[SDOI2013]项链

传送门求每个珠子的方案数即有序的求三元组 (x,y,z),x,y,z≤a(x,y,z),x,y,z\le a(x,y,z),x,y,z≤a 满足 gcd(x,y,z)=1gcd(x,y,z)=1gcd(x,y,z)=1设 GiG_iGi 表示 iii 个小于等于 aaa 的有序数字，满足 gcd=1gcd=1gcd=1 的方案数容斥得到要求的16(G3+2G2+3G1)\frac{1}...

2018-12-26 15:44:00 210

原创 BZOJ4659：lcm

传送门题目所给的不合法的条件可以转化为∃p,p2∣gcd(a,b)⇔μ(gcd(a,b))≠0\exists p,p^2|gcd(a,b) \Leftrightarrow \mu(gcd(a,b))\ne 0∃p,p2∣gcd(a,b)⇔μ(gcd(a,b))̸=0那么ans=∑a=1A∑b=1B[μ(gcd(i,j))≠0]abgcd(a,b)ans=\sum_{a=1}^{A}\su...

2018-12-26 11:40:59 203

原创 51NOD1965：奇怪的式子

传送门拆开变成∏i=1nσ0(i)μ(i)∏i=1nσ0(i)i\prod_{i=1}^{n}\sigma_0(i)^{\mu(i)}\prod_{i=1}^{n}\sigma_0(i)^{i}i=1∏nσ0(i)μ(i)i=1∏nσ0(i)i考虑 ∏i=1nσ0(i)μ(i)\prod_{i=1}^{n}\sigma_0(i)^{\mu(i)}∏i=1nσ0(i)μ(i)运用...

2018-12-26 10:30:26 188

原创 BZOJ4815: [Cqoi2017]小Q的表格

传送门重点 111f(a,a+b)a(a+b)=f(a,b)ab=f(a,b−a)a(b−a)=f(a,b&nbsp;mod&nbsp;a)a(b&nbsp;mod&nbsp;a)=f(d,d)d2\frac{f(a,a+b)}{a(a+b)}=\frac{f(a,b)}{ab}=\frac{f(a,b-a)}{a(b-a)}=\frac{f(a,b~mod~a)}{a(b~mod~a)}=\...

2018-12-24 16:30:57 237 1

原创 BZOJ3512：DZY Loves Math IV

传送门Sol好神仙的题目。。一开始就直接莫比乌斯反演然后就 GGGGGG 了orz 题解permui枚举 nnn，就是求 ∑i=1nS(i,m)\sum_{i=1}^{n}S(i,m)∑i=1nS(i,m)其中S(n,m)=∑i=1mφ(ni)S(n,m)=\sum _{i=1}^m\varphi (ni)S(n,m)=∑i=1mφ(ni)设 n=∏ipicin=\prod_...

2018-12-24 15:20:56 245

原创 Luogu 4240：毒瘤之神的考验

传送门Sol分开考虑 φ(ij)\varphi(ij)φ(ij) 中 ijijij 的质因子那么φ(ij)=φ(i)φ(j)gcd(i,j)φ(gcd(i,j))\varphi(ij)=\frac{\varphi(i)\varphi(j)gcd(i,j)}{\varphi(gcd(i,j))}φ(ij)=φ(gcd(i,j))φ(i)φ(j)gcd(i,j)直接莫比乌斯反演设 g(x...

2018-12-24 12:47:31 271

原创 BZOJ3193: [JLOI2013]地形生成

传送门Sol第一问可以考虑按照山的高度从大到小放但是这样如果遇到高度相同的就不好考虑，那么同时要求数量限制从小到大这样每次放的时候后面的一定不会影响前面，并且高度相同的时候前面能放的位置后面的也能放直接乘起来就好了对于第二问，此时高度相同的会有影响对于高度相同的一段，强制要求数量限制从小到大，并且后面的位置必须小于前面设 fi,jf_{i,j}fi,j 表示放了 iii 个到 j...

2018-12-23 19:35:26 170

原创 ARG102E：Stop. Otherwise...

传送门Sol对于每个 iii ，可以把 kkk 个数字分成 (x,i−x)(x,i-x)(x,i−x) 的若干组。那么就是求每组只能其中选择一个且可以重复的方案数。预处理 f[i][j]f[i][j]f[i][j] 表示从 jjj 个组内选 iii 个，每个组必须选的方案数。f[i][j]=(f[i−1][j]+f[i−1][j−1]×2)f[i][j]=(f[i-1][j]+f[i-1...

2018-12-22 22:26:55 234

原创 51NOD1847：奇怪的数学题

传送门Sol设 f(d)f(d)f(d) 表示 ddd 所有约数中第二大的，lowdlow_dlowd 表示 ddd 的最小质因子f(d)=dlowdf(d)=\frac{d}{low_d}f(d)=lowdd那么∑i=1n∑j=1nsgcdk(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n}sgcd^k(i,j)i=1∑nj=1∑nsgcdk(i,j)=∑i...

2018-12-22 19:28:26 385

原创 CF932G Palindrome Partition

传送门Sol首先 nnn 为奇数肯定无解当 nnn 为偶数时老套路，把串 SSS 变成 S1SnS2Sn−1S_1S_nS_2S_{n-1}S1SnS2Sn−1，设为 TTT那么满足条件的 SSS 的划分相当于 TTT 中的划分，使得每一段为长度为偶数的回文串下面就只考虑 TTT 的划分设 fif_ifi 表示前 iii 个字符合法划分的方案数，用 PAMPAMPAM 可以...

2018-12-21 18:48:21 270

原创 51nod1538：一道难题(常系数线性递推/Cayley-Hamilton定理)

传送门Sol考虑要求的东西的组合意义，问题转化为：有 nnn 种小球，每种的大小为 aia_iai，求选出大小总和为 mmm 的小球排成一排的排列数有递推 fi=∑j=1nfi−ajf_i=\sum_{j=1}^{n}f_{i-a_j}fi=∑j=1nfi−aj常系数线性递推求一个满足 kkk 阶齐次线性递推数列 fif_ifi 的第 nnn 项fn=∑i=1kai×fn...

2018-12-20 22:18:19 594

原创 LOJ572: Misaka Network 与求和

传送门假设 fk(i)f^k(i)fk(i) 就是 f(i)f(i)f(i)莫比乌斯反演得到ans=∑i=1N⌊Ni⌋2∑d∣if(d)μ(id)ans=\sum_{i=1}^{N}\lfloor\frac{N}{i}\rfloor^2\sum_{d|i}f(d)\mu(\frac{i}{d})ans=i=1∑N⌊iN⌋2d∣i∑f(d)μ(di)令 g(N)=∑i=1N(f×μ)...

2018-12-20 17:29:11 334

原创 UOJ188. 【UR #13】Sanrd

传送门Sol设 fif_ifi 表示 iii 的次大质因子题目就是要求∑i=lrfi\sum_{i=l}^{r}f_ii=l∑rfi考虑求 ∑i=1nfi\sum_{i=1}^{n}f_i∑i=1nfi所求的东西和质因子有关，考虑 min25min25min25 筛的那一套理论设 s(n,j)=∑i=1n[lowi≥pj]fis(n,j)=\sum_{i=1}^{n}[lo...

2018-12-20 15:36:14 296

原创 Catalan数的通项公式(母函数推导)

首先hn=∑ihihn−i−1h_n=\sum_{i}h_ih_{n-i-1}hn=i∑hihn−i−1写出 hhh 的母函数 H(x)H(x)H(x)那么H(x)=H2(x)x+1,H(x)=1−1−4x2xH(x)=H^2(x)x+1,H(x)=\frac{1-\sqrt{1-4x}}{2x}H(x)=H2(x)x+1,H(x)=2x1−1−4x(解二元一次方程取符号时候...

2018-12-19 20:06:07 1892

原创 BZOJ2173: 整数的lqp拆分

传送门Sol构造 fibfibfib 数列的母函数 F(x)F(x)F(x)那么答案就是[xn]∑i=1∞Fi(x)=[xn]F(x)1−F(x)[x^n]\sum_{i=1}^{\infty}F^i(x)=[x^n]\frac{F(x)}{1-F(x)}[xn]i=1∑∞Fi(x)=[xn]1−F(x)F(x)而F(x)=xF(x)+x2F(x)+x,F(x)=x1−x−x2F(...

2018-12-19 17:57:41 225

原创 LOJ2541:「PKUWC2018」猎人杀

传送门Sol第一步就不会问题转化杀人后将其打上标记，仍可以以他为目标重复选，直到选到一个未打标记的人。这和原问题等价，而且这样每轮选中每人的概率都不变，只是游戏变成了无穷轮数这样就好做多了考虑容斥，枚举在 111 后面被标记的猎人集合 SSS，设其 www 的和为 AAA，总的 www 的和为 BBB，那么ans=∑S(−1)∣S∣w1B∑i=0∞(1−A+w1B)i=∑S(−1...

2018-12-18 16:17:56 221

原创 CF954I Yet Another String Matching Problem

传送门Sol考虑如何计算两个等长串的距离相当于两个匹配的字符之间连边，同一个连通块内可以互相转化，答案就是并查集合并的次数本题的字符集大小只有 666，那么考虑枚举两种字符匹配连边匹配就是一个非常套路的反转 +++ FFTFFTFFT 了# include &lt;bits/stdc++.h&gt;using namespace std;typedef long long ll;...

2018-12-13 16:57:02 374

原创拆系数FFT(任意模数FFT)

拆系数FFT对于任意模数 modmodmod设m=modm=\sqrt {mod}m=mod把多项式A(x)A(x)A(x)和B(x)B(x)B(x)的系数都拆成a×m+ba\times m+ba×m+b的形式，时a,ba, ba,b都小于mmm提出，那么一个多项式就可以拆成两个多项式的加法一个是a∗ma*ma∗m的，一个是bbb的直接乘法分配律，aaaaaa一遍，ababab一遍，...

2018-12-13 16:12:35 657

原创 BZOJ2187：fraction

Sol分情况讨论⌊ab⌋+1≤⌈cd⌉−1\lfloor\frac{a}{b}\rfloor+1\le \lceil\frac{c}{d}\rceil-1⌊ba⌋+1≤⌈dc⌉−1直接取 q=1,p=⌊ab⌋+1q=1,p=\lfloor\frac{a}{b}\rfloor+1q=1,p=⌊ba⌋+1a=0a=0a=0那么 q&amp;gt;pdcq&amp;gt; \frac{...

2018-12-09 16:58:47 304

原创 BZOJ2987：Earthquake(类欧几里德算法)

Sol设 n=⌊ca⌋n=\lfloor\frac{c}{a}\rfloorn=⌊ac⌋问题转化为求∑i=0n⌊c−axb⌋+1=∑i=0n⌊−ax+b+cb⌋\sum_{i=0}^{n}\lfloor\frac{c-ax}{b}\rfloor+1=\sum_{i=0}^{n}\lfloor\frac{-ax+b+c}{b}\rfloori=0∑n⌊bc−ax⌋+1=i=0∑n⌊b−...

2018-12-09 16:04:01 218

原创 BZOJ2229: [Zjoi2011]最小割(最小割树)

传送门最小割树算法初始时把所有点放在一个集合从中任选两个点出来跑原图中的最小割然后按照 sss 集合与 ttt 集合的归属把当前集合划分成两个集合，递归处理这样一共跑了 n−1n − 1n−1 次最小割可以证明图中任意一对点之间的最小割的数值都包含在这 n−1n − 1n−1 个数值当中把每次求出的最小割看成是两个点之间的边，可以建出一棵树定理1任意三点之间的最小割一定是两个相...

2018-12-09 14:40:05 402 1

原创 UOJ42. 【清华集训2014】Sum

传送门Sol(−1)a=1−2(a mod 2)=1−2a+4⌊a2⌋(-1)^a=1-2(a~mod~2)=1-2a+4\lfloor\frac{a}{2}\rfloor(−1)a=1−2(a mod 2)=1−2a+4⌊2a⌋那么原式变成 n−2∑i=1n⌊dr⌋+4∑i=1n⌊dr2⌋n-2\sum_{i=1}^{n}\lfloor d\sqr...

2018-12-09 11:14:10 307

原创 UOJ46. 【清华集训2014】玄学

传送门Sol考虑对于操作时间建立线段树，二进制分组那么现在主要的问题就是怎么合并信息你发现一个性质，就是每个修改只会在整个区间内增加两个端点那么我们二进制分组可以得到每个区间内最多只有区间长度级别段，每一段的修改都是一样的那么可以直接一层层归并上来最后询问就是二分每一个线段树的节点的询问段即可修改复杂度 Θ(nlogn)\Theta(n log n)Θ(nlogn) 询问复杂度 Θ...

2018-12-08 19:14:55 388

原创一类划分关系和指数级生成函数，多项式exp的关系

划分关系姑且这么叫着设满足性质 AAA 的集合为 SAS_ASA，每个元素有标号如果 SBS_BSB 是由若干个 SAS_ASA 组成的一个大集合设 aia_iai 表示大小为 iii 的 SAS_ASA 的个数设 bib_ibi 表示大小为 iii 的 SBS_BSB 的个数构造指数级生成函数A(x)=∑i=0∞aixii!A(x)=\sum_{i=0}^{\inft...

2018-12-08 18:50:22 551

原创有标号的DAG计数系列问题

传送门II设 fif_ifi 表示 iii 个点的答案那么枚举至少 jjj 个点的出度为 000∑j=0i(−1)j(ij)fi−j2(i−j)j=0\sum_{j=0}^{i}(-1)^j\binom{i}{j}f_{i-j}2^{(i-j)j}=0j=0∑i(−1)j(ji)fi−j2(i−j)j=0所以fi=∑j=1i(−1)j+1(ij)fi−j2(i−j)jf_i=\...

2018-12-06 17:50:14 1243

原创 CF891E Lust

传送门题目大意你有 nnn 个数 a1,a2...ana_1,a_2...a_na1,a2...an要进行 kkk 次操作每次随机选择一个数 xxx，使得答案加上 ∏i≠xai\prod_{i \neq x}a_i∏i̸=xai ，并将 axa_xax 减去 111求最后答案的期望，对 1e9+71e9+71e9+7 取模Sol设 bib_ibi 表示 iii 选择了...

2018-12-06 15:44:47 268

原创 BZOJ3812: 主旋律

传送门Sol考虑容斥强联通图反过来就是一些缩点后的 DAGDAGDAG一个套路就是对出(入)度为 000 的点进行容斥设 gS,hSg_S,h_SgS,hS 分别表示选了奇数个 000 入度和偶数个的，集合为 SSS 的方案数那么通过钦定一个特殊的点 uuu 有gS=∑T⊂S,u∈TfThS−Tg_S=\sum_{T\subset S,u \in T}f_Th_{S-T}gS...

2018-12-04 11:02:53 311

原创 BZOJ5372: PKUSC2018神仙的游戏

传送门Sol自己还是太 naivenaivenaive 了，上来就构造多项式和通配符直接匹配，然后遇到 borderborderborder 相交的时候就 gggggg 了神仙的游戏蒟蒻还是玩不来一个小小的性质：存在长度为 lenlenlen 的 borderborderborder 的充要条件是 ∀i,si=sn−len+i\forall i,s_i=s_{n-len+i}∀i,si...

2018-12-01 22:44:10 232

空空如也

空空如也