:-D

[Leonidas]

So, welche sind die mittleren Ziffern?

Such Dir eben eine der Möglichkeiten aus

Ja, das war auch genau der Punkt der mir das meiste Kopfzerbrechen bei meiner Lösung gekostet hat, weil ich mich nicht entscheiden konnte und das wohl nicht definiert war. Ich finde das aber recht blöd, weil es damit quasi zwei Arten des Algorithmus gibt, die unterschiedliche Ergebnisse liefern.

[Karl]

Ja, deshalb ist das korrekterweise auch kein Algorithmus.
Jedenfalls ist er zu unpräzise beschrieben.
Ich hab gerade angefangen, das in Brainfuck zu schreiben, aber ich bin noch nicht wirklich weit :p das is ne richtige Qual.

[numerix]

Es scheint so zu sein, dass man bei der Quadratmittenmethoden nach Neumann, um die es hier ja geht, in solchen Fällen, die "Mitte" nach rechts verschiebt.

Beleg - leider ohne weitere Erläuterung, warum gerade dies die Mitte ist:
http://statistik.mathematik.uni-wuerzbu ... seudo.html

[BlackJack]

Eine Lösung in C, die `libgmp` verwendet. Mittelpunktberechnung ist aber bei meinen Lösungen auf die hypothetischen doppelt so vielen Ziffern der Ausgangszahl bezogen, die beim Quadrieren entstehen können. Alles andere war mir zu kompliziert, weil man sich Sonderfälle einhandelt, wenn der aktuelle Zustand weniger als die hälfte der Ziffern der Ausgangszahl hat.

Code: Alles auswählen

#include <alloca.h>
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>

typedef struct {
    mpz_t state;
    unsigned int length;
} QRandom;

void QRandom_init(QRandom *r, unsigned long n)
{
    mpz_init_set_ui(r->state, n);
    r->length = mpz_sizeinbase(r->state, 10);
}

void QRandom_destroy(QRandom *r)
{
    mpz_clear(r->state);
    r->length = 0;
}

unsigned long QRandom_next(QRandom *r)
{
    char *tmp;
    unsigned int offset;
    
    mpz_mul(r->state, r->state, r->state);
    tmp = alloca(r->length * 2 + 1);
    mpz_get_str(tmp, 10, r->state);
    offset = r->length / 2;
    tmp[offset + r->length] = '\0';
    mpz_set_str(r->state, tmp + offset, 10);
    return mpz_get_ui(r->state);
}

int main(int argc, char** argv) {
    QRandom random;
    unsigned int i;
    
    QRandom_init(&random, 11111);
    for (i = 0; i < 20; ++i) {
        printf("%lu\n", QRandom_next(&random));
    }
    QRandom_destroy(&random);
    
    return EXIT_SUCCESS;
}

Wobei mir gerade auffällt, dass vielleicht auch mal eine 0 oder eine 1 als Zustand auftreten kann, dann sieht's blöd aus mit den folgenden "Zufallszahlen".

[mitsuhiko]

Und hier in Vimscript: http://paste.pocoo.org/show/73064/

Muss mal wieder Schmerzen spüren

[dennda]

Hier noch ein bisschen sehr schlechtes Java.

Wer mich dafür kritisiert wird erschossen.

[audax]

Rekursiv in Java?
Mutig ;D

[Hyperion]

Rekursiv in Java?
Mutig ;D

Wieso?

@dennda:
Du könnest Dir den häßlichen Cast in Zeile 25 sparen:

Code: Alles auswählen

String subsequence = cast.substring(mid-2, mid+2);

[mitsuhiko]

Und weil's so schön ist nochmal in C:

Code: Alles auswählen

#include <stdio.h>
#include <string.h>
#include <stdlib.h>


int
nrnd(int n, int x)
{
        int l, nl, offset, i;
        char s[100], s2[100];
        sprintf(s, "%d", n);
        l = strlen(s);

        for (i = 0; i != x; ++i) {
                sprintf(s, "%d", n * n);
                nl = (strlen(s) - l) / 2;
                memcpy(s2, s + nl, l);
                s2[l] = '\0';
                n = atoi(s2);
        }

        return n;
}