Теорема Вильсона гласит:

- простое, тогда и только тогда, когда

делится на

.
Предлагаю следующее обобщение:
Введём обозначение:

,
Замечание: для

получается обычный факториал.
Тогда
Теорема :

- простое, тогда и только тогда, когда

делится на

.
Где

- натуральные числа,

- остатки от деления на

(

), а

- решение сравнения

такое что

.
Замечание: 
нельзя выразить формулой. Но

можно выразить алгоритмически, и показать что

всегда существует.
Примеры: 

Для

справедливо
соотношение:

При произвольных

оно может быть не верно.
Пример: 
-----------------------------------------------------------------------------
Замечание: определение

может казаться искусственным, но, на самом деле, оно встречается при разложении в ряды некоторых неэлементарных интегралов;
Можно попробовать:

т.е.

где

- матрицы,

- единичный вектор,

- скалярное произведение,

-

-тые строки соотв. матриц.
Теорема: Если есть целочисленные матрицы

, причём в каждой

есть

, то есть и матрица

и константы

, такие что (см. выше).
-----------------------------------------------------------------------------
Можно обобщать и дальше: добавить степени, целые числа итп.
(Не исключены опечатки)
Сейчас внимательно посмотрел Википедию, очень похоже что Гаусс меня опередил
