2017年11月_君诺

原创关于数的容斥定理的代码实现

前述：其实容斥定理的用法并非自己感悟出来的，还是从大量博客里学习的，感觉这种深邃的思想我学不来，看题目的时候都想不到为什么要那样用，但是容斥定理本身是很简单的，公式我就不再展示了。但是还好因为集合的交集的求法比较困难，使得这方面能用固定代码解决的问题的种类大大减少。其中比较广泛的一个应用就是求1-m里面和n互质的数的个数。思路是用m-所有和n非互质的数的和，那么，和n非互质又有一个说法就是和n有大

2017-11-30 21:16:30 437

原创组合数的代码实现

前述：又是忙碌的一周，基本前三天都在看资料，大概是星期三晚上之前看完了老师发的所有课件，真的只是了解了知识点，对怎么码题一无所知，晚上的英语课去老师的博客里偷来了组合计数的几个博客资料，一直读到十点多，一想还有场练习赛，挑了个a的人最多的题目，花十分钟签了个到就休息了。当然并不是白休息，周四的凌晨0点爬起来看组合数专题的题目，把题干短的题目还有从博客资料里面见过的题目都仔细读了一遍，建立了第一想

2017-11-30 20:46:15 3539 1

原创网络流之最大流

前述：这段时间一直在忙数论，图论这块因为知识点都看得差不多了（按dalao的博客还差几小块知识，大头已经差不多了），周六写完作业想忙里偷闲弄会图论（复习完板子就该从博客做题入手了），结果晚上没弄完那一点板子，还是拖到了周末。不仅仅是因为有一些东西遗忘了，更是第二遍看的时候更需要把里面一些比较细致的东西看明白，然后弄懂，不能留下像第一遍看的时候的那样的漏洞。网络流的重要概念（按自己理解写的）

2017-11-26 19:25:35 359

原创 HDU-2588 GCD（欧拉函数）

题目：HDU-2588 GCD 题目大意：给定两个数n和m，求所有大于m的gcd(i,n)(i=1,2....n-1)的个数。解题思路：这个题睡前看的，睡着之前一直辗转反侧的想（其实也就一小会），数据量不小，暴力一定超时，欧拉函数是求与n互素的数的个数，也就是gcd(i,n)=1的这样的i的个数，突然不会，可是转念一想，gcd(a*c,b*c)=gcd(a,b),那么和n的最大公约数

2017-11-26 13:44:57 274

原创 HDU-3123 GCC(水题)

题目：HDU-3123 GCC 题目大意：给定一个数n,求（0!+1!+.....n!）mod m。解题思路：乍一看n挺大挺吓人，可是m很小，大于等于m的数的阶乘mod m就是0，这样一看这个题目就非常的水了。 AC代码：#include#include#include#include#include#define maxn 1000010using

2017-11-26 13:22:39 173

原创 POJ-3641 Pseudoprime numbers(快速幂）

题目：POJ-3641 Pseudoprime numbers 题目大意：给了base-a Pseudoprime numbers一个定义，就是首先某个数p不是素数，然后这个数对于一个数a,满足a^p mod p=a。现在给定两个数p,a让我们判断p是否为base-a Pseudoprime numbers。解题思路：按照定义判断即可，p比较大，所以用到了快速幂。 AC代码：

2017-11-26 13:08:22 197

原创 POJ-2689 Prime Distance（筛素数的平移)

题目：POJ-2689 Prime Distance 题目大意：给定一个区间[l,r]，让求区间当中相邻素数组中距离最小的一组素数以及距离最大的一组素数。解题思路：分析一下数据发现，虽然区间的右端点的最大值可以达到2^32左右，但是，区间的跨度确实极小的，应用这个，可以通过筛两次素数，第一次正常筛出2^16以内的素数，第二次对筛法进行区间平移，从区间的左端点开始筛除即可。这里注意左端

2017-11-26 12:55:43 253

原创 POJ-1995 Raising Modulo Numbers（快速幂简单应用）

题目：POJ-1995 Raising Modulo Numbers 题目大意：说实话，并不知道题目大意，看样例就能过的一道题，就是让求(A1^B1+....An^Bn)MOD M 解题思路：一看数据量不大，直接快速幂过之，第一发tle了，因为用了快乘，去了就a了，看起来不溢出的题还是最好别用快乘，不是一个很快的东西。算是一个快速幂的模板题了吧。 AC代码：#include

