2016年04月_里阿奴摩西

原创 [可并堆] BZOJ 2809 [Apio2012]dispatching

考虑倒着想每次弹出最大值知道费用#include#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

2016-04-30 21:10:53 307

原创 [可并堆] BZOJ 2333 [SCOI2011]棘手的操作

诶，脑子不好使了一道裸题调了一晚上...#include#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

2016-04-30 20:42:02 478

原创 [可并堆] BZOJ 1078 [SCOI2008]斜堆

Orz mato大神的题解：http://www.cppblog.com/MatoNo1/archive/2012/10/07/192131.html实际上这道题并没有要写可并堆，只是考察了斜堆的性质斜堆是可合并堆的一种实现形式，其更稳定的实现是左偏树（斜堆只能做到均摊logN，而左偏树则可以严格做到每次操作O(logN)）。斜堆最典型的特点，上面已经说过了，如果一

2016-04-30 16:46:16 393

原创 [Random Heap 并查集] BZOJ 1455 罗马游戏

裸的可并堆 Random Heap 还是很兹瓷啊#include#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)+fre

2016-04-30 16:07:31 377

原创 [LCT 线段树 dfs序] BZOJ 3779 重组病毒

%%%PoPoQQQ http://blog.csdn.net/popoqqq/article/details/45919545”题目大意：给定一棵树，初始每个点都有一个颜色，支持三种操作： 1.将某个点到根的路径上所有点染上一种新的颜色 2.将某个点到根的路径上所有点染上一种新的颜色，然后把根设为这个点 3.定义一个点的代价为这个点到根路径上颜色的种类数，求某个点子树中所有点

2016-04-30 14:05:07 703

原创 [三进制倍增 || 可并堆] BZOJ 4003 [JLOI2015]城池攻占

蒟蒻想起来自己还不会可并堆诶出题人卡倍增，那么淡定的写了个三进制倍增#include#include#include#include#include#include#include#define V G[p].vusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,

2016-04-30 14:02:10 685

原创 [乱搞暴力并查集] BZOJ 3563 DZY Loves Chinese

好神奇的一道题不多说最后一次询问暴力并查集#include#include#includeusing namespace std;inline int read(int &x){ static char c=getchar(); for (;!(c>='0' && c<='9');c=getchar()) if (c=='\n' || c==EOF) return

2016-04-30 13:57:20 652

原创 [高斯消元物理] BZOJ 2419 电阻

首先我们设电流为1A 终点电势为零点i的电势为Ui由于电流是流显然对于每个点(点1和点n除外) 有总流入等于总流出即Σ(Ui-Uj)/Rij=0 (i!=1,i!=n)Σ(U1-Uj)/R1j=1Σ(Un-Uj)/Rnj=-1Un=0联立方程组高斯消元即可最后输出点1的电势就是答案#include#include#includ

2016-04-29 07:13:21 555

原创 [dfs序树状数组] BZOJ 4034 [HAOI2015]T2

两个log可以树链剖分day2听课时惊闻一个log的做法但我打low了，并不快...#include#include#include#define V G[p].vusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,

2016-04-27 19:52:30 545

原创 [SA] BZOJ 2251 [2010Beijing Wc]外星联络

见后缀数组后暴力就好了考前练版子，竟然WA了几发#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,std

2016-04-27 19:45:02 435

原创 [线段树期望] BZOJ 2752: [HAOI2012]高速公路(road)

分母很显然分子的话展开后发现只要维护vi，vi*i， vi*i*i的和就好了#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=

2016-04-27 19:42:52 537

原创 [LCT 组合期望] BZOJ 3091 城市旅行

%%%PoPoQQQ：http://blog.csdn.net/popoqqq/article/details/40823659第四个操作怎么办？！做过2752就好了#include#include#include#define V G[p].vusing namespace std;typedef long long ll;inline char nc(){ s

2016-04-27 19:40:52 425

原创 [Contest] Hnoi2016 题目汇总

把最近做的hnoi汇总一下按做的顺序排好了大数：莫队————————————线割分是我 >w———————————————– 序列：莫队单调栈— ———————————线割分是我 >w网络：树链剖分线段树堆————————————线割分是我 >w最小公倍数：分块并查集————————————线割分是我 >w

2016-04-23 10:29:19 538

原创 [Contest] Sdoi2016 题目汇总

把最近做的sdoi汇总一下按做的顺序排好了排列计数：组合错排————————————线割分是我 >w———————————————– 游戏：树链剖分线段树标记永久化— ———————————线割分是我 >w生成魔咒：后缀数组————————————线割分是我 >w征途：斜率优化DP————————————线割分

2016-04-23 10:18:16 570

原创 [二分图费用流] BZOJ 4514 [Sdoi2016]数字配对

不会有奇环，那么就建二分图然后跑最大费用流，直到费用小于零为止，处理下零头#include#include#include#include#define cl(x) memset(x,0,sizeof(x))#define oo 1LL<<60#define V G[p].vusing namespace std;typedef long long ll;inl

