Длинная арифметика. Факториал 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 **[ТС]**

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

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

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

Спасибо.

Новые блоги и статьи Все статьи Все блоги /
PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .
Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .

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

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

Решение