2015年04月_poursoul

原创【BNU】40719 Arithmetic Progressions【分块+FFT】

传送门：【BNU】40719 Arithmetic Progressions题目分析：用分块+FFT强行AC了这题…… 之前一直TLE……然后改了好久把姿势改的优美点了……终于过了……大概思路是：我们考虑分块，假设每一块的大小为S，一共分了B块然后我们分两种情况讨论： 1.第二个数在第i块，第一个数在(1~i-1)块内，第三个数在(i+1~B)块内。 2.至少两个数在同一块内。对于第一种情况，

2015-04-30 22:06:29 1049

原创【HDU】5197 DZY Loves Orzing 【FFT启发式合并】

传送门：【HDU】5197 DZY Loves Orzing题目分析：首先申明，我不会dpdp方程= =……这个东西给队友找出来了，然后我就是套这个方程做题的Qrz……对于这题，因为n2n^2个数互不相同，所以每一列都可以单独考虑。设dpnidp_ni表示长度为nn的排列，能恰好看见ii个人的方案数，根据队友的发现，dpnidp_ni就等于|sni||s_ni|，其中snis_ni是第一类Stirl

2015-04-30 15:56:11 1418

原创【ZOJ】3874 Permutation Graph 【FFT+CDQ分治】

传送门：【ZOJ】3874 Permutation Graph#include <stdio.h>#include <string.h>#include <math.h>#include <algorithm>using namespace std ;typedef long long LL ;#define clr( a , x ) memset ( a , x , sizeof a )

2015-04-29 22:22:50 1828

原创【SPOJ】Triple Sums【FFT】

传送门：【SPOJ】Triple Sums题目分析：首先我们不考虑i<j<ki<j<k这个条件，构造多项式： Y=∑xai\qquad\qquad\qquad Y = \sum x^{a_i} 那么ai+aj+ak=Sai+aj+ak=S的个数即xai+aj+ak=Sx^{a_i+a_j+a_k=S}的个数，等价于Y3中xSY^3中x^S的系数。然后我们考虑容斥： (∑x)3=∑x3+3∑

2015-04-29 16:49:15 1088

原创【BNU】33943 Super Rooks on Chessboard 【FFT】

【BNU】33943 Super Rooks on ChessboardUVA上的题，然而我怎么会蠢到去UVA呢！（其实是百度首先跳出来的是BNU→\to_→\to）题目分析：设numxnumx为NN个车没有覆盖的行数，numynumy为NN个车没有覆盖的列数。首先我们考虑没有主对角线覆盖这一条件时，总共的没有被覆盖的面积就是numx∗numynumx \ast numy。现在我们考虑主对角

2015-04-29 15:27:29 1243

原创【HDU】4609 3-idiots 【FFT】

传送门：【HDU】4609 3-idiots

2015-04-28 18:01:30 1122

原创【HDU】1402 A * B Problem Plus 【FFT】

传送门：【HDU】1402 A * B Problem Plus题目分析：这就是大数乘法题，问两个大数相乘的结果，由于O(n2)O(n^2)的算法复杂度太大，所以我们用FFT来优化他。关于FFT网上资料很多，我就不多说啦。这是我做的第一道FFT，FFT是看算法导论学来的，感觉算导讲的很不错，简单易懂~

2015-04-28 11:13:57 832

原创【HDU】5213 Lucky 【分块（在线算法）】

传送门：【HDU】5213 Lucky 【分块（在线算法）】题目分析：我来说一下这题的在线做法。首先我们将区间分成n√\sqrt n块，用f[x][y]表示第x块的数和第y块的数相加等于K的对数，这个可以O(nn√)O(n \sqrt n)的预处理。然后还有g[x][y]表示在第1~x块中有的大小为y的数的个数，这个的复杂度同样O(nn√)O(n \sqrt n)。接下来，对于每组询问，我们考虑

2015-04-26 23:17:09 934

原创【HDU】4348 To the moon 【可持久化线段树】

