C語言常用代碼

1、-WORD格式-可編輯-#includeusing namespace std;void fullxunhuan(int*a,int n)int t=0;int m1,m2;if(n%2!=0)n+；for(int i=0;in%2;i+)for(m1=i,m2=n-i-1;m1i;m2-)am1m2=t%n+1;t+; for(;m1i;m1-) am1m2=t%n+1; t+; for(;m2n-i-1;m2+) am1m2=t%n+1; t+;void print(int *a,int n)for(i nt i=0;i n ;i+)for(intj=O;jvn;j+)coutaijvvt

2、;coutn;void main()int *a;int n;coutvv請輸入Nvn;a=new int *n;申請一個N行N列的數(shù)組for(inti1= 0;i1vn;i1+)ai1=new intn;for(int i=0;ivn;i+)給數(shù)組初始化for(int j=O;jvn;j+) aij=O； fullxunhuan(a,n); print(a,n);乘大數(shù)）4. 精度計算一一加法5. 精度計算減法6. 任意進制轉換7. 最大公約數(shù)、最小公倍數(shù)8. 組合序列9. 快速傅立葉變換（FFT）10. Ronberg算法計算積分11. 行列式計算12. 求排列組合數(shù)字符串處理：1. 字符

3、串替換2. 字符串查找3. 字符串截取計算幾何：1. 叉乘法求任意多邊形面積2. 求三角形面積3. 兩矢量間角度4. 兩點距離（2D 3D）5. 射向法判斷點是否在多邊形內部6. 判斷點是否在線段上7. 判斷兩線段是否相交8. 判斷線段與直線是否相交9. 點到線段最短距離ACM、組內部預定函數(shù)數(shù)學問題：1.精度計算一一大數(shù)階乘2. 精度計算一一乘法（大數(shù)乘小數(shù)）3. 精度計算一一乘法（大數(shù)10. 求兩直線的交點11. 判斷一個封閉圖形是凹集還是凸集12. Graham掃描法尋找凸包數(shù)論：1. x的二進制長度2. 返回x的二進制表示中從 低到高的第i位3. 模取幕運算4. 求解模線性方程5. 求

4、解模線性方程組(中國余數(shù)定理)6. 篩法素數(shù)產生器7. 判斷一個數(shù)是否素數(shù)圖論：I. Prim算法求最小生成樹2. Dijkstra算法求單源最短路徑3. Bellman-ford算法求單源最短路徑4. Floyd算法求每對節(jié)點間最短路徑排序/查找：1. 快速排序2. 希爾排序3. 選擇法排序4. 二分查找數(shù)據(jù)結構：1. 順序隊列2. 順序棧3. 鏈表4. 鏈棧5. 二叉樹、數(shù)學問題1. 精度計算一一大數(shù)階乘語法：int result=factorial(int n);參數(shù)：n：n的階乘返回值：階乘結果的位數(shù)本程序直接輸出n!的結 果，需要返回結果請保留long a需要math.h源程序：in

5、t factorial(int n)long a10000;int i,j,l,c,m=0,w;a0=1;for(i=1;i=n ;i+)c=0;for(j=0;j0)m+;am=c;w=m*4+log10(am)+1; printf(n%ld,am);for(i=m-1;i=0;i-)printf(%4.4ld,ai);return w;2. 精度計算一一乘法(大數(shù)乘小數(shù))語法：mult(char c,chart,int m);參數(shù)：c:被乘數(shù)，用字符串表示，位數(shù)不限t:結果，用字符串表示 m乘數(shù)，限定10以內返回值： null汪意:需要 string.h源程序：void mult(char

6、 c,chart,int m)int i,l,k,flag,add=O; char s100; l=strlen(c);for (i=0;il;i+)sl-i-1=ci-O;for (i=0;i=10)si=k%10;add=k/10;flag=1; elsesi=k;flag=0;add=0;if (flag)l=i+1;si=add; elsefor (i=0;il;i+)tl-1-i=si+0;tl=0;3. 精度計算一一乘法(大數(shù)乘大數(shù))語法：mult(char a,charb,char s);參數(shù)：a:被乘數(shù)，用字符串表示，位數(shù)不限b:乘數(shù)，用字符串表示，位數(shù)不限t:結果，用字符串表

