INCINCIBLE的博客

原创 题目汇总２

NO.1 HDU - 3448分类:　枚举，meet-in-the-middle简要题解:暴力枚举每个物品是否被选的时间复杂度是O(2n)O(2n)O(2^n)枚举前一半物品，处理出n2n2\frac n2 个set,其中集合SiSiS_i里面存着不多于iii 个物品的所有可能重量和再枚举后一半物品，再对应集合中二分查找即可时间复杂度降为O(n∗2n2)O(n∗2n2)O(n*...

原创 题目汇总1

NO.1 CF Gym 100738A Fitting boxes分类: 数学简要题解:设大矩形(长宽为c,dc,dc,d)水平放置,小矩形(长宽为a,ba,ba,b),与水平方向夹角为θθ\theta,那么需要满足:b∗sinθ+a∗cosθ<=da∗sinθ+b∗cosθ<=cθ∈(0,π2)b∗sinθ+a∗cosθ<=da∗sinθ+b∗cosθ<...

原创 【CQ Training 2014 Day3】燃烧的天空

Description给出若干个三角形，被i个三角形覆盖的区域总面积记为Area(i)Area(i),输出每一个Area(i)(1<=i<=n)Area(i) (1<=i<=n).Solution扫描线。 先处理所有顶点和线段交点的横坐标的集合X。 考虑处于X[i]X[i]到X[i+1]X[i+1]之间的图形对答案的贡献，可以看成是一个一个的梯形。 可以用一个变量k记录“当前正在讨论对Area

原创 6月题目汇总1（6.6---6.13）

6.6上午NO.1 来源 : UVa 1318 Monster Trap分类： 计算几何，平面直线图简要题解：首先从所有线段的端点中出不在任何一条线段（端点不算）上的点的集合S，再加上(0,0),(inf,inf)(0,0),(inf,inf)这两个点。对于SS中的两个点a,ba,b，如果线段a,ba,b不和任何给出的线段相交，那么给a,ba,b连上边。但这样会在线段的公共端点处出现问题

原创 5-6月题目汇总（5.30---6.4）

5.30晚NO.1 来源 : BZOJ 2342 [Shoi2011]双倍回文分类： 回文，Manachar简要题解：先用Manachar求出f(i)f(i),令p(i)p(i) 表示以i,i+1i,i+1 中间为对称轴能延伸的最长长度，也就是 s[i−p(i)]=s[i+1+p(i)]s[i-p(i)]=s[i+1+p(i)] ,显然p(i)=f(i∗2+1)−12p(i)= \frac{

原创 HN2015集训 永远亭的竹笋采摘

第一反应肯定是dp啦，然后就会愉快的T成一头象拔蚌。 那么，说说正解吧。 显然，选取的区间一定可以是两头为差值最小的。那就好办了，只需要预处理出所有的这样的区间，再dp即可，由于数据是随机生成的，所以这样的区间不会太多，有人验证了，好像是O(n)级别的。 预处理怎么搞？分块优化，pre[i][j]表示在第i块内的数，与j差值的最小值。每一块内的数取出来排序，两个下标跟着动就行了。 然后把那些

原创 HN2015集训 月宫的符卡序列

可能看到这个题有的人第一反应是回文自动机之类的，但是注意到这个题是跟回文串中心相关，而回文自动机则是与回文串末尾相关，可能无法处理。 那么与回文串中心相关的算法还有什么，manacher(马拉车)。 学过的同学都知道吗，本质不同的回文串最多n个，也就是拉车到i的时候，本来已有f[i]所表示的回文串是之前已经出现过的，暴力新扩展出来的才是新出现的回文串。 那么就构造出了一套拓扑序，最后反拓扑序异

原创 5月题目汇总（5.22---5.28）

5.23晚NO.1 来源： BZOJ 4197 [NOI2015] 寿司晚宴分类： 状压 dp简要题解：注意到n≤500n \le 500，至有一个大于2020 的质因数，而小于20的质数只有8个，我们可以令f[S][T]f[S][T] 表示第一个人选的小于20的质数集合为S，第一个人选的小于20的质数集合为T的方案数。把2到500的数按照大于20的质因数（没有就是1）排序，依次讨论大质因

原创 PKUSC2017 酱油记(by spark)

Day 1报了到，领了盒饭之后就去逛校园啦。嗯，不会骑自行车的蒟蒻腿都走软了。Day 2首先是数学笔试，回忆1s题目：x4−3x3+3x2+ax+bx^4-3x^3+3x^2+ax+b 可以被 x2−3x+2x^2-3x+2 整除，求a,ba,b 的值。令 x2−3x+2=0x^2-3x+2=0 解得 x1=1,x2=2x_1=1,x_2=2x1,x2x_1,x_2 也是 x4−3x3+3x2+

原创 5月题目汇总2（5.7---5.13）

5.7晚NO.1 来源： BZOJ 4565 [Haoi2016]字符合并分类： 区间,状压 dp简要题解：注意到长度为lenlen 的01串最后会被合并成长为x=(len−1) mod (k−1)+1x=(len-1)\ mod\ (k-1)+1 的01串。设f(i,j,S)f(i,j,S) 表示区间[i,j][i,j] 的01串最后变成二进制状态SS 的最大得分。讨论SS 的第一个字符

原创 Codeforces, 587F: Duff is Mad

Duff is Mad给出nn个字符串S1,S2......SnS_1,S_2......S_n ，qq 个询问。 每次询问给定l,r,kl,r,k, 要求输出 ∑ri=loccur(Si,Sk)\sum_{i=l}^{r}occur(S_i,S_k),. 其中occur(S,T)occur(S,T)表示SS在TT中出现的次数。 n,q,∑ni=1|Si|≤105n,q ,\sum _{i=

原创 5月题目汇总1 （4.30-5.6）

4.30晚NO.1 来源：[ZJOI2015] 地震后的幻想乡分类： dp, 概率与期望简要题解：orzorzorz 题解看这里 5.1上午NO.2 来源： BZOJ 4456 [Zjoi2016]旅行者分类： 分治简要题解：离线处理询问，分治平面区域。每次分治处理某一维坐标（范围较小的递归下去），下面讨论递归处理yy，讨论xx的情形。假设现在平面区域是(x1,y1)→(x2,

原创 图片暂存

原创 [ZJOI2015] 诸神眷顾的幻想乡

P3256【ZJOI2015 Day1】诸神眷顾的幻想乡时间限制 : 20000 MS   空间限制 : 524288 KB问题描述        幽香是全幻想乡里最受人欢迎的萌妹子，这天，是幽香的2600岁生日，无数幽香的粉丝到了幽香家门前的太阳花田上来为幽香庆祝生日。        粉丝们非常热情，自发组织表演了一

原创 [ZJOI2015] 幻想乡战略游戏

P3254【ZJOI2015 Day1】幻想乡战略游戏时间限制 : 120000 MS   空间限制 : 512000 KB问题描述        傲娇少女幽香正在玩一个非常有趣的战略类游戏，本来这个游戏的地图其实还不算太大，幽香还能管得过来，但是不知道为什么现在的网游厂商把游戏的地图越做越大，以至于幽香一眼根本看不过来，

原创 [ZJOI2015] 地震后的幻想乡

P3255【ZJOI2015 Day1】地震后的幻想乡时间限制 : 20000 MS   空间限制 : 262144 KB问题描述        傲娇少女幽香是一个很萌很萌的妹子，而且她非常非常地有爱心，很喜欢为幻想乡的人们做一些自己力所能及的事情来帮助他们。        这不，幻想乡突然发生了地震，所有的道路都崩塌

原创 BZOJ 1095 [ZJOI2007]Hide 捉迷藏

传送门 BZOJ 1095大致题意有一棵树，每个点是黑色或者白色，最开始全是黑点，有两种操作。 1 Get : 输出最远两个黑点之间的距离。 2 Change(x): 改变x号点的颜色。动态点分治大法好，具体解法如下。1. 构建 重心树按照静态点分治的方法，记录每一个重心xx的上一层重心··par(x)par(x) ,显然会形成一棵树高为log(n)log(n) 级别的树。 对于一个重

原创 4月题目汇总2(4.23---4.29)

4.23下午NO.1 来源： GDOI 2014 拯救莫莉斯分类： 状压DP简要题解：设f(i,S,T)f(i,S,T),表示前i−1i-1 行已经满足要求，第i−1，ii-1，i行的选取状态分别为集合S，TS，T.令cover(S)=S∪(S<<1)∪(S>>1)cover(S)= S \cup (S<<1) \cup (S>>1) 表示 SS 选这个集合中的点能覆盖的点的集合。f(i

原创 [HNOI2016] 网络

P4054[Hnoi2016 day1]网络时间限制 : - MS   空间限制 : 165536 KB 评测说明 : 2s问题描述一个简单的网络系统可以被描述成一棵无根树。每个节点为一个服务器。连接服务器与服务器的数据线则看做一条树边。两个服务器进行数据的交互时，数据会经过连接这两个服务器的路径上的所有服

原创 [HNOI2016] 树

P4055[Hnoi2016 day1]树时间限制 : - MS   空间限制 : 165536 KB 评测说明 : 2s问题描述小A想做一棵很大的树，但是他手上的材料有限，只好用点小技巧了。开始，小A只有一棵结点数为N的树，结点的编号为1,2,…,N，其中结点1为根；我们称这颗树为模板树。小A决定通过这棵

原创 [HNOI2016] 最小公倍数

P4053 [Hnoi2016 day1]最小公倍数时间限制 : - MS   空间限制 : 565536 KB 评测说明 : 4s问题描述给定一张N个顶点M条边的无向图(顶点编号为1,2,…,n)，每条边上带有权值。所有权值都可以分解成2^a*3^b的形式。现在有q个询问，每次询问给定四个参数u、v、a和

原创 4月题目汇总1(4.16---4.22)

从4月17日开始4.17 晚来源： BZOJ 2811: [Apio2012] Guard分类：**贪心 区间问题 简要题解： 1: 去掉没有忍者的区间 2：离散化剩下的点和剩下的区间（只离散化点就可以了） 3： 区间之间有包含关系的 ，去掉大区间，只考虑小区间。具体实现用一个队列就好。 4: 贪心原则：从左到右没有忍者的区间在右端点放一个 5

原创 CQOI2017流水账(滚粗记) ---NKwxh (1879570236) 本博客为两个Oier共有，你们不要晕啊！

时间过得真快，从noip2016转眼间就到了CQOI2017。 好吧，废话不多说，谈谈CQOI2017的经历以及一些心得。    --NK wxh  （1879570236）4月6日，我们打完了最后一场信心赛，255rank1，但我一点信心都没有！4月7日，上午复习了一下各种数学板子以及计算几何板子，心里面是慌的，其实并没有怎么复习进去。下午去重邮试机，得知今年用的

原创 CQOI2017酱油记(spark)

CQOI2017在樱花盛开的重庆邮电大学举行。 靠着NOIP的一点点微弱优势，勉强没有滚粗，但是就这两天的比赛而言，我的表现明显是比较失败的。 先说DAY1 ： 50 分 T1： 树形DP，一下过了样例，非常兴奋，只粗略看了一下DP 数组就没在管，其实后来有足够的时间对拍，但是抱着侥幸心里觉得T1应该没有大问题，于是果断不管，50分没有了。 T2： 计算几何。 不会正解的算法，但本应该拿到

原创 大组合数取模

考虑C(n,m)%P 情况一：n,m很大，P为素数 处理小范围的阶乘和阶乘的逆元 用卢卡斯定理即可。卢卡斯定理: 情况二： 当P= p1∗p2∗p3∗...∗pn% MathType!MTEF!2!1!+-% feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9w

原创 BZOJ 1500 维修序列 Splay

题目链接： (http://www.lydsy.com/JudgeOnline/problem.php?id=1500) 各种区间操作，但是方法都是统一的 对于操作区间[L,R] 先把(L-1)旋转到根节点，再把(R+1)旋转到根的右儿子，那么需要操作的区间[L,R]就是根的 右儿子的左儿子。 但是L==1和R==n都需要特判，非常麻烦。所以不妨在前后都加一个数字，变成N+2长度的序列，就

原创 作诗--莫队

【分块·莫队】作诗Time Limit:30000MS Memory Limit:131072K Case Time Limit:3000MSDescription神犇SJY虐完HEOI之后给傻×LYD出了一题： SHY是T国的公主，平时的一大爱好是作诗。 由于时间紧迫，SHY作完诗之后还要虐OI，于是SHY找来一篇长度为N的文章，阅读M次，每次只阅读其中连续的一段[l,r]，从这一段中

原创 小Z的袜子--莫队

【分块·莫队】小Z的袜子Time Limit:10000MS Memory Limit:524288K Case Time Limit:1000MSDescription作为一个生活散漫的人，小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天，小Z再也无法忍受这恼人的找袜子过程，于是他决定听天由命…… 具体来说，小Z把这N只袜子从1到N编号，然后从编号L到R(L 尽管小Z

原创 Miller_Robin素数判定和Pollard_rho质因数分解模板

#include<cstdio>#include<iostream>#include<cstdlib>#include<cmath>#include<cstring>#include<queue>#include<vector>#include<algorithm>#define ULL unsigned long long using namespace std;const int

原创 CQOI2014 排序机械臂--splay膜版题

【CQOI2014】排序机械臂Time Limit:20000MS Memory Limit:565536K Case Time Limit:2000MSDescription为了把工厂中 高低不等的物品按从低到高排好序，工程师发明了一种排序机械臂。它遵循一个简单的排序规则，第一次操作找到最低的物品位置P1，并把从左起第1个至第P1个之间的物品反序；第二次找到第二低的物品的位置P2，并把左起第二

原创 KD树 CQOI 2016 K 远点对

Description 已知平面内N个点的坐标，求欧氏距离下的第K远点对。分析： 维护一个存有当前2k远距离的堆（最开始全部是0）。建一个KD树，每个点查询一次，更新堆即可。对于KD树的每个节点，存储Max[],Min[],存下子树中每一维的最值，用来查询时剪枝。 代码：#include <cstdio>#include <iostream>#include <algorithm

原创 后缀自动机+DP NOI2015 Day 2 品酒大会

很显然可以用后缀自动机来搞。 将输入的字符串翻转，构造SAM。 对每一个节点x，求出： (1)子树中 满足LCA(u,v)==x 的点对 的对数 (2)子树中 满足LCA(u,v)==x 的点对 的美味值乘积最大值。 注意最大值有可能由两个最小的负数相乘得到，所以最大、最小值都要记。 最后的答案为ans1[],ans2[],如果节点x表示的最长子串长度为Max, 那么x的答案可

原创 后缀自动机+DP BZOJ 3238 差异

传送门：BZOJ 3238 差异 分析： 其实题目就是要求任意两个后缀T[i]和T[j] (i#include<cstdio>#include<iostream>#include<cstdlib>#include<cmath>#include<cstring>#define LL long longusing namespace std;const int maxn=1000000

原创 dfs序+分块求众数

铅导体问题描述 何老板要求第三题要很简单，最好是铅导体的难度。 于是，nodgd把N个铅块用N−1根导线相连，就形成了一个铅导体。只要是在这个基础上出题，就符合何老板的要求。nodgd为了方便，就把其中的一个铅块固定在了墙上，其他铅块在导线的作用下自然下垂。每个铅块有个固定的纯度，若干个相同纯度的铅块可以聚变发电，发电的电压与铅块的数量成正比。每当nodgd需要电疗的时候，就在铅导体上选一个铅块

原创 FFT多项式快速幂，对于x^num取模，顺便再模一个998244353

多项式快速幂 时间限制 : 60000 MS 空间限制 : 524288 KB问题描述：给一个n次多项式，求它的k次方。没关系，随手模一个998244353就行了。没关系，再随手模一个xm就行了。输入格式：第一行n，意义如上。 第二行n+1个数，a0,a1,…,an，分别是0,1,…,n次项系数。 第三行k，意义如上。 第四行m，意义如上。输出格式一行，b0,b1,…,bm-1，分别是0

原创 CDQ分治——NKOJ3210 围牛群

题意： 维护一个点集，支持两种操作： (1) 插入点P（x,y）. (2)询问所有点是否在 直线Ax+By=C的同一侧。 分析： 显然可以用平衡树维护凸包，每次二分查找切线。但是代码比较复杂。 考虑cdq分治，按照输入的顺序，左区间的点对右区间的询问都会有影响。所以直接求出左区间点构成的凸包，静态查询右边区间的直线是否切割凸包即可。 对于直线Ax+By+C=0，令f(x)=Ax+By+

原创 CDQ分治+ DP BZOJ 1492 Cash

题目链接：[NOI2007]货币兑换Cash分析请见CDQ论文： 从《Cash》谈一类分治算法的应用代码如下：#include#include#include#include#includeusing namespace std;const int maxn=200000+5;const double inf= 999999999999.00, eps= 1e-9;

原创 CDQ 分治 NKOJ 3655 菊花的故事1

问题描述何老板很喜欢菊花，所以现在他买了n朵菊花.每朵花有三个属性：花形(s)、颜色(c)、气味(m)，又三个整数表示。现要对每朵花评级，一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽，当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然，两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。输入格式第一行为N,K

原创 高斯消元（xor方程组）模板

#include#include#include#include#includeusing namespace std;const int inf=0x3f3f3f3f;template inline void _read(T& x){ char ch=getchar();bool sign=true; while(!isdigit(ch)){if(ch=='-')

原创 动态树不完整模板---hdu4010 Query on The Trees

#include#include#include#include#include#include#includeusing namespace std;template inline void _read(T& x){ char ch=getchar();bool sign=true; while(!isdigit(ch)){if(ch=='-')sign=fal