传送门：【HDU】4348 To the moon题目分析：可以将标记留在节点上，查询的时候累加就行了，这样避免了可持久化线段树打标记内存吃紧的窘态。my code：#include #include #include using namespace std ;typedef long long LL ;#define rep( i , a , b ) f

2015-04-22 21:38:18 2014

原创【FZU】Problem 2137 奇异字符串【后缀数组】

传送门：【FZU】Problem 2137 奇异字符串题目分析：枚举x所在位置，向左右暴力扩展，lcp(L,x+1)>=x-L就累加ans。复杂度是O(nlogn+26*n)，复杂度可以这么证明，假设我们以字母a为中心，那么不存在字母a的位置最多n-1个，直接枚举过去就好了，其他字母同理。my code：#include #include #include #in

2015-04-22 19:40:31 641

原创【BZOJ】3289: Mato的文件管理【区间逆序对，莫队套树状数组】

传送门：【BZOJ】3289: Mato的文件管理my code：#include #include #include #include using namespace std ;typedef long long LL ;#define rep( i , a , b ) for ( int i = ( a ) ; i < ( b ) ; ++ i )#defi

2015-04-21 23:55:29 1451

原创【BZOJ】2286: [Sdoi2011消耗战【虚树DP】

传送门：【BZOJ】2286: [Sdoi2011消耗战my code：#include #include #include using namespace std ;typedef long long LL ;#define rep( i , a , b ) for ( int i = ( a ) ; i < ( b ) ; ++ i )#define For(

2015-04-21 21:17:10 808

原创【POJ】1509 Glass Beads 【后缀自动机】

传送门：【POJ】1509 Glass Beadsmy code：#include #include #include #include using namespace std ;typedef long long LL ;typedef unsigned long long ULL ;typedef pair pii ;#define rep( i , a

2015-04-21 16:05:38 1091

原创【BZOJ】3572: [Hnoi2014]世界树【虚树DP】

传送门：3572: [Hnoi2014]世界树my code：#include #include #include #include using namespace std ;typedef long long LL ;typedef unsigned long long ULL ;typedef pair pii ;#define rep( i , a ,

2015-04-21 12:48:13 1124

原创【HDU】5208 Where is Bob 【DP】

传送门：【HDU】5208 Where is Bob题目分析：设状态为dp[cur][l1][r1][l2][r2]，cur表示当前二进制位的位置，l1表示第一个人的数的下界，r1表示第一个人的数的上界，l2表示第二个人的数的下界，r2表示第二个人的数的上界。l1，r1，l2，r2都是01变量，对于第一个数，如果数卡着下界不能变小时l1置为0，否则l1置为1表示可以变小。第一个数的

2015-04-19 19:38:16 1153 1

原创【Codeforces】Codeforces Round #299 (Div. 1) E. Tavas on the Path 【树链剖分+区间合并】

传送门：【Codeforces】Codeforces Round #299 (Div. 1) E. Tavas on the Path大概题意：首先对于一个串s，我们可以提取m个只包含1的块，每个块都是s里的一个最长连续1子串。然后我们设x1,x2,x3...xm分别为这m个块中1的长度。现在我们有函数T(S)，定义，f(xi)是f关于xi的函数，下面会给出。现在我们

2015-04-17 23:51:42 990

原创【ZJU】3863 Paths on the Tree 【浙大2015年4月校赛D题】树分治

传送门：【ZJU】3863 Paths on the Tree题意：给一棵树，问树上有多少个路径对有不超过K个公共节点的，路径a->b和b->a等价，路径对（A，B）和（B，A）只有当A和B是同一条路径时相同。分析：反过来考虑，考虑有超过K+1个公共节点的路径对数。我们考虑重叠的路径部分，这个可以用树分治来搞，然后路径对的两端延伸出去的部分不重叠，我们要预处理出这个部分。最

2015-04-13 23:24:42 1133

原创【一道不知道名字的题】AC自动机+dp

太懒了，好久没写题解了。Time Limit : 3000/3000ms (Java/Other) Memory Limit : 65535/65536K (Java/Other)Total Submission(s) : 34 Accepted Submission(s) : 8Problem DescriptionLet Str be a long string

2015-04-07 12:52:01 868

poursoul