| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1352 人关注过本帖, 2 人收藏
标题:请第一个正确完成倒数和高精度计算的兄弟czz5242199接分
取消只看楼主 加入收藏
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
结帖率:100%
收藏(2)
已结贴  问题点数:100 回复次数:2 
请第一个正确完成倒数和高精度计算的兄弟czz5242199接分
呵呵,做人要讲诚信。既然czz5242199第一个正确完成任务,那一百分当之无愧。

关于算法就不多说什么了,这里我想聊聊关于为了得到小数点后K位的精确值,实际运算中至少应该保留多少位。

假设我们实际计算中使用了P位,P位之后的值舍去,那么它的误差将小于10^P。

用这样的值做加法时,误差会线性累加。对于前100个倒数的和来说,累积误差将小于100 * 10^P,即10^(p - 2)。

那么之前的P-3位应该是精确的。

由此可知,如果想要1000位的精确值,实际计算中至少应该使用1003位来计算。

下面是我关于这一任务的代码,供各位参考。
程序代码:
#include<stdio.h>
#define BASE    100000
#define LENGTH    202
void show(int *a, int n)
{
    int i;
    printf("%d.", a[0]);
    for(i = 1; i < n; printf("%05d ", a[i++])) if(i % 10 == 1) putchar('\n');
}
void getInv(int *a, int n, int b)
{
    int i, r = 1;
    for(i = 0; i < n; i++)
    {
        a[i] = r / b;
        r = r % b * BASE;
    }
}
void add(int *a, int *b, int n)
{
    int i, f = 0;
    for(i = n - 1; i >= 0; i--)
    {
        a[i] += b[i] + f;
        if(a[i] >= BASE)
        {
            a[i] -= BASE;
            f = 1;
        }
        else f = 0;
    }
}
int main()
{
    int a[LENGTH] = {1}, b[LENGTH], i;
    for(i = 2; i <= 100; i++)
    {
        getInv(b, LENGTH, i);
        add(a, b, LENGTH);
    }
    show(a, 201);
    return 0;
}
搜索更多相关主题的帖子: 小数点 倒数 兄弟 
2011-11-02 13:07
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
回复 6楼 DefineFeel
用数组的一个元素表示一个十进制位的方法你应该理解的吧?
一次计算一位效率有点低,所以我的代码一次计算5位。这样既便于运算,也便于按我想要的格式输出。
其实可以理解为我是将数以十万进制来运算的。当然也可以以其实进制来运算,这题以1千万进制计算最经济,但以5位输出时不太方便。
另外我对于数是以定点小数的方式来处理的。a[0]为小数点前十万进制的个位,其它元素为小数点后的位。
具体该使用定点数还是浮点数,就具体问题具体分析。本题非常适合用定点数来处理。

重剑无锋,大巧不工
2011-11-02 21:41
beyondyf
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:贵宾
威 望:103
帖 子:3282
专家分:12654
注 册:2008-1-21
收藏
得分:0 
怎么把这贴子给翻出来了。
czz5242199这兄弟的水平也很不错,是值得多交流的朋友。

重剑无锋,大巧不工
2011-12-04 12:10
快速回复:请第一个正确完成倒数和高精度计算的兄弟czz5242199接分
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.033028 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved