rpSebastian

摘要： var n:longint; nt:array[0..200008,0..26] of longint; size,l,fa:array[0..200008] of longint; s:array[0..200008] of char; cnt,last:longint; procedure in 阅读全文

摘要： 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... 阅读全文

摘要： #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... 阅读全文

摘要： 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... 阅读全文

摘要： 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... 阅读全文

摘要： 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... 阅读全文

摘要： 给定方程 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... 阅读全文

摘要： 求一个字符串的最小表示法。 将字符串S倍长，从根走length(s)步所走路径即为最小表示法。 记所到达位置为x，则这个最小表示法的起点为a[x]-len(s)+1 阅读全文

摘要： 用后缀自动机统计出出现1~n次的串的数量f[i] 对于ans[k]=sigma（f[i]*C(i,k)) i>=k 阅读全文

摘要： http://blog.sina.com.cn/s/blog_7812e98601012cim.html //后缀自动机建立的详细介绍 http://www.tuicool.com/articles/Mjuu2y //后缀自动机学习指南（习题列表） 阅读全文

摘要： 一个正整数的价值就是把这个数的十进制写出来之后，最长的等差子串的长度。求出在[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... 阅读全文

摘要： 由于方块最多涉及3行，于是考虑将每两行状压起来，dfs搜索每种状态之间的转移。这样一共有2^12种状态，显然进行矩阵快速幂优化时会超时，便考虑减少状态。进行两遍bfs，分别为初始状态可以到达的状态，和可以到达终止状态的状态。同时出现在两次bfs中的状态即为有效状态，一共有141种。这样就可以跑出来了... 阅读全文

摘要： 计算若干个圆与一个矩形的面积并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,... 阅读全文

摘要： 首先 1+x+x^2+x^3+...+x^∞=1/(1-x)对于题目中的几种食物写出生成函数 (对于a*x^b , a表示方案数 x表示食物，b表示该种食物的个数）f(1)=1+x^2+x^4+...+x^∞=1/(1-x^2)f(2)=1+xf(3)=1+x+x^2f(4)=x+x^3+x^5+.... 阅读全文

摘要： 1 var s,t:ansistring; 2 n,op:longint; 3 p:array[0..2000008] of longint; 4 procedure pre; 5 var i:longint; 6 begin 7 s:='$*'; 8 for i:=1 to length(t) do 9 begin 10 ... 阅读全文