文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

一篇文章带你了解论C语言中算法的重要性

2024-04-02 19:55

关注

一、问题一(打印阶乘)

问题描述:

打印出数字一到数字20的阶乘

一开始,我总会多打印出一个1,这令我十分苦恼,并且从n等于13开始,数据就开始溢出

问题分析:

让我们分析一下问题,这里面存在着两个问题:

1.多打印出一个1

2.数据溢出

解决方案:

1.让我们检查一下结果,发现问题很有可能是循环的时候没有循环本身

for (i = 1; i < num; i++)//这句话明显错了

改成


for (i = 1; i <= num; i++) {//i的值要乘以它本身!
			n = n * i;
		}

2.这里要引入C++中STL库的一个知识点

常规的32位整数只能够处理40亿以下的数。

如果遇到比40亿要大的数,就要用到C++的64位扩展。不同的编译器对64位整数的扩展有所不同。这个我也是听别人科普的,大家可以站内搜索一下。

优化后的代码如下:


#include <stdio.h>
void main() {
	__int64 fac(int num);
	int n = 1;
	int num;
	for (num = 0; num <= 20; ++num) {
		printf("%3d! = %I64d\n", num, fac(num));
	}
}
__int64 fac(int num) {
	register __int64 n = 1, i;     //寄存器变量
		for (i = 1; i <= num; i++) {//i的值要乘以它本身!
			n = n * i;
		}
		return n;
}

二、问题二(比较多项式计算时间)

问题描述:

这里先科普几个测试代码中的知识点:

这个表示本程序开始计时:

start = clock();

本程序结束计时:

stop = clock();

clock tick :时钟打点

CLK_TCK:机器时钟每秒所走的时钟打点数

问题分析:

首先这个问题有两种算法:

直接算

p1 += (pow(x, i)/i);

把x当成公因式提出计算(秦九韶法)

p2 = 1.0/(a[i - 1])+ (x*p2);

然后我发现了三个问题:

1.测量不出时间

2.程序重复性高

3.第一种结果和第二种结果不一致

解决方案:

1.让被测函数重复运行多次,使得测出的总时钟打点间隔充分长,最后计算被测函数除以运行次数即可得出平均每次的运行时间


duration = ((double)(stop - start)) / CLK_TCK / MAXK;

2.可以通过多设置几个函数,并调用函数解决问题

3.这是算法的问题

这个地方真的特别容易出错,我改了不知道多少遍。。。。。。


double f2(int n, double a[], double x) {
	int i;
	double p2 = 1.0/a[n];
	for (i = n; i > 0; i--) {
		p2 = 1.0/(a[i - 1])+ (x*p2);//算法思路出毛病了(数学)
	}
	return p2;
}

总体的代码:


#include <stdio.h>
#include <math.h>
#include <time.h>
clock_t start, stop;
double duration;
#define MAXN 101//数组里元素个数(多项式的系数),如果看n值需要减一,因为有a0
#define MAXK 1000//重复调用的次数
double f1(int n, double a[], double x)
{
	double p1 = a[0];//a[0]都已经算出来了,循环从1开始
	for (int i = 1; i <= n; i++) {
		p1 += (pow(x, i)/i);
	}
	return p1;	
}
double f2(int n, double a[], double x) {
	int i;
	double p2 = 1.0/a[n];
	for (i = n; i > 0; i--) {
		p2 = 1.0/(a[i - 1])+ (x*p2);//算法思路出毛病了(数学)
	}
	return p2;
}
double ceshijian()
{
	stop = clock();//停止计时
	duration = ((double)(stop - start)) / CLK_TCK / MAXK;//计算单次运行时间
	printf("ticks=%f\n", (double)(stop - start));
	printf("duration=%6.2e\n", duration);
	return 0;
}
int main()
{
	int i;
	double a[MAXN];
	for (i = 0; i < MAXN; i++) {
		a[i] = (double)i;
	}//输入的早就是i值了
	a[0] = 1;
	//不在测试范围内的准备工作写在clock()调用之前
	start = clock();//开始计时
	for (int i = 0; i < MAXK; i++)//重复调用
		f1(MAXN - 1, a, 1.1);//被测函数,这里如果写数组的话就越界了,而且要调用某个值是不确定的,只能写a,因为要定义的就是a值
	printf("第一种结果为%f\n", f1(MAXN - 1, a, 1.1));
	ceshijian();

	start = clock();//开始计时
	for (i = 0; i < MAXK; i++)
		f2(MAXN - 1, a, 1.1);//被测函数,这里如果写数组的话就越界了,而且要调用某个值是不确定的
	printf("第二种结果为%f\n", f2(MAXN - 1, a, 1.1));
	ceshijian();
	return 0;
}

结果如下

总结

本篇文章就到这里了,希望能给你带来帮助,也希望您能够多多关注编程网的更多内容!

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     813人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     354人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     318人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     435人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     224人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