2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу Пред.  1, 2, 3
 
 Re: Тест простоты
Сообщение01.02.2020, 06:55 
Заслуженный участник


20/12/10
9107
Dmitriy40
Если $n$ простое число, то указанное в гипотезе сравнение заведомо верно --- это банальный факт, и многочлены Фибоначчи здесь ни при чем (верно для любого нормированного многочлена степени $n-1$). Так что контрпример --- это составное $n$, для которого сравнение также было бы верно. Т.е. вычислять неприятную сумму нужно только для составных $n$.

Еще один упрощающий момент: сумму можно сократить вдвое, поскольку при нечетном $n$ многочлен $F_n(x)$ является четной функцией.

 Профиль  
                  
 
 Re: Тест простоты
Сообщение01.02.2020, 11:06 
Заслуженный участник


20/08/14
11867
Россия, Москва
nnosipov в сообщении #1437781 писал(а):
Т.е. вычислять неприятную сумму нужно только для составных $n$.
Ну собственно я и не гнался за скоростью. А это почти и не ускоряет, до $n<3000$ всего 430 простых, т.е. проверять надо более 85% чисел. Дальше простых ещё меньше. Так что потеря менее 15% скорости даёт зато и проверку простых тоже.
nnosipov в сообщении #1437781 писал(а):
Еще один упрощающий момент: сумму можно сократить вдвое, поскольку при нечетном $n$ многочлен $F_n(x)$ является четной функцией.
И тоже, ускорение навскидку примерно 25%. Если б гнаться за скоростью, то да, но для первого подхода сойдёт и так.

С практической точки зрения данный тест простоты интереса не вызывает из-за своей медлительности (перебор делителей и то на порядки быстрее), а с точки зрения поиска контрпримера ... Может кому и интересно будет.

 Профиль  
                  
 
 Re: Тест простоты
Сообщение01.02.2020, 11:20 
Заслуженный участник


20/12/10
9107
Dmitriy40 в сообщении #1437789 писал(а):
с точки зрения поиска контрпримера ...
По-моему, только это здесь и интересно.

 Профиль  
                  
 
 Re: Тест простоты
Сообщение02.02.2020, 13:25 
Заслуженный участник


08/04/08
8562
Пусть $f(x)=\sum\limits_{k=0}^{n-1}a_kx^k$ - произвольный, upd: fix: $a_0=1$.
Тогда $S=\sum\limits_{x=0}^{n-1}f(x)=\sum\limits_{x=0}^{n-1}\sum\limits_{k=0}^{n-1}a_kx^k=\sum\limits_{k=0}^{n-1}a_k\sum\limits_{x=0}^{n-1}x^k=n+\sum\limits_{k=1}^{n-1}a_k\sum\limits_{x=0}^{n-1}x^k$. (upd: fix: здесь для $k=0$ получается слагаемое $na_0$ - исправил здесь и далее)
Тогда $S\equiv -1\pmod n \Leftrightarrow \sum\limits_{k=1}^{n-1}a_k\sum\limits_{x=0}^{n-1}x^k \equiv -1 \pmod n$.
Пусть $n$ - свободно от квадратов, $p\mid n$ - любой простой делитель $n$. Тогда $S\equiv -1\pmod n \Leftrightarrow (\forall p) \sum\limits_{k=1}^{n-1}a_k\sum\limits_{x=0}^{n-1}x^k \equiv -1 \pmod p$
Имеем: $\sum\limits_{x=0}^{n-1}x^k \equiv \sum\limits_{0\leqslant x < n, p\nmid n}x^k \equiv \frac{n}{p}\sum\limits_{x=1}^{p-1}x^k \equiv -\frac{n}{p}[p-1|k] \pmod p$.
(Здесь использована нотация Айверсона: $|P(t)|=1 \Leftrightarrow P(t), |P(t)|=0 \Leftrightarrow \neg P(t)$).
Тогда $S\equiv -1\pmod n \Leftrightarrow (\forall p \mid n) \frac{n}{p}\sum\limits_{k=1}^{n-1}a_k[p-1|k]\equiv \frac{n}{p}\sum\limits_{l=1}^{\left[\frac{n-1}{p-1}\right]}a_{(p-1)l}\equiv 1\pmod p$.
Получили общий критерий для произвольного многочлена $f(x)$.
Для многочленов Фибоначчи $a_k=\binom{\frac{n+k-1}{2}}{k}[k\not\equiv n\pmod 2]$.
Возьмем нечетные $n$, тогда $a_{(p-1)l}=\binom{\frac{n+(p-1)l-1}{2}}{(p-1)l}$. Нам остается исследовать выполнимость сравнения $\sum\limits_{l=1}^{\left[\frac{n-1}{p-1}\right]}\binom{\frac{n+(p-1)l-1}{2}}{(p-1)l} \equiv 1\pmod p (\forall p\mid n)$.

