数论
文章平均质量分 55
数论
Stockholm_Sun
这个作者很懒,什么都没留下…
展开
-
编年史
编年史(chronicle.pas/c/cpp)题目描述公元1582 年罗马天主教教宗格里高利十三世把儒略历1582 年10 月4 日的第二日定为格里高利历1582 年10 月15 日,同时修改了儒略历的置闰法则。这就是格里高利历法,又称为西历,是现在的国际通用历法。现在 applepi 手上有一本十分古老的编年史,这本史书记录了很多著名的历史事件。于是applepi 有了一个奇怪的想法……他想知道原创 2017-03-25 21:44:14 · 503 阅读 · 0 评论 -
因子和(数论,算数基本定理,快速幂,逆元)
求逆元的方法,将在接下来讲到。代码#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<cmath>using namespace std;long long n,i,j,m,k;long long x,t;long long r(){ long long p=0,f=1;原创 2017-10-28 21:49:33 · 518 阅读 · 0 评论 -
P1114 “非常男女”计划(数论,线性结构)
P1114 “非常男女”计划(数论,线性结构) 题目描述 近来,初一年的XXX小朋友致力于研究班上同学的配对问题(别想太多,仅是舞伴),通过各种推理和实验,他掌握了大量的实战经验。例如,据他观察,身高相近的人似乎比较合得来。万圣节来临之际,XXX准备在学校策划一次大型的“非常男女”配对活动。对于这次活动的参与者,XXX有自己独特的选择方式。他希望能选择男女人数相等且身高都很接近的一些人。这种选择原创 2017-10-13 10:34:41 · 817 阅读 · 0 评论 -
P1170 兔八哥与猎人(简单数论,GCD)
P1170 兔八哥与猎人今天来看一道水题。题目描述兔八哥躲藏在树林旁边的果园里。果园有M × N棵树,组成一个M行N列的矩阵,水平或垂直相邻的两棵树的距离为1。兔八哥在一棵果树下。猎人背着猎枪走进了果园,他爬上一棵果树,准备杀死兔八哥。如果猎人与兔八哥之间没有其它的果树,猎人就可以看到兔八哥。现己知猎人和兔八哥的位置,编写程序判断兔子所在的位置是否安全.输入输出格式输入格式: 第一行为n,表示有n原创 2017-10-13 10:42:57 · 1103 阅读 · 0 评论 -
快速幂和慢速乘
今天看看数论方面的事情。快速幂 a^b mod c递推式一般我个人比较推荐递推式,分两种情况,一个是b为奇数,另一个是b为非0偶数,弹出条件为b==0。int power(int a,int b,int p){ if(b==0) return 1%p;//最好写上%p,更加严谨。 int c=power(a,b/2,p); c=c*c%p; if(b&1) c=c*原创 2017-10-05 19:12:32 · 797 阅读 · 0 评论 -
NOIP2016组合数问题
NOIP2016组合数问题原创 2017-03-19 10:50:36 · 2167 阅读 · 0 评论 -
[TJOI2010]中位数(优先队列,堆)
[TJOI2010]中位数(优先队列,堆)题目概述: 给定长度为N的初始序列和M个操作,操作有两种:1.插入一个元素;2.查询当前序列的中位数。 数据规模: N<=100000,M<=50000 思路: 这个题比较自然的思路就是先把初始队列排序后对半分开,然后用一个大根堆存1~N/2号元素,小根堆存N/2+1~N号元素,这样大根堆堆顶元素和小根堆堆顶元素就是当前的中位数的备选元素。然后执行原创 2017-08-24 16:48:49 · 810 阅读 · 0 评论 -
CQOI2007 余数求和
P2261 [CQOI2007]余数求和题目概述:给出正整数n和k,计算G(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。 数据规模:n,k <= 10^9 这个题暴力算法只能过30%的点。 思路:显然,这是一道数论题,要注意数据规模虽然在int范围内,但是不能保证ans在int范围内,所以要开长整型格原创 2017-08-09 17:39:14 · 445 阅读 · 0 评论 -
[SDOI2012]Longge的问题(数论,欧拉函数)
[SDOI2012]Longge的问题P2303题目描述 Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。输入输出格式 输入格式: 一个整数,为N。 输出格式: 一个整数,为所求的答案。输入输出样例 输入样例#1: 6 输出样例#1: 15说明 对于60%的数据,0 < N <原创 2017-11-04 21:22:27 · 539 阅读 · 0 评论 -
泰勒公式
泰勒公式O.绪言高考之后闲来无事,又搞了几天数学,感觉可以总结总结泰勒公式这块儿的内容。然后发现泰勒公式和二项式组合数学上的联系,所以说一说。 这篇博文适合学完导数的同学食用。I.神奇的泰勒公式先把公式挂上来, 我们先不管Rn(x)是怎么回事,只看前面这一部分的话就可以。然后我们就想套公式吧,举个栗子吧,我们令函数f(x)=e^x,式子中的a是一个已知函数值及其各阶...原创 2018-07-13 11:49:26 · 7895 阅读 · 0 评论 -
凸包入门
凸包入门凸包是信息学奥赛的计算几何里面一个比较重要的一块儿。今天写一下一些有关于凸包的入门内容。什么是凸包顾名思义,听见凸包的名字,大家大概对它的形状有一定的认识了。而凸包里面“包”的是什么呢?是一堆点。 就像这样。 就是一个凸凸的包把一堆点全包住啦 值得一提的是,凸包的各个顶点都是这些点里面的几个。 怎么解凸包呢今天主要讲讲葛立恒(Ronald Graham,...原创 2018-08-09 23:01:53 · 964 阅读 · 0 评论 -
公约数的和(数论,欧拉函数)
接下来我们来看看几道数论的题目。 代码#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<cmath>using namespace std;long long n,i,j,m,k;long long f[2000001];long long r(){ long lo原创 2017-10-28 21:15:23 · 572 阅读 · 0 评论 -
斐波那契数列(数论+矩阵加速)
斐波那契数列输入一个long long范围的一个数x,要求输出斐波那契数列第x项。思路显然这个题目,用O(N)的朴素做法是肯定不行的。所以我们要寻求一种新的方法。 接下来介绍一种计算:矩阵加法。这样的一种运算具有满足结合律的性质。对这个有兴趣的同学可以参阅人教版数学选修4-2的第一章和第二章(第一章带你入门,第二章教你矩阵乘法)。 大概矩阵加速就是和快速幂类似的,下面代码里面的操作可以教会我们其原创 2017-10-11 21:33:06 · 1115 阅读 · 0 评论 -
随机数生成器 (NOI 2012 第一试 第一题)
随机数生成器(NOI 2012 第一试第一题)题目描述栋栋最近迷上了随机算法,而随机数是生成随机算法的基础。栋栋准备使用线性同余法(Linear Congruential Method)来生成一个随机数列,这种方法需要设置四个非负整数参数m,a,c,X[0],按照下面的公式生成出一系列随机数{Xn}: X[n+1]=(aX[n]+c) mod m其中mod m表示前面的数除以m的余数。原创 2017-10-12 18:56:40 · 598 阅读 · 0 评论 -
三元上升子序列(二重顺序对与乘法原理的应用)
P1637 三元上升子序列 题目概述: 定义三元上升序列: THAIR{i,j,k}满足:i#include<iostream>#include<cstdio>#include<cstring>#define LL long longusing namespace std;LL i,j,k,m,n,temp;LL a[30001],b[30001],c[30001],ans;LL原创 2017-08-10 14:48:32 · 363 阅读 · 0 评论 -
联合权值(NOIP2014提高组 第一试 第二题)
联合权值(NOIP2014提高组 第一试 第二题)题目描述 无向连通图G 有n 个点,n - 1 条边。点从1 到n 依次编号,编号为 i 的点的权值为W i ,每条边的长度均为1 。图上两点( u , v ) 的距离定义为u 点到v 点的最短距离。对于图G 上的点对( u, v) ,若它们的距离为2 ,则它们之间会产生Wu×Wv 的联合权值。 请问图G 上所有可产生联合权值的有序点对中,联合权原创 2017-08-21 14:16:37 · 509 阅读 · 0 评论 -
无聊的函数(二分查找+数论)
无聊的函数 时空限制:1s,128MB 题目描述 定义一个函数: f(x)=x^x (k=1) f(x)=(x^x)^x (k=2) 其中 k 为给定的值(1 或 2)。 使得 f(x) 达到或超过 n 位数字的最小正整数 x 是多少? 此题由Luogu P2759 奇怪的函数 改装而来,题目难度上升了一点点。 思路: 对于函数f(x)的位数,可以用lg(x^x)表示,即:lg(原创 2017-08-15 10:44:26 · 368 阅读 · 0 评论 -
[USACO06FEB]数字三角形Backward Digit Su…
题目描述3 1 2 4 4 3 6 7 9 16有这么一个游戏: 写出一个1~N的排列a[i],然后每次将相邻两个数相加,构成新的序列,再对新序列进行这样的操作,显然每次构成的序列都比上一次的序列长度少1,直到只剩下一个数字位置。 现在想要倒着玩这样一个游戏,如果知道N,知道最后得到的数字的大小sum,请你求出最初序列a[i],为1~N的一个排列。原创 2017-08-23 10:38:11 · 453 阅读 · 0 评论 -
平方根(sqrt.pas/c/cpp)(数论)
平方根(sqrt.pas/c/cpp)【问题描述】 给出一个正整数 n (1<n≤2^31-1),求当 x,y 都为正整数时,方程 的解中, x 最小值为多 少? √n=√x-√y 【输入文件】 输入文件只有一行,一个正整数 n。 【输出文件】 输出文件只有一行,即满足条件的最小 x 的值。 【文件样例】 sqrt.in sqrt.out 4 9 【数据规模】 30%的数据满原创 2017-09-16 20:14:56 · 1111 阅读 · 0 评论 -
最长线段(几何证明题)
最长线段(chord.pas/chord.in/chord.out)(LYOI信息学综合模拟20090321Problem 1) 问题描述 给定两个圆各自的圆心坐标和半径长。过其中一个交点作直线,该直线与圆的另外两个交点分别为A、B。线段AB最长是多少? 输入数据 第一行有三个用空格隔开的整数x1,y1,r1,依次表示第一个圆的圆心坐标和半径; 第二行有三个用空格隔原创 2017-09-10 10:46:47 · 2263 阅读 · 0 评论 -
循环(NOIP2005普及组第四题)
循环(circle.pas/circle.in/circle.out)(NOIP2005普及组T4) (LYOI20090321信息学综合模拟Problem2)问题描述 乐乐是一个聪明而又勤奋好学的孩子。他总喜欢探求事物的规律。一天,他突然对数的正整数次幂产生了兴趣。 众所周知,2的正整数次幂最后一位数总是不断的在重复2,4,8,6,2,4,8,6……我们说2的正整数次幂最后一位的循环长度是4原创 2017-09-10 10:51:30 · 3471 阅读 · 0 评论 -
素数筛法
素数筛法这一次我们一步一步的来了解素数筛法。O(nlogn)筛法一般情况下我们可能比较熟悉的就是O(nlogn)的筛法。for(int i=2;i<=n;i++){ if(!vis[i]) { for(int j=2;j*i<=n;j++) vis[j*i]=1; }}时间复杂度计算: n/1+n/2+…+n/n=nlognO(nlog原创 2017-10-05 20:26:07 · 310 阅读 · 0 评论 -
费马小定理的思考与证明
费马小定理定义费马小定理(Fermat’s little theorem)是数论中的一个重要定理,在1636年提出,其内容为: 假如p是质数,且gcd(a,p)=1,那么 a(p-1)≡1(mod p),即:假如a是整数,p是质数,且a,p互质(即两者只有一个公约数1),那么a的(p-1)次方除以p的余数恒等于1。证明后序个人认为费马小定理证明的难点就是在于设出两个数列,这一点想起来比较困难。再一个原创 2017-10-06 20:41:02 · 769 阅读 · 0 评论 -
Num(数论(或打表))
Num(num.* 时间空间限制: 1S, 128M)) (LYOI20171011模拟赛第一试第二题) 题目描述: Tky 的 GF 也是一个大牛,现在她给 Tky 出了一道题, Tky 觉得这题太简单,扔给你了。 题目是这样的: 给定正整数 N, 令 f(i)表示 i 的所有正因子的和. 求 f(1) xor f(2) xor f(3) ….xor f(n) 输入格式: 一个正整原创 2017-10-11 20:15:37 · 2904 阅读 · 0 评论 -
欧拉公式
欧拉公式山东省聊城第一中学 Stockholm_Sun最近在玩一些niubility的算法和公式,学到傅里叶的时候有点蒙蔽,据说要用到欧拉公式,所以今天复习下欧拉公式吧。欧拉公式,你以为是R+V-E=2吗?不不不,今天我们讲讲数论上的欧拉公式。 推导下面给出这个公式的推导。(字丑勿喷) 先导知识:泰勒公式 ...原创 2018-08-15 21:21:22 · 3887 阅读 · 0 评论