2017-11-26 12:43:45 186

原创 POJ-1811 Prime Test（Miller_Rabin算法和Pollard_rho算法模板题）

题目：POJ-1811 Prime Test 题目大意：非常简单的题意，就是让判断一个大数是否为素数，若不是，求出他的最小素因子。解题思路：Miller_Rabin算法和Pollard_rho算法模板题-算法详解 AC代码：#include#include#include#include#include#include

2017-11-26 12:10:34 292

原创 POJ-1730 Perfect Pth Powers(思维：大数分解素因子)

题目：POJ-1730 Perfect Pth Powers 题目大意：判断一个数n，若有x^p=n,求p的最大值。解题思路：这个题一开始瞄了一眼题解，发现pow还可以开根号，但是并不想怎么做，晚上睡觉的时候想起来可以用分解素因子的做法做，因为对数a,b,有a^p*b^p=(a*b)^p,这是其一，然后有a^(p1*p2)=（a^p1)^p2,综上所述，也就是找到数n的所有素因子a1

2017-11-26 12:03:11 324

原创 POJ-1284 Primitive Roots（原根的个数定理)

题目:POJ-1284 Primitive Roots 题目大意:给定一个奇素数，求他的原根的个数，其中，一个数n的原根是指一个小于n的正整数x，若有x^1mod n,x^2 mod n...x^n-1 mod n刚好是n-1的一个全排列，那么

2017-11-26 11:49:50 518

原创 POJ - 1006 Biorhythms（CRT中国剩余定理）

前述：对做过的题目进行复习，顺便整合一下专题里面比较经典的数论题目。 POJ - 1006 Biorhythms 这是一道CRT的模板题，也是最经典的例题，题意就是人的体力，情感，智力分别有高峰周期m1=23,m2=28,m3=33,给定一个三者达到高峰的初始日期（从一年的第一天开始的日期），问重新达到高峰的日期与给定日期的差值，若求得的数小给定日期，则加上21252。因为三个周期满足

2017-11-26 11:16:25 252

原创一元线性同余方程&&一元线性同余方程组

