Blum-Blum-Shub-Generator

Der Blum-Blum-Shub-Generator (BBS-Generator; auch „s² mod n - Generator“) ist ein Pseudozufallszahlengenerator, entwickelt 1986 von Lenore Blum, Manuel Blum und Michael Shub. Anwendung findet das System u. a. in der Kryptologie im Entwurf komplexitätstheoretisch sicherer Kryptosysteme.

Definition

Der BBS-Generator ist definiert als Folge (s_i) durch die Iterationsvorschrift

$s_0 \!\, = s^2 \; \bmod \; n$,

$s_{i+1} \!\, = s_i^2 \; \bmod \; n$ (dabei bezeichnet $\mod$ den Divisionsrest; siehe Modulo).

Der Modul n = pq ist das Produkt zweier verschiedener Primzahlen p und q, die von der Form 4k + 3 sind, d. h. $p \equiv q \equiv 3 \; ( \bmod \; 4 )$. Eine Zahl n mit diesen Eigenschaften wird auch Blum-Zahl genannt. Der Startwert s ist zu n teilerfremd: $\operatorname{ggT}(s,n) = 1$.

Der Parameter n sollte außerdem folgenden Bedingungen genügen, damit n möglichst schwer zu faktorisieren ist und der Generator garantiert hochwertige Zufallszahlen erzeugt (siehe Kapitel „Sicherheit“):

• n sollte hinreichend groß sein; für kryptografische Anwendung mindestens etwa 200 Dezimalstellen.
• p und q sollten etwa gleichdimensioniert sein, aber nicht zu nah beieinander liegen, etwa 2 < p / q < 1000.
• p − 1 und p + 1 sowie q − 1 und q + 1 sollten jeweils einen großen Primfaktor haben, größer als ca. $\sqrt[4]{n}$.

Manchmal ist es praktisch, k Iterationsschritte des Generators auf einmal zu berechnen. Dies geht mit der Formel

$s_{i+k} = s_i^{2^k \, \bmod \, \lambda(n)} \; \bmod \; n$, wobei hier $\lambda(n) = \operatorname{kgV}(p-1,q-1)$, siehe kgV.

Periodenlänge

Sei Z die Menge der zu n teilerfremden Zahlen in $\{1,\dots,n-1\}$.

Der BBS-Generator arbeitet auf der Menge $Q \subset Z$ der Quadratreste modulo n. Die s_i sind trivialerweise in Q, da sie als Rest des Quadrats einer zu n teilerfremden Zahl berechnet werden. Q enthält genau ein Viertel der Zahlen in Z. Jedes $x \in Q$ hat genau eine Wurzel in Q: $W_x = \{w \in Z \; | \; w^2 \; \bmod \; n = x \} \; ; \; |W_x \cap Q| = 1$. Die Iterationsvorschrift des Generators bildet also Q bijektiv auf Q ab. Somit zerfällt Q in mehrere Teilmengen Q_j, die jeweils eine Periode des Generators bilden. Die Periodenlänge $\left| Q_j \right|$ ist immer ein Teiler von λ(λ(n)), wobei λ die Carmichael-Funktion ist.

Beispiel

Sei $n = 3 \cdot 7 = 21$. Dann ist $Z \!\, = \{1,2,4,5,8,10,11,13,16,17,19,20\}$ und die Quadratreste sind $Q \!\, = \{1,4,16\}$. Die Wurzeln der Quadratreste sind

$W_1 \!\, = \{ 1,8,13,20 \}$ mit $1 \!\, \in Q$,

$W_4 \!\, = \{ 2,5,16,19 \}$ mit $16 \!\, \in Q$ und

$W_{16} \!\, = \{ 4,10,11,17 \}$ mit $4 \!\, \in Q$.

Damit haben wir die beiden Perioden $Q_1 \!\, = \{ 1 \}$ und $Q_2 \!\, = \{ 4, 16 \}$.

Nachweis der Periodenlänge

Es ist schwierig, die Parameter s und n so zu bestimmen, dass eine ausreichende Periodenlänge garantiert ist. Bei der Verwendung des BBS-Generators in der Kryptographie wird dieses Problem oft vernachlässigt, denn die Wahrscheinlichkeit einer zu kurzen Periode ist sehr klein. Absolute Sicherheit kann ohnehin nicht erreicht werden, da ein Angreifer den Modul n mit Glück faktorisieren könnte, etwa indem er einige Millionen zufällig erzeugte Probeteiler ausprobiert.

