数论
我在浪里
I LOVE Coding,冲刺中。。。
展开
-
2020-09-07 因式分解定理
当然,答案肯定在[-10,10],这题按步长0.01枚举能过72%的数据。。。搜了一下,这个叫因式分解定理,上式f(x)=0,则x=-q/p是有理解,q是a0的因子,p是an的因子。def gcd(a,b): if a%b==0: return b return gcd(b,a%b)def fenjie(n): m = [] while n!=1: for i in range(2,int(n+1)): if n % i == 0: m.append(i..原创 2020-09-07 10:51:08 · 323 阅读 · 0 评论 -
poj1061——青蛙的约会(扩展欧几里得)
青蛙的约会Time Limit: 1000MS Memory Limit: 10000KTotal Submissions: 114186 Accepted: 23372Description两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰原创 2017-05-04 11:34:17 · 239 阅读 · 0 评论 -
扩展欧几里得最小整数解
我们先讨论a,b,c都大于0的情况 (都为整数)[cpp] view plain copy int exgcd(int a,int b,int &x,int &y) { int d=a; if(b!=0) { d=exgcd(b,a%b,y,x); y转载 2017-05-04 11:36:55 · 3609 阅读 · 2 评论 -
poj2891——Strange Way to Express Integers(扩展欧几里得解中国剩余定理)
Strange Way to Express IntegersTime Limit: 1000MS Memory Limit: 131072KTotal Submissions: 15859 Accepted: 5245DescriptionElina is reading a book written by Ruji原创 2017-05-04 14:37:25 · 347 阅读 · 0 评论 -
lightoj1236——Pairs Forming LCM(素因子分解)
B - Pairs Forming LCMTime Limit:2000MS Memory Limit:32768KB 64bit IO Format:%lld & %lluSubmit StatusDescriptionFind the result of the following code:long long pairsFormLCM(原创 2017-05-04 16:43:58 · 378 阅读 · 0 评论 -
uva10892——LCM Cardinality(最小公倍数质因子分解)
10892 - LCM CardinalityTime limit: 3.000 secondshttp://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=467&page=show_problem&problem=1833A pair of numbers ha原创 2017-05-04 16:49:01 · 354 阅读 · 0 评论 -
uva10791——Minimum Sum LCM(分解质因数)
题意:给一个数字n,范围在[1,2^23-1],这个n是一系列数字的最小公倍数,这一系列数字的个数至少为2 例如12,是1和12的最小公倍数,是3和4的最小公倍数,是1,2,3,4,6,12的最小公倍数,是12和12的最小公倍数……………… 那么找出一个序列,使他们的和最小,上面的例子中,他们的和分别为13,7,28,24……显然最小和为7 首先假设我们知道了一系列数字a1,a2原创 2017-05-04 16:57:20 · 320 阅读 · 0 评论 -
hdu4497——GCD and LCM(数论&容斥原理or排列组合)
GCD and LCM Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Submission(s): 2710 Accepted Submission(s): 1186Problem DescriptionGiven two positiv原创 2017-04-22 22:10:30 · 521 阅读 · 0 评论 -
hzau1202——GCD(斐波那契&矩阵快速幂)
1202: GCDTime Limit: 1 Sec Memory Limit: 1280 MBSubmit: 184 Solved: 26[Submit][Status][Web Board]DescriptionInput The first line is an positive integer T . (1Output原创 2017-04-24 19:18:25 · 344 阅读 · 0 评论 -
Torry的困惑(基本型)(判断素数)
问题描述 Torry从小喜爱数学。一天,老师告诉他,像2、3、5、7……这样的数叫做质数。Torry突然想到一个问题,前10、100、1000、10000……个质数的乘积是多少呢?他把这个问题告诉老师。老师愣住了,一时回答不出来。于是Torry求助于会编程的你,请你算出前n个质数的乘积。不过,考虑到你才接触编程不久,Torry只要你算出这个数模上50000的值。输入格式 仅包含一个...原创 2019-02-02 14:13:34 · 198 阅读 · 0 评论 -
素因子去重(分解素因子模板)
问题描述 给定一个正整数n,求一个正整数p,满足p仅包含n的所有素因子,且每个素因子的次数不大于1输入格式 一个整数,表示n输出格式 输出一行,包含一个整数p。样例输入1000样例输出10数据规模和约定 n<=10^12 样例解释:n=1000=2^3*5*3,p=2*5=10注意n应该是long long类型#include...原创 2019-02-02 16:05:15 · 607 阅读 · 0 评论 -
北航本科密码学实验一
# while True: # try: # except: # break# 厄拉多塞筛法O(nloglogn) # N = 1000000# isprime = [0 for i in range(N)]# prime = [0 for i in range(N)] # n = ...原创 2019-02-28 10:45:21 · 817 阅读 · 0 评论 -
密码学实验二
MILLER-RABIN检测(转)Miller-Rabin检测是目前应用比较广泛的一种二次探测定理:如果p是一个素数,且0<x<p,则方程x^2%p=1的解为:x=1或x=p-1费马小定理:a^(p-1) ≡ 1(mod p)这就是Miller-Rabin素性测试的方法。不断地提取指数n-1中的因子2,把n-1表示成d*2^r(其中d是一个奇数)。那么我们需要计算的东西就...原创 2019-03-13 11:24:05 · 881 阅读 · 0 评论 -
7.1一步之遥
从昏迷中醒来,小明发现自己被关在X星球的废矿车里。 矿车停在平直的废弃的轨道上。 他的面前是两个按钮,分别写着“F”和“B”。小明突然记起来,这两个按钮可以控制矿车在轨道上前进和后退。 按F,会前进97米。按B会后退127米。透过昏暗的灯光,小明看到自己前方1米远正好有个监控探头。 他必须设法使得矿车正好停在摄像头的下方,才有机会争取同伴的援助。或许,通过多次操作F和B可以办到。矿...原创 2019-04-09 23:08:44 · 266 阅读 · 0 评论 -
伪随机大素数生成
from random import randint#使用randint需要加上这句def growbin(w): #w是设定的二进制位数 list = [1] for i in range(w-2): list.append(randint(0,1)) list.append(1) s2 = [str(_) for _ in list] #['...原创 2019-05-22 11:37:51 · 861 阅读 · 0 评论 -
NOIP 2001-2——最大公约数和最小公倍数问题(简单推导/分解质因数)
问题描述输入二个正整数x0,y0(2条件: 1.P,A是正整数2.要求P,Q以x0为最大公约数,以y0为最小公倍数.试求:满足条件的所有可能的两个正整数的个数.样例输入:x0=3 yo=60输出:4说明(不用输出)此时的 P Q 分别为: 3 6015 1212 1560 3所以:满足条件的所原创 2017-05-03 16:56:50 · 1543 阅读 · 0 评论 -
乘法逆元小结
乘法逆元的定义貌似是基于群给出的,比较简单地理解,可以说是倒数的概念的推广。记a的关于模p的逆元为a^-1,则a^-1满足aa^-1≡ 1(mod p)加减乘与模运算的顺序交换不会影响结果,但是除法不行。有的题目要求结果mod一个大质数,如果原本的结果中有除法,比如除以a,那就可以乘以a的逆元替代。在mod p的运算中,a存在乘法逆元当且仅当a与p互质。一般题目给的是一个大质数,所原创 2017-05-03 10:54:59 · 443 阅读 · 0 评论 -
hdu 4990 矩阵快速幂
递推公式:f(n)=2*f(n-2)+f(n-1)+1数据量太大 矩阵快速幂#include #include #include #include using namespace std;long long int mod;const int m = 3;struct Matrix{ long long int mat[3][3];};原创 2017-03-24 14:48:27 · 191 阅读 · 0 评论 -
hdu4910 找规律+大素数检测
题意:给定一个数字,求出1 - n之间与他互质的数的乘积mod n 规律:对于 1,2,4,P^n, 2*P^n,答案为 N-1,其余情况都是1。也就是说,1,2,4,以及只有一个质因子(奇数)或者它的1/2只有一个 质因子(偶数),答案是N-1.其余情况答案均是 1然后由于这题的n很大,也没法直接判定,可以这样搞,先原创 2017-03-17 14:00:39 · 416 阅读 · 0 评论 -
sdut3257——Square Number(打表)
Square NumberTime Limit: 1000MS Memory Limit: 65536KBSubmit Statistic DiscussProblem DescriptionIn mathematics, a square number is an integer that is the square of an integer. In other原创 2017-04-28 20:17:22 · 489 阅读 · 0 评论 -
sdut3257——Cube Number(打表)
Cube NumberTime Limit: 2000MS Memory Limit: 65536KBSubmit Statistic DiscussProblem DescriptionIn mathematics, a cube number is an integer that is the cube of an integer. In other words原创 2017-04-29 13:00:32 · 446 阅读 · 0 评论 -
判断互质的五种方法
转自:http://res.tongyi.com/resources/old_article/student/1380.html判断互质数的五种方法一. 概念判断法公约数只有1的两个数叫做互质数。根据互质数的概念可以对一组数是否互质进行判断。如:9和11的公约数只有1,则它们是互质数。二. 规律判断法根据互质数的定义,可总结出一些规律,利用这些规律能转载 2017-05-01 16:25:30 · 17309 阅读 · 0 评论 -
sdut2605——A^X mod P
A^X mod PTime Limit: 5000ms Memory limit: 65536K 有疑问?点这里^_^题目描述It's easy for ACMer to calculate A^X mod P. Now given seven integers n, A, K, a, b, m, P, and a function f(x) whic原创 2017-04-19 21:59:31 · 273 阅读 · 0 评论 -
hdu5108——Alexandra and Prime Numbers(素数)
Alexandra and Prime NumbersTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2208 Accepted Submission(s): 723Problem DescriptionAle原创 2017-04-16 20:49:18 · 327 阅读 · 0 评论 -
hdu5104——Primes Problem(素数)
Primes Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3384 Accepted Submission(s): 1515 Problem DescriptionGiven a num...原创 2017-04-16 20:42:55 · 323 阅读 · 0 评论 -
hdu1787——GCD Again(欧拉函数入门)
GCD AgainTime Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 3583 Accepted Submission(s): 1575Problem DescriptionDo you have spent so原创 2017-05-02 19:32:02 · 427 阅读 · 0 评论 -
HDU1695——GCD(莫比乌斯反演)
三个启发性连接:http://blog.csdn.net/lixuepeng_001/article/details/50577932http://www.cnblogs.com/femsub/p/5723780.htmlhttp://blog.csdn.net/acdreamers/article/details/8542292GCDTime Limit:原创 2017-05-02 19:25:16 · 651 阅读 · 0 评论 -
hdu2504——又见GCD(水)
思路gcd(a/b,c/b)==1 #include#include#includeusing namespace std;typedef long long ll;int n,m;int gcd(int a,int b){ if(b==0)return a; return gcd(b,a%b);}int main(){ in原创 2017-05-02 19:56:26 · 516 阅读 · 0 评论 -
hdu2588——GCD(欧拉函数)
GCDTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 2136 Accepted Submission(s): 1087Problem DescriptionThe greatest common diviso原创 2017-05-02 20:22:03 · 549 阅读 · 0 评论 -
poj1222——高斯消元基础题
题意:给出一个5*6的图,每个灯泡有一个初始状态,1表示亮,0表示灭。每对一个灯泡操作时,会影响周围的灯泡改变亮灭,问如何操作可以使得所有灯泡都关掉。思路:因为每盏灯,如果操作两次就相当于没有操作,所以相当于(操作次数)%2,即异或操作。考虑一个2*3的图,最后需要的状态是 :,如果初始状态为:。对这两个矩阵的每个数字做异或操作可以得到线性方程组每个方程的答案。总共6盏灯,原创 2017-04-05 12:45:47 · 300 阅读 · 0 评论 -
hdu4675——GCD of Sequence(莫比乌斯反演+组合数取模+乘法逆元+快速幂)
GCD of SequenceTime Limit: 6000/3000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Submission(s): 1595 Accepted Submission(s): 540Problem DescriptionAlice is playin原创 2017-05-03 10:53:16 · 567 阅读 · 0 评论 -
10^11内素数统计
#include #define MAXN 100 // pre-calc max n for phi(m, n)#define MAXM 100010 // pre-calc max m for phi(m, n)#define MAXP 666666 // max primes counter#define MAX 10000010 // max prime原创 2016-09-19 19:09:45 · 1292 阅读 · 0 评论