2012年03月_huyuncong

原创【线性筛法求解积性函数】Archer

计算sigma（a=1~n）sigma（b=1~n）lcm(a,b). oimaster的好题。考虑a与b枚举是同域的，对于一个a只需统计比他小的的b，然后乘2即可，则构造函数f[a]=-a+2*lcm(a,b)【b=1~a】；设gcd（a,b）=d f[a]=-a+2a*sigma(d|a)sigma(b 令b/d=k, f[a]=-a+2a*sigma(d|a)sigma(

2012-03-10 15:55:48 935

原创【树链剖分+维护凸线】Saber

秋哥是saber控，每次saber都很难，但是还有两种颜色的saber没出。有一棵 n 个点的树，每个点有四个权值 x,y,p,q，给出 m 个询问(a,b)，假设 i,j 为 a 到 b 的路径上的可以重合的两个点，求（yi+qj）/(xi+pj)的最大值。二分答案，就化出了（x,y）与（p，q）分开的式子，分别对此两式求最大值判断，用斜率的经典分析来看（x,y）若有一k比i优

2012-03-10 15:28:53 933 4

原创【数位统计】Tickets

有一位售票员给乘客售票。对于每位乘客，他会卖出多张连续的票，直到已卖出的票的编号的数位之和不小于给定的正数K。然后他会按照相同的规则给下一位乘客售票。初始时，售票员持有的票的编号是从L到R的连续整数。请你求出，售票员可以售票给多少位乘客。数据规模：1 ≤ L ≤ R≤ 1018，1 ≤ K ≤1000。刘聪论问题，以前做此类题都是按位考虑，可是偏偏此题用树考虑极为方便（居然不满足区

2012-03-10 15:17:12 639

原创【构造】后缀数组求逆

根据sa[]数组反求原数组，只能26小写字母。有了sa[]，自然有了rank[]，从原数组最后的往前考虑，我们尽量要将原数组用小的字符集表示，因为只有26个字母可以表示，那么排名比当前枚举的小一位的的字符串，如果第二位比其小的话，那么排名比当前枚举的小一位的的字符串的首位可与其首位相同，同样的向比它大的考虑一遍则可使字符集尽量小。 #include #include #include

2012-03-10 15:05:08 694

原创【矩阵优化连通性状压dp】哈密尔顿路径

以前觉得插头dp好晕啊，结果稍微推了一下转移发现并不难推，只是转移多了点罢了，可是noip模拟的时候270行的恶心转移dp都打过，100+的状压也就还好了。本题有一维特别大，于是我们压缩小的那一维，然后裸转一行后，用矩阵加速即可。转移第一次写会觉得麻烦，最后自己推一边在跟别人推的对一下检查有没有漏转或多转，至于要一模一样则没必要，我完全是按自己的喜好来的。 ural有道题是有障碍的哈密尔

2012-03-10 14:54:21 932

原创【mt noip模拟题】最大公倍树

求一个无向带权连通图的所有极小连通子图的所有边权的最大公约数的最小公倍数。更清楚的说就是，令带权连通图的一个极小连通子图的权值为所有边权的最大公约数，带权连通图的权值为所有极小连通子图的最小公倍数。答案可能很大，请对 2^31-1 取模。就如mt的xor题是分位考虑一般，此题要分质因子考虑，对于一个质因子最后出现在答案里是最大瓶颈树上的最大瓶颈边，最后求出lcm即可。

2012-03-10 14:43:33 756

原创【树链剖分+线段树维护最小费用流】大灾变

恶心题要写两个线段树还得是区间覆盖的，其中一个是树链剖分的，另一个是维护括号序列的。有两种边，一种是免费的，一种是有代价的，每次取最短路，树形结构大约就不要退流了吧。每次树链剖分找最短边，括号序列维护距离。 #include #include #include const long long oo=1073741819; const int maxn=200000,maxr=5

2012-03-09 22:13:02 998

原创【动态树】弹飞绵羊

简洁暴力动态树类似于树链剖分，只不过重边不以子树大小唯一确定，而是随着询问适时改变，从而摆脱静态的限制，可以动态添边删边。边分为实边与虚边，实边连的称为路径，用splay维护（以保证可以动态加删边），从一节点往根走时将路径上的虚边改为实边，通过splay加速上升，而两点lca为最后一个虚变实的节点。具体实现时，每个节点要存父亲，而根节点的父亲为此链的最浅处连出去的虚边父亲，虚变实时将新

2012-03-09 22:08:30 1761

huyuncong的专栏