自以为心若顽石，却终究人非草木。

转载 中国剩余定理 证明

中国剩余定理问题解同余方程组 其中m1,m2,m3...mk为两两互质的整数求x的最小非负整数解定理令即M是所有 mi 的最小公倍数ti为同余方程的最小非负整数解则有一个解为通解为特别的，最小非负整数解为（x m...

原创 初等数论四大定理

威尔逊定理 当且仅当 p 为素数时：（p-1)！≡ -1 （mod p）或者这么写：（p-1）！≡ p-1 （mod p）或者说若 p 为素数，则 p | （p-1）!+1 欧拉定理 若n,a为正整数，且n,a互质，即gcd(a,n) = 1，则a^φ(n) ≡ 1 (mod n)欧拉定理拓展：（A,C不互质）A^B mod C = A^（B mod φ...

转载 欧拉定理 证明及推论

欧拉定理： 若正整数a , n互质，则aφ(n)≡1(mod n)其中φ(n)是欧拉函数（1~n)与n互质的数。证明如下： 不妨设X1，X2...... Xφn是1~n与n互质的数。　　首先我们先来考虑一些数：aX1，aX2...... aXφn　　这些数有如下两个性质：（1）任意两个数模n余数一定不同：（反证） ...

原创 费马小定理 （证明）

费马小定理费马小定理：a^(p-1) ≡ 1(mod p)前提：p为质数，且a，p互质 //互质：a和p相同的因数为1.先来看一下≡是什么：a ≡ b (mod p) <=> a mod p = b mod p //a和b在模p的意义下同余注释：<=> 两边相等在证明之前，先给出引理：(1)：如果p,c互质,并且a*c≡b*c(m...

原创 浅谈 逆元&乘法逆元

逆元除以一个数再取模 = 乘以这个数的逆元再取模设 inv[b] 是b的逆元，则（a/b） mod p = （a * inv[b]） mod p（1）一个数的逆元是什么 一个数x在模p的条件下的逆元是多少（2）一个数的逆元有无穷个，只求最小正整数即可（3）一个数x在模p的条件下不一定有逆元，x关于p的逆元存在当且仅当x和p互质推导：设 a 为 x 的逆元，b 为任意整数...

原创 浅谈拓展欧几里得算法（辗转相除法）

前言前天讲了一堆的数论，还没什么时间整理，拓欧我也只是第一次接触（果然还是太蒟），搞了半个晚上才弄懂。避免遗忘或弄丢，写篇博文。拓展欧几里得算法了解一下：扩展欧几里得算法，简称 exgcd，一般用来求解不定方程，求解线性同余方程，求解模的逆元等引理：存在 x , y 使得 gcd(a,b)=ax+by，求 x，y证明：part 1：（1）由gcd（a，b）= gc...

原创 浅谈欧几里得算法（辗转相除法）证明

前言之前只知道一味的背辗转相除法的biao，连它是欧几里得算法都不知道，所以决定写一篇证明百度百科了解一下欧几里德算法又称辗转相除法，是指用于计算两个正整数a，b的最大公约数。应用领域有数学和计算机两个方面。计算公式gcd(a,b) = gcd(b,a mod b)。//这就是我接下来要证明的东东证明假设1.求n、m的最小公倍数2、a、b分别是 m div n 和 ...

原创 筛素数方法浅谈

前言很久之前学了欧拉筛，可现在都差不多忘了，希望写一篇博文记录下来，也希望对各位初学线性筛的oier有一定的帮助。筛素数普遍四大方法1.二重循环暴力找素数 O(n^2)2.在第一种的基础下，只判断到 O(n·sqrt（n）)3.埃氏筛法 ...

转载 浅谈树状数组