Wenn n und s so gewählt werden, dass gilt:

$\operatorname{ord}_{\lambda(n)/2}(2) = \lambda(\lambda(n))$, und

$\operatorname{ord}_n(s_0) = \lambda(n)/2$,

dann beträgt die Periodenlänge λ(λ(n)).

$\operatorname{ord}_b(a)$ bezeichnet dabei die Ordnung des Elements a der primen Restklassengruppe $(\Z / b \Z)^\times$:

$\operatorname{ord}_b(a) = \min \{ m \in \N^+ \; | \; a^m \equiv 1 \; (\bmod \; b) \}$.

Zur effizienten Berechnung kann ausgenutzt werden, dass die Elementordnung laut dem Satz von Lagrange ein Teiler der Gruppenordnung sein muss:

$\operatorname{ord}_b(a) \; | \; \phi(b)$.

Dafür muss die Faktorisierung der Gruppenordnung ϕ(b) bekannt sein (siehe Eulersche φ-Funktion).

Mann muss n also so konstruieren, dass die Faktorisierungen von p − 1 und q − 1 bekannt sind oder mit vertretbarem Aufwand berechnet werden können, und ebenso die Faktorisierungen der um 1 verminderten Primfaktoren von p − 1 und q − 1. Damit können die benötigten Größen und die Faktoren der Gruppenordnungen effizient bestimmt werden. Mit der binären Exponentiation kann man anschließend jeweils $a^m \; (\bmod \; b)$ für alle Teiler m von ϕ(b) effizient berechnen.

Anwendung

Erzeugung von Zufallsbits

Aus jedem s_i werden ein oder mehrere Zufallsbits gewonnen. Im einfachsten Fall nimmt man das niederwertigste Bit, also

$b_i = s_i \; \bmod \; 2$,

oder man berechnet das Paritätsbit zu s_i:

$b_i = \operatorname{bitnum}(s_i) \; \bmod \; 2$.

Die Funktion $\operatorname{bitnum}(x)$ liefert die Zahl der Bits mit dem Wert 1 in der Binärdarstellung von x.

Eine weitere Möglichkeit ist die Bestimmung des Positionsbits, das von der Position von s_i im Intervall [1,n − 1] abhängt:

$b_i = \begin{cases} 1, & \text{wenn } s_i > \lfloor n/2 \rfloor\\ 0, & \text{wenn } s_i \leq \lfloor n/2 \rfloor \end{cases}$.

Am besten ist es jedoch, wenn das Paritätsbit von einigen fest gewählten Bits aus s_i bestimmt wird. Dazu wählt man vorab eine Konstante z als Maske, die etwa so groß wie n ist und eine unregelmäßige, „zufällige“ Binärdarstellung aufweist, und berechnet

$b_i = \operatorname{bitnum}(s_i \wedge z) \; \bmod \; 2$.

Dabei bezeichnet $\wedge$ die bitweise UND-Verknüpfung.

Aus einem s_i kann man mehrere Zufallsbits erhalten. Die Erfinder Blum, Blum und Shub haben schon früh vorgeschlagen, das niederwertigste Bit und das Positionsbit zugleich zu nutzen:

$b_{2i} = s_i \; \bmod \; 2, \quad b_{2i+1} = \begin{cases} 1, & \text{wenn } s_i > \lfloor n/2 \rfloor\\ 0, & \text{wenn } s_i \leq \lfloor n/2 \rfloor \end{cases}$.

Man kann zeigen, dass der BBS-Generator kryptografisch auch dann noch sicher ist, wenn bis zu $a = \lfloor \log_2 \log_2 n \rfloor$ Bits aus jedem s_i extrahiert werden. Meist werden einfach die a niederwertigsten Bits genommen:

$b_{ai+j} = \lfloor s_i \, / \, 2^j \rfloor \; \bmod \; 2 \quad \text{mit} \; j = 0, 1, \ldots , a-1$,

oder etwas elaborierter, mit „disjunkten“ Masken z_j:

