![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
[Number Theory]
一刀不二
这个作者很懒,什么都没留下…
展开
-
POJ_2992_divisors
//懒得解释#include #include #include using namespace std;#define Type int#define Max 450bool prime[Max];Type primeTable[Max];Type primeRadiation(const Type& size){ prime[1]=fals原创 2013-10-14 14:38:29 · 726 阅读 · 0 评论 -
POJ_2407_relatives
//懒得解释#include #include #include #include using namespace std;#define Max 100000bool isPrime[Max];void primeRadiation(const int& size){ memset(isPrime,true,sizeof(isPrime));原创 2013-10-14 14:44:02 · 758 阅读 · 0 评论 -
素数判定 费马测试
;; Fermat's Little Theorem:;; If N is a prime number and A is any positive integer less ;; than N, then A raised to the N-th power is congruent to A modulo N;; Two numbers are said to be congruent原创 2014-07-05 01:54:50 · 2894 阅读 · 0 评论 -
埃拉托塞尼筛选法
#include #include #include #include using namespace std;#define Max 1000100000#define Type long longbool arr[Max];void getPrime(Type size){ for(Type i=2;i<=sqrt(size);++i){原创 2013-10-14 14:27:31 · 1900 阅读 · 0 评论 -
POJ_1809_Regetni
//懒得解释#include #include using namespace std;struct Node{ Node():x(0),y(0){}; int x; int y;};double area(Node a,Node b,Node c){ return fabs(a.x*b.y-a.y*b.x+b.x*c.y-b原创 2013-10-14 14:37:09 · 999 阅读 · 0 评论 -
素数辐射法(比埃拉托塞尼块)
//懒得解释#include #include #include #include #include #include using namespace std;#define Type int#define Max 1000000010bool prime[Max];bool prime1[Max];void primeRadiationMeth原创 2013-10-14 14:29:16 · 1594 阅读 · 0 评论 -
POJ_1730_Perfect Pth Powers
#include #include using namespace std;const double eps=1e-12;int main(){ double n; while(1){ cin>>n; if(n==0) break; if(n>0){ for(int i=31原创 2013-10-14 14:34:36 · 724 阅读 · 0 评论 -
POJ_2034_anti-prime_sequences
//懒得解释#include #include #include #include #include #include #include using namespace std;#define Max 3300#define Size 3200bool isPrime[Max];int n,m,d;void primeRadiation(){原创 2013-10-14 14:41:43 · 1107 阅读 · 0 评论 -
快速幂(二进方法)
#include using namespace std;// Right to left binaryint quickPower( int base, int exp ){ int res = 1; int product = base; while( exp ){ if( exp & 1 ){ res *= produ原创 2013-11-30 23:37:19 · 1063 阅读 · 0 评论 -
Cantor展开式
X=an*(n-1)!+an-1*(n-2)!+...+ai*(i-1)!+...+a2*1!+a1*0! 其中,a为整数,并且0康托展开的应用实例:{1,2,3,4,...,n}表示1,2,3,...,n的排列 如 {1,2,3} 按从小到大排列一共6个。123 132 213 231 312 321 。 代表的数字 1 2 3 4 5 6 也就是把10进制数与一个排列对应起来原创 2013-11-14 02:59:31 · 1484 阅读 · 0 评论 -
乌拉姆数列 Euler167
昨天 yep 问了一个问题:For two positive integers a and b, the Ulam sequence U(a,b) is defined by U(a,b)1 = a, U(a,b)2 = b and for k > 2, U(a,b)k is the smallest integer greater than U(a,b)(k-1) which can原创 2014-08-31 19:45:21 · 2172 阅读 · 0 评论 -
数论学习网站备份
http://www.math.uiuc.edu/Algebraic-Number-Theory/伊利诺伊州立大学数学系代数数论档案网,在英国还有镜像站点,包括代数数论方面的各种讨论课题。 http://www.maa.org/BLL/numtheory.htmMAA(美国数学协会)在线的数论主页,包括数论的各个分支领域研究、工作组、研究课题等内容。 htt转载 2014-07-26 18:45:49 · 1718 阅读 · 0 评论 -
POJ_2773_happy_2006(欧拉公式)
//懒得解释#include #include #include using namespace std;#define Max 1000000int prime[Max];int GCD(int a,int b){ if(a==0) return b; return GCD(b%a,a);}int main(){ int n;原创 2013-10-14 14:44:46 · 788 阅读 · 0 评论 -
POJ_1365_Prime_Land
//懒得解释#include #include #include #include using namespace std;#define Max 110000#define Size 100000#define Type intbool isPrime[Max];Type primeTable[Max];Type expCount[Max];v原创 2013-10-14 14:39:56 · 809 阅读 · 0 评论 -
POJ_2739_Sum_of_Consecutive_Prime
//懒得解释#include #include #include using namespace std;#define Max 10020#define Size 10010bool isPrime[Max];int primeTable[Max];void primeRadiation(const int& size){ memset(isPr原创 2013-10-14 14:43:01 · 675 阅读 · 0 评论