Длинная арифметика. Факториал 1 <= n <= 10000

@Dennis Ritchie · Регистрация: 27.07.2014

Студворк — интернет-сервис помощи студентам

Добрый вечер.
Подскажите, пожалуйста, возможно ли написать длинную арифметику для вычисления факториала до 10000 включительно на C++ без использования boost и т. д.? На python'e, например, это будет выглядеть так

:

Python
1
2
3
4
5
n = int(input())
a = 1
for i in range(2, n + 1):
    a *= i
print(a)

А что делать на C++?

@Dennis Ritchie · 04.03.2015, 01:46 **[ТС]**

Сообщение от Новичок

Не, вот так вроде работает...

Проверил в тестирующей системе (на сложение длинных чисел там тоже задача есть) на числах, состоящих из 10000 цифр. Программа прошла все тесты.

Сообщение от jvbel

изо всей длинной арифметики здесь реализовано только умножение длинного числа на обычное, но этого для вычисления факториала вполне хватает

Программа проходит только 2 теста из 23. Что-то не так.

@Ромаха · 04.03.2015, 08:33

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
#include <stdio.h>
#include <string.h>
#define radix 10000
#define digits 4
#define size (10000 / digits + 1)
   
typedef struct
{
    int d[size], len;
} bignum_t;
   
void read_bignum (bignum_t *num)
{
    char s[size*radix];
    int i, j;
 
    strcpy (s, "1");
    num->len = 1;
    num->d[0] = 0;
   
    for (i = 0; s[i]; i++) {
        int cr = (s[i] - '0');
        for (j = 0; j < num->len; j++) {
            cr = num->d[j]*10 + cr;
            num->d[j] = cr % radix;
            cr /= radix;
        }
        if (cr > 0)
            num->d[num->len++] = cr;
    }
}
   
void write_bignum (bignum_t *num)
{
    int i;
    printf ("%d", num->d[num->len-1]);
    for (i = num->len-2; i >= 0; i--)
        printf ("%04d", num->d[i]);
}
   
int greatere (bignum_t *a, bignum_t *b)
{
    int i;
    if (a->len > b->len)
        return 1;
    else if (a->len < b->len)
        return 0;
   
    for (i = a-> len-1; i>0; i--)
        if (a->d[i] <= b->d[i])
            break;
   
    return a->d[i] > b->d[i];
}
   
bignum_t pr_bignum(bignum_t *a, int b)
{
    int i, cr = 0;
    bignum_t t;
      
    
    for (i = 0; i < a->len; i++) {
        cr += a->d[i] * b;
        t.d[i] = cr % radix;
        cr /= radix;
    }
     
    t.len = a->len;
    if (cr)
        t.d[t.len++] = cr;
   
    return t;
}
   
int main ()
{
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout); 
      
    bignum_t a;
    int b, i;
 
     
    scanf ("%d", &b);
   
    if (b == 0) {
        printf ("1");
        return 0;
    }
     
    a.d[0] = b;
    a.len = 1;
    for (i = b-1; i > 1; i--)
        a = pr_bignum(&a, i);
    
    write_bignum (&a);
    return 0;
}

Вот.. Писалось давно и на Си

@Dennis Ritchie · 04.03.2015, 08:43 **[ТС]**

Ромаха, для 10000 не работает.

@Ромаха · 04.03.2015, 09:21

Нужно было всего-то размер увеличить..
Заодно потер не нужное..

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#include <stdio.h>
#include <string.h>
#define radix 100000
#define digits 4
#define size (100000 / digits + 1)
   
typedef struct
{
    int d[size], len;
} bignum_t;
   
void write_bignum (bignum_t *num)
{
    int i;
    printf ("%d", num->d[num->len-1]);
    for (i = num->len-2; i >= 0; i--)
        printf ("%04d", num->d[i]);
}
   
 
bignum_t pr_bignum(bignum_t *a, int b)
{
    int i, cr = 0;
    bignum_t t;
      
    
    for (i = 0; i < a->len; i++) {
        cr += a->d[i] * b;
        t.d[i] = cr % radix;
        cr /= radix;
    }
     
    t.len = a->len;
    if (cr)
        t.d[t.len++] = cr;
   
    return t;
}
   
int main ()
{
    bignum_t a;
    int b, i;
 
     
    scanf ("%d", &b);
   
    if (b == 0) {
        printf ("1");
        return 0;
    }
     
    a.d[0] = b;
    a.len = 1;
    for (i = b-1; i > 1; i--)
        a = pr_bignum(&a, i);
    
    write_bignum (&a);
    return 0;
}

@Dennis Ritchie · 04.03.2015, 10:00 **[ТС]**

Сообщение от Ромаха

Нужно было всего-то размер увеличить..

Спасибо. Но я не могу понять, почему тестирующая система засчитывает решение на Python'е на 100%, а на C/C++ проходит только 2 теста из 23.
Вот можете сами проверить:
Факториал
Что не так в C++?

@Ромаха · 04.03.2015, 10:13

Сообщение от Dennis Ritchie

Что не так в C++?

Память.. Вы ее слишком много брали.. Отправляйте изначальный вариант и он зайдет на 100..

Добавлено через 2 минуты
Это чудо тоже на полный балл..

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
#include <stdio.h>
#include <string.h>
#define radix 10000
#define digits 4
#define size (10000 / digits + 1)
   
typedef struct
{
    int d[size], len;
} bignum_t;
   
