Ваша задача - дать два целых числа aи bвычислить модульную мультипликативную инверсию по модулю b, если она существует.

Модульная обратная по aмодулю bэто число cтакое, что ac ≡ 1 (mod b). Это число является уникальным по модулю bдля любой пары aи b. Он существует, только если наибольшим общим делителем aи bявляется 1.

Страницу Wikipedia для модульного мультипликативного обратного можно обратиться , если вам необходима дополнительная информация о теме.

Вход и выход

Входные данные задаются либо двумя целыми числами, либо списком из двух целых чисел. Ваша программа должна выводить либо одно число, либо модульную мультипликативную инверсию, которая находится в интервале 0 < c < b, либо значение, указывающее, что инверсии нет. Значение может быть любым, кроме числа в диапазоне (0,b), а также может быть исключением. Однако значение должно быть одинаковым для случаев, когда нет обратного.

0 < a < b можно предположить

правила

• Программа должна завершиться в какой-то момент и решить каждый контрольный пример менее чем за 60 секунд.
• Применяются стандартные лазейки

Контрольные примеры

Тестовые случаи ниже приведены в формате, a, b -> output

1, 2 -> 1
3, 6 -> Does not exist
7, 87 -> 25
25, 87 -> 7
2, 91 -> 46
13, 91 -> Does not exist
19, 1212393831 -> 701912218
31, 73714876143 -> 45180085378
3, 73714876143 -> Does not exist

счет

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

Это и это схожие вопросы, но оба задают конкретные ситуации.

Mathematica, 14 байтов

Обязательный Mathematica встроен :

ModularInverse

Это функция, которая принимает два аргумента ( aи b) и возвращает инверсию мода b, если он существует. Если нет, возвращается ошибка ModularInverse: a is not invertible modulo b..

JavaScript (ES6), 79 73 62 61 байт

Возвращает, falseесли обратного не существует.

Он использует расширенный евклидов алгоритм и решает все тестовые случаи практически мгновенно.

f=(a,b,c=!(n=b),d=1)=>a?f(b%a,a,d,c-(b-b%a)/a*d):b<2&&(c+n)%n

Контрольные примеры

f=(a,b,c=!(n=b),d=1)=>a?f(b%a,a,d,c-(b-b%a)/a*d):b<2&&(c+n)%n

console.log(f(1, 2)) // -> 1
console.log(f(3, 6)) // -> Does not exist
console.log(f(7, 87)) // -> 25
console.log(f(25, 87)) // -> 7
console.log(f(2, 91)) // -> 46
console.log(f(13, 91)) // -> Does not exist
console.log(f(19, 1212393831)) // -> 701912218
console.log(f(31, 73714876143)) // -> 45180085378
console.log(f(3, 73714876143)) // -> Does not exist

Желе , 2 байта

æi

Попробуйте онлайн!

Он использует встроенную функцию для модульной инверсии и возвращает 0 для отсутствия модульной инверсии.

Желе , 7 байт

R×%⁸’¬T

Попробуйте онлайн!

Выводит пустой набор (представленный в виде пустой строки) без модульной инверсии. Запускается нехватка памяти на TIO для самых больших тестовых случаев, но должна работать при достаточном объеме памяти.

Как это работает

R×%⁸’¬T  
R        Generate range of b
 ×       Multiply each by a
  %⁸     Mod each by b
    ’    Decrement (Map 1 to 0 and all else to truthy)
     ¬   Logical NOT
      T  Get the index of the truthy element.

Если вы хотите работать с большими тестовыми примерами, попробуйте эту (относительно безвкусицу) версию, которая требует больше времени, чем памяти:

Желе, 9 байт

×⁴%³’¬ø1#

Попробуйте онлайн!

Как это работает

×⁴%³’¬ø1#
        #   Get the first
      ø1      one integer
            which meets:
×⁴            When multiplied by a
  %³          And modulo-d by b
    ’         Decrement
     ¬        Is falsy

Python 2 , 34 байта

f=lambda a,b:a==1or-~b*f(-b%a,a)/a

Попробуйте онлайн!

Рекурсивная функция , которая дает Trueдля print f(1,2), который я считаю приемлемыми, и ошибки для недействительных входов.

$x$$a\cdot x\equiv 1\pmod{b}$

$a\cdot x-1=k\cdot b$$k$

