2016年04月_里阿奴摩西

原创 [数学] BZOJ 4544 椭圆上的整点

和 bzoj1041 圆上整点很类似只要把A改为3*a*a枚举就行了这种暴力233跑得太慢#include#include#include#includeusing namespace std;typedef long long ll;inline ll Jud(ll x){ ll isqr=sqrt(x); if (isqr*isq

2016-04-19 12:29:44 1044

原创 [Contest] Cqoi2016 题目汇总

把最近做的cqoi汇总一下按做的顺序排好了密钥破解：Rho大整数分解手机号码：数位DP路由表：Trie树单调栈伪光滑数：堆不同的最小割：voidK远点对：voidK远点对：void

2016-04-18 20:45:17 436

原创 [堆思路题] BZOJ 4524 [Cqoi2016]伪光滑数

可持久化可并堆+DP：http://blog.csdn.net/liuguangzhe1999/article/details/51132255我也不想这个有趣的方法就这么绝迹了呢但是我不会打堆：http://blog.csdn.net/dropd/article/details/51138254#include#include#include#incl

2016-04-18 20:37:11 598

原创 [Trie树单调栈] BZOJ 4523 [Cqoi2016]路由表

建一棵字典树，记一下时间然后就是在字典树上匹配因为匹配长度是单调增的，所以维护一个时间的单调递增栈#include#include#include#include#define cl(x) memset(x,0,sizeof(x))using namespace std;typedef long long ll;inline char nc(){ stati

2016-04-18 20:32:07 599

原创 [Trie树] BZOJ 1819 [JSOI]Word Query电子字典

Trie树上暴力复杂度不虚#include#include#include#define cl(x) memset(x,0,sizeof(x))using namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(

2016-04-18 18:53:59 602

原创 [暴力乱搞] BZOJ 2338 [HNOI2011]数矩形

枚举对角线，两条对角线合法必定中点重合长度相同暴力两辆计算竟然没T#include#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)

2016-04-18 18:51:06 331

原创 [数位DP] BZOJ 4521 [Cqoi2016]手机号码

题解：http://www.cnblogs.com/ccz181078/p/5379967.html数位DP和喜闻乐见的前缀和#include#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p

2016-04-18 18:49:31 506

原创 [Rho大整数分解] BZOJ 4522 [Cqoi2016]密钥破解

这种裸题就是Rho的直接应用长点记性啊，linux下RAND_MAX=2^31 狂T啊啊#include#include#include#include#include#includeusing namespace std;typedef long long ll;typedef pair abcd;inline char nc(){ static cha

2016-04-18 18:47:13 678

原创 [数学拉格朗日四平方和定理 Rho大整数分解] BZOJ 2904 平方和

题解:JudgeOnline/upload/201604/sol.docPS linux下RAND_MAX是2^31！！TLE了快一页另一题4522跟着一起T#include#include#include#include#include#includeusing namespace std;typedef long long ll;inline

2016-04-18 18:45:42 1696

原创 [CDQ分治凸包] BZOJ 2961 共点圆

%%% PoPoQQQ ：http://blog.csdn.net/PoPoQQQ/article/details/42318247“题目大意：给定平面，多次插入点和圆，每次插入点时询问当前插入的点是否在之前插入的所有圆中并且至少在一个圆中直接用数据结构维护这些点和圆不是很好写，我们考虑CDQ分治对于每层分治，我们需要对于[mid+1,r]中的每个点求出[l,mid

2016-04-18 18:41:44 468

原创 [单调栈] BZOJ 1345 [Baltic2007]序列问题Sequence

单调栈维护自己体会一下应该问题不大#include#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,1000

2016-04-18 18:38:04 409

原创 [单调栈] BZOJ 1113 [Poi2008]海报PLA & BZOJ 1628 [Usaco2007 Demo]City skyline

双倍经验就是把下面能合并的合并成一张单调栈维护1113#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin

2016-04-18 18:36:47 413

原创 [欧拉函数] BZOJ 2226 [Spoj 5971] LCMSum

显然ans=Σ[d|n]d*F(n/d)/d=nΣF(d)F(n)为n以内与n互质的数之和F(n)=n*phi(n)/2 由互质的对称性[ gcd(i,n)==1 gcd(n-i,n) ] 可知#include#include#include#includeusing namespace std;typedef long long ll;inline char

2016-04-18 18:33:20 347

原创 [数论 CRT] BZOJ 2976 [Poi2002]出圈游戏

CA爷把这道题搬到了BC上，据说因为题目描述有误炸了，于是就来做做就是一同余方程组用中国剩余定理解但是自己的CRT打得好low#include#include#include#include#define cl(x) memset(x,0,sizeof(x))using namespace std;typedef long long ll;typedef pa

2016-04-18 18:27:23 597

原创 [分块随机化] BZOJ 2223 [Coci 2009]PATULJCI

不会主席树分块加随机化卡时过但是另一道双倍经验就过不去啦233这个坑以后再填#include#include#include#includeusing namespace std; inline char nc(){ static char buf[100005],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=bu

