06 2015 档案

斜率dp cdq 分治
摘要：f[i] = min { f[j] + sqr(a[i] - a[j]) }f[i]= min { -2 * a[i] * a[j] + a[j] * a[j] + f[j] } + a[i] * a[i]由于a[i]不是单调递增的，不能直接斜率dp。考虑有cdq分治来做，复杂度(nlog2n) 1... 阅读全文
posted @ 2015-06-27 17:10 rpSebastian 阅读(182) 评论(0) 推荐(0)
POJ2449 （k短路）
摘要：#include #include #include #include #include using namespace std; #define maxn 2008 #define maxm 2000008 #define INF 2000000000 int lt[maxn],LT[maxn],sum=1,SUM=1; int h[maxn]; int s,t,k,n,m; b... 阅读全文
posted @ 2015-06-26 20:45 rpSebastian 阅读(211) 评论(0) 推荐(0)
C++ 优先队列类实现
摘要：1 #include 2 #include 3 #include 4 #include 5 #include 6 #include 7 using namespace std; 8 9 int a[10000],b[10000]; 10 class t 11 { 12 public: 13 int x, y; 14 bool operatorp.x... 阅读全文
posted @ 2015-06-26 16:43 rpSebastian 阅读(250) 评论(0) 推荐(0)
BZOJ1576 （最短路+并查集）
摘要：1 #include 2 #include 3 #include 4 #include 5 #define inline __attribute__((optimize("O2"))) 6 using namespace std; 7 8 #define maxn 100008 9 #define maxm 1000008 10 #defin... 阅读全文
posted @ 2015-06-24 14:35 rpSebastian 阅读(342) 评论(0) 推荐(0)
SWUST0249 （凸包面积）
摘要：1 type node=record x,y:longint; end; 2 const maxn=3008; 3 var k,q,qq:longint; 4 sum:double; 5 f,g:array[0..maxn] of node; 6 m,i,j,a... 阅读全文
posted @ 2015-06-13 14:43 rpSebastian 阅读(234) 评论(0) 推荐(0)
道路修建 （网络流）
摘要：1 const INF=2000000000; 2 const maxn=4008; 3 var r:array[0..maxn] of longint; 4 eg:array[0..1000008] of record u,v,w,nt:longint; end; 5 el:longint; 6 lt:array[0..maxn] of long... 阅读全文
posted @ 2015-06-12 17:46 rpSebastian 阅读(244) 评论(0) 推荐(0)
HDU3930 （原根）
摘要：给定方程 X^A = B (mol C) ，求 在[0,C) 中所有的解 ， 并且C为质数。设 rt 为 C 的原根 , 则 X= rt^x （这里相当于求 A^x =B (mol C) 用大步小步算法即可）那么 ( rt^x ) ^ A = b (mol C) rt^Ax = b (mol C... 阅读全文
posted @ 2015-06-12 14:53 rpSebastian 阅读(260) 评论(0) 推荐(0)
ZOJ2006 (后缀自动机）
摘要：求一个字符串的最小表示法。 将字符串S倍长，从根走length(s)步所走路径即为最小表示法。 记所到达位置为x，则这个最小表示法的起点为a[x]-len(s)+1 阅读全文
posted @ 2015-06-10 16:16 rpSebastian 阅读(233) 评论(0) 推荐(0)
Codechef2015 May - Chef and Strings (后缀自动机）
摘要：用后缀自动机统计出出现1~n次的串的数量f[i] 对于ans[k]=sigma（f[i]*C(i,k)) i>=k 阅读全文
posted @ 2015-06-10 11:21 rpSebastian 阅读(314) 评论(0) 推荐(0)
后缀自动机
摘要：http://blog.sina.com.cn/s/blog_7812e98601012cim.html //后缀自动机建立的详细介绍 http://www.tuicool.com/articles/Mjuu2y //后缀自动机学习指南（习题列表） 阅读全文
posted @ 2015-06-10 08:44 rpSebastian 阅读(171) 评论(0) 推荐(0)
Digit (数位DP）
摘要：一个正整数的价值就是把这个数的十进制写出来之后，最长的等差子串的长度。求出在[l,r]范围内的数字的价值总和。(lb then exit(a) else exit(b); 8 end; 9 function dfs(now,ml,l,d,pre,st,lim:longint):int64;10 va... 阅读全文
posted @ 2015-06-07 18:15 rpSebastian 阅读(578) 评论(2) 推荐(0)
瓷砖铺放 （状压DP+矩阵快速幂）
摘要：由于方块最多涉及3行，于是考虑将每两行状压起来，dfs搜索每种状态之间的转移。这样一共有2^12种状态，显然进行矩阵快速幂优化时会超时，便考虑减少状态。进行两遍bfs，分别为初始状态可以到达的状态，和可以到达终止状态的状态。同时出现在两次bfs中的状态即为有效状态，一共有141种。这样就可以跑出来了... 阅读全文
posted @ 2015-06-03 21:04 rpSebastian 阅读(367) 评论(0) 推荐(0)
农场阳光 (simpson)
摘要：计算若干个圆与一个矩形的面积并simpson公式 ans = ( f[l] + f[r] + 4 * f[mid] ) * (r-l) / 6 1 uses math; 2 type arr=record x,y:double; end; 3 const eps=1e-8; 4 var a,b,n,... 阅读全文
posted @ 2015-06-03 20:36 rpSebastian 阅读(286) 评论(0) 推荐(0)