C语言常用数值计算算法(素数、公约数、级数、方程根和定积分)

素数判断

#include<stdio.h> #include<math.h> int main() { int n,min,max,isprime; scanf('%d %d',&min,&max); if(min<=2){ printf('%4d',2); min=3; } if(min%2==0) min++; for(n=min;n<=max;n+=2){ for(isprime=1,i=2;t&&i<=sqrt(n);i++) if(n%i==0) isprime=0; if(isprime) printf('%4d',n); } return 0; }

最大公约数

1.brute-force算法

#include<stdio.h>
int main()
{
    int x=30,y=45,z;
    z=x;
    while(!(x%z==0&&y%z==0))
        z--;
    printf('%d',z);
    return 0;
}

2.欧几里得算法

#include<stdio.h> int main() { int x=35,y=45,r; while((r=x%y)!=0){ x=y; y=r; } printf('%d',y); return 0; }

穷举法

例 解方程: ①x+y+z=100 ②5x+3y+z/3=100

#include<stdio.h>
int main()
{
    int x,y,z;
    for(x=0;x<=20;x++)
        for(y=0;y<=33;y++)
            for(z=0;z<=100;z++)
                if(x+y+z==100&&5*x+3*y+z/3==100&&z%3==0)
                    printf('x=%d,y=%d,z=%d\n');
    return 0;
}

#include<stdio.h> int main() { int x,y,z; for(x=0;x<=20;x++) for(y=0;y<=33;y++){ z=100-x-y; if(5*x+3*y+z/3==100&&z%3==0) printf('x=%d,y=%d,z=%d\n',x,y,z); } return 0; }

级数近似

#include<stdio.h>
#include<math.h>
int main()
{
    double s=1,a=1,x,eps;
    int n;
    scanf('%lf%lf',&x,&eps);
    for(n=2;fabs(a)>eps;n++){
        a=a*x/(n-1);
        s+=a;
    }
    printf('%f',s);
    return 0;
)

#include<stdio.h> #include<math.h> int main() { double sum,x,eps=1e-6,fn,tn; int s=1,n=2; scanf('%lf',&x); s=fn=x; do{ s=-s; fn=fn*(2*n-3)/(2*n-2)*x*x; tn=sign*fn/(2*n-1); sum=sum+tn; n++; }while(fabs(tn)>eps); printf('%f',sum);

一元非线性方程求根

一、牛顿迭代法

　　1.基本概念：如果函数连续，且待求零点孤立，那么在零点周围存在一个区域，当初值在这个邻域内时，牛顿法收敛。如果零点不为0, 那么牛顿法将具有平方收敛的性能。

　　2.基本公式：x_k+1=x_k-f(x_k)/f'(x_k)

　　3.判断条件：|f(x_n+1)|<ε或|x_n+1-x_n|<ε是否为真。若为真则x_n+1就是方程f(x)=0在x₀附近误差ε范围内的一个近似根。

　　4.实际应用：求cos(x)-x=0的近似解，精确到10^-6。

#include<stdio.h>
#include<math.h>

int main()
{
    double x1=0, x2 = 2;
    while (fabs(cos(x2)-x2)>1e-6&&fabs(x2-x1)>1e-6){
        x1 = x2;
        x2 = x1 + (cos(x1) - x1) / (sin(x1) + 1);
    }
    printf('x=%f\nf(x)=%.6f',x2,fabs(cos(x2)-x2));
    return 0;
}

二、二分法

　　1.基本概念：对于区间[a，b]上连续不断且 f(a)f(b)<0的函数y=f(x) ，通过不断地把函数 f(x) 的零点所在的区间一分为二，使区间的两个端点逐步逼近零点，进而得到零点近似值的方法叫二分法。

　　2.实际应用：求exp(x)+x=0在(-1, 0)的根，精确到10^-6。

#include<stdio.h> #include<math.h> double f(double x); int main() { double a=-1, b=0, c; c = (a+b)/2; do{ if(f(a)*f(c)>0) a = c; else b = c; c = (a+b)/2; } while (fabs(f(c)) > 1e-6&&fabs(a-b)>1e-6); printf('x=%.6f', c); return 0; } double f(double x) { return exp(x) + x; }

三、弦截法

　　1.基本概念：弦截法是求非线性方程近似根的一种线性近似方法。它是以与曲线弧AB对应的弦AB与x轴的交点横坐标作为曲线弧AB与x轴的交点横坐标的近似值μ来求出方程的近似解。

　　2.实际应用：求（(x+2)*x-2)*x-1=0在(-1, 0)的根，精确到10^-6。

#include <math.h>
#include <stdio.h>
float f(float x)
{
    float y;
    y = ((x + 2.0) * x - 2.0) * x - 1.0;
    return y;
}

int main()
{
    float x, x1, x2;
    x1 = -1, x2 = 0;
    do{
        x = (x1 * f(x2) - x2 * f(x1)) / (f(x2) - f(x1));
        if (f(x) * f(x1) > 0)
            x1 = x;
        else
            x2 = x;
    } while (fabs(f(x)) >= 1e-6);
    printf('the root is %f\n', x);
    printf('((x+2.0)*x-2.0)*x-1.0=%f\n', f(x));
    return 0;
}

定积分近似计算

1.梯形法

　　算法思想

　　程序实现

#include<stdio.h> #include<math.h> int main() { long n, i; double a=-1.57,b=1.57,h,s,x; scanf('%ld',&n); h=(b-a)/n; s=(cos(a)+cos(b))/2; x=a; for(i=1;i<n;i++){ x+=h; s+=cos(x); } printf('%f',s*h); return 0; }

2.矩形法

　　算法思想

　　程序实现

#include<stdio.h>
#include<math.h>
int main()
{
    long n,i;
    double a=-1.57,b=1.57,h,s=0,x;
    scanf('%ld',&n);
    h=(b-a)/n;
    x=a;
    for(i=0;i<n;i++){
        s+=cos(x);
        x+=h;
    }
    printf('%f',s*h);
    return 0;
}