2015年12月_Hillan_

原创 BZOJ1900: Uva10601 Cubes

burnside的定理。。。然后看到1s的时限整个人都炸了！！！！好吧那我就打表好了23333(1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(2 3 4 1) (6 7 8 5) (10 11 12 9)(3 1) (4 2) (7 5) (8 6) (11 9) (12 10)(4 3 2 1) (8 7 6 5)

2015-12-27 21:13:29 654

原创 Miller_Rabin 大素数判断

满足费马定理公式的非素数不怎么多就是了。。。然后这就是第一次筛选然后就根据二次探测定理二次探测定理如果p是一个素数，0然后就有一般步骤0、先计算出m、j，使得n-1=m*2^j，其中m是正奇数，j是非负整数1、随机取一个b，22、计算v=b^m mod n3、如果v==1，通过测试，返回4、令i=15、如果v=n-1，通过测试，返回

2015-12-24 20:33:25 446

原创 PAM模板

对着别人程序打的。。。好歹打出来了是吧#include#include#includeusing namespace std;struct Node{ int len; Node *ch[31]; Node *last; int num; Node(){for(len=0;len<=30;len++)ch[len]=NULL;last=NULL;len=num

2015-12-24 17:00:34 522

原创回文自动机 PAM

http://blog.csdn.net/lwfcgz/article/details/48739051这篇文章翻译地很详细了感觉现在概念差不多了晚上开始打代码PAM和SAM感觉差不多？ 0.0然后貌似有很神奇的用处？（虽然我现在SAM学的都不怎么扎实。。。

2015-12-24 14:17:37 774

原创 BZOJ2301: [HAOI2011]Problem b 莫比乌斯反演

同上一题然后就是拆操作了注意一下abcd除以k后区间的是闭还是开的#include #include #include #include #include #include #define LL long longusing namespace std;char c;inline void read(int&a){ a=0;do c=g

2015-12-23 21:13:46 534

原创 BZOJ1101: [POI2007]Zap 莫比乌斯反演

很裸的吧转换一下就好了然后统计相同a/d,b/d 的值#include #include #include #include #include #include #define LL long longusing namespace std;char c;inline void read(int&a){ a=0;do c=getchar();whil

2015-12-23 21:12:06 623

原创莫比乌斯反演

首先莫比乌斯函数有个性质∑d|nμ(d)={1 (n=1)0 (n>1)证明： ①n=1时，不做多余说明。。。 ②n>1,根据唯一分解定理，可以分解n=∏ki=1paii 对于那些含平方因子也就是存在ai不为1的数，它的函数值为0，对答案没有任何贡献。所以我们来看看那些是互异素数乘积的数，每一个成为它约数的数是什么样的情况。（1）若d中有0个质因子(d=1

2015-12-23 19:56:58 539

原创 BZOJ1803: Spoj1487 Query on a tree III

裸的主席树+DFS序然后我差错查了好长时间发现时查询的时候k没有减应该到右子树结果到了左子树。。。#include#include#include#include#include#include#includeusing namespace std; char c;inline void read(int&a){ a=0;do c=getchar();

2015-12-22 21:09:25 715

原创 12.22.2015

突然发现这个月学了好多数据结构。。。一坨一坨的指针调试的时候要疯了周末的时候都是一天几十k的代码量希望自己能够快点强大起来吧！毕竟生在zj 就要付出更大的努力！！

2015-12-22 17:48:55 413

原创 BZOJ3674: 可持久化并查集加强版&&BZOJ3673: 可持久化并查集 by zky

妈个鸡。。。。没加强的没有要求在线。。。我日了狗了查了半天错。。。。反正就是用可持久化线段树维护一个可持久化数组具体看代码/************************************************************** Problem: 3674 User: liutian Language: C++

2015-12-22 17:43:02 556

原创 BZOJ2738: 矩阵乘法

整体二分然后在二维树状数组上更新#include#include#include#include#include#include#includeusing namespace std;char c;inline void read(int &a){ a=0;do c=getchar();while(c'9'); while(c='0')a=(a<<3)+(a

2015-12-20 20:44:29 541

原创 BZOJ1975: [Sdoi2010]魔法猪学院

第一道A*。。。今天才会做k短路。。。。#include#include#include#include#includeusing namespace std;struct Node{ int x; double dis; inline friend bool operatorb.dis; }};struct node{ double f,g; int

2015-12-20 19:34:45 563

原创 BZOJ2132: 圈地计划

又是不得了的建图姿势。。。。#include#include#include#define T 10001const int INF=1<<29;using namespace std;int n,m,cnt=1,tot,ans,head[10001],h[10002];int a[101][101],b[101][101],c[101][101];int xx[4]={

2015-12-20 18:06:36 521

原创 BZOJ1078: [SCOI2008]斜堆

http://www.cppblog.com/MatoNo1/archive/2013/03/03/192131.html可怕的结论%%%#include#include#includeusing namespace std;char c;inline void read(int &a){a=0;do c=getchar();while(c'9');while(c='0')a=

2015-12-20 16:14:05 539

原创 BZOJ1978: [BeiJing2010]取数游戏 game

反正就是DP然后用因子更新减少复杂度#include#includeusing namespace std;inline int max(int a,int b){return a>b?a:b;}int pre[10000005],f[50005],n,l;char c;inline void read(int &a){ a=0;do c=getchar();while(c

2015-12-20 15:25:11 621

原创 BZOJ2223: [Coci 2009]PATULJCI&&BZOJ3524: [Poi2014]Couriers

自己的智商又不够用了。。。。。一直错查了一天发现时New_Node的时候没有打return。。。。。这个错误之前做的题目已经出现过很多次了感觉自己瞬间拉低整条街智商很裸的静态区间k大值用主席树2223：#include#include#include#include#includeusing namespace std;setD;char c;boo

2015-12-20 14:52:13 520

原创 BZOJ1492: [NOI2007]货币兑换Cash

CDQ的论文讲的很清楚了其实这一道题写成斜率优化的形式之后我们会发现一个很严重的问题斜率不是单调的然后能用斜率优化么？可以。。。用CDQ将前一半的决策点处理一边后一半就能线性更新了然后总的复杂度是CDQ分治的复杂度%%%CDQ代码的话是直接copy的昨天打代码手打残了不想打了#include#include#include#include#incl

2015-12-20 12:58:45 940

原创主席树（可持久化线段树）学习笔记

首先离散化然后当前版本的线段树可以有上个版本的线段树修改lgn个节点得到的还满足区间减法好神#include#include#include#include#includeusing namespace std;setD;char c;bool flag;inline void read(int&a){a=0,flag=false;do c=getchar

2015-12-20 09:45:10 996

原创 BZOJ2527: [Poi2011]Meteors

这个。。。一开始用的是longlong 然后改成int就wa了。。。。时间垫底。。。。。可怕全局分治然后用线段树维护的时候直接永久化标记不用下传然后这一题和上一道树套树一样，又是因为自己傻逼少了一倍的线段树节点然后一直OLE不知道怎么了。。。人傻没办法#include#include#include#include#includeusing

2015-12-19 19:25:19 896

原创 BZOJ1146: [CTSC2008]网络管理Network

。。。我居然调试了一个下午　　因为一个else　一个if。。。。果然自己太弱听Claris说这一题可以转化为DFS序然后容斥可以直接变成logn*logn*n 只能膜拜我的是logn*logn*logn*logn*n然后Jrmz说，这还不如暴力。。。。。反正就是很普通的一道树剖然后在链上找k大值用树套树由于只会外层线段树维护区间的写法于是就打了3个log的#in

2015-12-19 15:43:47 2233

原创 BZOJ3262: 陌上花开 CDQ

反正就是裸题。。。。先是去重然后按x排序降一维在根据y来插排最后z在树状数组上查询修改就好了#include#include#include#includeusing namespace std;char c;inline void read(int &a){ a=0;do c=getchar();while(c'9'); while

2015-12-17 20:50:54 1119

原创 BZOJ1176: [Balkan2007]Mokia CDQ分治

很裸的一道CDQ分治吧拆操作询问区间的时候拆成询问四个子区间然后加加减减根据时间顺序来二分然后用插排降一维最后一维树状数组维护前缀和#include#include#include#includeusing namespace std;char c;inline void read(int &a){ a=0;do c=getchar();while(

2015-12-15 19:44:55 409

原创 BZOJ1047: [HAOI2007]理想的正方形

恩数据水被我卡过去了2333333我的是n^2logn的。。每一次扫一遍multiset。。然后居然就真的过去了23333#include#include#include #includeusing namespace std;multiset Sum,Par[1001];char c;bool flag;inline void read(int &a){ fla

2015-12-15 17:05:55 459

原创 BZOJ1053: [HAOI2007]反素数ant

最多不会超过13个不同的素因数然后就每一个素因数去搜能乘几次 2^2*3^3的因数是和2^3*3^2一样的所以需要前一个素数的指数比后一个素数指数大每次更新最优值#include#include#includeusing namespace std;#define ll long longint prime[101]={1,2,3,5,7,11,13,17,

2015-12-15 16:24:05 381

原创 FFT bzoj2179: FFT快速傅立叶

一个晚上和FFT相爱相杀的故事。。。。。。反正就是裸的。。。。然后这个是递归式的自己还是太弱了。。。。话说发现自己的代码很长诶。。。。。#include#include#include#includeusing namespace std;const double PI=acos(-1);struct Complex{ double Real,fake; Com

2015-12-14 21:10:53 449

原创 bzoj1324: Exca王者之剑&&bzoj1475: 方格取数

双倍经验。。都是最大流学到一个很了不得的建图姿势呢代码就是单纯的dinic我就不贴了既然没有在我的模板库里找到dinic 我还得手打。。。。。可怕

2015-12-13 13:56:22 516

原创【bzoj1023】[SHOI2008]cactus仙人掌图

http://z55250825.blog.163.com/blog/static/150230809201412793151890/%%%%反正自己就是sb。。第一道仙人掌dp。。。感觉很神奇啊#include#include#includeusing namespace std;char c;inline void read(int &a){a=0;do c=ge

2015-12-13 10:13:26 558

原创 Bzoj 3409: [Usaco2009 Oct]Barn Echoes 牛棚回声

SAM练手题反正就是构造两个SAM然后匹配就好了。。。。#include居然会错..要#include才行。。。#include#include#includeusing namespace std;struct Node{ Node *ch[200]; Node *last; int len; bool end; Node(){for(len=0;len<

2015-12-10 20:45:40 709

原创 BZOJ3196: Tyvj 1730 二逼平衡树

#include#include#includeusing namespace std;inline int min(int a,int b){return (a<b?a:b);}inline int max(int a,int b){return (a>b?a:b);} inline int abs(int a){return (a<0?-a:a);}char c;boo

2015-12-10 16:55:34 620

原创 BZOJ 1500: [NOI2005]维修数列

我的AC率啊 QAQ一开始看到陈老师300行左右的代码我就淡淡的说了句怎么打的这么多然后我就打了450行。。。。。调试了三天重写了两边果然是毒瘤Splay。。。。提取区间打标记下传标记就好了么。。还是蛮裸的果然还是自己太弱了#include#include#include#includeusing namespace std;#define ll int#d

2015-12-08 17:52:55 530

原创组合数学

http://wenku.baidu.com/link?url=rZhh5B-ZAdNES26Q_qP25wqxnAuqXRer1IJetlHzhuZvzHh7ab5xrAq5Lj8tfuF_8P8GbGwUp5D3r7wLf5ezhr2BgKcpeBA37Q3y1ha9AlG&qq-pf-to=pcqq.c2c很清楚

2015-12-01 21:10:34 561

Don't fake it till you made it