AC，∑ndless

先预处理下：在每个字符的两边都插入一个特殊的符号，比如abba变成#a#b#b#a#，aba变成#a#b#a#（因为Manacher算法只能处理奇数长度的字符串）。同时，为了避免数组越界，在字符串开头添加另一特殊符号，比如$#a#b#a#。以字符串3212343219为例，处理后变成S[] = "$#3#2#1#2#3#4#3#2#1#9#"。然后用一个数组Len[i]来记录以字符S[i

#includeusing namespace std;int a[140];bool ok[140];/*输出从n个数中选m个数的组合字典序以n=7,m=3为例首先构造ok=1110000然后ok的前一个排列为1101000再前一个排列为1100100...1100001然后是1011000...直到0000111*/void printfC(int n, int m

http://acm.hdu.edu.cn/showproblem.php?pid=2602/*31ms,240KB*/#include#include#includeusing namespace std;const int mx = 1005;int w[mx], v[mx], dp[mx];void _01pack(int n, int maxw){ me

#includeusing namespace std;const int mx = 100005;int fa[mx], rk[mx];int find(int x) {return fa[x] = (fa[x] == x ? x : find(fa[x]));}void union(int x, int y){ x = find(x), y = find(y); if (

http://poj.org/problem?id=3468/*2391ms,4508KB*/#include #include using namespace std;#define lson l , m , rt << 1#define rson m + 1 , r , rt << 1 | 1#define root 1, N, 1#define LL long l

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=505&page=show_problem&problem=4147递推写法：L = strlen(text);for (i = L - 1; i >= 0; --i){ trie.find_prefixes(text +

此题相当于解方程组x+d≡p(mod 23)x+d≡e(mod 28)x+d≡i(mod 33)由于23,28,33两两互素，所以M=23*28*33=21252，M1=924，M2=759，M3=644.则有924y1≡1(mod 23)，由924/23=(40,5,1,3)计算得y1=6(mod 23)则有759y2≡1(mod 28)，由759/28=(27,9,2,1)计算得y2=19(mod 28)则有644y3≡1(mod 33)，由644/33=(19,1,1,16)计算得y3=

http://poj.org/problem?id=1276模板题。完整代码：/*110ms,4316KB*/#include#includeconst int mx = 1005;const int mxw = 1000005;int w[mx], m[mx], dp[mxw], deq[mxw], deqv[mxw];int solve(int m

#include#include#include#include#include#includeusing namespace std;const int mx = 10005;typedef pair P; ///first是最短距离，second是顶点编号struct edge{ int to, cost;};vector G[mx];int disTo[mx]

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=941由于n很小，floyd算法写起来方便，先用这个A了一下：/*0.162s*/#includeusing namespace std;const int mx = 105;i

测试数据：8 135 4 0.354 7 0.375 7 0.285 1 0.324 0 0.380 2 0.263 7 0.391 3 0.297 2 0.346 2 0.403 6 0.526 0 0.586 4 0.93测试结果：5 to 0 : 0.735 to 1 : 0.325 to 2 : 0.625 t

一个人的旅行http://acm.hdu.edu.cn/showproblem.php?pid=2066Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Problem Description虽然草儿是个路痴（就是在杭电待了一年多，居然还会在校园里

代码思想来自埃拉托斯特尼筛法，参考：http://blog.youkuaiyun.com/synapse7/article/details/18727405#include const int mx = 50001;int prime[mx], phi[mx];bool unprime[mx];///O(NloglogN)，推荐void phi_table(){ int i, j

一、简单介绍参加wiki页面：http://en.wikipedia.org/wiki/Partition_(number_theory)二、基本形式例题：POJ 3014 Cake Pieces and Plates题意：(m的n划分的总数)给你m个东西，放在n个相同的盒子中，且每个盒子可以放任意多，问有多少种放法。思路：记

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2733算是一种模板题吧。。两个发现：1. 重复初始化queue相当费时间！2. 能不用iterator就不用！完整代码：/*0.732s*/#includeusi

10305 - Ordering TasksTime limit: 3.000 seconds http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=24&page=show_problem&problem=1246模板：/*0.022s*/#includeusi

http://poj.org/problem?id=1113使用Graham-Scan算法。为什么扫描阶段的复杂度是O(n)的？因为在扫描的过程中，每个点至多进栈一次，至多出栈一次，也就是说对每个点我们至多进行两次叉积运算，所以复杂度是O(n)的。完整代码：/*16ms,592KB*/#include#include#includeusing name

复杂度均为O(VE)DFS实现，适用于稠密图：#includeusing namespace std;const int mx = 105;vector G[mx];int match[mx]; ///match表示匹配的对象编号bool vis[mx];bool dfs(int i){ vis[i] = true; for (int j = 0; j < G[i].

http://codeforces.com/problemset/problem/401/D解释全部在代码的注释中：/*78ms,205464KB*/#includeusing namespace std;const int mx = 1 << 18;long long dp[mx][100];///dp[mask][j]表示余数为j时的mask对应的x的个数///

【背景】如何求出在给定区间[A,B]内，符合条件P(i)的数i的个数？条件P(i)一般与数的大小无关，而与数的组成有关，有一下几种P(i)：数i是递增/递减的：1234, 2579,…双峰的：19280,26193,…含/不含某一数字的，比如含49：49, 149, 1492,… 被某一数m整除的，比如m=13：39,130,650...【思路】采用

题目：求一个一维数组arr[n]中的最长递增子序列的长度，如在序列1,5,8,3,6,7中，最长递增子序列长度为4 (即1,3,6,7)。方法一：一般的DP方法(O(N^2))像LCS一样，从后向前分析，很容易想到，第i个元素之前的最长递增子序列的长度要么是1（单独成一个序列），要么就是第i-1个元素之前的最长递增子序列加1，这样得到状态方程：

http://poj.org/problem?id=1141怎么记录路径？定义divide_pos[i][j]为i到j这一段的最佳添加括号位置，若这一段本身就不需要添加括号，则置-1完整代码：/*0ms,416KB*/#include#includeconst int mx = 100;char s[mx];int dp[mx][mx];///最小添加

一、类型转换char数组转int、__int64或double：#include#include#include#includechar str[100];int main(void){ char *stopstr; gets(str); if (strchr(str, '.') == NULL) { if (strlen(str) < 10) print

测试数据：8 164 5 .354 7 .375 7 .280 7 .161 5 .320 4 .382 3 .171 7 .190 2 .261 2 .361 3 .292 7 .346 2 .403 6 .526 0 .586 4 .93测试结果：1.81代码：#include#include#

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=24&page=show_problem&problem=241注意代码中的递推式并不能指出实际的位置，但是要让m是满足题意的，这个位置（在减一映射后）必须>=k当然你也可以打表实现。/*0.079s*/#include

/*请从main开始看*//*SampleInput:6 8 70 52 42 31 20 13 43 50 21 40 01 12 23 34 45 54 2 20 12 32 10 0*/#include#include#include#include#includeusing namespace std;const int

#include #include #include using namespace std;const int mx = 1000000 + 1; ///在(1,mx)的范围内寻找素数const int sqrt_mx = (int)sqrt((double)mx);bool vis[mx];int prime[mx / 10]; ///在mx>65000时建议写成 int pr

#include#include#include#include#include#includeusing namespace std;const int Size = 32;const double log_2 = log(2.0);char str[Size];int main(void){ int num; char *endstr; while (~scan

http://lightoj.com/volume_showproblem.php?problem=1255/*0.068s,7548KB*/#includeusing namespace std;const int mx = 1000005;char t[mx], p[mx];int f[mx];void getfail(){ f[0] = f[1] = 0

http://acm.hdu.edu.cn/showproblem.php?pid=4741模板题。理论知识见代码注释。这题背景居然是命运石之门。。/*218ms,276KB*/#include#includeconst double eps = 1e-9;struct P3{ double x, y, z; P3(double x = 0.0, do

http://poj.org/problem?id=2932先给圆的最左端和最右端的点排个序，当两点x相同时，左端点排在前面。然后就是扫描了，若扫描到的是圆的左端点，就判断圆心(y坐标)在其上方且离其最近的圆是否包含此圆，以及圆心(y坐标)在其下方且离其最近的圆是否包含此圆，若包含就continue，不包含就insert到set中；若扫描到的是圆的右端点，就从set中era

树的最近公共祖先(Lowest Common Ancestor)问题是树结构上最经典的问题之一。给一棵树T，每个询问形如：“点u和点v 的公共祖先是哪个点？”，此问题的答案被记为LCA(u, v)。LCA问题的算法分为在线和离线（Tarjan算法）两种，前者要求在回答后一个问题之前必须给出前一个问题的输出，而离线问题允许在读入所有询问之后一次性给出所有问题的答案。LCA问题的应用很多，例如

311 - PacketsTime limit: 3.000 seconds http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=113&page=show_problem&problem=247A factory produces products packed in sq

443 - Humble NumbersTime limit: 3.000 secondshttp://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=24&page=show_problem&problem=384A number whose only prime factors a

****转载请注明http://blog.youkuaiyun.com/synapse7/article/details/19096049****概述：使用getchar()和putchar()，以及math.h头文件中的一些函数，基本实现了以下函数scanf("%u",&x)scanf("%d",&x)scanf("%lf",&x)printf("%u",x)printf("%d",x

http://acm.hdu.edu.cn/showproblem.php?pid=1166此题考查RSQ(Range Sum Query)。模板如下：/*343ms,744KB*//*我们使用如下缩写：root--整棵树的根节点rt--当前区间的根节点l,r--区间左右端点m--区间的二等分点，注意m是落在左子区间的(右端点)lson,rson-

首先介绍扩展欧几里得定理：对于两个不全为0的整数a,b，必存在一组解x,y，使得ax+by=gcd(a,b)换句话说，形如ax+by的最小正整数等于gcd(a,b)。实现代码如下：（一般题目都要用64位）（复杂度：O(log max(a,b))）

零、先修课程首先，在开始理解DP的思想前，你需要1. 完成HDU里面的递推求解专题练习（For Beginner）那7道题（这些题很简单，题解请在博客中搜索），这对你理解DP有很大的帮助。2. 对递归搜索（比如深度优先搜索，DFS）有一定了解。一、递归与记忆化搜索我们从POJ 3176入手来学习这一思想。（题目很短，请快速读完）从上往下看，最大和值无非是往左走和往右

一、威尔逊定理若p为质数，则p|(p-1)!+1亦：(p-1)! ≡ p-1 ≡ -1(mod p)例题：HDU 2973 YAPTCHA (威尔逊定理及其逆定理)解题报告见http://blog.youkuaiyun.com/synapse7/article/details/18728157二、费马小定理假如p是质数，且gcd(a,p)=1，那么a^(

/*请从main开始看*//*SampleInput:6 8 30 52 42 31 20 13 43 50 21 52 42 54 2 10 12 32 1*/#include#include#include#include#includeusing namespace std;const int maxn = 1000002;li