[理学]C语言程序设计循环控制.ppt

例：搬砖问题:36块砖，36人搬。男搬4，女搬3，童子2搬1。问男、女、童各几何？ 解：man+woman+child=36 4* man+3*woman+child/2=36 main() //见c3_1_10.cpp { int man,woman,child; man=0; while (man=9) { woman=0; while (woman=12) { child=36-woman-man; If (4*man+3*woman+child/2.0==36.0) printf(“men=%d,women=%d,children=%d\n”); } woman++; } man++; } } 例 打印九九乘法表 算法： 打印表头； for (i=1,i=9;i++) { 打印i; for (j=1,j=9;j++) 打印i*j; 输出换行； }//见c3_1_14.cpp 例6.6 利用 求π 解：令ak=(-1)k-1/(2k-1),则以上求和到| ak |eps为止(eps=10-7) #define eps 1.0e-7 //tc3_10.c main() { int s; float n,t,pi; t=1.0;pi=0;n=1.0;s=1; while(fabs(t)=eps) { pi+=t; n+=2.0; s=-s; t=s/n; } pi=pi*4; printf(\n pi=%lf\n“pi); } 例 求正整数u,v的最大公约数 //见c3_1_12.cpp main() { int x,y,u,v,r; printf(input u,v:); scanf(%d,%d,x,y); u=x;v=y; if ((u=0)||(v=0)) printf(\n数据输入错！\n); else { while (v) { r=u%v; u=v; v=r; } printf(\ngcd(%d,%d)=%d\n,x,y,u); } } 例 求 sinx-x/2=0 在[π/2， π]中的根。 //见c3_1_11.cpp #define eps 0.0000001 #define pai 3.1415926 #define f(x) sin(x)-x/2.0 void main(void) //二分法求方程的根 {double x0, x1,f0,f1,f2,x2; int i=0; x0=pai/2.0; x1=pai; f0=f(x0);f1=f(x1); while (fabs(x0-x1)=eps) {i++; x2=(x0+x1)/2.0; f2=f(x2); if (fabs(f2)eps) x0=x1=x2; else if (f2*f00) {x1=x2;f1=f2; } else {x0=x2;f0=f2;} } printf(i=%d,root=%lf\n,i,x2); } 例6.9 求100~200之间的素数 #define N1 100 #define N2 200 main() //见c3_1_13.cpp { int p,k,c=0; if(N1%2)p=N1; //保证p从奇数开始 else p=N1+1; for(; p=N2;p+=2) { for(k=2;(k*k=p) (p%k!=0); k++); if(k*kp) //p为素数 { printf( %d,,p); c++; //素数个数计数 if (c%10==0) //每10个换行 { printf(\n); c=0; } } } } 例、求定积分的梯形法。 令 x0=a xi=a+i*h, i=1,..n; h=(b-a)/n; ti=(f(a+(i-1)*h)+f(a+i*h))*h/2; 例 求 main() //tc3_18.c { float x,s,t; int i,n; scanf(“%f,%d”,x,n); for (t=x,s=0.0,i=1;i=n;i++) { s+=t; t