算法代码大全.doc

数论算法： 一.求两数的最大公约数 int gcd(int a,int b) { if(b==0) return a; else return gcd(b,a%b); } 二.分解质因数 void fjzys(int n)//将质因数储存在a[]数组中 不包括* {long n,i,j=0; a[j++]=1; for(i=2;i=n;i++) for(;n%i==0;n/=i) a[j++]=i; } 三.求质数表 #includefstream using namespace std; bool h[50002]={false}; int main(void) { h[0]=true; h[1]=true; h[2]=false; int j; ofstream f1(z.txt); for(int i=2;i50001;i++) { j=i+i; if(!h[i]) { while(j50001) { h[j]=true; j+=i; } } } f1int prime={; for(int i=0;i50001;i++) { if(!h[i]) f1i,; } f1};; } 高精度系列 零.通用: int lc,x[100],y[100]; void chartoint(char a[],char b[])//字符串数组转整型数组 { int la=strlen(a),lb=strlen(b); int i; if(lalb) lc=la; else lc=lb; for(i=0;ila;i++) x[la-i-1]=a[i]-48; for(i=0;ilb;i++) y[lb-i-1]=b[i]-48; } void print(int c[])//输出高精度计算结果 稍微一改即可变为输出在数组中 { for(int i=lc-1;i=0;i--) coutc[i]; } 一.高精度加法： void gjjf(int a[],int b[],int c[])//a[]与b[]所记录的数相加得到c[] lc为a b最长位数 { int m=0,i; for(i=0;ilc;i++) { c[i]=(a[i]+b[i]+m)%10; m=(a[i]+b[i]+m)/10; } if(m0) c[lc++]=1; } 二.高精度减法: void gj_f(int a[],int b[],int c[])//a,b已经倒置的被减数和减数数组c倒置输出结果 lc为a b最长位数 { int i; for (i=0;ilc;i++) { c[i]=c[i]+(a[i]-b[i]); if(c[i]0) //借位问题 {c[i]=c[i]+10; c[i+1]--;} } while (lc0c[lc-1]==0) lc--; } 三.高精度乘法： ①.高精度乘低精度 void gjcdj(int a[],int b,int c[])//a为高精度数，b为单精度,c为计算结果的倒置，！转化函数不再适用 { int i,j; for(i=0;i=lc;i++) c[i]=a[i]*b; for(j=0;jlc||a[j]10;j++) { c[j+1]=c[j+1]+c[j]/10; c[j]=c[j]%10; } while(c[lc-1]==0lc0) lc--; } ②.高精度阶乘(引申) void gjdjc(int n,int c[])//求n! c数组倒置储存结果 lc初始值为0 { int i,j; //循环计数变量 c[0]=1; for(i=2;i=n;i++,lc=j) //算阶乘 { for(j