Обчисліть останні цифри числа Грема


13

Число Грема закінчується на 7. Це масове число, теоретично вимагає для зберігання більше інформації, ніж розмір самої Всесвіту. Однак можна обчислити останні кілька цифр числа Грема.

Останні кілька цифр:

02425950695064738395657479136519351798334535362521
43003540126026771622672160419810652263169355188780
38814483140652526168785095552646051071172000997092
91249544378887496062882911725063001303622934916080
25459461494578871427832350829242102091825896753560
43086993801689249889268099510169055919951195027887
17830837018340236474548882222161573228010132974509
27344594504343300901096928025352751833289884461508
94042482650181938515625357963996189939679054966380
03222348723967018485186439059104575627262464195387

Ваша програма може не містити цих (або подібних номерів), але повинна їх обчислити. Він повинен обчислити 200 або більше цифр.

Вихід у stdout. Час роботи максимум 2 хвилини на гідному обладнання. Найкоротша програма виграє.


Скільки цифр слід надрукувати?
Догберт

@Dogbert D'oh. Я пропустив це. 200 і більше було б добре.
Томас О

Рубі навіть не підрахує, 3**7625597484987тоді як Python це робить :)
gnibbler

@ gnibbler, гм як? результат отримав би більше 3 трлн цифр.
Догберт

1
@Dogbert, маючи достатню кількість пам'яті та часу, Python піде вперед, обчисливши його, використовуючи його longs. Рубі навіть не зробить 3 ** 5000000. Здається, там є якась межа
гніблер

Відповіді:


9

постійний струм - 21 символ

[3z202>xO200^|]dsxxrp

На моєму комп’ютері це займає близько хвилини, а значення, які перевищують 200, займуть більше часу. Це не виводить нульові нулі.

Ось трохи довша, але швидша версія (26 символів):

[3rAz^|dz205>x]dsxxAz6-^%p
3[3rAz^|dz202>x]dsxxAz3-^%p # or an extra character for a few less iterations

4

Хаскелл, 99

Продуктивність не зоряна, але їй вдається обчислити 500 цифр за хвилину на моєму десятилітньому старому обладнанні.

f a b|b==0=1|odd b=mod(a*f a(b-1))m|0<1=f(mod(a^2)m)$div b 2
main=print$iterate(f 3)3!!500
m=10^500

(btw, я хотів би почути про його ефективність на більш сучасному обладнання)


На моєму ПК потрібно близько 19 секунд. З боку примітки, це не надрукує провідний 0 перед виходом.
Догберт

Так, це помилка на всіх підрахунках цифр з провідними нулями. Просто обчисліть для 501 ;-) Дякую за тест. Ви запускали його інтерпретованим чи складеним?
JB

Я склав його ghc -o g.exe g.hs. Не впевнений, чи це найкращий спосіб компілювати.
Догберт

Я щойно запустив ghc -O3 graham.hs рекомендовані варіанти зловживань від інтернет-документа -O2 -fvia-C. (і схоже на те, що мій GHC вже на кілька випусків)
JB

Здається, він працює з однаковою швидкістю з обома, -O3і -O2 -fvia-Cприблизно за 18,3 секунди.
Догберт

3

Пітон - 41 символ

499 цифр

x=3;exec'x=pow(3,x,10**500);'*500;print x

500 цифр

x=3;exec'x=pow(3,x,10**500);'*500;print'0'+`x`

1
Ви використовуєте знання, що 500-та цифра ззаду - 0. Це дасть неправильну відповідь, скажімо, за 200.

1
@Tim Проблема вимагає "200 цифр або більше". Просто жорсткий код підрахунку, який працює і з ним можна зробити. (або залиште його як таке: він надруковує 499 цифр, і це досить добре для запитаного запитання)
JB

@JB: Звичайно, я був би задоволений 499, якщо 0 не залишиться. Однак тепер передбачається, що конкретна цифра дорівнює 0.