Взятие $\mod{a}$$-1\equiv k\cdot b\pmod{a}$$-k\cdot b\equiv1\pmod{a}$$k$

$k$$f(-b\%a,a)$

$a$$a$$b$

$k$$a\cdot x-1=k\cdot b$$x$$\frac{k\cdot b+1}{a}$

R + числа , 15 байтов

numbers::modinv

возвращается NAдля тех, кто aбез инверсий мод b.

R-Fiddle, чтобы попробовать это!

R , 33 байта (не конкурирует)

Это приведет к ошибке очень большого размера, bпоскольку фактически создает вектор 32*bбитов размера .

function(a,b)which((1:b*a)%%b==1)

Попробуйте онлайн!

Возвращает integer(0)(пустой список) для тех aбез инверсии мод b.

Mathematica, 18 байт

PowerMod[#,-1,#2]&

вход

[31, 73714876143]

Python 2 , 51 49 54 53 51 49 байт

^{-1 байт благодаря officialaimm
-1 байт благодаря Shaggy}

a,b=input()
i=a<2
while(a*i%b-1)*b%a:i+=1
print+i

Попробуйте онлайн!

Печатает, 0когда нет решения.

Japt , 9 8 байт

Принимает входные данные в обратном порядке. Выходы -1не совпадают. Вылетает, когда большее целое становится больше.

Ç*V%UÃb1

Попробуй это

• Сохранено 1 байт благодаря ETH, указывающему на ошибочное и очень очевидное пространство.

Python 3 + gmpy , 23 байта

Я не думаю, что это может стать немного короче в Python.

gmpy.invert
import gmpy

Попробуйте онлайн! (не будет работать, если у вас не установлен gmpy)

Python 3 , 49 байт

lambda a,b:[c for c in range(b)if-~c*a%b==1][0]+1

Попробуйте онлайн!

Python 3 , 50 байт

lambda a,b:[c for c in range(1,b+1)if c*a%b==1][0]

Попробуйте онлайн!

Это бросает IndexError: list index out of rangeв случае, если нет модульного мультипликативного обратного, как это разрешено правилами.

8 , 6 байтов

Код

invmod

объяснение

invmodявляется восьмым словом, которое вычисляет значение, обратное a, по модулю b. Возвращается nullпри переполнении или других ошибках.

Использование и тестовые случаи

ok> 1 2 invmod .
1
ok> 3 6 invmod .
null
ok> 7 87 invmod .
25
ok> 25 87 invmod .
7
ok> 2 91 invmod .
46
ok> 13 91 invmod .
null
ok> 19 1212393831 invmod .
701912218
ok> 31 73714876143 invmod .
45180085378
ok> 3 73714876143 invmod .
null

Пари / ГП , 22 байта

a->b->lift(1/Mod(a,b))

Выдает ошибку, когда нет обратного.

Попробуйте онлайн!

J , 28 байт

4 :'(1=x+.y)*x y&|@^<:5 p:y'

Попробуйте онлайн!

Использует теорему Эйлера . Возвращает 0, если обратное не существует.

объяснение

4 :'(1=x+.y)*x y&|@^<:5 p:y'  Input: a (LHS), b (RHS)
4 :'                       '  Define an explicit dyad - this is to use the special
                              form `m&|@^` to perform modular exponentiation
                          y   Get b
                      5 p:    Euler totient
                    <:        Decrement
             x                Get a
                   ^          Exponentiate
               y&|@             Modulo b
       x+.y                   GCD of a and b
     1=                       Equals 1
            *                 Multiply

Pyth , 10 байт

3 байта сохранены благодаря @Jakube .

xm%*szdQQ1

Попробуй это здесь!

Возвращает -1без мультипликативного обратного.

Разбивка кода

xm%*szdQQ1      Let Q be the first input.
 m      Q       This maps over [0 ... Q) with a variable d.
   *szd         Now d is multiplied by the evaluated second input.
  %    Q        Now the remained modulo Q is retrieved.
x        1      Then, the first index of 1 is retrieved from that mapping.

Pyth , 15 13 байт

KEhfq1%*QTKSK

Выдает исключение, если мультипликативного обратного не существует.

Попробуй это здесь!

Pyth , 15 байт

Iq1iQKEfq1%*QTK

Это добавляет много байтов для обработки случая, когда такого числа не существует. Программа может быть значительно сокращена, если этот случай не нужно обрабатывать:

fq1%*QTK

Попробуй это здесь!

C (gcc) , 115 байтов

#define L long long
L g(L a,L b,L c,L d){return a?g(b%a,a,d-b/a*c,c):b-1?0:d;}L f(L a,L b){return(g(a,b,1,0)+b)%b;}

Попробуйте онлайн!

Расширенный евклидов алгоритм, рекурсивная версия

C (gcc) , 119 байт

long long f(a,b,c,d,t,n)long long a,b,c,d,t,n;{for(c=1,d=0,n=b;a;a=t)t=d-b/a*c,d=c,c=t,t=b%a,b=a;return b-1?0:(d+n)%n;}

Попробуйте онлайн!

Расширенный евклидов алгоритм, итерационная версия

C (gcc) , 48 110 104 байта

#define f(a,b)g(a,b,!b,1,b)
long g(a,b,c,d,n)long a,b,c,d,n;{a=a?g(b%a,a,d,c-(b-b%a)/a*d):!--b*(c+n)%n;}

Попробуйте онлайн!

Это должно работать со всеми входами (которые подходят в течение длительного времени) в течение 60 секунд.

Редактировать. Я уже злоупотребляю nпеременной, поэтому могу предположить, что gcc вставляет первое присваивание %rax.

Python 2 + sympy, 74 байта

from sympy import*
def f(a,m):i,_,g=numbers.igcdex(a,m);return g==1and i%m

Попробуйте онлайн!

Взято из исходного кода Jelly.

Аксиома, 45 байт

f(x:PI,y:PI):NNI==(gcd(x,y)=1=>invmod(x,y);0)

0 для ошибки, иначе верните z с x * z Mod y = 1

Python 2 , 52 байта

-3 байта благодаря мистеру Xcoder.

f=lambda a,b,i=1:i*a%b==1and i or i<b and f(a,b,i+1)

Попробуйте онлайн!

Выходы Falseна нет решения и ошибок, как bстановится больше.

Встроенный TIO

Я просто тестирую фреймы в Stack Snippets, и они работают просто фантастически.

html,body{height:100%;}iframe{height:100%;width:100%;border:none;}

<iframe src="https://tio.run/##PY3BCsIwEETv/Yq5CEndy6ZotbhfIh4SanBB09Lm4tdHU8HTvIHhzfzOjym5UqI8/SuMHp4CqfCgrd8FEfZphGJaoJeAWqLZJnu2Jd/XvEJwNUxwlmA6wrFmTzj1FdzhT4QzV@DuR7emidULTdhMA@ZFU/4@tGrLBw"></iframe>

JavaScript (ES6), 42 41 39 38 байт

Выходы falseне совпадают. Выдает ошибку переполнения, поскольку второе число становится слишком большим.

x=>y=>(g=z=>x*z%y==1?z:z<y&&g(++z))(1)

Желе , 27 байт

²%³
⁴Ç⁹Сx⁸
ÆṪ’BṚçL$P%³×gỊ¥

Попробуйте онлайн!

Используется теорема Эйлера с модульным возведением в степень. Поскольку у Jelly нет встроенной функции для модульного возведения в степень, его необходимо было реализовать, и он занял большую часть байтов.

Аксиома, 99 байт

w(a,b,x,u)==(a=0=>(b*b=1=>b*x;0);w(b rem a,a,u,x-u*(b quo a)));h(a,b)==(b=0=>0;(b+w(a,b,0,1))rem b)

он использует функцию h (); h (a, b) возвращает 0, если ошибка [не существует обратного], в противном случае она возвращает z так, что a * z mod b = 1 Это было бы нормально, даже если аргументы отрицательны ...

это будет общая функция egcd (), которая возвращает список int (так что они тоже могут быть отрицательными)

egcd(aa:INT,bb:INT):List INT==
   x:=u:=-1   -- because the type is INT
   (a,b,x,u):=(aa,bb,0,1)
   repeat
      a=0=>break
      (q,r):=(b quo a, b rem a)
      (b,a,x,u):=(a,r,u,x-u*q)
   [b,x, (b-x*aa)quo bb]

вот как это использовать

(7) -> h(31,73714876143)
   (7)  45180085378
                                                    Type: PositiveInteger

Я нахожу базовый алгоритм в Интернете с https://pastebin.com/A13ybryc