博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU1012 Calculate e
阅读量:5888 次
发布时间:2019-06-19

本文共 1778 字,大约阅读时间需要 5 分钟。

题目地址:http://acm.hdu.edu.cn/showproblem.php?pid=1012
解法一:最直观的方法,就是递归计算来求各项之和
#include <iostream>
#include <iomanip>
using namespace std;
double curItem(int n)
{//当前子项
    if (n==0)
    {
        return 1;
    }
    else
        return curItem(n-1)/n;
}
double sum(int n)
{
    if(n==0)
        return 1;
    else 
        return sum(n-1)+curItem(n);
}
void caculateE(int n)
{//输出指定n下e的值
    double tmp1,tmp2=1.0f,result=0.0f;
    int i;
    result = sum(n);
    if(n==0||n==1)
    {
        cout<<n<<" "<<static_cast<int>(result)<<endl;
    }
    else if(n==2)
    {
        cout<<setiosflags(ios::fixed)<<setprecision(1);
        cout<<n<<" "<<result<<endl;
    }
    else
    {
        cout<<setiosflags(ios::fixed)<<setprecision(9);
        cout<<n<<" "<<result<<endl;
    }
}
int main(int argc, char *argv[])
{
    int n,i;
    cout<<"n e"<<endl;
    cout<<"- -----------"<<endl;
    for(i=0;i<=9;++i)
    {
        caculateE(i);
    }
    cin>>i;
    return 0;
}
解法二:题中给出的计算e的式子是由e^x的泰勒级数展开而得,在计算之前可以使用个技巧,就是把它们叠乘起来改写成:
e=(1+(1+1/2(1+1/3(1+1/4(1+…1/(n-1)(1+1/n))))),从最里面的括号往外算,共做n次除法和加法得一段结果,运算效率也是O(N*M),但是由于收敛速度快些,所以N项节省一些,
#include <iostream>
#include <iomanip>
using namespace std;
double doCaculate(int n)
{//实际的计算
    double tmp1,tmp2=1.0f,result=0.0f;
    int i;
    for(i=n;i>=1;--i)
    {
        tmp1 = static_cast<double>(1)/static_cast<double>(i);
        tmp2 = 1.0f+tmp1*tmp2;
    }
    return tmp2;
}
void caculateE(int n)
{//输出指定n下e的值,主要是输出格式的处理
    double result = 0.0f;
    result = doCaculate(n);
    if(n==0||n==1)
    {
        cout<<n<<" "<<static_cast<int>(result)<<endl;
    }
    else if(n==2)
    {
        cout<<setiosflags(ios::fixed)<<setprecision(1);
        cout<<n<<" "<<result<<endl;
    }
    else
    {
        cout<<setiosflags(ios::fixed)<<setprecision(9);
        cout<<n<<" "<<result<<endl;
    }
}
int main(int argc, char *argv[])
{
    int n,i;
    cout<<"n e"<<endl;
    cout<<"- -----------"<<endl;
    for(i=0;i<=9;++i)
    {
        caculateE(i);
    }
    return 0;
}
本文转自Phinecos(洞庭散人)博客园博客,原文链接:http://www.cnblogs.com/phinecos/archive/2007/12/25/1013925.html,如需转载请自行联系原作者
你可能感兴趣的文章
白话讲反射技术 --- 适合初学者入门引导
查看>>
css变形 transform
查看>>
win7家庭版添加组策略编辑器
查看>>
iOS 一个开发者账号 多台Mac 共用
查看>>
lnmp环境搭建
查看>>
3.JUC之volatile
查看>>
oracle:win7手工卸载oracle数据库11g
查看>>
自定义session扫描器精确控制session销毁时间--学习笔记
查看>>
基于busybox搭建功能完善的小型linux(一)
查看>>
android The project target (Android 2.3.3) was not properly loaded
查看>>
【转】EDK简单使用流程(3)
查看>>
loj#2538. 「PKUWC2018」Slay the Spire
查看>>
这是一篇被河蟹了的博客
查看>>
一个两年Java的面试总结
查看>>
软件工程项目组Z.XML会议记录 2013/11/27
查看>>
科学计算库学习报告
查看>>
软件测试 -- 软件测试的风险主要体现在哪里
查看>>
Ubuntu中无法update的解决办法
查看>>
仿射变换
查看>>
decltype类型指示符
查看>>