浅谈树状数组之前也看过了好多关于树状数组的博客，结合这几天做的题，我一直想好好总结一下树状数组，这篇文章就来浅谈一下树状数组。1.前言首先我们要明白树状数组是一种数据结构，利用树状数组可以以空间换取时间，这一点和之前的线段树一样，但是树状数组访问会更快，效率更高，树状数组不同于线段数的一点就是这棵树的构成。二叉树或者线段树是这样的：而树状数组是这样的：...

转载 数据结构Data structures

堆Heap•二叉堆是一个基础的数据结构•它满足父结点的键值总是大于等于（或小于等于）子节点的键值。•堆的可以支持插入，删除或查找最大（或最小）键值的点。•且操作的时间复杂度为log级别。•基础例题：合并果子简单应用：计算哈夫曼树、数据维护……NOI2015 荷马史诗•n<=100000,k<=9HELP•先考虑k=2怎么做•题目大意，给你...

转载 莫比乌斯反演

莫比乌斯反演（PS：在评论区中众多dalao的催促下，我认真的写了三天三夜写完了这篇杜教筛，保证是精品！）前言（这大概是我第一次写学习笔记吧OvO）可能每一个刚开始接触莫比乌斯反演的OIer，起初都会厌恶这个神奇的东西。(我也一样233)每一个人厌恶的原因有许多，可能是这个烦人的式子，也可能仅仅只是因为不理解μμ函数而感到不爽。当然，莫比乌斯反演有一个小小的预备知识：整除分块那么我...

转载 杜教筛

杜教筛（似乎有很多人在催我的杜教筛呢......）前言话说，我是不是在自己的莫比乌斯反演中挖了许多杜教筛的坑啊...... 本文完整的总结介绍杜教筛，也算是将莫比乌斯反演中的坑全部填满吧！ 真诚地希望来阅读这篇学习笔记的每一个人，仔仔细细的看完每一段。 我相信，只要认真的看完整篇文章并跟着一起思考的读者，一定能够有所收获！ 如果您之前不会杜教筛，那么我希望这篇文章能够作为您学习...

原创 jzoj 3517. 【NOIP2013模拟11.6B组】空间航行(warp)

Description你是一艘战列巡洋舰的引擎操作人员，这艘船的船员在空间中侦测到了一些无法辨识的异常信号。你的指挥官给你下达了命令，让你制定航线，驾驶战列巡洋舰到达那里。船上老旧的曲速引擎的速度是0.1AU/s。然而，在太空中分布着许多殖民星域，这些星域可以被看成一个球。在星域的内部，你可以在任何地方任意次跳跃到星域内部的任意一个点，不花费任何时间。你希望算出到达终点的最短时间。In...

原创 jzoj 3515. 【NOIP2013模拟11.6B组】软件公司(company)

Description一家软件开发公司有两个项目，并且这两个项目都由相同数量的m个子项目组成，对于同一个项目，每个子项目都是相互独立且工作量相当的，并且一个项目必须在m个子项目全部完成后才算整个项目完成。这家公司有n名程序员分配给这两个项目，每个子项目必须由一名程序员一次完成，多名程序员可以同时做同一个项目中的不同子项目。求最小的时间T使得公司能在T时间内完成两个项目。Input第...

原创 jzoj 1435. 开花

Description　　在遥远的火星上，上面的植物非常奇怪，都是长方形的，每个植物用三个数来描述：左边界L、右边界R以及高度H，如下图所示描述一个植物：L=2,R=5和H=4。　　　　每天都有一个新植物长出来，第一天的植物高度为1，后面每天长出的植物比前一天的高1。　　当一个新植物长出来的时候，跟其他植物的水平线段相交处会长出一朵小花（前提是之前没有长出花朵），如果线段交于端点，是不会...

原创 jzoj 1434. 灌水

Description　　学生都很喜欢灌水，第一天只有Alice给她的每个朋友灌了一次水，从第二天开始，所有学生（包括Alice）将会有规律地去灌水：　　•如果前一天被灌了奇数次的水，他们将会给每个朋友灌一次水；　　•如果前一天被灌了偶数次的水，他们将会给每个朋友灌两次水。　　学生编号为1到N，Alice为1号，学生之间的朋友关系会给出。　　计算H天后一共灌了几次水。Input...