2016-04-18 18:24:32 438

原创 [Contest] IOI2011 day1 题目汇总

跟着ZZY神犇刷ioi题garden不会做，只会做倍增69分以后再填坑#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,

2016-04-16 14:07:07 742

原创 [树上倍增] BZOJ 4281 [ONTAK2015]Związek Harcerstwa Bajtockiego

树上倍增裸题%%%用树链剖分拿rank3的Evan大爷#include#include#include#include#include#define V G[p].vusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p

2016-04-16 13:55:10 543

原创 [CDQ分治] BZOJ 2244 [SDOI2011]拦截导弹

%%% PoPoQQQ ：http://blog.csdn.net/popoqqq/article/details/44962711“第一问CDQ分治裸上第二问用每个元素所在的LIS个数/总LIS个数就是答案每个元素所在的LIS自己必选，然后统计前面的方案数和后面的方案数以前面的方案数为例，令f[x]为以x结尾的LIS长度，那么有DP方程：g[i

2016-04-16 13:27:06 551

原创 [Simpson积分] BZOJ 1502 [NOI2005]月下柠檬树

%%% PoPoQQQ ：http://blog.csdn.net/popoqqq/article/details/39252719题目大意：给定一棵由圆台和圆锥构成的柠檬树，月光以α的夹角平行射向地面，求阴影部分面积求公切线#include#include#include#include#define PI acos(-1.0)#define EPS (

2016-04-16 13:24:09 578

原创 [Simpson积分 || 圆的离散化几何] BZOJ 2178 圆的面积并

不知道正解为何物：http://mojijs.com/2015/09/206783/index.html直接套用Simpson 函数值是一些线段的并#include#include#include#include#define PI acos(-1.0)#define eps 1e-13using namespace std;typedef pair abcd;i

2016-04-16 13:20:41 846

原创 [模拟退火二分答案] BZOJ 1038 [ZJOI2008]瞭望塔

正解半平面交，%%% PoPoQQQ : http://blog.csdn.net/popoqqq/article/details/39340759"确定瞭望塔的高度的时候我们选择二分处理对于二分的每一个值我们把折线上的端点从左到右枚举瞭望塔的塔尖到每个端点的连线必须从左到右逆时针顺序否则就会被遮挡"#include#include#include#include#in

2016-04-16 13:17:12 656

原创 [二分答案单调队列] BZOJ 3316 JC loves Mkk

%%%PoPoQQQ：http://blog.csdn.net/popoqqq/article/details/43059527看到平均值最大果断二分答案看到长度[L,R]果断单调队列对数组维护一个前缀和，对前缀和维护单调递增的单调队列每扫过一个数sum[i]，将sum[i-L]加入单调队列，再把距离i超过R的点删掉长度为偶数？对奇数位置和偶数位

2016-04-13 17:31:02 767

原创 [组合数学思路题] BZOJ 4402 Claris的剑

好神的题：http://blog.csdn.net/lych_cys/article/details/50603389实际上，任意一种可行的方案的序列都可以转化为本质相同的下述序列：先是若干个(1,2)，然后是若干个(3,2)，然后是若干个(3,4)……以此类推，当然如果长度为奇数就在最后面加上一个数，举个例子： {1,2,3,2,3,

2016-04-13 17:28:10 466

原创 [DP] BZOJ 2958 序列染色

题解：http://blog.csdn.net/lych_cys/article/details/51078401令f[i][j][k]表示到第i位，状态为j，第i位为k的方案数。其中状态0表示之前没有k个B和k个W；1表示只有k个B；2表示之前有k个B和k个W。转移的时候只需要得到这一位填某一个字母的情况下会不会产生新的k个连续串，然后该加的地方加上，改去掉重复方案的去掉即可。#in

2016-04-13 17:25:48 758

原创 [SPFA] BZOJ 2622 [2012国家集训队测试]深入虎穴 & BZOJ 2617 [Ioi2011]crocodile

题解：http://blog.csdn.net/wzq_QwQ/article/details/47997617其实总体东西就是个SPFA，不过这里面我们需要做一些改动。这只老虎是DP老虎，他每次会占据你的最优方案。所以使得你只能用次优方案来更新最优方案以及次优方案。这就是SPFA的路子。但是这个路子并不够。让我们考虑一个点只有一个儿子的时候。他可能会扫一圈后第一次

2016-04-13 10:49:23 1008

原创 [倍增堆] BZOJ 4458 GTY的OJ

题解：JudgeOnline/upload/201604/Solution-4458.rar By jinlifu1999超级钢琴...都已经成经典了本来想打树链剖分和ST表的，后来看题解发现倍增更好打，膜打树链的Evan题目后跟个题解链接是SMG#include#include#include#includeusing namespace std;

2016-04-13 10:16:20 818

原创 [线段树] BZOJ 4499 线性函数

好裸的线段树#include#include#include#define P 1000000007using namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf

