摘要： 用普通的筛法，1000 0000 的数，需要用时1600ms左右。#include <stdio.h>#include <string.h>#include <math.h> char dp[10000010]; const int inf = 10000001;int t = sqrt(inf); inline void prime ( ){ int i, j, k = 0,p; dp[1] = -1; for (i = 2; i <= t; i++) if( dp[i] != -1) { //p = inf / i; for (j = 2; j & 阅读全文

摘要： 扩展欧几里德1.定义：扩展欧几里德定理 对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。2.c语言实现#include <stdio.h>int q,x, y;void swap ( int *x1, int *y1){ int temp; temp = *x1; *x1 = *y1; *y1 = temp;} void fun( int a, int b){ if (b == 0 ) { x = 1; y = 0; q = a; } else { fun ( b, a % b); in 阅读全文

摘要： 1.同余定义：设两个整数a, b,正整数m, m | a - b, 即 m * c = k ( a - b),则称a 模m同余于b, 或a 与b 模m 同余，记做 a = b (mod m).2.同余性质：1） 自反性2) 传递性3） 对称性4） 若 a = b (mod m), c = d (mod m)则 a+b = b+ d(mod m) a - b= b - d(mod m) ac = bd (mod m);3.中国剩余定理CRT （中国剩余定理CRT）设m1,m2,...,mk是两两互素的正整数，即gcd(mi, mj) =1, i≠j, i,j = 1,2,...,k 则同余方程组 阅读全文