$b_{ai+j} = \operatorname{bitnum}(s_i \wedge z_j) \; \bmod \; 2 \quad \text{mit} \; j = 0, 1, \ldots , a-1 \; ; \quad j \neq k \Rightarrow z_j \wedge z_k = 0$.

Symmetrisches Kryptosystem

Zunächst wird der BBS-Generator zur Umsetzung einer Stromchiffre verwendet. Als geheimer Schlüssel zwischen Sender und Empfänger dienen n und der Startwert s des Generators.

Z. B. generiert der Sender aus n=7$\cdot$11=77 und s=64 nach der oben angegebenen Vorschrift die Folge der s_i. Die zugehörige Pseudozufallszahl b_i ergibt sich beispielsweise aus dem letzten Bit des jeweiligen Wertes von s_i, d. h. b_i=s_i mod 2. Um den Schlüsseltext zu bestimmen, wird der Klartext (im Beispiel: 0011) XOR mit der Pseudozufallszahlenfolge verknüpft.

 Generierte Folge         15 71 36 64 …
 Pseudozufallszahlenfolge  1  1  0  0 …
 Klartext                  0  0  1  1
 Schlüsseltext             1  1  1  1

Der Empfänger bestimmt seinerseits aus den geheimen Werten n und s die Folgen s_i und b_i. Mit Hilfe des übersendeten Schlüsseltextes wird wiederum mittels XOR der Klartext berechnet.

 Generierte Folge         15 71 36 64 …
 Pseudozufallszahlenfolge  1  1  0  0 …
 Schlüsseltext             1  1  1  1
 Klartext                  0  0  1  1

Asymmetrisches Kryptosystem

Zur Umsetzung eines asymmetrischen Kryptosystems eignet sich der BBS-Generator ebenfalls. Dieses Verfahren wurde 1984 von Manuel Blum und Shafi Goldwasser vorgeschlagen und wird auch als Blum-Goldwasser-Kryptosystem bezeichnet. Der geheime Schlüssel auf Seiten des Empfängers sind die Primfaktoren p und q.

Senderseitig laufen die Berechnungen analog zum obigen symmetrischen Fall ab. Zusätzlich zum Schlüsseltext $x_0 \oplus b_0, \ldots, x_i \oplus b_i$ wird aber noch s_i+1 gesendet. Da der Empfänger den Startwert nicht kennt, bildet er mit Hilfe der geheimen Primzahlen p und q die Folge der Pseudozufallszahlen ausgehend vom versendeten s_i+1 bis zum Startwert s zurück. Für das Beispiel bedeutet das, der Empfänger erhält 1111 sowie s₄=15.

s_i-1 = (up$\cdot$s_i^(q+1)/4 mod q + vq$\cdot$s_i^(p+1)/4 mod p) mod n

Der Ansatz bedient sich des Chinesischen Restealgorithmus, einem Spezialfall des chinesischen Restsatzes. Die beiden Unbekannten u und v sind von den Primfaktoren p und q abhängig und werden zu Beginn mittels des erweiterten Euklidischen Algorithmus bestimmt. Dabei gilt up+vq=1, also 2$\cdot$11-3$\cdot$7=1 im Beispiel. Damit ergibt sich die folgende Abarbeitung.

s₃ = (22$\cdot$15² mod 7 - 21$\cdot$15³ mod 11) mod 77

s₃ = (22$\cdot$1 - 21$\cdot$9) mod 77 = 64

s₂ = (22$\cdot$64² mod 7 - 21$\cdot$64³ mod 11) mod 77

s₂ = (22$\cdot$1 - 21$\cdot$3) mod 77 = 36

s₁ = (22$\cdot$36² mod 7 - 21$\cdot$36³ mod 11) mod 77

s₁ = (22$\cdot$1 - 21$\cdot$5) mod 77 = 71

s₀ = (22$\cdot$71² mod 7 - 21$\cdot$71³ mod 11) mod 77

s₀ = (22$\cdot$1 - 21$\cdot$4) mod 77 = 15

s = (22$\cdot$15² mod 7 - 21$\cdot$15³ mod 11) mod 77

s = (22$\cdot$1 - 21$\cdot$5) mod 77 = 64

