информационная безопасность
без паники и всерьез
 подробно о проектеRambler's Top100
Атака на InternetГде водятся OGRыЗа кого нас держат?
BugTraq.Ru
Русский BugTraq
 Анализ криптографических сетевых... 
 Модель надежности двухузлового... 
 Специальные марковские модели надежности... 
 Microsoft обещает радикально усилить... 
 Ядро Linux избавляется от российских... 
 20 лет Ubuntu 
главная обзор RSN блог библиотека закон бред форум dnet о проекте
bugtraq.ru / форум / theory
Имя Пароль
ФОРУМ
все доски
FAQ
IRC
новые сообщения
site updates
guestbook
beginners
sysadmin
programming
operating systems
theory
web building
software
hardware
networking
law
hacking
gadgets
job
dnet
humor
miscellaneous
scrap
регистрация





Легенда:
  новое сообщение
  закрытая нитка
  новое сообщение
  в закрытой нитке
  старое сообщение
  • Напоминаю, что масса вопросов по функционированию форума снимается после прочтения его описания.
  • Новичкам также крайне полезно ознакомиться с данным документом.
Я уже потом, по дороге домой, понял, что ты неправильно... 28.03.05 11:22  Число просмотров: 4402
Автор: DPP <Dmitry P. Pimenov> Статус: The Elderman
<"чистая" ссылка>
> Извини, не правильно понял.

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

> Я думал про арифметическую прогрессию, где основное
> свойство ряда а
> а(x-1)-a(x)=a(x)-a(x+1)=dx.
> откуда а(x-1)+a(x+1)=2а(х) и а(x-y)+a(x+y)=2а(х).
> поэтому я и рассматривал проивзедение (x+y)(x-y)=x*x-y-y.
> Я просто хотел попробовать "свернуть" часть ряда, чтобы
> быстрее его вычислить. Надеялся, что если свернуть 8 членов
> можно будет выиграть еще операций. В теории,
> если бы удалось для них получить многочлен вида
> x^8+k1*x^4+k2*x^2+k3. То был бы дополнительный выигрыш. И
> если бы тенденция сохранилась, то 2^m членов свернулось бы
> к полиному m-ой степени. :)
> Не получается. :( :)
> ??? не понял, это про сворачивание, или про то, что я
> предложил?

Я там тоже ошибся, в програмке, которая проверяля предположение. При очередном "сворачивании" "обнуляется" каждая вторая "производная".
Короче все равно "не катит".
Была мысль поработать над быстрым вычислением "простого" факториала, то есть произведения вех простых чисел не более заданного. Но похоже там все еще сложнее.

> > Над приближенным вычислением подумаю, им тоже можно
> > воспользоваться.
> Один из вариантов - если иметь приближенное значение, то
> можно пытаться искать тот самый НОД в некотором диапазоне
> от числа. Главное, чтобы приближение было в разумных
> пределах, а не районе корня (N).:)

Это понятно, но проблема будет в самом поиске в этом диапазоне. Вычислять сам факториал бессмыслено, поскольку его разрядность космическая. Остается вопрос как перебрать НОДы от чисел диапазона, зная НОД от приближенного значения.
Ну и опять же вопрос точности - количества переборов.

> Кстати, ты ничего не слышал о
> "Дэниэл Бернстайн (Daniel Bernstein) опубликовал статью,
> где описал новый метод вычисления факториалов больших
> чисел"
> http://www.pcsamara.ru/news.phtml?n=3&day=13&m=4&am
> p;y=2002&mode=prev
>
> И о формуле Стирлинга
> http://ru.wikipedia.org/wiki/Формула_Стирлинга#.D0.A4.D0.BE
> .D1.80.D0.BC.D1.83.D0.BB.D0.B0_.D0.A1.D1.82.D0.B8.D1.80.D0.
> BB.D0.B8.D0.BD.D0.B3.D0.B0

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

> ее "неприменимость" определяется,как мне кажется ,кроме
> неточности,
> вычислением степени т.е. (n/e)^n, однако, если учесть, что
> нам нужно
> значение факториала по модулю, то
> можно было бы рассчитать (n/e)^n mod n. Как поделить
> (n/e)^n на n - ясно - это
> (n/e)^(n-log по основанию n/e от n). Неясно, как остаток
> получить. :)

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

http://primenumber.narod.ru/

http://www.gamedev.ru/forum/?action=showtopic&group=0&topic=1939

http://ssl.stu.neva.ru/psw/crypto/timing.html
unsigned long st( unsigned long x, unsigned long y, unsigned long n ){

    unsigned __int64 s, t, u;

    s = 1;
    t = x;
    u = y;
    while( u ){
        if( u & 1 )
            s = ( s * t ) % n;
        u >>= 1;
        t = ( t * t ) % n;
    }

    return( s );
} 

long powmod( long a, long k, long n ){

    long b = 1;

    while( k ){
        if( k % 2 == 0 ){
            k /= 2;
            a *= a; // [ a = (a*a)%n; ]
        }
        else{
            k--;
            b *= a; // [ b = (b*a)%n; ]
        }
    }

    return b;
}

---
<theory> Поиск 






Rambler's Top100
Рейтинг@Mail.ru


  Copyright © 2001-2024 Dmitry Leonov   Page build time: 0 s   Design: Vadim Derkach