![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
lethalboy
这个作者很懒,什么都没留下…
展开
-
bzoj1968: [Ahoi2005]COMMON 约数研究 ——by lethalboy
1968: [Ahoi2005]COMMON 约数研究Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 1964 Solved: 1485[Submit][Status][Discuss]DescriptionInput只有一行一个整数 N(0 Output只有一行输出,为整数M,即f(1)到f(N)的累加原创 2016-11-06 16:50:44 · 380 阅读 · 0 评论 -
bzoj 2445 最大团(阶乘取模+中国剩余定理CRT)
题意即求下式:m∑x|nn!(x!)nx(nx!)m∑x|nn!(x!)nx(nx!)m∑x|nn!(x!)nx(nx!)m∑x|nn!(x!)nx(nx!)根据欧拉定理: 对于互质的正整数 a 和 n ,有 aφ(n) ≡ 1 mod n 可知,求指数式子模1e9-402的值然后快速幂即可。O(sqrt(n))分解质因数求值,注意原创 2017-04-18 10:00:10 · 875 阅读 · 0 评论 -
bzoj2001 [Hnoi2010]City 城市建设 动态最小生成树
昨晚水冬令营课件看到这题,感觉蛮有意思的,学习了一波,抽象式理解,今天又看了大佬的代码,彻底弄懂了这个东西。WC2013顾昱洲在《浅谈一类分治算法》中提到了动态最小生成树的分治做法,我来梳理下我的理解。这个算法有两个重要的操作:①reduction:对于一张图,reduction操作的目的是删除一定不会出现在最小生成树中的边,以此减小图的规模流程:我们假设对于当前这张图,有原创 2017-04-13 17:31:40 · 2067 阅读 · 0 评论 -
简单数论知识梳理(省选复习)
(noip数论算法汇总)①扩展欧几里得int ex_gcd(int a,int b,int &x,int &y){ if(!b) { x=1,y=0; return a; } int g=ex_gcd(b,a%b,x,y); int t=x;x=y,y=t-a/b*y; return g;}应用及要点:one : 求形如ax+by=gcd(a,b)的一组解(原创 2017-03-23 21:14:01 · 1802 阅读 · 0 评论 -
bzoj1385: [Baltic2000]Division expression
1385: [Baltic2000]Division expressionTime Limit: 5 Sec Memory Limit: 64 MBSubmit: 255 Solved: 185[Submit][Status][Discuss]Description除法表达式有如下的形式: X1/X2/X3.../Xk 其中Xi是正整数且XiInput先给出原创 2017-02-22 20:57:03 · 407 阅读 · 0 评论 -
线性求逆元算法
摘自:http://blog.csdn.net/acdreamers/article/details/8220787其实有些题需要用到模的所有逆元,这里为奇质数。那么如果用快速幂求时间复杂度为,如果对于一个1000000级别的素数,这样做的时间复杂度是很高了。实际上有的算法,有一个递推式如下 它的推导过程如转载 2016-08-23 16:57:16 · 13508 阅读 · 4 评论 -
bzoj1013---------算法模板系列(gauss消元)
1013: [JSOI2008]球形空间产生器sphereTime Limit: 1 Sec Memory Limit: 162 MBSubmit: 4771 Solved: 2491[Submit][Status][Discuss]Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在,你被困在了这个n维球体中,你只知道球面上n+1原创 2016-12-23 15:36:10 · 408 阅读 · 0 评论 -
bzoj1770-------(算法模板系列之gauss消元异或方程组)
异或方程组就是形如这个样子的方程组:M[0][0]x[0]^M[0][1]x[1]^…^M[0][N-1]x[N-1]=B[0]M[1][0]x[0]^M[1][1]x[1]^…^M[1][N-1]x[N-1]=B[1]…M[N-1][0]x[0]^M[N-1][1]x[1]^…^M[N-1][N-1]x[N-1]=B[N-1]其中“^”表示异或(XOR, exclus原创 2016-12-23 21:01:27 · 2421 阅读 · 0 评论 -
bzoj4407 于神之怒加强版(莫比乌斯反演+线性筛)
4407: 于神之怒加强版Time Limit: 80 Sec Memory Limit: 512 MBSubmit: 355 Solved: 174[Submit][Status][Discuss]Description给下N,M,K.求Input输入有多组数据,输入数据的第一行两个正整数T,K,代表有T组数据,K的意义如上所示,下面第二行到第原创 2016-12-15 17:59:31 · 349 阅读 · 0 评论 -
2440: [中山市选2011]完全平方数 (容斥原理+莫比乌斯函数)
2440: [中山市选2011]完全平方数Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2984 Solved: 1432[Submit][Status][Discuss]Description小 X 自幼就很喜欢数。但奇怪的是,他十分讨厌完全平方数。他觉得这些数看起来很令人难受。由此,他也讨厌所有是完全平方数的正整数原创 2016-12-11 17:18:41 · 303 阅读 · 0 评论 -
算法随笔(扩欧)
定理1 gcd(a,b)是ax+by的线性组合的最小正整数,x,y∈z;定理2 如果ax+by=c,x,y∈z;则c%gcd==0;定理3 如果a,b是互质的正整数,c是整数,且方程ax+by=c(1)有一组整数解x0,y0则此方程的一切整数解可以表示为x=x0+bt;y=y0-at;t∈z;(2)原创 2016-11-11 22:51:50 · 450 阅读 · 0 评论 -
bzoj1041 圆上的整点(一种新奇的思路)
用了一种新奇的方法重温了这道题目。学弟发来的定理很妙呦orz:有上述定理,则问题转化为求r^2的 %4余1因子数 和 %4余3因子数。用约数个数定理,去掉偶因子后可以求出两者之和因此,求其中一种即可。从%4余1的因子数下手吧易知%4余1的因子定是由任意多个%4余1的质数和偶数个%4余3的质数相乘得来的任意多个%4余1的质数实际上就是最大的%4余1的质数的因子数,原创 2017-06-15 10:05:58 · 1620 阅读 · 1 评论