
数学
silentsaber~
晨意微寒秋渐深,侧畔无事俏佳人
展开
-
bzoj1007[HNOI2008]水平可见直线
还记得做过一道简单的题:给出很多直线(给定斜率和在y轴上的截距),求可见直线的数量。那道题比较简单,n^2的枚举就可以过,就是对每条直线维护一个区间[l,r],只要[l,r]不为空,这条直线就能被看见。这里差不多也是求可见直线这里数据是50000条,n^2是肯定不能过的.考虑把直线按k升序排列对于一条直线i-1和一条直线i,因为k是升序排列,它们交点原创 2016-03-02 09:23:06 · 419 阅读 · 0 评论 -
bzoj2820&&YY的GCD
马丹,又是一道权限题。问我怎么搞到题面:vjudge:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=37167HOMEPROBLEMSTATUSCONTESTAdd ContestStatisticLOGOUTnikelongUPDATE原创 2016-03-28 17:35:22 · 563 阅读 · 0 评论 -
bzoj3122 线性同余方程
http://www.lydsy.com/JudgeOnline/problem.php?id=3122看了一眼,以为是脑残题,可以先找出一个循环,然后再搞搞,发现找出循环代价太大,或者说有可能有些数多次出现了,但并没有循环。然后推公式吧:(卧槽,就几个月没碰数列就tm不会做了~ 然后搞出一个通项)xn= a^(n-1) x1 +b* (a^(n-1)-1)/(a原创 2016-03-28 16:04:10 · 424 阅读 · 0 评论 -
bzoj2466高斯消元求解XOR方程
http://www.lydsy.com/JudgeOnline/problem.php?id=2466不会做,暴力- - 所以T掉不解释(n正解是高斯消元。预备知识:矩阵乘法,行列式的基本变换(其实不需要,只是掌握了之后可以把消元的过程看成是行列式转成上三角的过程),XOR操作.对于一个不会高斯消元的蒟蒻,看网上的大神的题解都是各种被虐,所以慢慢补一下高斯原创 2016-03-26 11:35:15 · 1370 阅读 · 1 评论 -
codevs1482【路线统计】
1482 路线统计 时间限制: 1 s 空间限制: 256000 KB 题目等级 : 钻石 Diamond题解 查看运行结果题目描述 DescriptionN个节点的有向图, 求从start到finish刚好经过时间time的总方案数 mod 502630原创 2016-02-28 21:23:35 · 436 阅读 · 0 评论 -
codevs1281-线性递推
1281 Xn数列 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 大师 Master题解 查看运行结果题目描述 Description给你6个数,m, a, c, x0, n, gXn+1 = ( aXn + c ) mod m原创 2016-02-26 21:56:07 · 431 阅读 · 0 评论 -
POJ3233[矩阵连续幂之和]
http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=10959题意,给一个n*n的矩阵A,求A+A^2+A^3+A^4+A^5+A^6+```+A^k, 矩阵对m取mod.有种东西叫二分- -------------------------------------------------------原创 2016-02-26 16:17:13 · 469 阅读 · 0 评论 -
矩阵快速幂 -codevs斐波拉契数列2
1732 Fibonacci数列 2 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond题解 查看运行结果在“1250 Fibonacci数列”中,我们求出了第n个Fibonacci数列的值。但是1250中,n9。现在,你的任务原创 2016-02-26 11:38:20 · 367 阅读 · 0 评论 -
矩阵快速幂入门-斐波拉契数列
cdoevs1250一道简单的矩阵快速幂。求斐波拉契第n项modq, f0=f1=1;我对于矩阵还是初学状态。------------------------------------------------------------------------------------------------------------我们有这样一个矩阵:[ fn-1原创 2016-02-26 09:01:49 · 513 阅读 · 0 评论 -
gcdlcm[组合数学]
gcdlcm[组合数学]:描述:问有多少个 k 元组(a1,a2,a3,...,ak)(ai>=1)满足 gcd(a1,a2,...,ak)=D,并且 lcm(a1,a2,a3,...,ak)=L。输入:第一行有三个整数 k,D,L 分别为题面所描述的。输出:输出答案 mod(1e9+7)输入输出样例:gcdlcm.in gcdlcm.out2 1 6 4原创 2016-02-22 21:35:15 · 525 阅读 · 0 评论 -
bzoj1008[简单组合数学]
一看还以为很难的样纸~~~比较水了。还OI题- -要求任意两个相邻,貌似很难算,还要容斥什么的- -考虑设答案为集合A总共的状态m^n,为集合U,A关于U的补集为B.B就是任意两个不相邻。第一个数选有m种,剩下n-1个数每个数和之前那个不一样就好了,每一个有m-1种。B=m*(m-1)^(n-1)然后A=U-B=m^n-M*(m-1)^(n-1)原创 2016-03-07 15:02:27 · 440 阅读 · 0 评论 -
莫比乌斯3题小结 hdu4746,1695 && poj2154
我们约定,[x]表示x取下整hdu4746: n m求∑ ∑ f(gcd(a,b)) i=1j=1 min(n,m) min(n/d,m/d)我们可以看成是:∑ ∑ u(t)* [n/dt]*[m/dt] 其中d的因子个数原创 2016-03-29 17:21:23 · 450 阅读 · 0 评论 -
生成树计数
补充:度数矩阵,邻接矩阵,切尔霍夫矩阵http://www.cnblogs.com/noip/archive/2013/03/03/2941853.html论文:http://wenku.baidu.com/view/0c086741be1e650e52ea990e.html花个把小时看看,理解清楚原创 2016-03-02 17:24:05 · 481 阅读 · 0 评论 -
bzoj4407于神之怒
看到好多大神留下的题解能看得懂,但是对于那个线性筛我表示我是蒟蒻我不会QAQ~于是我决定来水一发(顺便留下这道权限题,什么时候有权限了再去交一(交易)发)目前找了个过了的代码对拍发现并无错误&& 测试时间大致和那份程序差不多(应该能过吧~~)题意: n m求 ∑ ∑ gcd(i,j)^k mod 1e9+7原创 2016-03-29 09:20:07 · 676 阅读 · 1 评论