huyuncong的专栏

同1830，有关联系数为1，最后一列是是否要改变。由于题目要使元为1的尽量少，所以还要通过枚举解进行比较。我是用搜索进行枚举，将之前结果回带后，不定元行为（0,0,0,0....0）则有该元为1和0两种情况，为1时要改为（0,0,...1,....0）方便之后进行回带，否则为0；若为（0,0,0...1....0）则确定该元为1，改为（0,0,0...1...1）。最终解为最后一列为1的个数。var ans,n:longint; a,a1:array[0..22]of longint;proc

在n个数中选任意多个（不为0）进行xor，求结果中的第k小高斯消元后就可以确定自由元，对自由元0\1选取就可以快速确定第k小。有几个较为特殊的处理，首先是不能不选，因此0的选取要根据自由元的数量是否达到上限来判断，然后是对于非自由元必须跟着自由元变化，因此消元的时候要从头到尾消完，最后是必须从高位开始消元，否则的话，最高位如果不是自由元的话，就无法从高到低位来取0\1来逼近第k小数#i

互补松弛定理sigma(aij*yi | i=1~m)==cj or xi=0 j=1,2...,nsigma(aij*xj | j=1~n)==bi or yi=0 i=1,2...m题意：给一个n*n的矩阵，每次可以攻击一列或者一行，攻击一行的代价为Ai,攻击一列的代价为Bi,要求每个格子i,j至少被攻击F[i][j]次，要使得总代价最小，同时输出方案第一次接触对偶定理就是

之前用对偶定理解决了个网络流，这次又碰到了个半平面交...可惜漏了个约束条件，对偶后的变量少了一个首先把至少要跑的d先跑完，得到新的L和W设每个人跑的长度是xi,同时令xn=L-sigma(xi)我们的目标min(sigma(ti*xi)+tn*(L-sigma(xi)))变形得(-max(sigma((tn-ti)*xi)))+tn*L常数不去管它，只考虑max的式子再来

这个blog解释的挺不错:http://blog.youkuaiyun.com/acdreamers/article/details/16966369给定不相交的两个圆以及圆外一点，找一个经过给定点的圆与其他两个圆相切首先来看反演变换，首先是给定一个圆圆心为O，半径为R1、圆外一点P与圆内一点P‘会一一对应的反演OP*OP'=R*R2、经过O的圆，反演后成为不经过O的一条直线3、不经过O的

上学期写的一篇解题报告，还是传上来备份一下...算法：单纯形or对偶定理化为最小费用流对偶定理：max{cx | Ax≦b, x≧0}min{yb | yA≧c, y≧0}原题大意：给出一个n个点m条边的有向拓扑正权图，使得每个点都至少被一条由1到n的路径经过，现要求尽可能增加边权，使得由1到n的最长路不会增大。 抽象模型——>线性规划设修改图后1到i的最长距离

http://www.tsinsen.com/ViewGProblem.page?gpid=A1381

梭哈游戏题意：给定一系列比较规则，并给出A，B当前手牌，讨论牌抽完后，A能赢B的概率比较规则比较复杂，可以先判牌形，再将花色大小压成一个数，排序后我们就只要统计，在某副手牌之前的手牌中，有多少个与其无重复手牌，判重工作可以由容斥原理和hash表完成，对于每副手牌需要2^5判断复杂度为o（手牌数*2^5），手牌数最多只有c（51，5）5s时限还算充裕，对于hash，我们考虑用53进制压位，因

http://www.cnblogs.com/skyivben/archive/2008/07/23/1248413.html整理一下模板hdu1402sincos需要手写，因为hdu没有.

hdu4623通过打表发现很多数字在考虑互质的情况是dengjia

可以看出这些操作是可以用

上一篇：http://blog.youkuaiyun.com/huyuncong/article/details/18184873

找出不超过n次

hdu4304一段n*n的墙，用一个刷子来刷，要求刷子的第一根毛在第一个

matrix tree矩阵aij,表示i和j的lca的点权值

叉姐论文：http://www.docin.com/p-724323397.htmla[i]=sigma(a[j]*b[k-j]) 求第n项标准的常系数线性递推，用矩阵乘法可以做到o(k^3*logn)但是如果用特征多项式优化的话，可以把矩阵乘法换成多项式乘法，从而做到o(k^2logn)首先根据hamilton-cayley定理，一个矩阵的特征多项式是这个矩阵的化零多项式，因

整数拆分后状态最多38000+，由于已经连的边不能用，所以把还剩多少边压进状态，也就700+条边，然后记g[s][k]为s状态,省k条边的数学期望，直接列公式转移，n==40显然是要打表的，话说这道题大丰中学出过一次，也依稀记得见过，结果最后还是没想出来...把表全部打出来大概2~3分钟#include #include #include #include #include #d

用更相减损术的思想化为维护差分就很容易做了，lct与树链剖分均可，边界要注意，中间直接弄成开区间codechef常数也卡得很厉害...#include #include #include #include #include #include #include using namespace std;int a[200000],b[200000],gcc[200000],ad

