用高精加、高精乘求阶乘之和

1.题目描述

用高精度计算出 S = 1! + 2! + 3! + ....+ n!(n ≤ 50)。

其中 ! 表示阶乘,定义为n! = n × (n−1) × (n−2) × ⋯ ×1。例如,5! = 5 × 4 × 3 × 2 × 1 = 120。

2.输入格式

一个正整数。 

3.输出格式

一个正整数x,表示计算结果。 

4.输入输出样例

输入:3

输出:9 

输入:22

输出:1177652997443428940313

说明/提示【数据范围】:对于100% 的数据,1 ≤ n ≤ 50。

 5.关键点

1.数组a[]作为每次接受(i)的值的数组

2.数组c[]赋值给b[],b[]作为上一个(i-1)阶乘的值,在新的一轮中再与a[]做高精乘。

 (做到了保留(i-1)!的结果,省事儿了)

3.将赋值后的c[]用memset函数清空数组,接收下一次a[]*b[]的结果

3.f[]作为最后的结果接受数组,与每个阶乘的结果b[]做高精加,最后输出。

6.代码实现 

#include<bits/stdc++.h>
using namespace std;
const int maxx = 100;
int main()
{
    int a[maxx]= {},b[maxx]= {},c[maxx]= {},f[maxx]= {};
    int len_a,len_b = 1,len_c = 1,t,m=1;
    int n,i,j,k;
    cin >> n;
    b[0] = 1;//b[0]要赋初值
    for(i = 1; i <= n; i++)
    {
        len_a = 0;
        int p = i;
        while(p > 0)
        {
            a[len_a++] = p % 10;//存储i到数组a[]
            p = p / 10;
        }
        int jw;
        for(k = 0; k < len_a; k++)
        {
            jw=0;
            for(j = 0; j < len_b; j++)
            {
                c[k+j] = a[k]*b[j] + jw + c[k+j];
                jw=c[k+j]/10;
                c[j+k] = c[k+j]%10;
            }
            c[k+len_b] = jw;
        }
        len_c = len_a+len_b;

        //错在了len_c的长度啊啊啊啊啊啊啊啊啊,检查好长时间呜呜呜呜
        for(int test =len_c-1; test >=0; test--)
        {
            if(0 == c[i] && test > 1)
                test--;
            else
                break;
        }

        t = len_b;
        //c[]求出的两数之积赋给b[]继续运算。
        for(k = len_c - 1; k >=0; k--)
        {
            b[k] = c[k];
        }
        len_b = len_c;

        //高精加,m是和的长度
        m = max(m,len_c)+1;//+1以防进位

        len_c = len_a+t;//以防万一,取最长清零
        memset(c,0,sizeof(c));
        //加法运算
        for(j = 0; j < m; j++)
        {
            f[j] = f[j]+ b[j];
            if(f[j] > 9)
            {
                f[j+1] += f[j]/10;
                f[j] = f[j]%10;
            }
        }
    }

    //删除前导0
    for(i = m - 1; i >=0; i--)
    {
        if(0 == f[i] && m > 1)
            m--;
        else
            break;
    }
    //倒序输出
    for(i = m-1; i >=0; i--)
        cout << f[i];
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
当然可以,以下是一个简单的高精减乘除的 C 代码示例: ```c #include <stdio.h> #include <string.h> #define MAX_LEN 1000 typedef struct { int len; int num[MAX_LEN]; } BigInt; void init(BigInt *a) { memset(a->num, 0, sizeof(a->num)); a->len = 0; } void assign(BigInt *a, char *s) { init(a); int len = strlen(s); for (int i = len - 1, j = 0, k = 1; i >= 0; i--) { a->num[j] += (s[i] - '0') * k; if (k == 10000) { j++; k = 1; } else { k *= 10; } } a->len = len / 4 + 1; while (a->len > 1 && a->num[a->len - 1] == 0) { a->len--; } } void print(BigInt *a) { printf("%d", a->num[a->len - 1]); for (int i = a->len - 2; i >= 0; i--) { printf("%04d", a->num[i]); } printf("\n"); } int compare(BigInt *a, BigInt *b) { if (a->len != b->len) { return a->len - b->len; } for (int i = a->len - 1; i >= 0; i--) { if (a->num[i] != b->num[i]) { return a->num[i] - b->num[i]; } } return 0; } void add(BigInt *a, BigInt *b, BigInt *c) { init(c); int carry = 0; for (int i = 0; i < a->len || i < b->len; i++) { int sum = a->num[i] + b->num[i] + carry; c->num[c->len++] = sum % 10000; carry = sum / 10000; } if (carry > 0) { c->num[c->len++] = carry; } } void sub(BigInt *a, BigInt *b, BigInt *c) { init(c); int borrow = 0; for (int i = 0; i < a->len || i < b->len; i++) { int diff = a->num[i] - b->num[i] - borrow; if (diff < 0) { diff += 10000; borrow = 1; } else { borrow = 0; } c->num[c->len++] = diff; } while (c->len > 1 && c->num[c->len - 1] == 0) { c->len--; } } void mul(BigInt *a, BigInt *b, BigInt *c) { init(c); for (int i = 0; i < a->len; i++) { int carry = 0; for (int j = 0; j < b->len; j++) { int sum = a->num[i] * b->num[j] + c->num[i + j] + carry; c->num[i + j] = sum % 10000; carry = sum / 10000; } if (carry > 0) { c->num[i + b->len] += carry; } } c->len = a->len + b->len; while (c->len > 1 && c->num[c->len - 1] == 0) { c->len--; } } void div(BigInt *a, BigInt *b, BigInt *c, BigInt *d) { init(c); init(d); for (int i = a->len - 1; i >= 0; i--) { d->num[0] = a->num[i]; d->len = 1; int left = 0, right = 9999, mid; while (left <= right) { mid = (left + right) / 2; BigInt tmp, mul; init(&tmp); init(&mul); tmp.num[0] = mid; tmp.len = 1; mul(&tmp, b, &mul); if (compare(&mul, d) <= 0) { left = mid + 1; } else { right = mid - 1; } } c->num[i] = right; BigInt tmp, mul, sub; init(&tmp); init(&mul); init(&sub); tmp.num[0] = right; tmp.len = 1; mul(&tmp, b, &mul); sub(a, &mul, &sub); assign(a, ""); add(&sub, &tmp, a); } while (c->len > 1 && c->num[c->len - 1] == 0) { c->len--; } } int main() { BigInt a, b, c, d; char s1[MAX_LEN], s2[MAX_LEN]; scanf("%s%s", s1, s2); assign(&a, s1); assign(&b, s2); add(&a, &b, &c); sub(&a, &b, &c); mul(&a, &b, &c); div(&a, &b, &c, &d); print(&c); print(&d); return 0; } ``` 这个代码实现了高精减乘除四个基本运算,可以处理超过 long long 范围的整数。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值