7、示返回值： null空間復雜度為o(nA2) 需要 string.h 源程序：void mult(char a,charb,char s)inti,j,k=O,alen,blen,sum=O,res6565=0,flag=0;char result65;alen=strlen(a);blen=strl en(b);for (i=0;iale n;i+)for (j=0;j=0;i-)for(j=ble n-1;j=O;j-)sum=sum+resi+blen-j-1j;resultk=sum%10;k=k+1;sum=sum/10;for(i=blen-2;i=0;i-)for(j=0;j=i

8、;j+)sum=sum+resi-jj;if(strlen(s)!=strlen(a )&s0=0)strcpy(s,s+1);elsebreak;4. 精度計算一一加法語法：add(char a,charb,char s);參數(shù)：a:被乘數(shù)，用字符串表示，位數(shù)不限resultk=sum%10;k=k+1;sum=sum/10;if (sum!=0)resultk=sum;k=k+1;for (i=0;i=0;i-)si=resultk-1-i;sk=0;while(1)空間復雜度為o(nT) 需要 string.h 源程序：void add(char a,charb,char back)in

9、t i,j,k,up,x,y,z,l;char *c;if(strlen(a)strlen(b)l=strlen(a)+2; elsel=strlen(b)+2;c=(char *)b:乘數(shù)，用字符串表示，位數(shù)不限t:結果，用字符串表示返回值： nullmalloc(l*sizeof(char); i=strlen(a)-1; j=strlen(b)-1; k=O;up=O; while(i=O|j=O) if(i0)x=0;else x=ai; if(j9)up=1;z%=10; else up=0;ck+=z+0;i-;j-;if(up) ck+=1;i=0;ck=O;for(k-=1;k

10、=0;k-)backi+=ck;示，位數(shù)不限t:結果，用字符串表示返回值： null默認s1=s2，程序未處理負數(shù)情況需要 string.h源程序：void sub(char s1,chars2,char t)int i,l2,l1,k;I2=strlen(s2);l1=strlen(s1);tl1=0;l1 -;for(i=l2-1;i=0;i-,l1-)if(s1l1-s2i=0)tl1=s1l1-s2i+O;elsebacki=0;tl1=1O+s1l1-s2i+05. 精度計算減法語法：sub(char s1,chars2,char t);參數(shù)：s1:被減數(shù)，用字符串表示，位數(shù)不限s2

11、:減數(shù)，用字符串表s1l1- 1=s1l1 -1-1;k=l1;while(s1k=0)tl1=s1l1;l1-;loop:if (t0=0)I1=strlen(s1);for(i=0;ivl1 -1;i+)ti=ti+1;tl1 -1=0;goto loop;if (strlen(t)=0)t0=0;t1=O;6. 任意進制轉換語法：conversion(chars1,char s2,longd1,long d2);參數(shù)：s:原進制數(shù)字，用字符串表示s2:轉換結果，用字符串表示d1：原進制數(shù)d2：需要轉換到的進制數(shù)返回值：null高于9的位數(shù)用大寫AZ表示，216位進制通過驗證源程序：voi

12、d conversion(char s,char s2,long d1,long d2) long i,j,t,num; char c;num=0;for(i=0;si!=0;i+)if(siv=9&si=0)t=si-O; elset=si-A+10;num=num*d1+t;i=0;while(1)t=num%d2;if (t=9)s2i=t+O; elses2i=t+A-10;num/=d2;if (num=0)break;i+;for (j=0;ji/2;j+)C=s2j;s2j=si-j;s2i-j=c;s2i+1=0;7. 最大公約數(shù)、最小公倍數(shù)語法：resulet=hcf(int

13、a,int b) 、result=lcd(inta,int b)參數(shù)：a:int a，求最大公約數(shù)或最小公倍數(shù)b:int b，求最大公約數(shù)或最小公倍數(shù)返回值：返回最大公約數(shù)(hcf)或最小公倍數(shù)(led )lcd需要連同hcf 使用源程序：int hcf(int a,int b)int r=0;while(b!=0)r=a%b;a=b;b=r;return(a);lcd(int u,int v,int h)return(u*v/h);8. 組合序列 int ml, int* a, int head)參數(shù)：m 組合數(shù)C的上參數(shù)n1: 組合數(shù)C的下參數(shù) m1:組合數(shù)C的上參數(shù)，遞歸之用*a :1n

14、的整數(shù)序列數(shù)組head:頭指針返回值： nulla需要自行產生初始調用時，m=m1 head=0 調用例子：求C(m,n)序列： m_of_n(m,n,m,a,0);-源程序：void m_of_n(int m, intn1, int ml, int* a, inthead)int i,t;if(m1n1)return;if(m仁=n1)for(i=0;im;i+)coutvvaivv ; /輸出序列coutvn;return;m_of_n(m,n1-1,m1,a,head);語法：m_of_n(int m, int n1,/遞歸調用t=ahead;ahead=an1-1+head;an1-1

15、+head=t;m_of_n(m,n1-1,m1-1,a,head+1); /再次遞歸調用t=ahead;ahead=an1-1 +head;an1-1+head=t;9. 快速傅立葉變換(FFT)語法：kkfft(doublepr,double pi,intn,int k,doublefr,double fi,intl,int il);參數(shù)：prn:輸入的實部 pin:數(shù)入的虛部 n， k: 滿足 n=2Ak frn:輸出的實部 fin:輸出的虛部l :邏輯開關，0 FFT，1ifFTil :邏輯開關，0輸出按實 部/虛部；1輸出按模/幅角返回值： null需要math.h源程序：voidk

16、kfft(pr,pi,n,k,fr,fi,l,il)int n,k,l,il;doublepr,pi,fr,fi；intit,m,is,i,j,nv,l0;double p,q,s,vr,vi,poddr,poddi;for (it=0; it=n-1;it+)m=it; is=0;for (i=0; i=k-1;i+)j=m/2;is=2*is+(m-2*j); m=j; frit=pris; fiit=piis;pr0=1.0;pi0=0.0;p=6.283185306/(1.0*n);pr1=cos(p); pi1=-sin(p);if (l!=0)pi1=-pi1;for (i=2;

17、i=n-1;i+) p=pri-1*pr1; q=pii-1*pi1;s=(pri-1+pii-1)*(pr1+pi1);pri=p-q;pii=s-p-q;for (it=0; it=0;10-)m=m/2; nv=2*nv;for (it=0;it=(m-1)*nv; it=it+nv) for (j=0;jv=(nv/2)-1; j+)p=prm*j*frit+j+nv/2;q=pim*j*fiit+j+nv/2;s=prm*j+pim*j;s=s*(frit+j+nv/2+fiit +j+nv/2);poddr=p-q; poddi=s-p-q;frit+j+nv/2=frit+j 卜

18、p oddr;fiit+j+nv/2=fiit+j-p oddi;frit+j=frit+j+poddr;fiit+j=fiit+j+poddi; if (l!=0) for (i=0; i=n-1; i+)fri=fri/(1.0*n);fii=fii/(1.0*n); if (il!=0) for (i=0; i=n-1; i+)pri=sqrt(fri*fri+fii*fii)；if(fabs(fri)0)pii=90.0; else pii=-90.0;elsedouble integral(doublea,double b)double h=b-a;double10.Ronberg算

19、法計算積分語法：result=integral(doublea,double b);參數(shù)：a:積分上限b:積分下限function f : 積分函數(shù)返回值：f 在( a,b )之間的積分值pii=atan(fii/fri)* 360.0/6.283185306;return;t1=(1+f(b)*h/2.0;int k=1;doubler1,r2,s1,s2,c1,c2,t2;loop:double s=0.0;double x=a+h/2.0;while(xvb)s+=f(x);x+=h; t2=(t1+h*s)/2.0; s2=t2+(t2-t1)/3.0;if(k=1) function

20、 f(x)需要自行修改，程序中用的是 sina(x)/x 需要math.h默認精度要求是1e-5 源程序：double f(double x)return sin(x)/x; /在這里插入被積函數(shù)k+;h/=2.0;t1=t2;s1=s2;goto loop;c2=s2+(s2-s1)/15.0;if(k=2)c仁 c2;k+;h/=2.0;t1= t2;s1=s2;goto loop;r2=c2+(c2-c1)/63.0;if(k=3)r1=r2;c1=c2;k+;h/=2.0;t1= t2;s1=s2;goto loop;while(fabs(1-r1/r2)1e-5)r仁 r2;c1=c

21、2;k+;h/=2.0;t1= t2;s1=s2;goto loop;return r2;11.行列式計算語法：result=js(ints,int n)參數(shù)：s:行列式存儲數(shù)組n:行列式維數(shù)，遞歸用返回值：行列式值函數(shù)中常數(shù)N為行列式維度，需自行定義源程序：intjs(s,n)int sN,n;int z,j,k,r,total=0;int bNN;/*bNN用于存放，在矩陣sNN中元素s0的余子式*/if(n2)for(z=0;zn;z+)for(j=0;jn-1;j+)for(k=0;k=z)bjk=sj+1k+1;else bjk=sj+1k;if(z%2=0)r=s0 z *js(b

22、,n-1);歸調用*/elser=( -1)*s0z*js(b,n-1)/*遞total=total+r;else if(n=2)total=s00*s11-s 01*s10； return total;12.求排列組合數(shù)語法：result=P(long n,longm); / result=long C(longn,long m);參數(shù)：m 排列組合的上系數(shù) n： 排列組合的下系數(shù) 返回值： 排列組合數(shù)不能為空串swap:替換的字符串，可以為空串，為空串表示在源字符中刪除key返回值： null符合數(shù)學規(guī)則：mJ n源程序：long P(long n,long m)long p=1;whil

23、e(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;二、字符串處理1. 字符串替換語法：replace(charstr,char key,charswap);參數(shù)：str:在此源字符串進行替換操作key:被替換的字符串，默認str長度小于1000，如否，重新設定設定tmp大小需要 string.h 源程序：void replace(charstr,char key,charswap)intI1,l2,l3,i,j,f

24、lag; char tmp1000; l1=strlen(str); l2=strlen(key); l3=strlen(swap);for(i=0;i=l1 -I2;i+)flag=1;for(j=0;jl)return 0;for(i=start;istart+len;i+)strbackk+=stri;strbackk=O;return 1;三、計算幾何1. 叉乘法求任意多邊形面積語法：result=polygonarea(Point*polygon,int N);參數(shù)：polygon :多變形頂點數(shù)組 N:多邊形頂點數(shù)目返回值：多邊形面積支持任意多邊形，凹、凸皆可多邊形頂點輸入時按順時

25、 針順序排列 源程序：typedef struct double x,y; Point;double polygonarea(Point*polygon,int N)int i,j; double area = 0;for (i=0;iN;i+) j = (i + 1) % N;area +=polyg on i.x *polygonj.y;area -=polygoni.y *polygonj.x;area /= 2;return(area PI)dtheta -= PI*2;while (dtheta -PI)dtheta += PI*2;return(dtheta);4. 兩點距離(2D

26、 3D)語法：result=distance_2d(floatx1,float x2,floaty1,float y2);參數(shù)：x/y/z12： 各點的x、y、z坐標返回值： 兩點之間的距離返回角度為弧度制，并且 以逆時針方向為正方向需要math.h源程序：需要math.h源程序：floatdistance_2d(floatx1,float x2,floaty1,float y2)return(sqrt(x1-x2)*(x1- x2)+(y1-y2)*(y1-y2);float distance_3d(floatx1,float x2,floaty1,float y2,floatz1,floa

27、t z2)return(sqrt(x1-x2)*(x1- x2)+(y1 -y2)*(y1 -y2)+(z1- z2)*(z1-z2);5. 射向法判斷點是否在多邊 形內部 語法：result=insidepolygon(Poi nt *polygon,int N,PointP);參數(shù)：*polygon :多邊形頂點數(shù)組N:多邊形頂點個數(shù)p：被判斷點返回值：0:點在多邊形內部；1:點在多邊形外部若p點在多邊形頂點或者邊上，返回值不確定，需另行判斷需要math.h源程序：#define MIN(x,y) (x y ?x : y)typedef struct double x,y; Point;i

28、nt insidepolygon(Point *polygon,int N,Point p) int counter = 0;int i;double xinters;Point p1,p2;pl = polygon0;for (i=1;iMIN(p1.y,p2.y) if (p.y =MAX(p1.y,p2.y) if (p.x=MAX(p1.x,p2.x) if(p.y-p1.y)*(p2.x-p1.x)/(pi.y != p2.y) xinters =p2.y-p1.y)+p1.x;if (pl.x = p2.x | p.x =xinters)counter+;pl = p2；if (c

29、ounter % 2 = 0) return(OUTSIDE); else return(INSIDE);p1、 p2: p：返回值：6. 判斷點是否在線段上語法：result=Pointonline(Pointp1,Point p2,Point p);參數(shù)：線段的兩個端點 被判斷點0：點在不在線段上；1:點在線段上汪意：若p線段端點上返回需要math.h源程序：#define MIN(x,y) (x y ?x : y)typedef struct double x,y; Point;int FC(double x1,doublex2)if(x1- x2- 0.000002) return 1

30、; else return 0;int Pointonline(Pointp1,Point p2,Point p)double x1,y1,x2,y2;x1=p.x-p1.x;x2=p2.x-p1.x;y1=p.y-p1.y;y2=p2.y-p1.y;if(FC(x1*y2 -x2*y1,0)=0)return 0;if(MIN(p1.x,p2.x)v=p.x&p .xv=MAX(p1.x,p2.x)&x : y)return 0;(MIN(p1.y,p2.y)v=p. y&p.yv=MAX(p1.y,p2.y)return 1; else7. 判斷兩線段是否相交語法：result=secti

31、ntersect(Point p1,Point p2,Pointp3,Point p4);參數(shù)：p14:兩條線段的四個端 占八、返回值：0:兩線段不相交；1:兩線段相交；2兩線段首尾相接注意：p1!=p2;p3!=p4;源程序：#define MIN(x,y) (x y ?x : y)typedef struct double x,y; Point;int lineintersect(Pointp1,Point p2,Pointp3,Point p4)Point tp1,tp2,tp3;if(p1.x=p3.x&p1.y=p3.y )|(p1.x=p4.x&p1.y=p4 .y)|(p2.x=

32、p3.x&p2.y= p3.y)|(p2.x=p4.x&p2.y=p4.y)return 2;/快速排斥試驗if(MIN(p1.x,p2.x)vp3.x&p 3.xMAX(p1.x,p2.x)&MIN( p1.y,p2.y)vp3.yMAX(p1.y ,p2.y)|(MIN(p1.x,p2.x)vp4.x&p3 .xMAX(p1.x,p2.x)&MIN(p 1.y,p2.y)vp3.y=0) return 1; elsereturn 0;8. 判斷線段與直線是否相交語法：result=lineintersect(Poi nt p1,Point p2,Point p3,Point p4);參數(shù)：

33、pl、p2：線段的兩個端點p3、p4:直線上的兩個點返回值：0:線段直線不相交；1:線段和直線相交返回值： 點q到線段p1p2的距離如線段在直線上，返回1源程序：typedef struct double x,y; Point;int lineintersect(Point p1,Point p2,Point p3,Point p4)Point tp1,tp2,tp3; tp1.x=p1.x-p3.x; tp1.y=p1.y-p3.y; tp2.x=p4.x-p3.x; tp2.y=p4.y-p3.y; tp3.x=p2.x-p3.x; tp3.y=p2.y-p3.y;if(tp1.x*tp2

34、.y-tp1.y*tp2.x)*(tp2.x*tp3.y-tp2.y*tp3.x)=0) return 1; elsereturn 0;9點到線段最短距離語法：result=mindistance(Pointp1,Point p2,Point q);參數(shù)：pl、p2： 線段的兩個端點q： 判斷點需要math.h源程序：#define MIN(x,y) (x y ?x : y)typedef struct double x,y; Point;double mindistance(Point p1,Point p2,Point q) int flag=1;double k;Point s;if (

35、p1.x=p2.x)s.x=p1.x;s.y=q.y;flag=0;if (p1.y=p2.y)s.x=q.x;s.y=p1.y;flag=0;if (flag)k=(p2.y-p1.y)/(p2.x-p1.x);s.x=(k*k*p1.x+k*(q.y-p1. y)+q.x)/(k*k+1);s.y=k*(s.x-p1.x)+p1.y;if(MIN(p1.x,p2.x)v=s.x&s. xv=MAX(p1.x,p2.x) returnsqrt(q.x-s.x)*(q.x-s.x) +(q.y-s.y)*(q.y-s.y); else returnMIN(sqrt(q.x-p1.x)*(q.x

36、-p1.x)+(q.y-p1.y)*(q.y-p1. y),sqrt(q.x-p2.x)*(q .x-p2.x)+(q.y-p2.y)*(q.y -p2.y);10. 求兩直線的交點語法：result=mindistance(Pointp1,Point p2,Point q);參數(shù)：p1p4:直線上不相同的兩點*p :通過指針返回結果返回值：1:兩直線相交;2:兩直線平行如需要判斷兩線段交點，檢驗k和對應k1 (注釋中)的值是否在01之間，用在01之間的那個求交點源程序：typedef struct double x,y; Point;int linecorss(Pointp1,Point p

37、2,Pointp3,Point p4,Point *p)double k;/同一直線if(p4.x-p3.x)*(p1.y-p3.y) -(p4.y-p3.y)*(p1.x-p3.x) =0& (p2.x-p1.x)*(p1.y-p3.y)- (p2.y-p1.y)*(p1.x-p3.x)=0) return 2;/平行，不同一直線if(p4.y-p3.y)*(p2.x-p1.x)-(p4.x-p3.x)*(p2.y-p1.y)=0) return 0;k=(p4.x-p3.x)*(p1.y-p3.y)-(p4.y-p3.y)*(p1.x-p3.x)/(p4.y-p3.y)*(p2.x-p1.

38、x)-(p4.x-p3.x)*(p2.y-p1.y);k1=(p2.x-p1.x)*(p1.y-p3.y)-(p2.y-p1.y)*(p1.x-p3.x)/(p4.y-p3.y)*(p2.x-p1.x)-(p4.x-p3.x)*(p2.y-p1.y);(*p).x=p1.x+k*(p2.x-p1.x)；(*p).y=p1.y+k*(p2.y-p1.y)；return 1;/ 有交點11. 判斷一個封閉圖形是凹集還是凸集語法：result=convex(Point*p :n：返回值:集；0：*p,int n);參數(shù)：封閉曲線頂點數(shù)組封閉曲線頂點個數(shù)1：凸集；-1 :凹曲線不符合要求無法計算默認曲

39、線為簡單曲線：無交叉、無圈源程序：typedef struct double x,y; Point;int convex(Point *p,int n)int i,j,k;int flag = 0;double z;if (n 3) return(0);for (i=0;i n ;i+) j = (i + 1) % n;k = (i + 2) % n;z = (pj.x -pi.x) * (pk.y - pj.y); z -= (pj.y - pi.y) * (pk.x - pj.x);if (z 0)flag |= 2; if (flag = 3) return/CONCAVEif (fla

40、g != 0)return 1;/CONVEXelsereturn 0;-1;12. Graham掃描法尋找凸包 語法：Graham_scan(Point PointSet,Point ch,int n,int &len);參數(shù)：PointSet: 輸入的點集 ch:輸出的凸包上的點集，按照逆時針方向排列n:PointSet中的點的數(shù)目len :輸出的凸包上的點的個數(shù)返回值：null源程序：struct Point float x,y;；float multiply(Pointp1,Point p2,Point p0)return(p1.x-p0.x)*(p2.y -p0.y)-(p2.x-p0.x)*(p1.y -p0.y);float distance(Point p1,Point p2) return(s