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

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

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

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

Спасибо.

Новые блоги и статьи Все статьи Все блоги /
Саморегулирующийся социальный контракт для сервера cross-section. Hrethgir 14.08.2026 С кодом конечно таких глубоких размышлений пока не было, впрочем я уже привык к алгоритмизации. Суть предмета записи: снова в диалоге с нейросетью (я взял пока себе ник для учётки админа - Rector). . . .	Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .	Установка MinGW GCC 16.2 и CMake 8Observer8 10.08.2026 VK Видео: https:/ / vkvideo. ru/ video-240781534_456239017 YouTube: eY5-5PyI9NM Текстовая версия
Неделя из жизни имитационной модели склада: мои кривые руки растут, откуда надо anaschu 10.08.2026 Неделя из жизни имитационной модели склада: как я почти написал неправильную логику и что с этим делать Работаю сейчас над учебно-рабочим проектом: строю в AnyLogic имитационную модель процессов. . .	Калькулятор для расчета родства russiannick 07.08.2026 1. Задача: Создать калькулятор для расчета родства. Родственных связей существует 8 ступеней, такие как: p - отец P - мать q - муж Q - жена b - брат B - сестра s - сын S - дочь	Мир по моей воле kumehtar 07.08.2026 Когда-то кажется, что всё просто. Ты весь такой светлый. Причиняешь добро. Борешься за справедливость в этом тёмном мире. Потом начинаешь замечать одну неприятную вещь. Почти каждый хороший. . .	Кредитный калькулятор Maks 05.08.2026 Решение задачи по прикладной информатике средствами 1С. Задача: Напишите приложение-калькулятор, которое помогает рассчитывать параметры кредита для аннуитетного и дифференцированного видов. . .

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

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

Решение