鱼C论坛

 找回密码
 立即注册
查看: 1502|回复: 13

[已解决]数据存储问题求助

[复制链接]
发表于 2021-10-29 14:51:35 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
各位大佬们,突然有一个问题,
C语言如何存储超大的数据,
就比如计算100!(100的阶乘);

怎样存储结果,并把它打印出来?
最佳答案
2021-10-29 15:24:31
用其他库函数,例如 gnu 的 mpfr 库
The Multiple Precision Floating-Point Reliable Library

google 翻译
多精度浮点可靠库 

官方文档: https://www.mpfr.org/mpfr-4.0.1/mpfr.pdf

下面计算 10^1000000
^_^
#include <stdio.h>
#include <mpfr.h>

int main(void) {
    mpfr_t a, b, c;
    mpfr_inits2(10000000, a, b, c, NULL);
    mpfr_set_str(a, "10", 10, MPFR_RNDD);
    mpfr_set_str(b, "1000000", 10, MPFR_RNDD);
    mpfr_pow(c, a, b, MPFR_RNDD);
    mpfr_printf("%.0Rf\n", c);
    mpfr_clears(a, b, c, NULL);
    mpfr_free_cache();
    return 0;
}
gcc -g -Wall -o main main.c $(pkg-config --cflags --libs mpfr)
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2021-10-29 15:07:53 | 显示全部楼层
double
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 15:11:50 | 显示全部楼层
本帖最后由 村里小黑 于 2021-10-29 15:12 编辑
#include<stdio.h>
#define Ma 10000   //a的最大容量,必须大于na
int pa=0;//指向数组a的有效末端
int p=2;//求阶乘时的当前乘数
unsigned int carry=0;
int memory_over=0;
 
union data
{   //表示范围0--2^32-1
        unsigned long int b;//这里也可以用数组usigned long int ,但是用公用体可以存储不同类型的数据
}a[Ma];
 
void main()
{
        unsigned int n;//求n的阶
        void facto(unsigned int n);
        printf("Input n:");
        scanf("%u",&n);  
        
        a[0].b=1;//初始化
        facto(n);//测试后至少可以计算10000!
        
        if(memory_over==0)
        {
                printf("the result include %dNO:\n",pa+1);//高位输出
                printf("%u",a[pa--].b);
                for(;pa>=0;pa--)
                {
                        printf("%04u",a[pa].b);//04%d用来补充0限制每个int型中存放4位数
                }
                printf("\n");
        }
        getchar();
}
 
 
 
 
void facto(unsigned int n)
{
        void multiple();
        pa=0;        
        while(pa< Ma-1 && p<=n)//容量限制
        {
                multiple();
                p++;//每一轮乘一个阶数p
        }
        if(p<=n)//此时pa>=Ma-1退出while循环说明超出数组所能存储的范围
        {
                printf("memory out!\n");
                memory_over=1;
        }//如果当前的存储结果的数组a[Ma]不够用!应提高Ma
        
}
 
void multiple()
{
        int i=0;
        carry=0;
        while(i<=pa)//i指向当前处理的元素a[i],每一轮用一个位与阶数p相乘
        {
                a[i].b=a[i].b*p+carry;//计算结果,要考虑来自低位的进位
                carry=a[i].b/10000;//计算进位
                a[i].b=a[i].b%10000;//计算余数
                i++;
        }
        //储代码块后carry=0  
        if(carry>0)
                a[++pa].b=carry;
}
 
原文链接
https://blog.csdn.net/kai8wei/article/details/45110859
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 15:12:22 | 显示全部楼层
#include <stdio.h>

double f(int n){
    double res = 1;
    for(int i = n; i > 0; i--){
        res *= i;
    }
    return res;
}

int main()
{
    printf("%.0lf", f(100));
    return 0;
}
93326215443944175354307254139643190247129328132295862491935879110669343325734178368282822618707234467717279847537548956702435362278960753539491860335688679424
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 15:24:31 | 显示全部楼层    本楼为最佳答案   
用其他库函数,例如 gnu 的 mpfr 库
The Multiple Precision Floating-Point Reliable Library

google 翻译
多精度浮点可靠库 

官方文档: https://www.mpfr.org/mpfr-4.0.1/mpfr.pdf

下面计算 10^1000000
^_^
#include <stdio.h>
#include <mpfr.h>

int main(void) {
    mpfr_t a, b, c;
    mpfr_inits2(10000000, a, b, c, NULL);
    mpfr_set_str(a, "10", 10, MPFR_RNDD);
    mpfr_set_str(b, "1000000", 10, MPFR_RNDD);
    mpfr_pow(c, a, b, MPFR_RNDD);
    mpfr_printf("%.0Rf\n", c);
    mpfr_clears(a, b, c, NULL);
    mpfr_free_cache();
    return 0;
}
gcc -g -Wall -o main main.c $(pkg-config --cflags --libs mpfr)
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 15:25:33 | 显示全部楼层
本帖最后由 jackz007 于 2021-10-29 15:27 编辑

        通过 Python 计算出来的结果是:
    100 ! = 93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000
        C 语言如果涉及大数运算,一般是要借助于大数运算库的,这些超大数据一般也都是直接用字符串来进行存储和表达的。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-10-29 16:11:51 | 显示全部楼层
村里小黑 发表于 2021-10-29 15:11
原文链接
https://blog.csdn.net/kai8wei/article/details/45110859

这篇文章我刚才还看了
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-10-29 16:17:28 | 显示全部楼层

大佬,double类型不是8个字节吗?为啥能存储这么长的数据,不会溢出吗?
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-10-29 16:18:03 | 显示全部楼层
人造人 发表于 2021-10-29 15:24
用其他库函数,例如 gnu 的 mpfr 库

谢谢大佬,我去看看
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 16:18:46 | 显示全部楼层
番杰 发表于 2021-10-29 16:11
这篇文章我刚才还看了

感觉这个方法不错,相当于把数据给割开存放
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 16:42:38 From FishC Mobile | 显示全部楼层
番杰 发表于 2021-10-29 16:17
大佬,double类型不是8个字节吗?为啥能存储这么长的数据,不会溢出吗?

很大的数用double存的话精度都不够,打个比喻,999536存入,变成9.99e5,低位数没了。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 17:02:54 | 显示全部楼层
番杰 发表于 2021-10-29 16:17
大佬,double类型不是8个字节吗?为啥能存储这么长的数据,不会溢出吗?


不会溢出,不过精度会下降,如 11 楼大佬说的一样。
而且确实很多程序员为了计算超大数值都会用 MPSR 库如 5 楼大佬所说的一样。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-29 23:44:37 | 显示全部楼层
傻眼貓咪 发表于 2021-10-29 17:02
不会溢出,不过精度会下降,如 11 楼大佬说的一样。
而且确实很多程序员为了计算超大数值都会用 MPSR  ...

是 mpfr
不是 mpsr
^_^
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2021-10-30 08:33:06 | 显示全部楼层
人造人 发表于 2021-10-29 23:44
是 mpfr
不是 mpsr
^_^

哈哈,手误
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-9-22 19:22

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表