2016-04-23 10:06:35 591

原创 [斜率优化 DP] BZOJ 4518 [Sdoi2016]征途

就是求平方和的最小值然后就是裸的斜率优化了#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,

2016-04-23 10:02:03 354

原创 [进制贪心] BZOJ 1110 [POI2007]砝码Odw

因为都是倍数，可以考虑以最小的数为进制然后从小到大装，不够的向高位借，类似小学减法，但感觉可卡？#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(b

2016-04-23 09:59:37 412

原创 [找规律递归子问题 || 数位DP] BZOJ 4513 [Sdoi2016]储能表

看不懂数位ＤＰ，蒟蒻就是蒟蒻：http://fancypei.github.io/2016/04/16/SDOI2016%20Round1/然后ＺＺＹ　Ｇｏｏｇｌｅ到了一个神奇的做法，很ｈａｐｐｙ：https://blog.menci.moe/sdoi2016-table/强力膜拜考虑异或的性质：性质一：对于任意 x性质二：对于任意 x≠yx≠y，必

2016-04-23 09:51:49 1302

原创 [后缀数组 set] BZOJ 4516 [Sdoi2016]生成魔咒

将串反过来后缀数组然后用set维护一个序列插入时修改答案即可#include#include#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p

2016-04-22 10:29:35 487

原创 [树链剖分线段树标记永久化] BZOJ 4515 [Sdoi2016]游戏

想到了维护 ax+b 但不知道怎么修改然后就听说了标记永久化对于两个直线方程把决定范围小的下放到子树中树链剖分两个log，标记一次最多下推log层三个log复杂度好虚#include#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],

2016-04-22 10:26:43 835

原创 [分块并查集] BZOJ 4537 [Hnoi2016]最小公倍数

题解：http://www.cnblogs.com/y7070/p/5412793.html暴力分块”对于每一组，将符合一组a的询问选出来，将这些询问和这一块之前的边(a一定小于这些询问)按b排序，然后交替插入，询问，对于一个询问，在当前块也有可能有满足的边，我们将其加入，考虑后并撤销“#include#include#include#include

2016-04-22 10:23:24 614

原创 [二分答案树链剖分] BZOJ 4326 NOIP2015 运输计划

考虑二分把所有超过mid的链取交集，然后取交集上最长的边权赋0，判断#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)+frea

2016-04-22 10:19:45 624

原创 [概率树形DP] BZOJ 3566 [SHOI2014]概率充电器

题解：http://www.cnblogs.com/ctlchild/p/5105390.html考虑补集转化。先dfs一遍，令f[u]=(1-p[u])*∏(1-(1-f[v])*w) f[u]表示u这个点通过其子树并不能联通的概率。然后考虑v从其父亲连过来的情况，设x=1-f[u]/(1-(1-f[v])*w)表示除掉v，u联通的概率，那么f[v]*=(1-x*w)然后加起

2016-04-22 10:16:38 433

原创 [概率DP] BZOJ 3036 绿豆蛙的归宿

好水的题#include#include#include#include#define V G[p].vusing namespace std;typedef pair abcd;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(bu

2016-04-22 10:14:45 432

原创 [生成函数 FFT 分块] BZOJ 3509 [CodeChef] COUNTARI

学着大神优化了FFT常数：http://blog.csdn.net/geotcbrl/article/details/50636401可以看出来这是个生成函数暴力的话是N次FFT 然后考虑分块 B次FFT + 块内统计#include#include#include#include#include#define cl(x) memset(x,0,sizeo

2016-04-21 13:55:29 773

原创 [三分贪心] BZOJ 3874 [Ahoi2014]宅男计划

三分答案然后贪心单峰的话，可以发现，如果我们外卖的次数过少，那么就会出现一些食品性价比不高的情况；如果次数过多，那么就会浪费外卖运费。#include#include#includeusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p

2016-04-21 13:35:26 621

原创 [树链剖分线段树堆] BZOJ 4538 [Hnoi2016]网络

裸的树链剖分线段树里套个堆最暴力的打法还有一种做法是考虑二分答案，如果大于答案的边都经过询问点，则往小的二分否则往大的二分。这样我们只要维护路径的交集。用rmq求lca的话可以做到O(1)求出两条路径的交集。可以离散后用线段树维护二分结构，并在对应的节点上记录路径交集O(logn)内回答询问了。

2016-04-21 13:32:40 489

原创 [数学找规律] BZOJ 3823 定情信物

%%% zky http://blog.csdn.net/iamzky/article/details/42211705加强数据，被卡出屎，不知道WA了几发对比下列算式：(x+2)^0=1(x+2)^1=x+2x+2)^2=x²+4x+4(x+2)^3=x^3+6x^2+12x+8可以归纳出：一个n维立方形(n-cube)所包含的k维元素个数等于(x+2)^n展开式的k

2016-04-21 13:18:33 553

原创 [组合数学卡特兰数] BZOJ 1856 [Scoi2010]字符串