Empfängerseitig wird nun analog zum symmetrischen Fall aus der eben rückwärts berechneten BBS-Generatorfolge die Folge der Pseudozufallszahlen bestimmt und letztlich durch XOR-Verknüpfung mit dem Schlüsseltext der Klartext generiert.

Ein so konstruiertes asymmetrisches Kryptosystem ist jedoch nicht sicher gegen aktive Angreifer, z. B. durch einen Angriff mit frei wählbarem Geheimtext (englisch: chosen-ciphertext attack).

Sicherheit

Die Sicherheit des BBS-Generators basiert auf der Faktorisierungsannahme (FA).^[1] Jeder, der BBS brechen kann, kann auch faktorisieren, was aber als praktisch unmöglich gilt. Folglich ist BBS sicher.

Faktorisierungsannahme (FA): Die Wahrscheinlichkeit, dass ein schnelles Faktorisierungsverfahren eine ganze Zahl n=pq mit Erfolg faktorisiert, sinkt rapide mit zunehmender Länge der Faktoren p und q.

Zur Zeit kann keine sichere Aussage getroffen werden, wie schwer Faktorisierung ist. Mit anderen Worten, die Frage nach einem Algorithmus, der in annehmbarer Zeit bei Eingabe beliebiger n die Primfaktorzerlegung in p und q durchführt, bleibt unbeantwortet. Somit kann die Problematik lediglich mit Hilfe einer Annahme abgeschätzt werden.

Für konkrete praktische Anwendungen fordert man dann, dass bei gegebener Länge der Primfaktoren nur ein bestimmter Teil in einer bestimmten Zeit mit maximal verfügbarer Rechnerkapazität und den besten bekannten Faktorisierungsverfahren faktorisiert werden kann, also z. B. bei einer Länge von 1024 Bit werden 2^-50 % aller n in einem Jahr faktorisiert.

Wer faktorisieren kann, kann auch BBS brechen. Faktorisieren ermöglicht das Brechen der Quadratischen-Reste-Annahme, was es erlaubt, die Pseudozufallsfolge vorherzusagen.

Quadratische-Reste-Annahme (QRA) (englisch: quadratic residuosity assumption): Es ist schwierig (im Sinne von aufwändig), von einer gegebenen Zahl a zu entscheiden, ob sie ein Quadratischer Rest in einem Restklassenring $\mathbb{Z}/n\mathbb{Z}$ ist, d.h. ob es eine Zahl b gibt, so dass $a \equiv b^2 \pmod n$ ist. Die QRA ist wie die FA nicht bewiesen.

Zwei Punkte erschweren diesen Test. Erstens gibt es in einem Restklassenring mehrere Wurzeln zu einer gegebenen Zahl. So haben z. B. im $\mathbb{Z}/4\mathbb{Z}$ die Zahlen 1 und 3 die gleichen Quadrate: $1^2 \equiv 3^2 \equiv 1 \pmod 4$. Zweitens interessiert man sich nur für solche Quadrate, die selbst Quadrate sind. Diesen Umstand kann man sich mittels der Definition der BBS-Generatorfolge verdeutlichen.

Zusammenfassend gilt daher: Die Sicherheit des BBS-Generators ist äquivalent zur Faktorisierungsannahme.

Einzelnachweise

1. ↑ Andrey Sidorenko und Berry Schoenmakers: Concrete Security of the Blum-Blum-Shub Pseudorandom Generator. In: Cryptography and Coding 2005. 2005, S. 355–375 (http://www.win.tue.nl/~asidoren/bbs.pdf).

Weblinks

• Lenore Blum, Manuel Blum, und Michael Shub. "A Simple Unpredictable Pseudo-Random Number Generator", SIAM Journal on Computing, Band 15, Seiten 364–383, Mai 1986.
• Lenore Blum, Manuel Blum, und Michael Shub. "Comparison of two pseudo-random number generators", Advances in Cryptology: Proceedings of Crypto '82.
• Martin Geisler, Mikkel Krøigård, und Andreas Danielsen. "About Random Bits", Dezember 2004. Als PDF und Gzipped Postscript.

Arithmetische Zufallszahlengeneratoren

Kongruenzgenerator (linearer, multiplikativer, gemischt linearer, Fibonacci-Generator) | Inverser Kongruenzgenerator | Blum-Blum-Shub-Generator | Mersenne-Twister | WELL | Xorshift