原创 jzoj 1433. 数码问题

Description　　Alice有一个N*N的格子，把1-N^2按照从上到下从左到右的顺序填进表格中，允许在表格上进行两种操作：　　(1) 旋转行——这一行的数向右移动一个位置，而最后一列的数会移到第一列；　　(2) 旋转列——这一列的数向下移动一个位置，最后一行的数会移到第一行。　　Alice想把数X移到（R,C）处可以采用以下方法：　　•如果X不在C这一列，通过旋转行操作把X移...

原创 jzoj 1432. 输油管道

Description　　请你帮忙设计一个从城市M到城市Z的输油管道，现在已经把整个区域划分为R行C列，每个单元格可能是空的也可能是以下7种基本管道之一：　　　　油从城市M流向Z，‘+’型管道比较特殊，因为石油必须在两个方向（垂直和水平）上传输，如下图所示：　　　　　 　　　　　现在恐怖分子弄到了输油管道的设计图，并把其中一个单元格中的管道偷走了，请你帮忙找...

原创 jzoj 3403. 【NOIP2013模拟】数列变换

Description小X 看到堆成山的数列作业十分头疼，希望聪明的你来帮帮他。考虑数列A=[A1,A2,...,An],定义变换f(A,k)=[A2,A3,,,,.Ak,A1,Ak+2,Ak+3,,,,A2k,Ak+1,...]，也就是把a 分段，每段k 个（最后如果不足k 个，全部分到新的一段里，见样例），然后将每段的第一个移动到该段的最后一个。现在，小 X想知道 f (f (f (f ...

原创 一只蒟蒻，大家多多关照

原创 jzoj 1308. 取数游戏

Description　　Alice想让Bob陪他去看《唐山大地震》，但由于Bob是个很感性的人，怕流泪不想去，但又不好意思以这个作为拒绝的理由，便提出玩一个游戏。　　N个正整数围成一圈，规则如下：　　•两个玩家轮流取数；　　•最开始先手的玩家可以取任意一个数x；　　•从第二步开始当前玩家只能取x（上一玩家刚刚取的数）左右两边相邻的数；　　•直到取完所有的数，游戏结束；　　•取得较...

原创 jzoj 1382. 区间

Description　　Alice收到一些很特别的生日礼物：区间。即使很无聊，Alice还是能想出关于区间的很多游戏，其中一个是，Alice从中选出最长的不同区间的序列，其中满足每个区间必须在礼物中，另序列中每个区间必须包含下一个区间。　　编程计算最长序列的长度。Input　　输入文件第一行包含一个整数N(1<=N<=100000)，表示区间的个数。　　接下来N行，每...

原创 jzoj 1381. 删除

Description　　Alice上化学课时又分心了，他首先画了一个3行N列的表格，然后把数字1到N填入表格的第一行，保证每个数只出现一次，另外两行他也填入数字1到N，但不限制每个数字的出现次数。　　Alice现在想删除若干列使得每一行排完序后完全一样，编程计算最少需要删除多少列。Input　　第一行包含一个整数N(1<=N<=100000)，表示表格的列数。　　接下...

原创 jzoj 1420. 佳肴

Description　　佳肴就是非常美味的菜的意思，佳肴最关键的是选择好原料。　　现在有N种原料，每种原料都有酸度S和苦度B两个属性，当选择多种原料时，总酸度为每种原料的酸度之积，总苦度为每种原料的苦度之和。　　正如大家所知，佳肴是既不酸也不苦的，因为要保证所选的原料使得总酸度和总苦度差的绝对值最小。　　由于佳肴不能只有水，所以必须至少选择一种佳肴。Input　　输入第一行包含...

原创 jzoj 2643. 【NOIP2011模拟10.31】六边形

Description棋盘是由许多个六边形构成的，共有5种不同的六边形编号为1到5，棋盘的生成规则如下：1.从中心的一个六边形开始，逆时针向外生成一个个六边形。2.对于刚生成的一个六边形，我们要确定它的种类，它的种类必须满足与已生成的相邻的六边形不同。3.如果有多个种类可以选，我们选择出现次数最少的种类。4.情况3下还有多个种类可以选，我们选择数字编号最小的。现在要你求第N个生成的六...

原创 jzoj 2644. 【NOIP2011模拟10.31】数列

Description给你一个长度为N的正整数序列，如果一个连续的子序列，子序列的和能够被K整除，那么就视此子序列合法，求原序列包括多少个合法的连续子序列？对于一个长度为8的序列,K=4的情况：2, 1, 2, 1, 1, 2, 1, 2 。它的答案为6，子序列是位置1->位置8,2->4,2->7,3->5,4->6,5->7。Input第一行：T...

原创 jzoj2642. 【NOIP2011模拟10.31】游戏

DescriptionAlice和Bob在玩一个游戏，游戏是在一个N*N的矩阵上进行的，每个格子上都有一个正整数。当轮到Alice/Bob时，他/她可以选择最后一列或最后一行，并将其删除，但必须保证选择的这一行或这一列所有数的和为偶数。如果他/她不能删除最后一行或最后一列，那么他/她就输了。两人都用最优策略来玩游戏，Alice先手，问Alice是否可以必胜？Input第一行：T，表示数...

原创 jzoj 2019.7.12【NOIP2012模拟10.26】比赛总结

https://blog.csdn.net/qq_42887171/article/details/95619073火炬手https://blog.csdn.net/qq_42887171/article/details/95625219雕塑前言有了昨天的教训，晚上再也不敢熬夜了，乖乖的早睡晚起（也就7:00起而已）今天精神倍儿棒<(￣︶￣)>今天的题终于又回到了正...

原创 jzoj 3059. 【NOIP2012模拟10.26】雕塑

DescriptionWcyz为了迎接百年校庆，美化校园，请了校友笨笨将n座雕塑，准备安置在校园内，整个校园可以抽象成一个n*n的大网格，每个1*1网格最多只能安置一座雕塑，但是某些1*1的网格上恰好是一个食堂或湖泊，这些网格是不能安置雕塑的，每个雕塑的造型相同，这样同一种安置方案中交换排列都算一种。任意雕塑在同一行或同一列是不合法的方案。学校想知道有多少种安置方案，笨笨想从中选择最好的一...

原创 jzoj 3057. 【NOIP2012模拟10.26】电影票

Description笨笨当了很久的道路调度员，笨笨也开始想体验生活，从生活中发现数学问题，锻炼自己思维。最近《变形金刚3》，《哈利波特7》同步放映，明显是决战雌雄，已知王府井中一共有n人买了《变形金刚3》的票，m人买了《哈利波特7》的票，并且n>=m，并且电影院中现在只有两种票，每次只有一个人买，(共有n+m次)，这n+m次组成一个排列，为了保证每一个人买票时，《变形金刚3》票...

原创 jzoj 3058. 【NOIP2012模拟10.26】火炬手

Description全运会就要开始了，笨笨想成为湖南地区的火炬手，经过层层选拔，最终到了最后一关，这一关给出了一个正整数n（N<=100000),求一个最小的正整数m，使得n*m的十进制表示形式中只含1和0.Input一行一个数nOutput输出一行，如果有解，输出最小的m，否则输出‘no solution’.Sample Input...

原创 jzoj 3235. 数字八

DescriptionFarmerJohn的奶牛最近收到一块大理石。但不幸的是，这块石头有些不完整。为了说明这块石头的状况，我们就可以用一个N*N正方形网格（5<=N<=300）来描述，其中字符'*'代表石头的缺损部分，'.'表示石头完美无瑕的部分。奶牛要在这一块大理石上雕刻数字“8”。然而，牛也需要FJ的帮助，以确定在这块大理石上最佳的雕刻位置。这里有几个要求来定义一个...

原创 jzoj 3234. 阴阳

DescriptionFarmerJohn正在在计划自己的农场漫步。他的农场的结构就像一棵树：农场有N个谷仓（1<=N<=100,000），分别由N-1条路链接。这样，他便可以通过这些谷仓间的道路遍及各个谷仓。FarmerJohn想要选择一条路线：这条路线的起点和终点分别为农场中两个不同的谷仓，这条路线不能重复经过一条边两次。FarmerJohn担心这条路径可能会偏长，...

原创 2019.07.11比赛总结

前言昨晚是真的没睡好，今天的比赛貌似都是睡过去的，呵呵。。。前两天的题都太淼了(*^▽^*)（可能是XC想给我们放松一下），今天直接扔了个 省选- 的比赛过来，把我们全轰了(*Φ皿Φ*)Farmer John此生最大的敌人！！！T1（照片）赛时：刚看到题，jio得还挺简单的，于是一波疯狂打码，结果打着打着，发现事情并不简单（嘿嘿嘿）于是开始思考正解，但是。。。一头栽倒...

原创 jzoj 3233. 照片

DescriptionFarmerJohn决定为他的N头排列好的奶牛（1<=N<=200,000）做一张全景合照。这N头奶牛分别以1..N进行编号。他一共拍了M（1<=M<=100,000）张相片，每张相片都只包含有一部分位置连续的奶牛：第i张照片涵盖着编号从a_i到b_i的所有奶牛。当然，这些照片合起来并不保证包含所有的牛。FarmerJohn拍摄完所...

原创 jzoj 1782. Travel

Description　　给出一个有n个顶点m条边的有向图，对于一条边长度为len的边有两种走法。　　1、如果a和b可以互达，则走过这条边的时间为len　　2、如果a和b不可以互达，则走过这条边的时间为2*len　　现在给出一个k，问，从顶点1到顶点n，满足第二种走法不超过k次的最短时间是多少。Input　　第一行有3个整数n,m,k(1<=n<=100,1<=...

原创 jzoj 1765. Circle

Description　　给定三个点（不共线）的坐标，要求以这三个点为圆心做三个圆，圆两两不相交，不包含，问三个圆的直径和最大为多少。Input　　第一行为测试数据组数t。接下来t行每行6个数为坐标。Output　　T行，每行一个整数表示直径和（取下整）。Sample Input10 0 0 1 1 0Sample Output3Data Const...

原创 jzoj 1781. Number

Description　　给出一个整数 ，你可以对 进行两种操作。　　1、将x变成4x+3　　2、将x变成8x+7　　问，最少通过多少次操作，使得x是1000000007的倍数？Input　　一行，一个整数x(1<=x<=1000000006)。Output　　一行，表示最少的操作步数。保证答案不超过10^5。Sample Input / Sample O...

原创 jzoj 1763. Gift

Description无Input　　输入的第一行为一个整数t。　　接下来t行，每行包含九个自然数。Output　　输出t行　　每行一个整数，表示2^a+2^b+2^c+2^d+2^e+2^f+2^g+2^h+i。Sample Input121 30 0 0 0 0 0 0 2147483647Sample Output3223322629Dat...

原创 jzoj2019.7.9比赛总结 NOIP2010四校联考模拟五

前言今天的状态还不错，一眼全是淼题，想着能拿个大几百，结果╥﹏╥…T1（Gift）赛时：刚看到题目时瞬间炸了，没有题目描述？？？（一脸懵逼*1）再往下看，好吧，题目全在输入输出里边（吓死宝宝了 ）截图给你们看呵呵。。。老师教我们要看数据范围，一瞟，妈呀！100% t<=100000 a,b,c,d,e,f,g,h<=60 i<=92233720368547758...