
数论
文章平均质量分 79
Hacker_vision
一件事,你可以不去做,要做,就要全力以赴,做到最完美。
展开
-
快速幂(整数+矩阵):hdu 2817+poj3233+poj1995+poj3070
快速幂的思想:hdu 2817(整数快速幂) A sequence of numbers 解题思路:分别写出等差、等比数列通项an的表达式: 快速幂取模+基本的模运算就做了,附代码:{优快云:CODE:5888原创 2015-01-25 21:23:33 · 946 阅读 · 0 评论 -
nefu117素数个数的位数(素数定理+大数位数公式)
素数个数的位数Problem : 117Time Limit : 1000msMemory Limit : 65536Kdescription小明是一个聪明的孩子,对数论有着很浓烈的兴趣。他发现求1到正整数10n 之间有多少个素数是一个很难的问题,该问题的难以决定于n原创 2015-01-19 14:29:27 · 711 阅读 · 0 评论 -
同余定理:poj2769 Reduced ID Numbers(标记数组)+hdu1021 Fibonacci Again(应用于递推公式)
Reduced ID NumbersDescriptionT. Chur teaches various groups of students at university U. Every U-student has a unique Student Identification Number (SIN). A SIN s is an integer in the range 0原创 2015-01-21 21:25:04 · 1159 阅读 · 0 评论 -
poj 2891 Strange Way To Express Integers(线性同余方程组)
Strange Way to Express IntegersDescriptionElina is reading a book written by Rujia Liu, which introduces a strange way to express non-negative integers. The way is described as following:原创 2015-01-22 11:11:01 · 1290 阅读 · 0 评论 -
hdu 1573: X问题(线性同余方程组求正整数解的个数)
X问题Problem Description求在小于等于N的正整数中有多少个X满足:X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 Input输入数据的第一行为一个正整数T,表示有T组测试数据。每组测试数据的第一行为两个正整数N,M (0原创 2015-01-22 21:34:43 · 1493 阅读 · 0 评论 -
梅森素数:nefu 120 梅森素数(卢卡斯—莱默判别法+大数乘积取模)
梅森素数description由于梅森学识渊博,才华横溢,为人热情以及最早系统而深入地研究2p-1 型的数(其中p为素数),为了纪念他,数学界就把这种数称为“梅森数”;并以Mp 记之(其中M为梅森姓名的首字母),即Mp=2p-1 。如果梅森数为素数,则称之为“梅森素数”。 比如p=2,3,5,7时,Mp都是素数,但211-1原创 2015-01-21 16:25:33 · 3404 阅读 · 1 评论 -
线性同余方程:poj 2115 C Looooops
C LooooopsDescriptionA Compiler Mystery: We are given a C-language style for loop of type for (variable = A; variable != B; variable += C) statement;I.e., a loop which starts by setti原创 2015-01-22 13:43:25 · 643 阅读 · 0 评论 -
【n!素因子p的幂】swjtuOJ 2090
【n!素因子p的幂】swjtuOJ 2090【注:交大的看到这篇文章要学会自己写,不要为了比赛而比赛!~】题目大意数论一道好题:给以两个大整数n,s(n<=10^18,s<=10^12),试找到最大的整数k使得n! % s^k ==0 数论一道不错的题目,很容易想到思路,但是数据会大一点,有可能爆long long ,笔者由于n!素因子p的幂采用累乘法,在10^12左右的一个素数爆掉了,QAQ原创 2015-08-07 10:10:27 · 1187 阅读 · 0 评论 -
【A/B%m+扩展欧几里得】hdu 1576 A/B
hdu 1576 A/B http://acm.hdu.edu.cn/showproblem.php?pid=1576思路A/B%m经典问题,在(B,m)=1的情况下,转换为Bx+my = n ,扩展欧几里得得到最小正整数解即可,数论水题!参考代码 #include<iostream>#include<cstdio>#include<cstdlib>using namespace std;t原创 2015-09-16 22:14:46 · 778 阅读 · 0 评论 -
【矩阵快速幂】ZOJ 2974 Just Pour the Water
Source : ZOJ 2974 Just Pour the Water http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1973 Problem Description Shirly is a very clever girl. Now she has two containers (A and B), each wit原创 2015-09-16 20:39:06 · 485 阅读 · 0 评论 -
【高斯消元】poj 1222
http://poj.org/problem?id=1222高斯消元,关键是初始化构造系数矩阵,有了增广矩阵套模板即可/* poj 1222 高斯消元题意: 类似于翻转棋,按一下上下中左右翻转,给定初始状态,问回到全0的状态的按键方案思路: 对于每个位置,倒着想,初始为0,上下中左右加起来为初始状态,每一个位置皆如此。 30个位置,30个变量,30个方原创 2015-10-11 20:12:36 · 438 阅读 · 0 评论 -
【高斯消元】poj 1830 开关问题
http://poj.org/problem?id=1830高斯消元—求方程组解的个数,注意是0-1方程,若有解个数为2^(自由变元数),唯一解是2^0=1满足/* poj 1830 高斯消元-求方程组解个数题意: 一堆开关,一个开关的变化会改变其他开关的状态,给定他们的制约关系, 以及初始和终了状态,判断可否操作实现思路: 每个开关看成0-1变量,n个开关有n个变元,原创 2015-10-11 21:27:35 · 537 阅读 · 2 评论 -
【高斯消元】poj 1681 Painter's Problem
http://poj.org/problem?id=1681高斯消元求最小步数,二进制枚举0-1变元,统计和的最小值/* poj 1222 高斯消元-枚举自由变元题意: 0-1开关,n*n个方程,n*n个未知量,求最小操作的步数思路: 充分理解高斯消元的过程,如果有变元,二进制枚举 利用初等行变换最终的上三角阵,自下向上求解剩余变量,统计1的个数原创 2015-10-12 20:56:17 · 446 阅读 · 0 评论 -
素因子分解:nefu118(n!后面有多少个0)+ nefu119 (组合素数)
n!后面有多少个0description从输入中读取一个数n,求出n!中末尾0的个数。input输入有若干行。第一行上有一个整数m,指明接下来的数字的个数。然后是m行,每一行包含一个确定的正整数n,1<=n<=1000000000。output原创 2015-01-21 09:51:40 · 616 阅读 · 0 评论 -
素数筛选(2次):poj2689 Prime Distance
Prime DistanceDescriptionThe branch of mathematics called number theory is about properties of numbers. One of the areas that has captured the interest of number theoreticians for thousands of原创 2015-01-20 22:52:30 · 712 阅读 · 0 评论 -
poj1061 青蛙的约会(扩展欧几里得算法求解同余方程)
题目分析:假设走了t次相遇,则有等式(x+mt)-(y+nt)=pL成立,等价于求解同余方程(n-m)t≡(x-y) (mod L)的最小整数解(a)对于一般同余方程ax=d mod b,方程有解,则有(a,d)| b ,所以问题第一步判断解的情况(b)有(n-m)t+pL=x-y,t、p均为未知变量,将问题转化为求解ax+by=d的最小整数x,扩展欧几里得算法:brief原创 2015-01-19 09:53:13 · 8869 阅读 · 2 评论 -
高次同余方程:poj 3243+poj 2417+hdu 2815 (Baby Step Giant Step 算法)
个人水平有限,套用的模板,学习资料如下:题目链接:poj 3243 http://poj.org/problem?id=3243 poj 2417 http://poj.org/problem?id=2417 hdu 2815 http://acm.hdu.edu.cn/showproblem.php?pid=原创 2015-01-23 15:36:43 · 918 阅读 · 0 评论 -
hdu 3579:Hello Kiki (线性同余方程组求正整数解)
Hello KikiProblem DescriptionHello Kiki is such a lovely girl that she loves doing counting in a different way. For example, when she is counting X coins, she count them N times. Each time she div原创 2015-01-22 23:16:27 · 711 阅读 · 0 评论 -
毕达哥拉斯三元组:poj 1305+佩尔方程:poj 1320,hdu 3292(特殊不定方程)
毕达哥拉斯三元组:正整数x,y,z构成一个本原毕达哥拉斯三元组且y为偶数,当且仅当存在互素的正整数m,n(m>n),其中m为奇数n为偶数,或者n为偶数m为奇数,并且满足:原创 2015-01-28 10:27:51 · 1054 阅读 · 0 评论 -
poj 1006+hdu 1788(中国剩余定理求解同余方程组)
中国剩余定理:求解此类同余方程组最小非负整数解的算法实现:ll China(){ ll M1,ans=0; M=1; for(int i=1;i<=3;i++) M*=m[i]; for(int i=1;i<=3;i++){ M1=M/m[i]; t=extended_gcd(M1,m[i],x,y); ans=(ans原创 2015-01-26 23:38:13 · 1441 阅读 · 0 评论 -
不定方程:poj 2142+poj 1091(数论+容斥)
二元一次不定方程ax+by=d与同余方程ax=d(mod b)等价,求解线性同余方程可参考链接:http://blog.youkuaiyun.com/u012717411/article/details/42869291poj 2142 The Balance题目链接:http://poj.org/problem?id=2142解题思路:扩展欧几里得很容易得到通解x=x0+b0*t,y=y0-原创 2015-01-27 15:30:40 · 630 阅读 · 0 评论 -
poj 3641 Pseudoprime numbers(Miller-Rabin素数测试)
题目链接:http://poj.org/problem?id=3641题目大意:根据已a为基的伪素数的定义,判断p是否是已a为基的伪素数。(p为合数且满足a^p mod p = a), p解题思路:p参考代码:#include #include #include #include #include #include #include using namespa原创 2015-01-28 20:17:30 · 1075 阅读 · 0 评论 -
整数分解(3种算法比较):hdu 1164 Eddy's research I+poj 1811 Prime Test
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1164算法一:试除法算法思路:i = 1 to n^(1/2) ,找到能整除的素因子。该方法容易想到,但进行了很多不必要的运算,时间复杂度比较高。参考代码:#include #include #include #include #include #include #in原创 2015-01-29 13:58:56 · 935 阅读 · 0 评论 -
因子和与因子个数的应用:poj 2992+poj 1845+hdu 1452+poj 2480
求解因子和与因子个数的本质是——素因子分解poj 2992: Divisors解题思路:把组合数写成阶乘的形式,以n!为例,小于等于n的prime都是素因子,利用数论的知识:n!素因子分解中素数p的幂为[n/p]+[n/p^2]+[n/p^3]+[n/p^4]+……可以求出每个素因子在n!中的幂,+1乘起来就是因子的总个数。注意:这道题很容易超时,先打表,少用lon原创 2015-01-31 15:39:54 · 738 阅读 · 0 评论 -
数论及应用之总结提纲
(一)GCD、LCM1.GCD、LCMll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}ll lcm(ll a,ll b){return a*b/gcd(a,b);}2.扩展欧几里得算法:求等式ax+by=gcd(a,b)中的x,y;返回d=gcd(a,b)ll extended_gcd(ll a,ll b,ll &x,l原创 2015-02-02 18:15:09 · 1175 阅读 · 0 评论 -
完全数+莫比乌斯反演hdu 1695+伪随机数 poj2183
完全数: 完全数,又称完美数或完备数,是一些特殊的自然数,它所有的真因子(即除了自身以外的约数)之和,恰好等于它本身。或者可以理解为,如果n是一个正整数,且所有正因子之和等于2n,那么n称为完全数。完全数有6个重要性质:性质1 完全数都能写成连续自然数之和;性质2 每个完全数的全部因数倒数之和都是2,因此每个完全数都是调和数;性质3 除了6以外的完全数,每个数都可以写成连续原创 2015-02-01 23:13:42 · 918 阅读 · 0 评论 -
欧拉函数:poj 2407, poj 1284,poj 2478,poj 3090,poj 3696,poj 3358
定义:欧拉函数 定理:求解欧拉函数的值可用如下定理实现(通过n的素因子分解) 算法实现:(一)直接实现:直接套用定理求解欧拉函数值ll phi(ll n){ ll res=n; for(int i=2;i*i<=n;i++){ if(n%i==0){ res=res-res/i;原创 2015-01-29 18:50:33 · 655 阅读 · 0 评论 -
nefu84 五指山(扩展欧几里得算法解同余方程)
五指山Problem : 84Time Limit : 1000msMemory Limit : 65536Kdescription西游记中孙吾空大闹天宫,如来佛祖前来降伏他,说道:“我与你打个赌赛;你若有本事,一筋斗打出我这右手掌中,算你赢,再不用动刀兵苦争战,就请原创 2015-01-19 12:54:49 · 858 阅读 · 0 评论 -
[tzc3012:Fibnacci Numbers]降幂公式+矩阵快速幂+斐波那契求和公式
题目链接:http://acm.tzc.edu.cn/acmhome/problemdetail.do?&method=showdetail&id=3012DescriptionThe fibnacci numbers are as follows:f[1] = 1; f[2] = 2;f[n] = f[n - 1] + f[n - 2];And s[n原创 2016-10-31 11:26:59 · 1175 阅读 · 0 评论