2016-04-12 19:59:15 624

原创 [DP 组合数学] BZOJ 4498 魔法的碰撞

传送门：http://blog.csdn.net/visit_world/article/details/51090964dp其中一维解释一下，表示可以放且必须放的一个位置#include#include#include#include#define P 1000000007using namespace std;typedef long long ll;inlin

2016-04-12 19:57:31 768

原创 [树链剖分线段树] BZOJ 3531 [Sdoi2014]旅行

树链剖分给每一个信仰开一棵线段树然后就是动态开点的打码问题了#include#include#include#include#define V G[p].vusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fr

2016-04-12 19:54:43 364

原创 [分类讨论思路题] BZOJ 3135 [Baltic2013] pipes

官方题解：http://boi2013.informatik-olympiade.de/wp-content/uploads/2013/05/pipes-spoiler.pdf输入构成一个n个m元方程组当m=n-1所以只需考虑 m==n-1 m==nm==n-1 这是一棵树从叶子往上推就好了在根节点判一下是否矛盾m==n这是一个环套树先处理到只剩一个

2016-04-12 14:00:39 512

原创 [Contest] BC Round #79

可以打div1了，感谢带我飞的神犇ZZY不熟悉莫比乌斯，T4没打出来，一直在给暴力调常数A懒得推了，乱搞，暴力枚举好了#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2

2016-04-11 17:08:25 494

原创 [数学] BZOJ 4459 [Jsoi2013]丢番图

就是个类似求n^2因子个数的东西#include #include using namespace std;typedef long long ll; ll n,ans; int main(){ freopen("t.in","r",stdin); freopen("t.out","w",stdout); ll ans=1; scanf(

2016-04-07 20:21:41 669

原创 [DP LIS] BZOJ 4282 慎二的随机数列

显然随机的数全部属于上升序列肯定能构造出最优解，那么直接认为随机的数都在最优解中。如果是求最长不上升子序列的话就直接去掉所有的随机数求一遍然后加上随机数的个数。现在是让求最长上升子序列，只需把每个数都减去前面随机的数的个数然后求一遍再加上随机数的个数即为答案！#include#include#includeusing namespace std;inline char n

2016-04-07 19:35:05 603

原创 [生成函数 FFT] BZOJ 3771 Triple

一个生成函数，用FFT做乘法重复的部分我们考虑用容斥原理来解决。为了方便描述我们不妨设三个多项式。第一个是仅取一个而构成的多项式。->x第二个是仅取相同的两个而构成的多项式。->y第三个是仅取相同的三个而构成的多项式。->z对于本题有三种情况。第一种是取一个，显然直接将x加到答案就好。第二种是取两个，则需要一小步容斥，即(x*x-y)/2第三种是取三个，则需

2016-04-07 18:50:03 642

原创 [乱搞] BZOJ 2296 [POJ Challenge] 随机种子

乱搞题既然要0~9 那么就摆在前面，后面六位补齐使得是倍数就行了10+6 刚好等于 16#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(

2016-04-07 18:46:58 658

原创 [贪心] BZOJ 3671 [Noi2014]随机数生成器

这题因为要求是字典序最小所以就可以贪心。一位位取看能不能取到。取的时候把取的这个数左下右上标记成不可取就行了。均摊复杂度O（NM）类似筛法的复杂度#include#include#include#include#include#define cl(x) memset(x,0,sizeof(x))using namespace std;ty

2016-04-07 18:44:44 325

原创 [线段树] BZOJ 4491 我也不知道题目名字是什么

用个线段树维护下#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF

2016-04-07 18:41:46 820

原创 [LCT 线性模方程] BZOJ 2759 一个动态树好题 & 4266 小强的动态方程

%%%PoPoQQQ：http://blog.csdn.net/popoqqq/article/details/40436165“题目大意：给定n个形如xi=ki*x_pi+bi mod p的同余方程组支持修改操作和求解操作确实好题感谢此题作者顺便吐槽一下作者的Splay不加空节点太蛋疼了0.0将每个点i的父亲设为pi 我们将会得到一座基环树林将环上的一条边

2016-04-07 18:40:35 647

原创 [莫队算法线段树斐波那契暴力] Codeforces 633H Fibonacci-ish II

题目大意：给出一个长度为n的数列a。对于一个询问lj和rj，将a[lj]到a[rj]从小到大排序后并去重，设得到的新数列为b，长度为k，求F1*b1+F2*b2+F3*b3+...+Fk*bk，其中F为斐波那契数列，F1=F2=1。对每个询问输出答案模m。区间查询离线用莫队算法开棵权值线段树，然后用斐波那契的性质updateF(n+m)=F(n+1)*F(m)+F(n)*F(

2016-04-05 19:38:01 1705

原创 [动态树 LCT] BZOJ 2157 旅游

动态树裸题把边建成一个点就好了取反，打个标记#include#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin);

2016-04-04 13:24:12 627

空空如也

空空如也