@ user475 - За властивостями енергетичних веж, якщо ви обчислюєте останні (d) цифри, а результат менше (d) цифр, то відсутні цифри (зліва) повинні бути "0". Тож добре додати пропущену цифру "0", але це слід зробити, вивчивши довжину результату та додавши відповідне число "0".
Кевін Феган

3

Пітон - 62 59 55 символів

x=3
for i in range(500):x=pow(3,x,10**500)
print"0%d"%x

На моєму ПК займає близько 12 секунд.

sh-3.1$ time python cg_graham.py
02425950695064738395657479136519351798334535362521430035401260267716226721604198
10652263169355188780388144831406525261687850955526460510711720009970929124954437
88874960628829117250630013036229349160802545946149457887142783235082924210209182
58967535604308699380168924988926809951016905591995119502788717830837018340236474
54888222216157322801013297450927344594504343300901096928025352751833289884461508
94042482650181938515625357963996189939679054966380032223487239670184851864390591
04575627262464195387

real    0m11.807s
user    0m0.000s
sys     0m0.015s
sh-3.1$

3
Рідний порох - вбивця :-)
JB

Ви можете використовувати10**500
гніблер

@JB, це єдина причина, що я використовував Python для цього запису :)
Dogbert

@ gnibbler, оновлено, дякую! Я новачок у Python :)
Dogbert

0

Аксіома, 63 байти

f()==(r:=3;d:=10^203;for i in 1..203 repeat r:=powmod(3,r,d);r)

ungolf і результат

--This find the Graham's number follow the algo found in wiki
--http://en.wikipedia.org/wiki/Graham%27s_number
ff()==
   r:=3; d:=10^203
   for i in 1..203 repeat r:=powmod(3,r,d)
   r

(3) -> a:=f()::String
   (3)
  "8871783083701834023647454888222216157322801013297450927344594504343300901096
  92802535275183328988446150894042482650181938515625357963996189939679054966380
  03222348723967018485186439059104575627262464195387"
                                                             Type: String
(4) -> #a
   (4)  203
                                                    Type: PositiveInteger

# a = 203 означає, що кількість len становить> 200, це також означає, що спочатку у нього немає 0 ...


0

Навушники, 602 байти

(h@HP0&Y+h8h (hx0RWc@4vYcx#@#PJ"B?[#CPx (h Lx$2(pl2YL;KD:T{9$2j<
 LSSh,ZT l2I<Pp,@4SX`,:xtc@T",($)<cKT\lbBAy44,dQl[yL"l+i,;9<*j0P
|)lD[+`\RBi!< LaD(LHPLyt{{@\iADRQdHTZQIT3[X`DB*`X$Cxh$*(T0$| ,[;
4:bit0DqAqi!lCYQ)<Ad(|1<$R4l+#tZrLPDatC[d*@0pDclJbh0|#S9<JRy!TP0
D+!|qiTXp<r$##Atj,B1ts;HLJ"Xp44I4cK4@|Q,4JI$|hp$Zyd+yl:y<s#\pD:9
4RDK,A!<X \cqLZ" h,kHp|qLRQIDh,+StZbL+{(|jqqL;9L"(xd"<s$8\:x,CY\
z0T[,(XdcxlbaD*D;+tDj\JIi4k[)LPDLBzP@DSc$jL $s4BjQ19|;!)|9t;TaQA
dLzit[0@l2!)I$,0P@$y<L4pLPLStQT"!a| $+8(DZ`00t ,RtX4C@$yQ11|KI\"
`|2<k3|R(Dyi<)LshTrzQ$sp D+DRbH|Q$CqT0D;AA\jdXd"ppdK3LzZl#\Bl`@t
k$*,11qTK+Xp|rqDZXp4{C!<Y4

Друкує останні 200 цифр.

Перед запуском видаліть нові рядки.


Як ми повинні це запустити?
caird coinheringaahing

Абсолютно ніякої ідеї (я щойно переклав це з BF). Але я шукав "навушники" на github, і, схоже, є кілька реалізацій (хоча посилання на виконання посилань здається мертвим).
Esolanging Fruit
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.