【正文】 { r=a%b。 a=b。 b=r。 } return(a)。} lcd(int u,int v,int h){ return(u*v/h)。}語法：m_of_n(int m, int n1, int m1, int* a, int head)參數(shù)：m： 組合數(shù)C的上參數(shù)n1： 組合數(shù)C的下參數(shù)m1： 組合數(shù)C的上參數(shù)，遞歸之用*a： 1～n的整數(shù)序列數(shù)組head： 頭指針返回值： null注意： *a需要自行產(chǎn)生 初始調(diào)用時，m=mhead=0 調(diào)用例子：求C(m,n)序列：m_of_n(m,n,m,a,0)。源程序： void m_of_n(int m, int n1, int m1, int* a, int head) { int i,t。 if(m10 || m1n1) return。 if(m1==n1) { for(i=0。im。i++) couta[i]39。 39。 // 輸出序列 cout39。\n39。 return。 } m_of_n(m,n11,m1,a,head)。 // 遞歸調(diào)用 t=a[head]。a[head]=a[n11+head]。a[n11+head]=t。 m_of_n(m,n11,m11,a,head+1)。 // 再次遞歸調(diào)用 t=a[head]。a[head]=a[n11+head]。a[n11+head]=t。} （FFT）語法：kkfft(double pr[],double pi[],int n,int k,double fr[],double fi[],int l,int il)。參數(shù)：pr[n]： 輸入的實部 pi[n]： 數(shù)入的虛部n，k： 滿足n=2^kfr[n]： 輸出的實部fi[n]： 輸出的虛部l： 邏輯開關(guān)，0 FFT，1 ifFTil： 邏輯開關(guān)，0 輸出按實部/虛部；1 輸出按模/幅角 返回值： null注意： 需要 源程序： void kkfft(pr,pi,n,k,fr,fi,l,il) int n,k,l,il。 double pr[],pi[],fr[],fi[]。 { int it,m,is,i,j,nv,l0。 double p,q,s,vr,vi,poddr,poddi。 for (it=0。 it=n1。 it++) { m=it。 is=0。 for (i=0。 i=k1。 i++) {j=m/2。 is=2*is+(m2*j)。 m=j。} fr[it]=pr[is]。 fi[it]=pi[is]。 } pr[0]=。 pi[0]=。 p=(*n)。 pr[1]=cos(p)。 pi[1]=sin(p)。 if (l!=0) pi[1]=pi[1]。 for (i=2。 i=n1。 i++) { p=pr[i1]*pr[1]。 q=pi[i1]*pi[1]。 s=(pr[i1]+pi[i1])*(pr[1]+pi[1])。 pr[i]=pq。 pi[i]=spq。 } for (it=0。 it=n2。 it=it+2) { vr=fr[it]。 vi=fi[it]。 fr[it]=vr+fr[it+1]。 fi[it]=vi+fi[it+1]。 fr[it+1]=vrfr[it+1]。 fi[it+1]=vifi[it+1]。 } m=n/2。 nv=2。 for (l0=k2。 l0=0。 l0) { m=m/2。 nv=2*nv。 for (it=0。 it=(m1)*nv。 it=it+nv) for (j=0。 j=(nv/2)1。 j++) { p=pr[m*j]*fr[it+j+nv/2]。 q=pi[m*j]*fi[it+j+nv/2]。 s=pr[m*j]+pi[m*j]。 s=s*(fr[it+j+nv/2]+fi[it+j+nv/2])。 poddr=pq。 poddi=spq。 fr[it+j+nv/2]=fr[it+j]poddr。 fi[it+j+nv/2]=fi[it+j]poddi。 fr[it+j]=fr[it+j]+poddr。 fi[it+j]=fi[it+j]+poddi。 } } if (l!=0) for (i=0。 i=n1。 i++) { fr[i]=fr[i]/(*n)。 fi[i]=fi[i]/(*n)。 } if (il!=0) for (i=0。 i=n1。 i++) { pr[i]=sqrt(fr[i]*fr[i]+fi[i]*fi[i])。 if (fabs(fr[i])*fabs(fi[i])) { if ((fi[i]*fr[i])0) pi[i]=。 else pi[i]=。 } else pi[i]=atan(fi[i]/fr[i])*。 } return。 } 語法：result=integral(double a,double b)。參數(shù)：a： 積分上限b： 積分下限function f： 積分函數(shù)返回值： f在（a,b）之間的積分值注意： function f(x)需要自行修改，程序中用的是sina(x)/x 需要 默認精度要求是1e5源程序： double f(double x){ return sin(x)/x。 //在這里插入被積函數(shù) }double integral(double a,double b) { double h=ba。 double t1=(1+f(b))*h/。 int k=1。 double r1,r2,s1,s2,c1,c2,t2。 loop: double s=。 double x=a+h/。 while(xb) { s+=f(x)。 x+=h。 } t2=(t1+h*s)/。 s2=t2+(t2t1)/。 if(k==1) { k++。h/=。t1=t2。s1=s2。 goto loop。 } c2=s2+(s2s1)/。 if(k==2){ c1=c2。k++。h/=。 t1=t2。s1=s2。 goto loop。 } r2=c2+(c2c1)/。 if(k==3){ r1=r2。 c1=c2。k++。 h/=。 t1=t2。s1=s2。 goto loop。 } while(fabs(1r1/r2)1e5){ r1=r2。c1=c2。k++。 h/=。 t1=t2。s1=s2。 goto loop。 } return r2。} 語法：result=js(int s[][],int n)參數(shù)：s[][]： 行列式存儲數(shù)組n： 行列式維數(shù)，遞歸用返回值： 行列式值注意： 函數(shù)中常數(shù)N為行列式維度，需自行定義源程序： int js(s,n) int s[][N],n。 { int z,j,k,r,total=0。 int b[N][N]。/*b[N][N]用于存放，在矩陣s[N][N]中元素s[0]的余子式*/ if(n2) { for(z=0。zn。z++) { for(j=0。jn1。j++) for(k=0。kn1。k++) if(k=z) b[j][k]=s[j+1][k+1]。 else b[j][k]=s[j+1][k]。 if(z%2==0) r=s[0][z]*js(b,n1)。 /*遞歸調(diào)用*/ else r=(1)*s[0][z]*js(b,n1)。 total=total+r。 } } else if(n==2) total=s[0][0]*s[1][1]s[0][1]*s[1][0]。 return total。 } 語法：result=P(long n,long m)。 / result=long C(long n,long m)。參數(shù)：m： 排列組合的上系數(shù)n： 排列組合的下系數(shù)返回值： 排列組合數(shù)注意： 符合數(shù)學(xué)規(guī)則：m＝n源程序： long P(long n,long m){ long p=1。 while(m!=0) {p*=n。n。m。} return p。} long C(long n,long m){ long i,c=1。 i=m。 while(i!=0) {c*=n。n。i。} while(m!=0) {c/=m。m。} return c。} 二、字符串處理語法：replace(char str[],char key[],char swap[])。參數(shù)：str[]： 在此源字符串進行替換操作key[]： 被替換的字符串，不能為空串swap[]： 替換的字符串，可以為空串，為空串表示在源字符中刪除key[]返回值： null注意： 默認str[]長度小于1000，如否，重新設(shè)定設(shè)定tmp大小 需要 源程序： void replace(char str[],char key[],char swap[]){ int l1,l2,l3,i,j,flag。 char tmp[1000]。 l1=strlen(str)。 l2=strlen(key)。 l3=strlen(swap)。 for (i=0。i=l1l2。i++) { flag=1。 for (j=0。jl2。j++) if (str[i+j]!=key[j]) {flag=0。break。} if (flag) { strcpy(tmp,str)。 strcpy(amp。tmp[i],swap)。 strcpy(amp。tmp[i+l3],amp。str[i+l2])。 strcpy(str,tmp)。 i+=l31。 l1=strlen(str)。 } }}語法：result=strfind(char str[],char key[])。參數(shù)：str[]： 在此源字符串進行查找操作key[]： 被查找的字符串，不能為空串返