void write_bignum (bignum_t *num)
{
    int i;
    printf ("%d", num->d[num->len-1]);
    for (i = num->len-2; i >= 0; i--)
        printf ("%04d", num->d[i]);
}
   
 
bignum_t pr_bignum(bignum_t *a, int b)
{
    int i, cr = 0;
    bignum_t t;
      
    
    for (i = 0; i < a->len; i++) {
        cr += a->d[i] * b;
        t.d[i] = cr % radix;
        cr /= radix;
    }
     
    t.len = a->len;
    if (cr)
        t.d[t.len++] = cr;
   
    return t;
}
   
int main ()
{
    bignum_t a;
    int b, i;
 
    scanf ("%d", &b);
   
    if (b == 0) {
        printf ("1");
        return 0;
    }
     
    a.d[0] = b;
    a.len = 1;
    for (i = b-1; i > 1; i--)
        a = pr_bignum(&a, i);
    
    write_bignum (&a);
    return 0;
}

@Dennis Ritchie · 04.03.2015, 10:25 **[ТС]**

Сообщение от Ромаха

Это чудо тоже на полный балл..

А почему мой компилятор факториал для 10000 не считает?

@Ромаха · 04.03.2015, 11:08

Был такой косяк.. Там малясь неверная оценка size..

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
#include <stdio.h>
#include <string.h>
#define radix 100000
#define digits 5
#define size 7132
 
typedef struct
{
    int d[size], len;
} bignum_t;
 
void write_bignum (bignum_t *num)
{
    int i;
    printf ("%d", num->d[num->len-1]);
    for (i = num->len-2; i >= 0; i--)
        printf ("%05d", num->d[i]);
}
 
 
bignum_t pr_bignum(bignum_t *a, int b)
{
    int i, cr = 0;
    bignum_t t;
 
 
    for (i = 0; i < a->len; i++) {
        cr += a->d[i] * b;
        t.d[i] = cr % radix;
        cr /= radix;
    }
 
    t.len = a->len;
    if (cr)
        t.d[t.len++] = cr;
 
    return t;
}
 
int main ()
{
    bignum_t a;
    int b, i;
 
    scanf ("%d", &b);
 
    if (b == 0) {
        printf ("1");
        return 0;
    }
 
    a.d[0] = b;
    a.len = 1;
    for (i = b-1; i > 1; i--)
        a = pr_bignum(&a, i);
 
    write_bignum (&a);
    return 0;
}

Такой код и правильно работает, и полный балл получает

@Dennis Ritchie · 04.03.2015, 11:17 **[ТС]**

Сообщение от Ромаха

Такой код и правильно работает, и полный балл получает

Ну вот теперь я полностью доволен.

Спасибо.

Новые блоги и статьи Все статьи Все блоги /
Алиса нашла кучу ошибок компиляции и запуска в проекте, который без проблем компилировался и запускался))) anaschu 30.06.2026 Я пока посмеюся, но завтра проверю. А вообще интерсно. Дал алисе файл, в котором точно нет ошибок компиляции и запуска, и попросил их найти. Нашла кучу))) Критические ошибки, мешающие компиляции и. . .	сукцессия 16. Общий обзор, в основном что бы другие ии поняли anaschu 29.06.2026 # Передаточный документ: модель микоризной сукцессии (для нового чата) Этот документ предназначен для того, чтобы новый чат Claude мог продолжить работу без необходимости заново разбираться в. . .	сукцессия 15 неявная схема anaschu 29.06.2026 Алиса Калибровка параметров симбиотической модели: технический обзор Содержание: Введение Постановка проблемы Технические аспекты реализации Процесс внедрения изменений	сукцессия 14. Обновленная схема модели anaschu 28.06.2026 ГЛОБАЛЬНАЯ ОПИСАТЕЛЬНАЯ СПЕЦИФИКАЦИЯ ЭКОСИСТЕМНОЙ МОДЕЛИ «SOIL CHEMISTRY & MYCORRHIZA 2. 0» https:/ / ibb. co/ NnkGpfMd Представленная интегрированная схема описывает непрерывную нелинейную. . .
сукцессия 13. Питон модель трехзонного мицелия, пока что в основном арбускулярного anaschu 28.06.2026 ## Разработка агентной модели микоризной сукцессии: от выявления артефактов к созданию комплексной системы ### Аннотация Представлено исследование по разработке агентной модели микоризной. . .	сукцессия 12. краткий список проверок модели перед запуском. anaschu 27.06.2026 Скрытые отказы в моделях систем динамики (SD-models) экологических систем: два случая из практики Контекст Разбирался прототип модели систем динамики (SD-модели) микоризной сукцессии: пять. . .	Сукцессия 11. Проверка орудий перед войной: разработка через тестирование anaschu 27.06.2026 Как не дать модели соврать самой себе: проверки для симуляции микоризной сукцессии Введение Когда вы строите математическую модель живой системы — грибов, растений, почвы — главная опасность. . .	10 сукцессия. Питон код войны грибов и растений anaschu 27.06.2026 import numpy as np class PlantAgent: def __init__(self, name, strategy, initial_biomass): self. name = name self. strategy = strategy # "greedy" (широколиственные) или. . .

@Dennis Ritchie 555 / 148 / 58 Регистрация: 27.07.2014 Сообщений: 2,446
	04.03.2015, 08:43 [ТС]
	Ромаха, для 10000 не работает. 0

Длинная арифметика. Факториал 1 <= n <= 10000

Решение