Побитовые операции с bigz в gmp - PullRequest
Новая
       74

Побитовые операции с bigz в gmp

0 голосов
/ 28 апреля 2018

Я перевожу некоторые скриптовые скрипты с Python на R. Кажется, что Python обрабатывает очень большие целые числа намного лучше, чем R может изначально: 

10593080468914978578954316149578855170502344604886137564370015851276669104055 >> 1
# 5296540234457489289477158074789427585251172302443068782185007925638334552027

Но мне известна библиотека gmp для R, которая хорошо с ними справляется (в основном): 

as.bigz("10593080468914978578954316149578855170502344604886137564370015851276669104055")

Для контекста, чтобы перевести эти сценарии, мне нужно использовать побитовые операции. Проблема в том, что эти bigz объекты кодируются как необработанные значения, и поэтому я не могу использовать для них базовые побитовые функции, поскольку они несовместимы.

Найти обходной путь для сдвига битов влево и вправо - прямо вперед, но мне нужно кое-что, что будет:

  • Выполните эквивалент bitwAnd и bitwOr
  • Вкл. bigz значения
  • БЕЗ потери точности.

Есть идеи?

Бонус: если вы можете дать интерпретацию bitwAnd и bitwOr в терминах базы 10, тогда это может сработать. Предпочтительно с некоторым примером кода на R, если нет, я могу обойти это.

1 Ответ

0 голосов
/ 28 апреля 2018

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

library(gmp)

z <- as.bigz("10593080468914978578954316149578855170502344604886137564370015851276669104055")  
w <- as.bigz("1234874454654321549879876546351546654456432132321654987584654321321")

#express as numeric vectors of 0s and 1s
z1 <- as.numeric(charToRaw(as.character(z, b=2)))-48
w1 <- as.numeric(charToRaw(as.character(w, b=2)))-48

#normalise the lengths
mx <- max(length(z1), length(w1))
z1 <- c(rep(0, mx-length(z1)), z1)
w1 <- c(rep(0, mx-length(w1)), w1)

#then do & or | and convert back to bigz
zandw <- as.bigz(paste0("0b", rawToChar(as.raw(1*(z1 & w1) + 48))))
zorw <- as.bigz(paste0("0b", rawToChar(as.raw(1*(z1 | w1) + 48))))

zandw
Big Integer ('bigz') :
[1] 905773543034890641004226585015137324621885921615658881499355162273

zorw
Big Integer ('bigz') :
[1] 10593080469244079490573747058454505131838753934720683775076011957361968263103
...