- 博客(101)
- 收藏
- 关注
RSS订阅原创 IOI2018退役记 + NOIP2018游记
图先鸽着。Day1在空中感受了一下日本的特别高的森林覆盖率。认识了向导，据说878787个国家至少有13\frac{1}{3}31的向导是中国人(向导都是筑波大学的留学生)，没有英文水平的选手甚至觉得中文比英文好用(捂脸)。官方发了两件T恤，一件雨衣，一个made in China\texttt{made in China}made&
2018-11-22 00:12:23
7047
15
原创 NOI2018同步赛垫底记
NOI 2018 Day 1打出GGGG\texttt{GG}看完题感觉AAA和CCC都挺可做的，稍微想了想感觉都能编出来做法，就先去写AAA了。AAA是个普及组并查集，几乎全场都过了，但是是多组数据。写完之后直接命令行测了测每个样例，都没什么问题就扔了。看了看BBB，感觉应该先打个表。打完表发现所有满足条件的序列的最长下降子序列都不超过222，然后感觉就是这个结论了。然后开始逐位...
2018-07-18 21:00:11
2145
2
原创 质因数分解: Pollard's Rho Algorithm and Quadratic Sieve Algorithm
好像很久没写博客了…随便写点东西吧。在很多数论题里都要用到整数分解，大家好像正常使用的都是 Pollard RhoPollard Rho\texttt{Pollard Rho} 算法，但是大家好像都不太会这个算法的正确写法(包括之前的我)，所以在这也算是普及一下吧。Pollard Rho Algorithm首先先考虑一下我们要做的事：我们想要找到一个 nn n 的非平凡因子 1&...
2018-07-16 18:03:24
2161
3
原创 AHOI & HNOI 2018 游记
紧张刺激的胡策+省选就从这一天开始了Day 0出发去合肥，感觉这条路走过很多次了，快到合肥的时候下雨了。 吐槽：为啥没跟大巴车走大巴车就快如闪电连服务区都不进？？？？今年的键盘比 noip 更差了，ahoi 主办方提供了 3 种优秀的键盘可供选择： 1. 正常键盘，仅限高水平选手 Marco_L_T 拥有。 2. 大 Enter 键盘 A，键盘敲不动，被 Enter 挤掉的 OR...
2018-04-27 15:14:56
1497
原创 AtCoder Peterzavodsk Contest 001 题解
A. Two integers求是 xx x 的倍数且不是 yy y 的倍数的一个数 zz z ，要求 zz z 不超过 10181018 10^{18} 。#include<bits/stdc++.h>#pragma GCC optimize("no-stack-protector")#pragma GCC optimize("-funroll-loops...
2018-02-13 14:35:43
1299
3
原创 清华集训2017题解
ioi赛制好评QAQ好像啥都记不得了就不写游记了Day 11.1 生成树计数(tree.cpp)1.1.1. 题意 有 n≤30000n≤30000 n \le 30000 个点，第 ii i 个点有点权 aiai a_i ，设第 ii i 个点的点度为 didi d_i ，那么求求所有生成树的权值和。一个生成树的权值 w(T)w(T) w(T)...
2018-02-12 21:43:44
2242
原创 [hdu4656] Evaluation - FFT,CRT
式子真难推…还是看了叉姐的课件#include"bits/stdc++.h"using namespace std;typedef long long LL;const int Mod=1000003,m1=998244353,m2=1004535809,pwMod=Mod-1;const LL P=1002772198720536577LL;inline int CR
2017-04-27 21:38:48
1112
2
原创 [BZOJ3512] DZY Loves Math IV -- 杜教筛,搜索
#include"bits/stdc++.h"#define last last_using namespace std;typedef long long LL;const int P=1000000007;struct node{ int n,m; inline LL hash(){return (LL)n<<30|m;} node(int _n=0,int _m=0){
2017-04-19 23:11:23
617
原创 [BZOJ1187][HNOI2007] 神奇游乐园 - 插头DP
这么裸的题调得我一脸懵逼 - - 我太弱啦#include"bits/stdc++.h"#define F(i,l,r) for(register int i=l;i<=r;i++)#define readi(x) scanf("%d",&x)using namespace std;const int N=105,M=7,S=3005,inf=0x3f3f3f3e;int dp[N]
2016-11-04 00:26:27
519
原创 [BZOJ1435] [ZJOI2009] 多米诺骨牌 - 插头DP/状压DP，容斥原理
题意简述：求n*m的带障碍棋盘中，摆放1*2或2*1的多米诺骨牌且相邻行和相邻列均存在骨牌跨立的方案数。n,m 本蒟蒻的辣鸡题解表示校内OJ有毒……n=m=15就完美地因为常数大了一倍而T掉了。。所以有组特判不要在意……这里不应该写插头的应该写个状压，状态就能咔嚓掉一半。。时间复杂度是大玄学 = = 真难调 弱智+2#include"bits/stdc++.h"#define
2016-11-02 17:55:52
2739
原创 [BZOJ2671] Calc和与积 - 容斥原理,莫比乌斯反演定理
本蒟蒻的题解#include"bits/stdc++.h"using namespace std;typedef long long ll;#define regi register intconst int N=50005;int pri[N],cnt,mu[N],n;bool f[N];void Linear_Shaker(){ f[1]=true;mu[1]=1; for(
2016-10-16 00:05:30
566
原创 [BZOJ4032] [HEOI2015] 最短不公共子串 - 后缀自动机,序列自动机,BFS
4032: [HEOI2015]最短不公共子串Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 398 Solved: 189[Submit][Status][Discuss]Description 在虐各种最长公共子串、子序列的题虐的不耐烦了之后，你决定反其道而行之。一个串的“子串”指的是它的连续的一段，例如bcd是ab
2016-09-18 13:53:00
883
原创 [BZOJ4028] [HEOI2015] 公约数数列 - 分块
大分块是nsqrt(n)logn的qwq#include"bits/stdc++.h"#define mmap map::iterator#define mpar make_pair#define fge getchar()using namespace std;typedef long long ll;templatevoid read(MyInt&x){ x=0;cha
2016-08-25 21:28:02
475
原创 [BZOJ2049] [SDOI2008] 洞穴勘测
LCT模板题。#include"bits/stdc++.h"using namespace std;const int N=10005; struct node{ node *fa,*c[2];bool rev;node(); inline bool go(){return fa->c[0]==this||fa->c[1]==this;} inline void pushdown
2016-08-22 19:24:22
232
原创 [BZOJ3262] 陌上花开 - CDQ分治
对第一关键字排序，分治每朵花，合并的时候通过归并左右的花的第二关键字，将值插入树状数组中求值。对于相同的花，先进行预处理即可。#include"bits/stdc++.h"#define lowbit(x) ((x)&-(x))using namespace std;const int N=100005,M=200005;int n,m;int sum[M];inline v
2016-08-21 21:42:57
444
原创 循环卷积 -- 第K大C - FFT
虽然这题我并没有在其他题站上找到，但是由于实在理解不能我为什么会奇奇怪怪地理解了就是不是很懂，外加一些奇奇怪怪的原因，所以就把贴在这了……题意：给定数组a[N][M] 和b[N]，定义：其中N是质数首先我们显然可以对b进行分类，将a拆成m个数组，那么我们就可以将问题转化为如下问题：求x[i]=∑y[j]z[i*j mod n]那么由于n是素数，那么我们求出n的原根g，则对于
2016-08-20 20:34:07
3216
原创 整体二分&CDQ分治：[BZOJ2527][POI2011] meteors [BZOJ3295][CQOI2011] 动态逆序对
首先我们先说说3295题的树套树做法：显然只需要求出每次修改x左边比自身大的数的个数和右边比自身小的个数，就是x的贡献，每次减去贡献再删除x即可，这样就相当于带修改的可持久化线段树就能解决问题。显然，我们用树状数组套权值线段树来解决这个问题，每次修改就相当于修改棵可持久化线段树上的一个权值，所以在修改上总的时间复杂度是，下面来考虑查询。我们每次查询显然可以分成段，这样我们显然可以来完成所有的查
2016-08-15 23:00:55
577
1
原创 [BZOJ4523] [CQOI2016] 路由表 - Trie - 单调栈
对于加入就直接跑一遍在最后的结点标记时间就行了，询问就是将整个段维护一个单调栈，保证每个时间所在的位置之后没有比其小的时间即可，具体看代码。#include"bits/stdc++.h"using namespace std;const int L=40000005;char _buff[L];int _pos;void read(){fread(_buff,L,1,stdin);}
2016-08-12 22:45:31
403
原创 [BZOJ4522] [CQOI2016] 密钥破解 - Pollard-Rho - exgcd
4522: [Cqoi2016]密钥破解Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 459 Solved: 223[Submit][Status][Discuss]Description 一种非对称加密算法的密钥生成过程如下：1.任选两个不同的质数p,q2.计算N=pq,r=(p−1)(q−1)3.选取小于r
2016-08-11 22:34:50
962
原创 [BZOJ4521][CQOI2016] 手机号码 - 数位DP
数位DP小题目，f[i][j][k][p][t][q]表示所有i位以内,以至少连续k+1个j的数开头,4和8的出现情况为p和t,满足条件的情况为q的方案数。然后水一水就好了ovo#include"bits/stdc++.h"using namespace std;typedef long long ll;ll f[12][10][3][2][2][2],l,r;void pre
2016-08-11 14:28:44
263
原创 [BZOJ4519] 不同的最小割 - 分治,最小割(Gomory-Hu Tree)
题目大意：求一张图的不同的最小割个数。这里我们用Gomory-Hu Tree的思想但却并不需要建立那样的一棵最小割树。首先我们要知道最小割不会相互跨立（ZZT不会证233），然后就可以每次随机选取两点，求出最小割后分治两边的点集即可，根据最大流最小割定理，直接跑一遍dinic就可以了。时间复杂度上界O(n^2mlog n)，实际复杂度完全可以接受。#include"bits/stdc++.
2016-08-10 17:42:00
1303
原创 [BZOJ4517] [SDOI2016] 征途 - 斜率优化DP
题目大意：将n个数分成m段，每段求和，求这m段和的方差的最小值*m^2由方差公式 可进一步推导到 （否则没法优化 = =），所以我们只要求出F[i][j]代表前j个数分成i段，这i段的和的平方的最小值即可，这样这就是一个愉快的经典问题了。显然f(i,j)=max{f(i-1,k)+(s[j]-s[k])^2} 令ans(k)=f(i-1,k)+(s[j]-
2016-08-09 21:36:11
394
原创 [BZOJ4516] [SDOI2016] 生成魔咒 - 后缀数组/后缀自动机
感觉还是SAM比较好写QAQ。显然ans就是∑len[i]-len[fail[i]]#include"bits/stdc++.h"#define mmap map::iteratorusing namespace std;const int N=200005;struct SAM{ mape[N];long long res; int len[N],fail[N],cnt,
2016-08-08 20:52:18
310
原创 [BZOJ4515][SDOI2016] 游戏 - 树链剖分 - 半平面交 - 标记永久化
4515: [Sdoi2016]游戏Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 390 Solved: 177[Submit][Status][Discuss]DescriptionAlice 和 Bob 在玩一个游戏。游戏在一棵有 n 个点的树上进行。最初，每个点上都只有一个数字，那个数字是 12345678912
2016-08-08 17:55:34
535
原创 [BZOJ4514] [SDOI2016] 数字配对 - 费用流
显然可以根据质因子个数建二分图，再连源点和汇点，连边就可以跑费用流了。每次跑最大费用，直到费用小于0就结束过程，记录总流量即可 - -#include"bits/stdc++.h"using namespace std;typedef long long ll; const int N=40005,M=50000;const int inf=(int)1e9;const ll ll
2016-08-04 10:53:00
397
原创 [BZOJ4513] [SDOI2016] 储能表 - 数位DP
数位DP博大精深阿弥陀佛。（回归正题）woc这取模要死，乘法到处爆炸，一不小心就炸longlong，考试肯定fst惨挂只有暴力分。对于子问题F(n,m,k)求∑[0首先将n和m自减后同k一起二进制分解，统一长度后按从高位到低位排，然后我们就可以DP了。这里DP[t][q][k][y][p]表示：t表示当前进行到的位数，即这个子问题求解到了分解出的二进制的前t位。当p=0，这
2016-08-03 00:14:32
699
原创 [BZOJ4036] set - 子集和变换 - 容斥原理/莫比乌斯反演定理
虽然有点长还是建议看这篇博客，大意就是通过子集和变换构造出函数F使得可以快速计算∑[0,+∞)F（当其收敛），然后通过消无穷，最后得出式子就可以辣。题解 #include"bits/stdc++.h"using namespace std; const double eps=1e-6;const int N=1048576+5;double a[N],ans;int n,m,o[N]
2016-08-01 22:36:53
2060
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人