- 博客(23)
- 收藏
- 关注
RSS订阅原创 【CSP 2019游记】
还是把游记补一下吧Day - 2考了场模拟赛三道暴力rk1rk1rk1，感觉rprprp要被耗完(flag∗1)(flag*1)(flag∗1)Day - 1打了一天板子心态稳健(flag∗2)(flag*2)(flag∗2)Day 0到酒店开电脑颓废发现网不错就开网盘...
2019-11-22 18:48:56
236
原创 【CSP-S 2019模拟题解】
T1：对每个TTT预处理一下到其他TTT的距离然后状压dpdpdp即可#include<bits/stdc++.h>using namespace std;#define cs const#define pb push_back#define pii pair<int,int>#define fi first#define se second#define...
2019-11-14 18:28:59
283
原创 【CSP-S 2019模板复习2】
虽然Stargazer\mathrm{Stargazer}Stargazer很大可能要退役了但还是把复习的模板汇总一个吧虽然显然有些cspcspcsp根本不可能考可持久化平衡树#include<bits/stdc++.h>using namespace std;#define cs const#define re register#define pb push_back...
2019-11-14 14:17:26
127
原创 【CSP-S 2019模板复习 1】
虽然Stargazer\mathrm{Stargazer}Stargazer很大可能要退役了但还是把复习的模板汇总一个吧虽然显然有些cspcspcsp根本不可能考割点割边注意对rootrootroot的特判#include<bits/stdc++.h>using namespace std;#define cs const#define re register#def...
2019-11-13 19:14:30
182
原创 【洛谷P4245】【模板】 任意模数NTT（MTT）
传送门首先是一种只用两次DFTDFTDFT的卷积写法若求A(x)∗B(x)A(x)*B(x)A(x)∗B(x)设P(x)=A(x)+iB(x)P(x)=A(x)+iB(x)P(x)=A(x)+iB(x)Q(x)=A(x)−iB(x)Q(x)=A(x)-iB(x)Q(x)=A(x)−iB(x)那么如果求出来P(x)=∑P(x)=\sum_{}P(x)=∑...
2019-11-09 16:38:00
260
原创 【LOJ #2271】「SDOI2017」遗忘的集合（多项式Ln+任意模数NTT）
这里写自定义目录标题欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入欢迎使用Ma...
2019-11-08 23:38:34
207
原创 【LOJ #2320】「清华集训 2017」生成树计数（生成函数+分治NTT+多项式Exp）
传送门参考的就是lojlojloj讨论里面jokerjokerjoker写的做法首先答案和度数有关，考虑树的pruferpruferprufer序列先不考虑那个∑idim\sum_{i}d_i^m∑idim对每个连通块构造生成函数ft(x)=∑iati+1(i+1)mxii!f_t(x)=\sum_{i}a_t^{i+1}(i+1)^m\frac{x^i}{i!}ft(x)=∑i...
2019-11-07 23:38:15
279
原创 【LOJ#2249】「NOI2014」购票（点分治+斜率优化DP）
传送门显然的斜率优化DPDPDP第一眼想到的是用可持久化数组维护单调栈但是觉得懒得写就放弃了点分治后用重心向上的链更新重心子树感觉也可以看做树上CDQCDQCDQ一样的东西代码写的我自己很丑#include<bits/stdc++.h>using namespace std;#define cs const#define re register#define pb...
2019-11-07 22:12:03
165
原创 【LOJ #2540】「PKUWC2018」随机算法（状压DP）
传送门设f[i][s]f[i][s]f[i][s]为考虑排列里第iii个点，当前独立集为sss的方案数暴力枚举下一个点转移即可可以加入独立集的一定还没考虑过复杂度O(n22n)O(n^22^n)O(n22n)但是转移次数很少跑的飞快#include<bits/stdc++.h>using namespace std;#define cs const#define re...
2019-11-07 22:08:50
188
原创 【LOJ #2290】「THUWC 2017」随机二分图（状压DP）
传送门暴力的做法是设f[s1][s2]f[s1][s2]f[s1][s2]为左右匹配状态转移考虑对于第二类拆成两条独立的，出现概率为50%50\%50%的边但是一起出现的概率少了25%25\%25%可以再加一个一起出现的25%25\%25%的情况第三类类似加一个一起出现为−25%-25\%−25%的情况然后暴力dpdpdp，用mapmapmap存状态每次匹配标号最小的一个点这样...
2019-11-07 22:05:44
181
原创 【LOJ #2136】「ZJOI2015」地震后的幻想乡（状压DP）
传送门考虑对于一个边集∣S∣|S|∣S∣如果刚好加入第iii条边的时候联通那么贡献就是∣S∣m+1∗P∣S∣\frac{|S|}{m+1}*P_{|S|}m+1∣S∣∗P∣S∣P∣S∣P_{|S|}P∣S∣表示SSS中的边恰好是最小的∣S∣|S|∣S∣条的概率设fs,if_{s,i}fs,i表示点集SSS联通，连了iii条边的方案数ggg为不连通的方案数那么有fs,i+gs...
2019-11-07 21:59:34
132
原创 【CSP-S 2019模拟】题解
noipnoipnoip组：T1：简单的莫比乌斯反演之后得到ans=∑T=1nnTmTf(T)ans=\sum_{T=1}^{n}\frac{n}{T}\frac m Tf(T)ans=∑T=1nTnTmf(T)其中fi=∑d∣i(dB)μ(id)f_i=\sum_{d|i}{d\choose B}\mu(\frac i d)fi=∑d∣i(Bd)μ(di)设gi=(iB...
2019-11-06 17:16:18
165
原创 【UOJ #50】【UR #3】—链式反应（生成函数+分治NTT/多项式Exp+常微分方程）
传送门首先这个题面就很胃疼而且感觉讲的不是很清楚实际上是要求满足如下条件的树的个数：对于每个非叶节点，有c+2c+2c+2个儿子，其中有c∈Ac\in Ac∈A个叶子节点和2个非叶节点且点的编号满足父亲小于儿子按照套路设f[i]f[i]f[i]为iii个点的答案可以列出dpdpdp式f[i]=12∑j∑k[i−j−k−1∈A](i−1j)(i−j−1k)fjfkf[i]=\frac...
2019-11-05 19:48:54
431
原创 【CSP-S 2019模拟】题解
T1：显然可以用之前有次cspcspcsp模拟的方法，在转移出去的里面选方案最大的作为重儿子倍增由于T3T3T3写太久了根本没时间写#include<bits/stdc++.h>using namespace std;#define cs const#define pb push_back#define bg begin#define pii pair<int,...
2019-11-05 18:30:10
199
原创 【洛谷P4233】— 射命丸文的笔记（竞赛图+多项式求逆）
传送门考虑强连通竞赛图哈密顿回路总数很好求为(n−1)!2(n2)−n(n-1)!2^{{n\choose 2}-n}(n−1)!2(2n)−n现在问题是求强连通竞赛图数设fif_ifi为iii个点的图数fn=2(n2)−∑i=1n−1fi(ni)2(n−i2)f_n=2^{{n\choose 2}}-\sum_{i=1}^{n-1}f_i{n\choose i}2^{{n-i\ch...
2019-11-02 17:55:21
175
原创 【LOJ #2264】「CTSC2017」吉夫特（Lucas定理+Dp）
传送门用LucasLucasLucas定理展开后发现就是不能有二进制位前111后000的情况可以直接从后往前枚举之前的子集dpdpdp3log3^{log}3log也跑得过去不过可以分块优化f[i][j]f[i][j]f[i][j]表示前292^929位为iii，后292^929维为jjj的子集的方案数#include<bits/stdc++.h>using namespa...
2019-11-02 17:47:55
156
原创 【LOJ #3120】「CTS2019 | CTSC2019」珍珠（生成函数+NTT）
传送门奇数的EGFEGFEGF为ex−e−x2\frac{e^x-e^{-x}}{2}2ex−e−x，偶数为ex+e−x2\frac{e^x+e^{-x}}{2}2ex+e−x考虑枚举奇数的个数ans=n!∑i=0n−2m(ex−e−x2y+ex+e−x2)D[xn][yi]ans=n!\sum_{i=0}^{n-2m}(\frac{e^x-e^{-x}}{2}y+\frac{e^x+e...
2019-11-02 17:41:20
185
原创 【LOJ #2541】「PKUWC2018」猎人杀（容斥+分治NTT）
传送门能发现对于一个死的人计算的时候不忽略他的概率而看做遇到他继续射下一个不会影响答案考虑活、死的人的概率之和为w1,w2w1,w2w1,w2那么下一次射杀iii的概率为wiw1\frac{w_i}{w1}w1wi如果把死的人也算进去那么x=w2w1+w2x+wiw1+w2x=\frac{w2}{w1+w2}x+\frac{w_i}{w1+w2}x=w1+w2w2x+w1+w2w...
2019-11-02 17:26:29
118
原创 【LOJ #2346】「JOI 2016 Final」断层（线段树）
传送门把坐标(x,y)→(x−y,x+y)(x,y)\rightarrow(x-y,x+y)(x,y)→(x−y,x+y)后发现就是往下和右移2l2l2l长度对于x,yx,yx,y分别维护一颗线段树线段树上二分开始移的位置区间加即可#include<bits/stdc++.h>using namespace std;#define re register#define ...
2019-11-02 17:14:11
238
原创 【LOJ #2106】「JLOI2015」有意义的字符串（矩阵快速幂）
传送门由于有个奇妙的限制b2≤d<(b+1)2b^2\le d<(b+1)^2b2≤d<(b+1)2所以(b−d)/2)n(b-\sqrt d)/2)^n(b−d)/2)n一定是一个小于111的小数加上这个东西((b+d)/2)n+(b−d)/2)n((b+\sqrt d)/2)^n+(b-\sqrt d)/2)^n((b+d)/2)n+(b−d)/2)n一定是...
2019-11-02 17:10:06
202
原创 【LOJ #2325】「清华集训 2017」小 Y 和恐怖的奴隶主（矩阵快速幂）
传送门沙雕题打一下发现状态很少可以矩乘优化但是复杂度不对考虑一个行向量乘矩阵是n2n^2n2的于是把2k2^k2k的矩阵预处理出来每次乘即可有些卡常#include<bits/stdc++.h>using namespace std;#define re register#define pb push_back#define cs const#define ...
2019-11-02 17:07:00
256
原创 【CSP-S 2019模拟】题解
T1：O(n2)O(n^2)O(n2)枚举后用逆元算第三个即可用hash tablehash\ tablehash table存即可脑残写了个奇奇妙妙的n2logn^2logn2log wa3\ wa3 wa3个T3T3T3个#include<bits/stdc++.h>using namespace std;#define cs c...
2019-11-02 14:26:08
154
原创 【CSP-S 2019模拟】题解
T1：相当于是一个dp of dpdp\ of\ dpdp of dp，第二维背包只用存前kkk项即可看起来像2knk2^knk2knk的但实际上有用状态很少#include<bits/stdc++.h>using namespace std;#define re register#define pb push_back#def...
2019-11-01 16:11:14
191
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人