- 博客(109)
- 收藏
- 关注
RSS订阅原创 表达式求值
#include#include#include#include#include#includeusing namespace std;templateclass expCal{ public: expCal(){ expMid.clear();expSuffix.clear();} expCal(string s){ expCal();
2013-10-06 18:10:55
751
原创 hdu 4686 （等比解法）
这题的正解是矩阵快速幂，但是在比赛的时候没有构造出矩阵，然后直接用数学的方法把这题ac了，方法如下：可以得出an和bn都是类等比数列，类等比数列暂且定义为等比数列之间的加减乘除组合后的数列，an,bn如下：an=A0*Ax^n+Ay*(Ax^(n-1)+Ax^(n-2)+......+Ax^1+Ax^0).(n>=1)bn=B0*Bx^n+By*(Bx^(n-1)+Bx^(n-
2013-08-22 12:53:38
812
原创 hdu 4680(splay tree)
解题思路： 第一种操作，用小集合暴力插入到大集合，可以证明总的时间复杂度O(n*logn); 第二种和第三种操作，每次操作O(logn) 第四种操作，因为任意3个数不能组成一个三角型，则就是长边大于剩下的两边之和，和斐波那契数列类似， 所以对于10^9的数据范围来说只有几十个，用贪心的思路，从最小的开始，暴力二分查找即可。 第五种操作，因为一个区间的gcd值
2013-08-17 19:29:24
1006
原创 hdu 4667（几何）
解题思路：求出所有可能在外围的点。也就是将两圆两两求外公切线，得到所有的切点；将圆和三角形的三个点依次做切线，也得到切点，再加上所有三角形的三个点。最后对这些点求凸包，对于凸包上的每条边，如果它们在同一个圆上，用相应的圆弧替代。#include#include#include#include#include#include#include#include#inclu
2013-08-15 13:15:39
1059
2
原创 hdu 3547(Polya计数)
//根据Polya计数的原理，我们先要求出每个置换群中的循环的个数//这里有24个置换群，为什么是24个了，很容易想到，把每一个面当作正面并按一个方向（顺时针或逆时针）旋转4次，总共6个面，所以就是24个置换群了；//我们可以先对立方体的8个顶点编号，然后就可以暴力预处理出24置换群的循环个数了，最后就是用Polya的计数公式计数了,代码如下：#include#include
2013-08-13 13:20:34
787
原创 自己实现用FFT加速多项式计算
/* * FFT之所以能加速DFT的计算，得益于n次单位复数根的几个性质： * 1.消去引理 w(d*n,d*k)=w(n,k); * 1.折半引理 w(n,k)^2=w(n/2,k)； * 2.求和引理，即：sum{w(n,k)^j}=0 (0 * 详见：算法导论第30章，p527～p535*/#include#include#include#include
2013-08-12 14:45:42
1613
原创 hdu 4622(后缀自动机｜后缀数组)
//昨天晚上一直用后缀数组来写，怎么写都超时，今天看了多校3的官方解题，然来是卡常数；//今天上午就去学习了缀自动机（以前没有主动去学后缀自动机，以为后缀自动机能解的，后缀数组都能解，这次尝到了苦头）；//这题官方解题说可以O(n^2)+O(1),可是我写的后缀自动机时间复杂度还是O(n*q),因为每次计算字串的个数的还是把所有的step遍历了的一遍，但是没事重新建树，建树的常数小了一些；
2013-07-31 16:37:03
1876
4
转载 zoj 题型分类
简单题#1001 -____-b A+B。#1110 Dick and Jane 胡乱枚举收缩一下情况就可以了。#1115 a[i+1] = a[i] 的全部数位上的加起来，直到剩下一个，直接模拟。#1414 太弱太弱，按照模 4 分类讨论一下即可。#1713 简单的字符串截取和字符计数。#1716 简单的二维数组区间求和，作累加，然后容斥一下；预处理 O(W*H)，查询
2013-07-26 10:17:58
2141
原创 图论知识点记录
//此知识点记录一下学习的思路：比较经典的题目了，首先按照边权分离边，从小到大添加进树中，可以知道相同边权的边形成的联通分支是一样的，所以在上一次操作存留下来的图中对每种边权先求一次生成树，记录哪些联通分支被连起来了，再按照这些联通分支进行构图，求这个图中不同的生成树的个数，这一步按照Kirchhoff定理求幺模矩阵的秩就可以了原博客：http://hi.b
2013-07-23 11:30:49
857
原创 hdu 湫秋系列故事——安排座位(组合dp)
/* * 以前比赛的时候做这题，一点思路都没有，后面就放了放； * * 今天突然想来弄下这题，于是去拜访了大牛的博客，了解到了组合dp的解法 * * 解法如下： 首先我们令F[i][j]表示前i个系的所有学生坐在座位上，出现了j对相邻学生是同一个系的 * * 那么对于接下来一个系的a个学生，我们可以把这个系的学生插入前面出现的j对相邻学生的中间，从而使他们变得不相
2013-07-22 16:29:15
1098
原创 hdu 3939 (Sticks and Right Triangle)
//数论（勾股数组相关），在acdream 群里看到的实现代码如下：#include#include#include#includeusing namespace std;vector fac[1000001];int prime[100000],phi[1000001],m;void init(){ //预处理phi[] int i,j; m=0
2013-07-22 13:29:06
893
原创 几何代码总结
struct Point{ int x,y; Point(int a=0,int b=0):x(a),y(b){}};Point operator-(Point a,Point b){ return Point(a.x-b.x,a.y-b.y); }Point operator+(Point a,Point b){ return Point(a.x+b.x,a.y+b.y)
2013-07-12 20:21:26
980
原创 zoj 1030 (dfs+几何)
//从每一个点开始，向某一个几何方向dfs搜索几何圈，然后判断几何圈里有没有其他的点，再把几何圈去重就是答案，O(n^3).#include#include#include#include#include#includeusing namespace std;struct Point{ int x,y; Point(int a=0,int b=0):x(a),
2013-07-12 20:04:15
904
原创 hdu 1625(flyd判正环)
//自己能到自己，则表示这个点在环上；#include#include#includeusing namespace std;int mat[31][31],n,m;int main(){ int a,b,i,j,k,cas=0; while(scanf("%d",&m)!=EOF) { memset(mat,0,sizeof(ma
2013-07-11 21:05:23
871
原创 hdu 3666 THE MATRIX PROBLEM(判负环)
解题思路：(U,L,mat[i][j])>0;因为题目要判定 L因为减号具有传递性（比如a-b=5,b-c=5由传递性得a-c=10），因为图即具有传递性（可以求出一颗最短路径树的图），又具有矛盾性（负环，下面会详细介绍），我们就可以把它抽象成图来判定判定负环log(ai)-log(bj)log(L/mat[i][j])=-log(U/mat[i][j])；如果
2013-07-10 19:29:28
732
原创 CodeforcesRound#191(Div.2)
A题：直接枚举#include#includeint A[110],c[110];int main(){ int n,i,j,ans,sum,x; while(scanf("%d",&n)!=EOF){ sum=x=0; memset(c,0,sizeof(c)); for(i=1;i<=n;i++) scanf("%d",&A[i]); f
2013-07-05 11:10:00
865
2
转载 ACM-ICPC 2013 World Finals 简单题解
原博客地址：http://jiazhipeng.byethost7.com/archives/268
2013-07-04 23:23:20
1984
原创 hdu 4328 (悬线法+dp)
悬线法用于求解最大子矩阵问题，资料详见：这篇论文对于red-and-blue crisscross 方阵可以用dp，既可以用合并两个小方阵得到一个大方阵，合并过程如下图： 我们可以把两个绿色的3*3小方阵合成一个4*4大方阵，只要符合color[i][j]!=color[i-1][j]&&color[i][j]!=color[i][j-1]&&color[i][j]==color[
2013-07-02 14:10:46
939
原创 hdu 4483(数论)
/* 这是一个经典的题目，以前遇到n范围都比较就小，可以直接枚举每一个长方形的两个对角点，最后答案就是 C(n*n,3)-2*n*C(n,3)-2*sum{sum{(n+1-i)*(n+1-j)*(gcd(i,j)-1)}}(1 然而这题的n范围比较大，直接枚举肯定不行，但是我们可以优化上面的这个公式； 第一次看到这题的时候，不会优化，到后面我做了下uva1
2013-06-30 23:06:18
944
原创 2sat 题目练习
2sat 问题的练习这个博客收集很多：http://www.cnblogs.com/ambition/archive/2011/07/30/2-sat.html hdu 3062：裸的，（基本的2选一，限制的2选1，如果能在限制的条件下实现基本的2选1表示题目有解），代码如下：#include#includestruct E{ int t,next;}edge
2013-06-30 14:48:57
742
原创 ZOJ 3717 Balloon (2sat+二分)
这题是一个二分半径+判定的题，但是开始不知道怎么去判定，以前听别人提起过sat问题，感觉这题是一个2sat模型，因为有一些限制比如：同一组两个里面只能选一个，两个球重叠时只能选其中的一个等，这很符合2sat的模型，于是去翻了下白书上2sat问题的解法实现，建图后直接用搜索判定，我用这种方法写了写，开始提交wa了几次，wa的原因是：题目说在round之后也要保证不重叠，那就只有舍掉小数点3位以后的数
2013-06-30 13:45:33
1488
原创 Codeforces Round #190 (Div. 2)
A题：直接输出就好了，除了最后一个人不能跳外，其他的每个人first time dancing一次，总共可以跳 n+m-1次，代码如下：#include#includeint main(){ int n,m,i; while(scanf("%d%d",&n,&m)!=EOF) { printf("%d\n",m+n-1); for
2013-06-29 19:16:45
756
原创 函数式线段树(主席树)学习实践spoj10628
//count the tree/*思路：每个点建立一颗线段树（增量建立），以遍历的时间为序，充分利用上一颗线段树的信息，在这题上一颗线段树就是父节点的线段树，因为我们每次更新的信息只有一个节点，一个节点被更新了，那么它的所有祖先节点也要相应的被更新，又因为在线段树中一个节点的祖先节点数不会超过（logN）个，所以这颗线段树和上一颗线段树大部分节点是一样的，只有刚刚说的logN个节点被改变
2013-06-04 21:40:30
1158
原创 以前写的AVL树实现map练习
AVL树的资料见这：http://blog.csdn.net/laziercs/article/details/8280074实现代码（范型模版实现方式，若有不妥的地方，请大家广言）：#include#include#includetemplate struct Type{ Tk key; Tv val;};template Type makepair(T
2013-06-04 19:38:31
1078
原创 2013金山西山居初赛第一场（hdu4545-hdu4547)
A题:用短串去长串里里面找自己（找的时候把字符转换一起考虑），找到就是happy；代码：#include#includeint mat[30][30];char s[1100],str[1200],a[3],b[3];int main(){ int cas,m,i,cass; scanf("%d",&cass); for(cas=1;cas<
2013-05-17 22:17:00
919
原创 hdu 4486 Pen Counts
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4486 题目大意： 把一个长度为n的绳子分成三段，用这3段绳子组成三角形，问有多少个不同的三角形 如果一个三角形的三边都不相等，可以算两次，其它的情况算一次； 解题思路： 假设一个三角形的三条边长为a,b,c，且a 这样我们可以枚举三角形的a边，再求在a边确定情况
2013-05-05 22:03:51
1130
2
原创 xmu网宿科技杯 厦门大学第四届程序设计积分赛 第四场
A题：题目链接：http://acm.xmu.edu.cn/JudgeOnline/problem.php?id=1465 题目描述： 将含有n个整数的数列A中的所有整数重新从小到大排序以后得到的新的数列B, 满足, 对于任意正整数i, j(1 个整数列C, 你需要求出将其变为连续数列所需花费的最小代价和。(将整数x变为整数 y需要花费|x - y|的代价, 所
2013-05-05 21:26:19
2475
2
原创 第十一届北京师范大学程序设计竞赛解题报告
A题：模拟#include#include#include#include#includeusing namespace std;int main(){ int i,n,cas; vector month[15]; month[11].push_back("Basic Training"); month[12].push_back("Basic
2013-04-30 21:42:38
1068
原创 2013 福州大学第十届程序设计竞赛
A题：中学物理题，直接求就好了；#include#includeint main(){ __int64 cas,s,h,vx,vy; scanf("%I64d",&cas); while(cas--) { scanf("%I64d%I64d%I64d%I64d",&s,&h,&vx,&vy); vx=s
2013-04-29 22:12:40
1480
9
原创 自己写的splay模板
#include#include#includeusing namespace std;#define N 300010const int INF=1<<30;int tree[N][2],pa[N],cnt[N],rev[N],val[N];int Max(int a,int b){ return a>b?a:b; }class splaytree{ public: i
2013-04-28 09:26:05
814
原创 DLX 模板
//以hust1017为例#include#include#define N 1100#define M 1100#define V N*M#define INF 100000000int D[V+1],U[V+1],L[V+1],R[V+1],C[V+1],mark[V+1];int S[M+1],H[N+1];void remove(int c){ int i,j
2013-04-18 21:52:40
772
转载 最小 k 度限制生成树
原文地址：http://blog.163.com/surgy_han/blog/static/2079912702012112884113774/所谓最小 k 度限制生成树，就是指有特殊的某一点的度不能超过k时的最小生成树。简单思想：设特殊的那点为v0，先把v0，删除，求出剩下连通图的所有最小生成树，假如有m棵最小生成树，那么这些生成树必定要跟v0点相连，也就是说这棵生成树的v0点
2013-04-06 09:58:42
930
原创 ZOJ3690Choosing number
/*题目：Choosing number题目连接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3690题目大意：n个人站成一行，这儿有m个数字，1，2...m，每个人选择一个数字，但是如果相邻的两个人选择了同一个数字这个数字不能超过k，请问他们有多少种选择数字的方式？ 解题思
2013-04-02 20:00:11
764
2
原创 ZOJ 3691 Flower
/*题目： Flower题目连接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4974题目大意： 三维空间中有n个点，序号为1-n，每个点有一个坐标，每个点上有F[i]朵花，每个点只允许L[i]朵花从这个点转移出到别的点，每两个点之间的距离为欧几里得距离，现在Gao的女朋友要把序号2-n点
2013-04-02 19:13:05
1180
原创 hdu4533-威威猫系列故事——晒被子
/*把(0,0),(t,t)看成是一个大矩形的话，那么这个大矩形的右上坐标x是等于y的，有了这个就好办了，我们可以维持一颗关于t的线段树，比如现在对一个X矩形（x1，y1）,（x2,y2）来说如果t>=Max（x2,y2），那么这个面积直接加上；在关于t的这颗线段树上操作也就是相当于更新(Max(x2,y2)~Max(t))这个区间，而对于（0~MAx(x1,y1)）这个区间是无影响的
2013-03-29 21:54:29
3515
4
转载 左偏树的讲解
//本文转自：这里2.1.2 可并堆的定义可并堆(Mergeable Heap)也是一种抽象数据类型，它除了支持优先队列的三个基本操作(Insert, Minimum, Delete-Min)，还支持一个额外的操作——合并操作：H ← Merge(H1,H2)Merge( ) 构造并返回一个包含H1和H2所有元素的新堆H。 O(n)，用它来实现可
2013-03-27 17:20:05
1201
转载 base64 编码解码（来自我同学）
base64是一种基于64个可打印字符来表示二进制数据的表示方法。由于26=64，所以每6位为一个单位，对应某个可打印字符。三个字节共24位，对应于4个base64单位，即3个字节需要用4个可打印字符来表示。它常用来作为电子邮件的传输编码。在base64中的可打印字符包括大写英文字母A-Z，小写英文字母a-z、阿拉伯数字0-9，这样共有62个字符，此外两个可打印符号在不同的系统中而不同，通常用加号
2013-03-25 12:11:01
1326
原创 矩阵快速幂
void mult(int a[2][2],int b[2][2]){ int i,j,k,c[2][2]={0,0,0,0}; for(i=0;i<2;i++) for(j=0;j<2;j++) for(k=0;k<2;k++) c[i][j]=(c[i][j]+(a[i][k]*b[k][j])%mod)%mod; for(i=0;i<2;i++)for(j=0;j<
2013-03-23 23:03:41
665
原创 hdu 3308 LCIS
/*题型：线段树；解题思路： 在线段树的节点上保存lmax,rmax,max三个值，lmax表示从该节点左端点数递增的最长长度，rmax表示终点是该节点的右端点的最长长度，max表示该区间的最长长度更新就是：如果左孩子的右端点小于右孩子的左端点，那么就合并左孩子rmax和右孩子的lmax；查询时和更新相似； 代码如下：*/ #include#inclu
2013-03-17 18:05:29
579
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人