终于弄懂了Catalan数的另类通项的实际意义了：http://www.cnblogs.com/jianglangcaijin/p/3443689.html#include#include#include#include#define P 20100403using namespace std;typedef long long ll;const int N=20000

2016-04-20 19:32:42 448

原创 [莫队单调栈] BZOJ 4540 [Hnoi2016]序列

题解：http://www.cnblogs.com/Ngshily/p/5409337.html给出一个序列，求某段区间的位置不同子串的最小值之和对于前4040，我们可以n2n2预处然后O(nn√)O(nn)预处理中zzz[i][j]zzz[i][j]表示以ii开头，结尾最长到jj的前缀子串的最小值之和同理fff[i][j]fff[i][j]表示倒过来的情况

2016-04-20 16:23:09 916

原创 [莫队] BZOJ 4542 [Hnoi2016]大数

题解：http://www.cnblogs.com/Ngshily/p/5409409.html如何判断一个炒鸡大的数nn能不能被另一个数PP整除，，，我们有如下结论若xmodP=axmodP=a，且(n∗10k+x)modP=a(n∗10k+x)modP=a，kk为xx的长度，gcd(P,10k)=1gcd(P,10k)=1那么nmodP=0nmodP=0胡乱证明分割

2016-04-20 16:21:10 457

原创 [组合容斥 DP] BZOJ 2863 愤怒的元首

题解：http://blog.csdn.net/werkeytom_ftd/article/details/51182412我们设f[i]表示i个带编号结点组成的合法DAG的个数。如果删去DAG中所有入度为0的结点，接下来的DAG也是一个合法DAG。但我们不知道有多少入度为0的结点，所以我们想要得到至少有0个入度为0的结点的DAG个数（对就是这么神奇）。那就容斥呗。

2016-04-20 16:18:17 637

原创 [数学二项式定理快速幂] BZOJ 3157 && BZOJ 3516 && BZOJ 4126 国王奇遇记

膜拜神犇n+e：http://trinkle.is-programmer.com/2015/6/30/the-adventures-of-the-king.100058.html本题按时间复杂度的不同共有三种解法。只会O(m2)的方法#include#include#include#define P 1000000007using namespace s

2016-04-20 16:13:03 805

原创 [二分答案+对偶图 || 三角剖分] BZOJ 3007 拯救小云公主 && BZOJ 4219 跑得比谁都快

3007是小数据可二分答案：http://blog.csdn.net/popoqqq/article/details/442248154219是大数据三角剖分平面欧几里德最小生成树：http://trinkle.is-programmer.com/2015/7/13/bzoj-3007.102710.html3007的代码 4219不会打... 先挖个坑#incl

2016-04-19 19:50:40 1268

原创 [树链剖分线段树] BZOJ 2908 又是nand

树剖然后开32棵线段树记录1/0 从左/右经过变成的值然后就是一些区间合并成一条有向路径#include#include#include#define V G[p].v#define nand(x,y) (!((x)&(y)))#define digit(x,k) (((x)>>((k)-1))&1)using namespace std;typedef lo

2016-04-19 18:41:43 617

原创 [乱搞暴力] BZOJ 1142 [POI2009]Tab

题解：http://blog.csdn.net/vmurder/article/details/45688219无论如何变换，该在一行的还是会在一行，该在一列的还是会在一列。#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=bu

2016-04-19 18:39:20 503

原创 [组合错排] BZOJ 4517 [Sdoi2016]排列计数

n个里选m个然后n-m个里错排***错排相关公式Dn=(n-1)*(Dn-1+Dn-2) D1=0,D2=1Dn=n*Dn-1+(-1)^(n-2)=n*Dn-1+(-1)^nDn=n!(1-1/1!+1/2!-1/3!+……±1/n!)#include#include#include#define P 1000000007using namespace

2016-04-19 18:36:59 407

原创 [分治最小割] BZOJ 4519 [Cqoi2016]不同的最小割

最多有N-1种zuixiaoge分治最小割裸题#include#include#include#include#include#define V G[p].v#define cl(x) memset(x,0,sizeof(x))using namespace std;namespace hash_map{ struct node{ int x; int ne

2016-04-19 18:32:43 538

原创 [并查集思路题] BZOJ 2503 相框

%%% PoPoQQQ http://blog.csdn.net/popoqqq/article/details/48031135#include#include#include#includeusing namespace std;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (p1==p

2016-04-19 18:30:49 529

原创 [线段树标记永久化单调队列] BZOJ 1171 大sz的游戏 && BZOJ 2892 强袭作战

很好的题解：http://blog.sina.com.cn/s/blog_76f6777d0101dizp.html“考虑裸的n^2暴力dp. dp[i]=min(dp[j])+1 ( d[i]-d[j]由于第二个相对较难处理，首先考虑它. 把x,y离散化后搞个线段树，支持插入删除一个区间的答案，并且询问区间的答案即可.关于如何维护区间里的答案，虽然答案是

2016-04-19 12:47:00 1026

USACO全部译题

空空如也