Модульная Арифметика: Удобная и простая работа

Правка ru2, от dmitryAdams, 2024-11-19 23:29:29

Всем от новичка до профессионала приходится иметь дело с модульной арифметикой в очень широком спектре задач, от хеширования строк, до подсчета дп на деревьях (да-да, сейчас кто-то скажет, что все сводится к хэшам и комбе, но тема от этого не менее часто встречается).

Мне же еще с самых первых дней в олимпиадном программировании не нравились задачи с модульной арифметикой, забудешь один раз ans %= MOD; И сиди с WA, и, либо минусом по задаче, либо просто трать время на дебаг, поэтому всегда приходилось внимательно расписывать формулы. Однажды меня познакомили с хэшами и удобным способом их написания: функциями необходимыми, чтобы не забыть о модулях по типу add(x, y), sub(x, y) и mul(x, y). Код, хоть и стал выглядеть гораздо лучше, но все еще оставался нечитаемым, да и задумываться о прикладных мелочах, пока считаешь математику, не хочется, поэтому логическим продолжением данной темы с функциями (имхо) будет класс модульных интов. Я предлагаю следующую реализацию: #include "iostream" template<class T, T MOD> class ModInt{ public: ModInt(const T &x = T(0)) : x_(x % MOD) {} friend ModInt& operator+=(ModInt &a, const ModInt &b){ a.x_ += b.x_; if (a.x_ >= MOD) a.x_ -= MOD; return a; } friend ModInt operator+(ModInt a, const ModInt &b){ return a += b; } friend ModInt& operator-=(ModInt &a, const ModInt &b){ a.x_ -= b.x_; if (a.x_ < 0) a.x_ += MOD; return a; } friend ModInt operator-(ModInt a, const ModInt &b){ return a -= b; } friend ModInt& operator*=(ModInt &a, const ModInt &b){ a.x_ = a.x_ * b.x_ % MOD; return a; } friend ModInt& operator*(ModInt a, const ModInt &b){ return a *= b; } ModInt& operator=(const ModInt &b){ x_ = b.x_; return *this; } friend std::istream &operator>>(std::istream &in, ModInt &a){ in >> a.x_; a.x_ %= MOD; return in; } friend std::ostream &operator<<(std::ostream &out, const ModInt &a){ out << a.x_; return out; } ModInt& operator++(){ ++x_; if (x_ == MOD) x_ = 0; return *this; } ModInt& operator--(){ --x_; if (x_ == -1) x_ = MOD &mdash; 1; return *this; } ModInt operator++(int ){ ModInt old = *this; ++x_; if (x_ == MOD) x_ = 0; return old; } ModInt operator--(int ){ ModInt old = *this; --x_; if (x_ == -1) x_ = MOD &mdash; 1; return old; } private: T x_; };

Теперь все, что вам нужно сделать в задаче, чтобы ваш ответ считался по модулю — это написать: #define int ModInt<int64_t, 998244353>

WARNING: Константы по типу INF или подобные могут испортиться, думайте, когда используете этот define

Да, такой #define использовать плохо, но кто мне запретит. Конечно, теперь решение работает незначительно медленнее, например, потому что счетчики в циклах for теперь тоже работают по модулю (а вы их 100% часто используете), но, если вам не нужно пихать задачу, то жизнь вам это хуже не сделает.

Плюсы: 1) Хэши теперь можно писать почти как в 2009, когда еще не был известен anti-hash test против модулей 2^k 2) Любую комбинаторику из тетрадки можно просто перенести в код без лишней мороки

Минусы: 1) Незначительно ухудшает скорость программы

Лично для меня + сильно перевешивают -. Если придется загонять задачу, то, скорее всего, так и так придется хардкодить. Здесь же сильно упрощена работа с любой арифметикой и код становится гораздо более читаемым.

Далее можно развить класс, добавив модульное деление, т.е. умножение на обратное по модулю, но это я уже оставлю в качестве тривиального упражнения читателю)

В заключение: Конечно, на туре такое писать вряд ли быстрее, чем написать пару функций или, вообще, все считать в одной строке, но не всегда цель — решить задачу как можно быстрее, иногда этим процессом хочется насладиться, или для поиска бага повысить читаемость кода, да и не всегда пишешь задачи в какой-то ограниченный по времени тур, иногда закинуть готовый шаблон гораздо приятнее, чем переписывать все формулы по модулю на фул)

Теги целочисленная арифметика, модульная арифметика

История

 
 
 
 
Правки
 
 
  Rev. Язык Кто Когда Δ Комментарий
ru7 Русский dmitryAdams 2024-11-19 23:51:38 12 Мелкая правка: '\nПлюсы:\n1) Хэши ' -> '\nПлюсы:\n\n\n1) Хэши ' (опубликовано)
ru6 Русский dmitryAdams 2024-11-19 23:47:40 12 Мелкая правка: 'изацию:\n```\n#includ' -> 'изацию:\n`\n#includ'
ru5 Русский dmitryAdams 2024-11-19 23:36:16 1 Мелкая правка: 'ModInt{\n public:\n' -> 'ModInt{\n public:\n'
ru4 Русский dmitryAdams 2024-11-19 23:35:38 19
ru3 Русский dmitryAdams 2024-11-19 23:35:04 7
ru2 Русский dmitryAdams 2024-11-19 23:29:29 24
ru1 Русский dmitryAdams 2024-11-19 23:28:22 4180 Первая редакция (сохранено в черновиках)