[seqfan] Re: help needed with sequences related to Kaprekar map
zak seidov
zakseidov at yahoo.com
Wed Aug 19 21:21:25 CEST 2009
Yes,
62964 is OK:
62964->71973->83952->74943->62964
but 61974 is smaller (and the smallest!)
61974->82962->75933->63954->61974
I used Mathematica.
Zak
--- On Wed, 8/19/09, Richard Mathar <mathar at strw.leidenuniv.nl> wrote:
> From: Richard Mathar <mathar at strw.leidenuniv.nl>
> Subject: [seqfan] Re: help needed with sequences related to Kaprekar map
> To: seqfan at seqfan.eu
> Date: Wednesday, August 19, 2009, 3:14 PM
>
> A151959(4) is 61974 (or smaller), not 62964 as currently
> listed:
> from the map 10026 -> 61974
> -> 82962 -> 75933 -> 63954
> The 62964 is in a cycle of length 4 as observed, but 61974
> is just smaller.
>
> The offset in A151959 should be changed to 1.
>
>
> In Maple 9:
>
> A151949 := proc(n)
> tdigs := convert(n,base,10) ;
> tup := sort(tdigs) ;
> add( (op(i,tup)-op(-i,tup))
> *10^(i-1),i=1..nops(tup)) :
> end:
>
> N := 10 :
> allc := [seq(10000000000000000000000000000000,i=1..N)] :
>
> cyK := proc(n,)
> global allc,N ;
> L := [n] ;
> k := n ;
> while true do
> k := A151949(k) ;
> if k in L then
>
> for s from 1 do
>
> if op(-s,L) = k then
>
> cur := op(s,allc) ;
>
> if k < cur then
>
>
> printf("%a %d %d\n",L,s,k) ;
>
>
>
>
> allc := subsop(s=k,allc) ;
>
> fi;
>
> RETURN(s) ;
>
> fi;
>
> od:
>
> break;
> fi;
> L := [op(L),k] :
> od:
> end:
>
> for n from 1 do
> cyK(n) ;
> od:
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
More information about the SeqFan
mailing list