Otrzymujesz funkcję Rand5 (). Ta funkcja zwraca idealnie losowe (równomiernie rozmieszczone) liczby całkowite od 1 do 5.
Podaj funkcję Rand7 (), która używa Rand5 () do tworzenia idealnie losowych liczb całkowitych od 1 do 7.
Otrzymujesz funkcję Rand5 (). Ta funkcja zwraca idealnie losowe (równomiernie rozmieszczone) liczby całkowite od 1 do 5.
Podaj funkcję Rand7 (), która używa Rand5 () do tworzenia idealnie losowych liczb całkowitych od 1 do 7.
Odpowiedzi:
Java - 61 znaków
int rand7(){int s=0,c=7;while(c-->0)s+=rand5();return s%7+1;}
Przetestuj sterownik do sprawdzania poprawności:
class Rand {
public static void main(String[] args) {
int[] nums = new int[7];
// get a lot of numbers
for(int i = 0; i < 10000000; i++) nums[rand7()-1]++;
// print the results
for(int i = 0; i < 7; i++) System.out.println((i+1) + ": " + nums[i]);
}
// just for rand5()
static java.util.Random r = new java.util.Random();
static int rand5() {
return r.nextInt(5)+1; // Random.nextInt(n) returns 0..n-1, so add 1
}
static int rand7(){int s=0,c=7;while(c-->0)s+=rand5();return s%7+1;}
}
Wyniki
C:\Documents and Settings\glowcoder\My Documents>java Rand
1: 1429828
2: 1429347
3: 1428328
4: 1426486
5: 1426784
6: 1429853
7: 1429374
C:\Documents and Settings\glowcoder\My Documents>
rand5
. Przeliczyłem je w Maple przy użyciu prostej algebry macierzowej, ale możesz to zrobić ołówkiem i papierem w ciągu kilku minut, jeśli chcesz. W każdym razie okazuje się, że Omar opublikował już te same dane (bez współczynnika normalizującego) w komentarzu do innej odpowiedzi kilka dni wcześniej. (Również ps., Możesz @ tylko powiadomić jednego użytkownika na komentarz, chociaż autor wpisu jest zawsze powiadamiany zawsze.)
sub rand7{($x=5*&rand5+&rand5-3)<24?int($x/3):&rand7}
Plus mogę użyć trójskładnikowego operatora ORAZ rekurencji. Najlepszy… dzień… kiedykolwiek!
OK, 47 znaków, jeśli używasz mod zamiast div:
sub rand7{($x=5*&rand5+&rand5)<27?$x%7+1:&rand7}
&
znaki, aby sprowadzić go do 46 znaków (w tym spacji, co oznacza, że twoja obecna wersja to 48).
Ruby - 54 znaki (na podstawie rozwiązania Dana McGratha, przy użyciu pętli)
def rand7;x=8;while x>7 do x=rand5+5*rand5-5 end;x;end
Ruby - 45 znaków (to samo rozwiązanie, przy użyciu rekurencji)
def rand7;x=rand5+5*rand5-5;x>7 ?rand7: x;end
(x=rand5+5*rand5-5)>7?
.
W Common Lisp 70 znaków:
(defun rand7()(let((n(-(+(rand5)(* 5(rand5)))5)))(if(> n 7)(rand7)n)))
Nawiasy zajmują więcej miejsca, niż bym chciał.
(defun rand7()(setq n(-(+(rand5)(* 5(rand5)))5))(if(> n 7)(rand7)n))
(defun rand7()(if(>(setq n(-(+(rand5)(* 5(rand5)))5))7)(rand7)n))
W c / c ++ przy użyciu próbkowania odrzucania
int rand7(){int x=8;while(x>7)x=rand5()+5*rand5()-5;return x;}
62 znaków.
while(x>7)
taki, że byłyby spełnione tylko liczby w prawidłowym zakresie.
Tłumaczenie na PHP, z odpowiedzi opublikowanej przez Dana McGratha.
function Rand7(){$x=8;while($x>7)$x=rand5()+5*rand5()-5;return $x;}
67 znaków.
W języku R (język zbudowany do obliczeń statystycznych) celowo oszustwo:
# Construct a Rand5 function
Rand5 <- function() sample(seq(5),1)
# And the golf
Rand7=function(r=Rand5())sample(1:(r/r+6),1)
# Or (same character count)
Rand7=function(r=Rand5())sample.int(r/r+6,1)
# Or even shorter(thanks to @Spacedman)
Rand7=function()sample(7)[Rand5()]
Dzięki leniwej ocenie argumentów wyeliminowałem średnik i nawiasy klamrowe.
Wyjście ponad 10 ^ 6 powtórzeń:
> test <- replicate(10^6,Rand7())
> table(test)
test
1 2 3 4 5 6 7
142987 142547 143133 142719 142897 142869 142848
library(ggplot2)
qplot(test)
Rand7=function(){r=Rand5();sample(7)[r]}
Rand7=function(){sample(7)[Rand5()]}
def rand7:Int={val r=5*(rand5-1)+rand5
if(r<8)r else rand7}
z 2 wejściami z rand5:
\ 1 2 3 4 5
1 1 2 3 4 5
2 6 7 8 ..
3 11 ..
4 ..
5
Pomnożę pierwszy-1 przez 5 i dodam drugi. Większość wyników jest ignorowana i prowadzi do nowych obliczeń. Wynikiem powinna być równa dystrybucja wartości od 1 do 25, z których wybieram tylko 7 pierwszych. Mogłem zaakceptować pierwsze 21 z budowaniem modulo, ale doprowadziłoby to do dłuższego kodu.
historyczny kod, który zawiódł, ale niezbyt oczywiście. Dzięki Ilmari Karonen za zwrócenie na to uwagi:
def rand7=(1 to 7).map(_=>rand5).sum%7+1
Dzięki Yoshiteru Takeshicie za to podejście scala-2.8.0, dzięki któremu „suma” była tak łatwa. Moje rozwiązanie wcześniej:
def rand7=((0/:(1 to 7))((a,_)=>a+rand5-1))%7+1
rand5:
val rnd = util.Random
def rand5 = rnd.nextInt (5) + 1
def rand7=(1 to 7).map(_=>rand5).sum%7+1
int Rand4()
{
int r = Rand5();
return r > 4 ? Rand4() : r;
}
inline int Rand8()
{
return (Rand4() - 1) << 2 + Rand4();
}
int Rand7()
{
int r = Rand8();
return r > 7 ? Rand7() : r;
}
int Rand4(){int r=Rand5();return r>4?Rand4():r;}int Rand7(){int r=Rand4()-1<<2+Rand4();return r>7?Rand7():r;}
Tłumaczenie na Javascript, z odpowiedzi zamieszczonej przez Dana McGratha.
function Rand7(){x=8;while(x>7)x=rand5()+5*rand5()-5;return x}
62 znaki
function Rand7(){for(x=8;x>7;x=rand5()+5*rand5()-5);return x}
jest nieco krótszy: P
function Rand7(){for(x=0,i=1;i<8;x^=i*((k=Rand5())%2),i*=1+(k<5));return x?x:Rand7()}
Wiem, że odpowiedź jest krótsza, ale chciałem pokazać test tej układanki. Okazuje się, że tylko odpowiedź Clyde Lobo przy użyciu próbkowania odrzucenia Dana McGratha jest poprawna (między odpowiedziami JS).
int Rand7()
{
int r = Rand5();
int n = 5;
do {
r = (r - 1) * 5 + Rand5();
int m = n * 5 / 7 * 7;
if (r <= m) {
return r % 7 + 1;
}
r -= m;
n = n * 5 - m;
} while (1);
}
Rozkład liczb (1000000 liczb całkowitych):
142935 142751 142652 143299 142969 142691 142703
Średnia liczba wywołań do Rand5 () na każdą wygenerowaną liczbę całkowitą wynosi około 2,2 (od 2 do 10+).
1 2 3 4 5 6 7 8 9 10
0 840180 112222 44433 2212 886 0 60 6 1
Inne rozwiązanie, które może być niepoprawne, w Pythonie:
rand7 = lambda: sum(rand5() for i in range(7)) % 7 + 1
To wydaje się zbyt proste, ale kiedy próbuję:
counter = [0] * 7
for i in range(100000):
counter[rand7()] += 1
Dostaję dość równomierną dystrybucję (wszystkie między 14000 a 14500).
Dobra, teraz, gdy ktoś głosował na ten post: Czy to rozwiązanie jest rzeczywiście poprawne? Bardziej opublikowałem to tutaj, aby ludzie go krytykowali. Cóż, jeśli jest to poprawne, moja wersja w golfa byłaby:
rand7=lambda:eval("+rand5()"*7)%7+1
który wychodzi na 37 znaków.
Java, 65 znaków:
int rand7(){int r;do{r=rand5()+5*rand5()-5;}while(r>7);return r;}
def rand7():
while True:
n=5*(rand5()-1)+(rand5()-1)
if n<21:return n%7+1
ale całkowicie poprawne na podstawie uzasadnienia tutaj .
sub rand7{1while($_=5*&rand5-rand5)>6;$_+1}
Daje to ostrzeżenie Ambiguous use of -rand5 resolved as -&rand5()
, ale działa poprawnie. Przygotowanie &
również do drugiego rand5
połączenia naprawia go kosztem jednego pociągnięcia. (I odwrotnie, drugi &
można również usunąć, jeśli rand5
został zdefiniowany za pomocą()
prototypu.)
Ps. Następująca 46-znakowa wersja jest około trzy razy szybsza:
sub rand7{1while($_=5*&rand5-rand5)>20;$_%7+1}
int rand7(){int s;while((s=rand5()*5+rand5())<10);return(s%7+1);}
Dłuższa niż poprzednia procedura, ale myślę, że ta zwraca równomiernie rozłożone liczby w krótszym czasie.
PostScript (46)
Używa binarnego kodowania tokenów, dlatego oto zrzut heksowy:
00000000 2f 72 61 6e 64 37 7b 38 7b 92 38 37 92 61 7b 92 |/rand7{8{.87.a{.|
00000010 40 7d 69 66 92 75 32 7b 72 61 6e 64 35 7d 92 83 |@}if.u2{rand5}..|
00000020 35 92 6c 92 01 35 92 a9 7d 92 65 7d 92 33 |5.l..5..}.e}.3|
0000002e
Aby go wypróbować, możesz go również pobrać .
Oto niepoznany i skomentowany kod wraz z kodem testowym.
% This is the actual rand7 procedure.
/rand7{
8{ % potentialResult
% only if the random number is less than or equal to 7, we're done
dup 7 le{ % result
exit % result
}if % potentialResult
pop % -/-
2{rand5}repeat % randomNumber1 randomNumber2
5 mul add 5 sub % randomNumber1 + 5*randomNumber2 - 5 = potentialResult
}loop
}def
%Now, some testing code.
% For testing, we use the built-in rand operator;
% Doesn't really give a 100% even distribution as it returns numbers
% from 0 to 2^31-1, which is of course not divisible by 5.
/rand5 {
rand 5 mod 1 add
}def
% For testing, we initialize a dict that counts the number of times any number
% has been returned. Of course, we start the count at 0 for every number.
<<1 1 7{0}for>>begin
% Now we're calling the function quite a number of times
% and increment the counters accordingly.
1000000 {
rand7 dup load 1 add def
}repeat
% Print the results
currentdict{
2 array astore ==
}forall
int result = 0;
for (int i = 0; i++; i<7)
if (((rand(5) + rand(5)) % 2) //check if odd
result += 1;
return result + 1;
Zdefiniuj rand7:
rand7=function(n)sample(7,n,T)
Ponieważ R zostało napisane z myślą o analizie statystycznej, to zadanie jest banalne, a ja korzystam z wbudowanej funkcji sample
z zamiennikiem ustawionym na PRAWDA.
Przykładowe dane wyjściowe:
> rand7(20)
[1] 4 3 6 1 2 4 3 2 3 2 5 1 4 6 4 2 4 6 6 1
> rand7(20)
[1] 1 2 5 2 6 4 6 1 7 1 1 3 7 6 4 7 4 2 1 2
> rand7(20)
[1] 6 7 1 3 3 1 5 4 3 4 2 1 5 4 4 4 7 7 1 5
Co powiesz na to?
int Rand7()
{
return Rand5()+ Rand5()/2;
}
/
operator wykonuje matematykę całkowitą? Co stanie się z wynikami, jeśli wykona matematy dziesiętny, zmiennoprzecinkowy lub całkowity?
[2/25, 4/25, 5/25, 5/25, 5/25, 3/25, 1/25]
. Niezupełnie jednolity.
int m=0;int rand7(){return(m=m*5&-1>>>1|rand5())%7+1;}
Test dystrybucji:
[1000915, 999689, 999169, 998227, 1001653, 1000419, 999928]
Algorytm:
> Liczby nie są już wzajemnie nieskorelowane, ale indywidualnie całkowicie przypadkowe.
Kod C / C ++ kod podstawowy ma tylko jedną linię!
static unsigned int gi = 0;
int rand7()
{
return (((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1;
}
//call this seed before rand7
//maybe it's not best seed, if yo have any good idea tell me please
//and thanks JiminP again, he remind me to do this
void srand7()
{
int i, n = time(0);
for (i = 0; i < n % 7; i++)
rand7();
}
Srand7 () jest zalążkiem rand7, musi wywołać tę funkcję przed rand7, podobnie jak wywołanie srand przed randem w C.
Jest to bardzo dobre, ponieważ wywołuje rand () tylko jeden raz, i nie ma pętli, nie zużywa dodatkowych wspomnień.
Pozwól mi to wyjaśnić: rozważ tablicę liczb całkowitych o rozmiarze 5:
1st get one number from 1 2 3 4 5 by rand5
2nd get one number from 2 3 4 5 6
3rd get one number from 3 4 5 6 7
4th get one number from 4 5 6 7 1
5th get one number from 5 6 7 1 2
5th get one number from 6 7 1 2 3
7th get one number from 7 1 2 3 4
Mamy więc TABELĘ, każda z 1-7 pojawia się w niej 5 razy i ma wszystkie 35 liczb, więc prawdopodobieństwo każdej liczby wynosi 5/35 = 1/7. I następnym razem
8th get one number from 1 2 3 4 5
9th get one number from 2 3 4 5 6
......
Po wystarczającej liczbie razy możemy uzyskać równomierny rozkład 1-7.
Możemy więc przydzielić tablicę, aby przywrócić pięć elementów 1-7 za pomocą przesunięcia w lewo-w pętli i uzyskać za każdym razem jedną liczbę z rand5. Zamiast tego możemy wygenerować wszystkie siedem tablic wcześniej i używać ich cyklicznie. Kod jest również prosty, ma wiele krótkich kodów, które mogą to zrobić.
Ale możemy użyć właściwości operacji%, więc wiersze tabeli 1-7 są równoważne z (rand5 + i)% 7, to znaczy: a = rand ()% 5 + 1 to rand5 w języku C, b = gi ++ % 7 generuje wszystkie permutacje w powyższej tabeli, a 0–6 zastępuje 1–7 c = (a + b)% 7 + 1, generuje 1–7 równomiernie. Wreszcie mamy ten kod:
(((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1
Ale nie możemy uzyskać 6 i 7 przy pierwszym wywołaniu, więc potrzebujemy zarodka, takiego jak srand dla rand w C / C ++, aby zmienić rozkład permutacji dla pierwszego oficjalnego połączenia.
Oto pełny kod do testowania:
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
static unsigned int gi = 0;
//a = rand() % 5 + 1 is rand5 in C language,
//b = gi++ % 7 generates all permutations,
//c = (a + b) % 7 + 1, generates 1 - 7 uniformly.
//Dont forget call srand7 before rand7
int rand7()
{
return (((rand() % 5 + 1) + (gi++ % 7)) % 7) + 1;
}
//call this seed before rand7
//maybe it's not best seed, if yo have any good idea tell me please
//and thanks JiminP again, he remind me to do this
void srand7()
{
int i, n = time(0);
for (i = 0; i < n % 7; i++)
rand7();
}
void main(void)
{
unsigned int result[10] = {0};
int k;
srand((unsigned int)time(0)); //initialize the seed for rand
srand7() //initialize the rand7
for (k = 0; k < 100000; k++)
result[rand7() - 1]++;
for (k = 0; k < 7; k++)
printf("%d : %.05f\n", k + 1, (float)result[k]/100000);
}
6
lub 7
zadzwonić raz ?
int main(){if(rand7()==6) printf("Hello, world!");}
, przybliżenie za pomocą pętli wyświetli „Witaj, świecie!” 1 na 7 razy, ale twój kod nie.