摘要： /* 对10000以内的素数打表 */ #include #include #include #include using namespace std; const int maxs = 10000+1;//只要改变这个的范围就能对多少范围内打素数表 int prime[maxs];//从1开始，prime[0]置为0 bool vis[maxs];//对数字做标记vis[100]... 阅读全文

摘要： /* 题目大意：对于m%a1=r1,m%a2=r2...m%ak=rk,求最小的非负的m的值 联立前面两个方程组则有a1*x-a2*y=r2-r1; 可利用欧几里得算法求出最小的非负x 那么满足前两个方程的一个特解m=a1*x+r1; 所有解M=m+x*LCD(a1,a2);---LCD(a1,a2)最小公倍数 在联立第3个方程，另a1 = LCD(... 阅读全文

摘要： 1 #include 2 #include 3 using namespace std; 4 5 int main() 6 { 7 int p,e,i,d,n,counts = 0,ans; 8 //(x+d)%23=p,(x+d)%28=e,(x+d)%33=i,使用中国剩余定理 9 //因为23,28,33，可直接使用定理,令x+d=n 10 /... 阅读全文

摘要： 通过poj1811整理这些算法的模板 所有代码如下： 阅读全文

摘要： 对于一个整数n,求小于n且和n互质的数的个数，可用欧拉函数求解。 例如eular(10)=4,互质的数有1,3,7,9. Euler函数表达通式：euler(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…(1-1/pn),其中p1,p2……pn为x的所有 素因数(素因数> 阅读全文

摘要： 对于整数a,b,x,y,c 有a*x+b*y=c,如果c不是a与b的最大公约数的倍数，那么此方程无解 证明：设gcd(a,b)=d，即最大公约数，那么a*x%d=0，b*y%d=0 则(a*x+b*y)%d=0,说明c是一个d的倍数，相反的，如果c不是d的倍数，那么次方程无解 对于欧几里得算法 in 阅读全文

摘要： 一.乘法快速取余 算a*b%n 二.乘方取模 算a^n%m 定理：要计算只包含加减乘的整数表达式除以整数m的余数时，可以在每步计算时对m取余 容易想到的代码如下： int ans = 1;for(int i=1;i<=n;i++) ans=ans*a%m; 这样做时间发杂度为O(n),如果n=10^ 阅读全文

摘要： Sorting It All Out Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 32966 Accepted: 11458 Description An ascending sorted sequence of distin 阅读全文

摘要： Borg Maze Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 12952 Accepted: 4227 Description The Borg is an immensely powerful race of enhanc 阅读全文

摘要： FDNY to the Rescue! Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 2917 Accepted: 896 Description The Fire Department of New York (FDNY) h 阅读全文