方案平均贡献*方案总数=总和 总方案*符合要求的方案概率=符合要求的方案数 定义NB数字为这个数字十进制的每一位之和为S，累乘的乘积为2p2*3p3*5p5*7p7，现在告诉你S、p2、p3、p5、p7求所有满足条件的NB数之和对500500573取模。（莫名其妙三个小点没过，所以不贴代码）设a[1]...a[9]为每种数字出现次数，则只需枚举四个数字即可推出全体数字，根据可重排列

xor方程组

省选前几天就看到过矩阵快速幂的题，结果没引起注意，也就没去做。可是，省选第一天第一题竟然就是矩阵乘法（就是雅礼下至3岁幼童都会的那一道）。今天去看了xqz的博文，虽然xqz不能进入省队令人惋惜，雅礼所作所为也令人不齿，但是，从博文中可以看出，xqz在考试时信念还是有所动摇，oi有我们的梦，有我们的伙伴，有我们的坚守，不论在什么情况下我们永不言败，希望xqz尽快恢复过来，在我们面前的，是整个

原题，结果还是没做出来题意：定义n为正整数，S(n)为 n的各位数字之和，令d[n]:={s[n] n=10} 求给定l到r中有多少n满足n=x*d[x]，x任意。不难求出d[n]:=(n-1) mod 9+1,通过这个式子：（x-1） mod 9=d[x]-1x

数学没学好还真吃亏。首先介绍下辛普森公式（没法帖图片，只能口述），对于一段函数，我们要求面积，除了求反导之外，有一种近似方法，即用辛普森公式。（来自维基）a,b分别指函数最左端和最右端，这个公式推导可以去维基上看，他有个直观的应用便是求拟柱体（比如棱柱）体积，f函

ac自动机，令人神往的名字。。。ac自动机：http://www.cppblog.com/mythit/archive/2009/04/21/80633.htmlpoj2778：http://hi.baidu.com/%D2%D5%C1%D6010/blog/item/6

备战noip ing背包问题其实算np问题，一般的动规算是伪多项式。按时间顺序来吧、、、、1、lmd的搜索匹配在n个数中选若干个数，使其和等于某个数。（n据他自己说，数据范围和有特点、、、、正解的数据范围是暴力的两倍、、、、结果几乎没人看出来。容易看出

在n*n平面上有n*n-m-1（除原点与另m个点），可以从原点引k条线，询问最多能穿过多少点。设法最简表示每种方案（i,j）->i,j互质，容易想到最多有（phi[2]+phi[3]+...+phi[n]）*2]种方案，而每种方案能射掉的点为n/max（i,j），接下来便好统计。#include #include #include struct point{ int x,y

神奇的口袋，考数学甚于考算法，【考虑无自由取球情况，那么，设A(X,Y)为第x次，取出颜色为Y的球，P{A(1,Y1),A(2，Y2)…}=P{A(1,Y1),A(2,Y1),A(3,Y1)}因为首先分母是确定的，这毋庸置疑的，而分母取值只与在它之前取过多少同颜色球有关，与其他颜色球无关，所以可以进行此交换，现在回到可自由取球情况，举个例子也就是P{A(1，K),A(2，Y1),A(3，Y2)}=

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

秋哥是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优

计算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）=df[a]=-a+2a*sigma(d|a)sigma(b令b/d=k,f[a]=-a+2a*sigma(d|a)sigma(

本来想一次性大批量更新的，每次都发现东西太多，又懒得写了干脆改成休息的时候，写博文放松算了，好吧，本博客处于半荒废状态梭哈游戏，参见ctsc2012d1t1（ctsc太颓了，酱油记都不想写）其实这道题，不是容斥恶心，而是牌局的大小判断恶心。分别搜出两人的牌局，然后统一排序，那么我们只需知道一个牌局前面有多少没有相同重复手牌的牌局即可，这部分就用容斥统计，然后比较大小各种恶心#i

题目描述： 有些代数式张得很像一张图，比如说 D*X % M。. 这些图就是传说中的数形图，有时候我们要求最小值就变成了最小数形图问题。 现在请求出满足 L=转化为L=0如果L/D!=R/D,那么l/d（+1）就是答案，如果L/D==R/D,那么L必与R处于相同循环节，我们考虑变换不等式D*X-R让X最小，也就是使Y最小，同时除以D，得（-R）%D很像拓展欧

题意：进行翻硬币实验，若k次向上则结束，进行第n次实验需花费2*n-1的费用，询问期望结束次数及期望结束费用设f[i]为结束概率f[i]=c(i-1,k-1)*p^k*(1-p)^(i-k)sigma(f[i])==1p^k*sigma(c(i-1,k-1)*(1-p)^(i-k))=1sigma(c(i-1,k-1)*(1-p)^(i-k))=1/(p^k)ans

我表示我完全是看题解莫名其妙对了，题目和代码暂存在这里。如有大牛路过，欢迎留评论【a^b≡a^(b mod phi(c)+phi(c))  （mod c）(b>=phi(c))】【公式的条件很苛刻，原以为（b>=phi(c)）没有意义】其中phi（c） 是欧拉函数原题目： 一、超级数【题目描述】：众所周知，多次加法运算是乘法，多次乘法运算是乘方，那么多次乘方运算呢？

1400：略#include #include #include #include #include #include const int dx[4]={-1,0,1,0};const int dy[4]={0,1,0,-1};const double oo=1e300;using namespace std;struct point{int x,y;};queue