Здесь скорее всего почти все слагаемые $a_{(p-1)l}\equiv 0\pmod p$, ненулевыми могут быть только несколько слагаемых сначала или с конца, м.б. даже только с конца. Но здесь у меня руки несколько опустились из-за довольно корявой структуры суммы и нехватки ресурсов. Здесь можно дальше ковырять, или расчехлить железку в поиска контрпримера - расчеты здесь будут сильно быстрее, чем вышеприведенная проверка в лоб. Если кто может и хочет - попробуйте.
Ну и выкладки нужно проверить - я легко мог наврать.
(очень жаль, что нет автоматических средств выполнения подобных преобразований - хорошо было бы поручить проверку ему, да и все)

 Профиль  
                  
 
 Re: Тест простоты
Сообщение02.02.2020, 14:00 
Заслуженный участник


20/12/10
9107
Sonic86 в сообщении #1437920 писал(а):
Пусть $f(x)=\sum\limits_{k=0}^{n-1}a_kx^k$ - произвольный.
Обнаружил у себя в одном из файлов критерий для случая $f(x)=x^{n-1}$. А именно, пусть $n \geqslant 2$. Сравнение $$1^{n-1}+2^{n-1}+\ldots+(n-1)^{n-1} \equiv -1 \pmod{n}$$ имеет место тогда и только тогда, когда $n$ --- простое число или $n$ --- число Кармайкла со свойством: $n/p-1$ делится на $p$ для любого простого делителя $p$ числа $n$. Таким образом, вопрос о контрпримере в данном случае сводится к поиску специальных чисел Кармайкла (существуют ли они, мне неизвестно).

Вообще, хотелось бы послушать начальника транспортного цеха ТС, почему все это может быть интересно.

 Профиль  
                  
 
 Re: Тест простоты
Сообщение02.02.2020, 19:55 
Аватара пользователя


12/10/16
637
Almaty, Kazakhstan
Вычислил некоторое количество полиномов Фибоначчи, для наглядности. $0 \leqslant n, x \leqslant13$

(Картинки)

Изображение

Изображение

Изображение


Всё ещё не до конца понял структуру сумм $f(n)=a \cdot f(n-1)+f(n-2)$ , а именно то - как накапливаются остатки по модулю $n$ .

nnosipov в сообщении #1437781 писал(а):
поскольку при нечетном $n$ многочлен $F_n(x)$ является четной функцией.

Как показала практика, при чётном $x$ и чётном $n$ - функция чётная. При $x$ чётное и $n$ нечётное - функция нечётная.
При $x$ нечётная - значения функции чередуются как для $n$ - чётная, так и для $n$ - нечётная. (через два нечётных одна чётная)

 Профиль  
                  
 
 Re: Тест простоты
Сообщение13.02.2020, 07:37 
Модератор
Аватара пользователя


11/01/06
5710
Pedja в сообщении #1437583 писал(а):
Dmitriy40 в сообщении #1437578 писал(а):
Обычно ищут более быстрые тесты

Это правда. Посмотрите тест с временем выполнения $O (\log ^ 3 n)$, описанным в гипотезе 5.

По поводу гипотезы 5 я уже высказал сомнения на MO.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 37 ]  На страницу Пред.  1, 2, 3

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: Stratim


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group