算法竞赛入门经典-第四章源代码

// 程序4-1 组合数
#include <stdio.h>

int f(int n)
{
    int i, m = 1;
    for(i = 1; i <= n; i++)
        m *= i;
    return m;
}

int main(void)
{
    int m, n;
    scanf("%d%d", &m, &n);
    printf("%d\n", f(n)/(f(m)*f(n-m)));
    return 0; // 如何避免中间结果溢出？
}

// 程序4-3 孪生素数
#include <stdio.h>
#include <math.h>
#include <assert.h>
int is_prime(int x)	// 惯用命名
{
    int i, m;
    assert(x >= 0); //当括号内不成立时，程序将异常终止，并给出提示信息
    if(x == 1)
        return 0; //特殊情况单独处理
    m = floor(sqrt(x) + 0.5); //浮点计算四舍五入误差，避免重复计算
    for(i = 2; i <= m; i++)
        if(x % i == 0)
            return 0;
    return 1;
}
int main(void)
{
    int i, m;
    scanf("%d", &m);
    for(i = m-2; i >= 3; i--)	// 不能为i=m，考虑边界
        if(is_prime(i) && is_prime(i+2)){
            printf("%d %d\n", i, i+2);
            break;
        }
    return 0;
}

// 程序4-6 用递归法计算阶乘
#include <stdio.h>

int f(int n)
{
    return n==0 ? 1:f(n-1)*n;
}

int main(void)
{
    printf("%d\n", f(10));
    return 0;
}

//(1)ax+by=c, dx+ey=f
#include <stdio.h>
#include <assert.h>

double fucx(double a,double b,double c,double d,double e,double f)
{
	assert((a*e-b*d)!=0);
	return (c*e-b*f)/(a*e-b*d);
}

double fucy(double a,double b,double c,double d,double e,double f)
{
	assert((a*e-b*d)!=0);
	return (a*f-c*d)/(a*e-b*d);
}

int main(void)
{
	double a,b,c,d,e,f;
	double x,y;
	scanf("%f%f%f%f%f%f",&a,&b,&c,&d,&e,&f);
	x=fucx(a,b,c,d,e,f);
	y=fucy(a,b,c,d,e,f);
	printf("%f %f\n",x,y);
	return 0;
}

//(2)f() g() h()
#include <stdio.h>

int f()
{
	static int i = 0;
	return ++i;
}
int g(){return 2;}
int h(){return 3;}

int main(void)
{
	int a, b;
	a = f();
	b = f();
	printf("%d %d\n", a, b);
	return 0;
}

//(3)
#include <stdio.h>
int a = 1;

int* f()
{
	int x = 1;
	return &x;
}

int main(void)
{
	int a = 2;
	printf("%d\n", a);
	a++;
	printf("%d\n", a);
	int* p = f();
	printf("%d\n", *p);
	return 0;
}