2. 你会求多项式吗

问题

在这里插入图片描述
如何编程求这个多项式?

code

很简单,我+再+再+不就出来了。

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

#define MAXK 10000000 //被测函数最大调用次数
#define MAXN 101 //多项式的最大项数,即阶数 + 1

clock_t start, stop; // clock_t是clock()函数返回的变量类型
double duration; //记录被测函数运行事件,以秒为单位

double f1(int n, double a[], double x)
{
    double p = a[0];

    for (int i = 1; i <= n; i++)
    {
        p += a[i] * pow(x, i);
    }
    return p;
}

double f2(int n, double a[], double x)
{
    double p = a[n];

    for (int i = n - 1; i >= 0; i--)
    {
        p = a[i] + p * x;
    }

    return p;
}

void printTime(int n, double a[], double x, bool flag)
{
    start = clock();

    for (int i = 0; i < MAXK; i++)
    {
        if (flag)
        {
            f1(n, a, x);
        }
        else
        {
            f2(n, a, x);
        }
    }

    stop = clock();
    duration = ((double)(stop - start)) / CLK_TCK / MAXK;       /* 计算单词执行用时 */

    printf("ticks = %lf\n", (double)(stop - start));
    printf("duration = %6.2e\n", duration);
}

int main()
{
    double a[MAXN];

    for (int i = 0; i < MAXN; i++)
    {
        if (i == 0)
        {
            a[i] = 1.0;
        }
        else
        {
            a[i] = (double)(1 / i);
        }
    }

    printf("f1:\n");
    printTime(MAXN - 1, a, 1.1, true);

    printf("f2:\n");
    printTime(MAXN - 1, a, 1.1, false);
}

结果

看下运行时间
在这里插入图片描述
差这么多?是不是惊呆了。

再来一个多项式

在这里插入图片描述

code

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

#define MAXK 10000000 //被测函数最大调用次数
#define MAXN 101 //多项式的最大项数,即阶数 + 1

clock_t start, stop; // clock_t是clock()函数返回的变量类型
double duration; //记录被测函数运行事件,以秒为单位

double f1(int n, double x)
{
    double p = (double)1;
    for (int i = 1; i <= n; i++)
    {
        p += pow(x, i) / i;
    }
    return p;
}

double f2(int n, double x)
{
    double p = (double)(1 / n);
    for (int i = n; i > 1; i--)
    {
        p += (1 / (i - 1)) + (x * p);
    }
    return p;
}

void printTime(int n, double x, bool flag)
{
    start = clock();

    for (int i = 0; i < MAXK; i++)
    {
        if (flag)
        {
            f1(n, x);
        }
        else
        {
            f2(n, x);
        }
    }

    stop = clock();
    duration = ((double)(stop - start)) / CLK_TCK / MAXK;       /* 计算单词执行用时 */

    printf("ticks = %lf\n", (double)(stop - start));
    printf("duration = %6.2e\n", duration);
}

int main()
{
    printf("f1:\n");
    printTime(MAXN - 1, 1.1, true);

    printf("f2:\n");
    printTime(MAXN - 1, 1.1, false);
}

结果

在这里插入图片描述
差这么多就离谱

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值