概述：求解一元线性同余方程以及一元线性同余方程组是拓展欧几里得算法的一个重要应用，以下是对问题所包含的几种情况向着拓展欧几里得算法转化的方法。一元线性同余方程(形如：a*x = b(mod m）)：这种方程的求法可以转化成为一个不定方程：a*x-m*y=b

2017-11-21 23:03:03 637

原创欧几里得算法&&拓展欧几里得算法

前述：欧几里得算法是求算两个数的最小公倍数的一个非常有效的算法，而由其经过推理改造得到的拓展欧几里得算法更是一个应用广泛的算法，由于我比较菜，做的题目很少，对这块知识点的理解一直不是特别的深刻，完全停留在有这样一个知识点，大概知道他的推理过程，以及有怎么个模板而已。于是遇到相关的变式题目屡屡受挫，于是今天花了一定时间，就算法的推导以及应用方面自己通了一遍，以下是相关总结。欧几

2017-11-21 22:22:26 297

原创素数相关

前述：有关于数论这一部分也看了不少了，正好今天去看病，路上总结一下有关于素数的一些相关知识，汇总一下我所见过的素数相关判别法，当然有可能不怎么官方的说法，按照自己的理解写的。素数的定义：对于一个数n，若n除了1和它本身没有其他因子，那么则称n是一个素数。素数的判别法： 1.定义法：这也就是为什么要再提一遍这个简单定义的原因，定义在数学中是极为重要的，所有的定理都是由公理和定义

2017-11-19 10:45:03 388

原创 Uva1434 YAPTCHA（威尔逊定理应用）

题目链接：http://uva.onlinejudge.org/index.phpoption=com_onlinejudge&Itemid=8&category=516&page=show_problem&problem=4180 这个题目就是让计算的值，分析一下可以看的出来，就是让计算1/3k+7加上（3k+6)!/3k+7的小数部分，是一个威尔逊定理的应用的题目。威尔逊定理：当

2017-11-17 22:16:44 917

原创 Miller_Rabin算法&&Pollard-Rho算法

前述：近期在uva上刷到一个题，prime test，题意非常简单，就是让判断一个2^54以内的数是否为素数，年轻的我满脑子都是相关的暴力加优化，做之前我还去百度搜了一下判断一个数是否为素数的最佳方法，博主给的是素数的线性筛法，很自信的敲完代码，然后很自信的tle。做不了了，搜了一下题解，原来并不是一道简单题，涉及到本文提及的两个算法，Miller_Rabin算法&&Pollard-Rho算法，当

2017-11-17 21:32:17 295

原创快速幂&&快速计算（a*b）mod m

概述：起初在dalao的博客里看到这个（a*b）mod m就觉得很不解，为什么（a*b）对m求余还需要用一个函数来实现，非常的不解。想想之后才发现自己真的蠢，要是m很大的话，虽然有a,b 简要思路：。由于二者思路类似，所以这里只说快速幂，记得以前知道这样一个概念，用1，2，4，8...2^n的线性组合可以表示2^n以内任意的数，快速幂正是利用了这个原理。这是基于二进制的一个原理，比如11的二

2017-11-17 19:09:50 1981

原创威尔逊定理

概述：威尔逊定理是判断一个数p是否是素数的又一个充分必要条件，但由于这个定理计算起来并不方便，所以并不作为判别素数的一个方式，但是却有着它独道的用处。定理内容：当且仅当p为素数时：( p -1 )! ≡ -1 ( mod p )。简要证明（摘自百度）：充分性如果“p”不是素数，那么它的正因数必然包含在整数1, 2, 3, 4, … ,p− 1 中，因此

2017-11-17 18:30:04 970

原创 2017开学训练第十二周周中总结

这周周一基本就是满课加补作业的状态，抄大物实验报告抄到手爆炸的那种。周二开完会回去之后开始做题，除去一道做过的CRT的题目加一道原根性质的题（忘了有这个性质了，差了题解之后更是一点不懂他的证明啥意思，暂且先记下来结论，就是euler（n-1））然后就是四道简单的思维题,G题据猜测是一个等比数列，后来发现只有素数满足这个性质，改了之后顺利ac。Q题有点类似于以前那个求2的3的4的5的次方次方

2017-11-16 22:09:27 206

原创 2017开学训练第十一周周末总结

这周周五回来之后去洗澡，然后回来补了补作业，写完常微分发现已然快断电了，又打开那个M题看了看，总感觉思路没错，可是为什么错始终想不出来，但是自己想了怎么久又不舍得放弃。周六补作业加继续想题，直到吃完午饭发现有人a了m题，连忙问了问怎么a的，不出所料是看的题解...但总归是拓展欧几里得算法，发现思路和我除了欧几里得模板部分其他完全不一样，我也懒得看他的想法，用他的代码试了几组数据，和我的完全

2017-11-12 17:55:34 246

原创 2017开学训练第十一周周中总结

这周周一满课，晚上回来写了写作业就到断电了，拿着手机当手电筒复习了一下周末看的图论二分图的基础知识，主要是定义方面的。二分图的典型问题就是结婚找对象这类的分配问题，简而言之就是两组之间配对，给了一组里面的想要的对象的范围，让尽可能多的人满意，就是二分图的匹配，要是都能心仪，那就是完美二分图。其中比较重要的一个概念就是增广路径（一个特殊的交错路径），匹配M的交错路径顾名思义就是交替出现在匹配M中和M

2017-11-09 21:59:56 253

原创 2017开学训练第十周周末总结

这周周末并没有什么很大的进展，想在去看病的路上看点东西，结果不能深入了去想，想的东西多了就有些晕，也就只看了些有关二分图的概念性的问题，回想了那个始终做不出来的M题，感觉要么就是题意都错了，怎么看代码不存在解决问题方面的错误，奇葩数据也测试过了，可就是wr。看来还需要再深入读题，建了个群发现还是好冷清，大家都不愿意交流有关做题方面的心得，要么埋头做自己的，要么等着我发出来思路，在凑合着写写，还是得

2017-11-05 21:05:51 291

原创 2017开学训练第十周周中总结

这周前几天主要在做数论的题目，因为周末要出去看病，还不知道折腾成什么样子，所以周末就只安排看图论的知识点了，这几天都没怎么碰图论，不知所云老师这周一直在加作业，主要她的两门课都不怎么能听得懂，作业就变得吃力的很。怎么说呢，越往后做越吃力了，其实还是差在读题上，百度翻译再也不敢开了，一到比赛就知道不熟练是多么的惨了，加上有个不懂英语的，长题都不敢动一下，只好自己慢慢加强了。J题，一看数据量怎

2017-11-02 16:56:09 299

m0_37772713的博客