С numba что-то не так

@Matrix3007 · Регистрация: 13.05.2019

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

Python
1
2
3
4
5
6
7
def fact(x):
    f=1
    while x != 0 and x != 1:
        f*=x*(x-1)
        x-=2
    return f
print(fact(100))

Вывод: 9332621544394415268169923885626670049071 5968264381621468592963895217599993229915 6089414639761565182862536979208272237582 51185210916864000000000000000000000000

Вывод при использовании numba: 0 (Код ниже)

Python
1
2
3
4
5
6
7
8
9
10
from numba import njit
 
@njit
def fact(x):
    f=1
    while x != 0 and x != 1:
        f*=x*(x-1)
        x-=2
    return f
print(fact(100))

На fact(50) вообще отрицательное значение выдаёт. Отчего такое может быть?

@Garry Galler · 21.01.2020, 20:19

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

Отчего такое может быть?

В Python числа безразмерные. А numba как я понимаю в С++ транслирует код и потом компилирует своим JIT.
И с длинной арифметикой там не так просто как в Python.

@Matrix3007 · 21.01.2020, 20:24 **[ТС]**

Garry Galler, Хм, не знал. В си я ещё не лез, не в курсе. Получается, мне надо будет как-то вручную прописывать эти даблы и лонги (или как они там называются)?

@dondublon · 22.01.2020, 16:24

Matrix3007, дайте-ка угадаю. Вы решаете какую-то задачку, типа олимпиадной, с факториалом. И не получается, потому что даже питон дохнет на таких больших числах. Но вы где-то краем уха слышали, что numba ускоряет работу с числами, и надеетесь, что он поможет. Так?

@Matrix3007 · 22.01.2020, 17:36 **[ТС]**

dondublon, Не совсем) Я вычисляю пи, факториал просто в пример. Питон не дохнет, но хотелось бы быстрее)

@__py__ · 26.01.2020, 13:20

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

Питон не дохнет, но хотелось бы быстрее)

через алгоритм побитового сдвига.

Добавлено через 8 минут

Python
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
# coding utf-8
 
 
def naive(a, b):  # алгоритм перемножения чисел
    xa = a
    z = 0
    while xa > 0:
        z = z + b
        xa = xa - 1
    return z
print(naive(2**2, 2**2))
 
 
def ru(a, b):
    x = a
    y = b
    z = 0
    while x > 0:
        if x % 2 == 1:
            z = z + y
        y = y << 1
        x = x >> 1
    return z
print(ru(77**7, 5**9))
# рус крестьянин
s = 0
for i in range(10):
    s = s + i
print(s)
 
 
def rec_rus(a, b):
    if a == 0:
        return 0
    if a % 2 == 0:
        return 2*rec_rus(a/2, b)
    return b + 2*rec_rus((a-1)/2, b)
print(rec_rus(678, 315))
# рекурс рус крестьянин (не перемножает в степенях ** - ошибка)

чтобы было понимание, что это
дальше всё просто

Добавлено через 13 минут
а ещё быстрее будет написать на си и вызвать си код из python

Добавлено через 17 минут
хотя на чистом си ещё быстрее.

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
    #include<stdio.h>  
      
    long factorial(int n)  
    {  
      if (n == 0)  
        return 1;  
      else  
        return(n * factorial(n-1));  
    }  
       
    void main()  
    {  
      int number;  
      long fact;  
      printf("Enter a number: ");  
      scanf("%d", &number);   
       
      fact = factorial(number);  
      printf("Factorial of %d is %ld\n", number, fact);  
      return 0;  
    }

Добавлено через 1 минуту
в общем я показал код осталось только вызвать его с помощью ctype

Добавлено через 14 секунд
сам напишешь надеюсь

Новые блоги и статьи Все статьи Все блоги /
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .

@Matrix3007 198 / 164 / 41 Регистрация: 13.05.2019 Сообщений: 844
	21.01.2020, 20:24 [ТС]
	Garry Galler, Хм, не знал. В си я ещё не лез, не в курсе. Получается, мне надо будет как-то вручную прописывать эти даблы и лонги (или как они там называются)? 0

@dondublon 4652 / 2072 / 366 Регистрация: 17.03.2012 Сообщений: 10,182 Записей в блоге: 6
	22.01.2020, 16:24
	Matrix3007, дайте-ка угадаю. Вы решаете какую-то задачку, типа олимпиадной, с факториалом. И не получается, потому что даже питон дохнет на таких больших числах. Но вы где-то краем уха слышали, что numba ускоряет работу с числами, и надеетесь, что он поможет. Так? 0

@Matrix3007 198 / 164 / 41 Регистрация: 13.05.2019 Сообщений: 844
	22.01.2020, 17:36 [ТС]
	dondublon, Не совсем) Я вычисляю пи, факториал просто в пример. Питон не дохнет, но хотелось бы быстрее) 0

Опции темы