use strict;
use warnings;
use feature 'say';
use PDL;

use Time::HiRes 'time';
my $t = time;

use constant MAX => 1e6;
use constant TOP => MAX < 20 ? MAX : 20;

my $seqs = 1 + sequence( longlong, MAX );
$seqs-> setbadat( 0 );
$seqs-> badvalue( 2 );

my $lengths = ones( short, MAX );
$lengths <<= 1;
$lengths-> set( 0, 1 );

while ( any my $good_mask = $seqs-> isgood ) {

    my ( $seqs_odd, $lengths_odd_masked )
        = where( $seqs, $lengths, $seqs & 1 );

    $lengths_odd_masked ++;
    $lengths-> where( $good_mask ) ++;
    ( $seqs_odd *= 3 ) ++;
    $seqs >>= 1; 
}

my $sorted_i = $lengths-> qsorti;
my $sorted   = $lengths-> index( $sorted_i );
my $value    = $sorted-> at( MAX - TOP );
my $pos      = vsearch_insert_leftmost( $value, $sorted );
my $top_i    = $sorted_i-> slice([ MAX - 1 , $pos ]);

( my $result = $lengths
    -> index( $top_i )
    -> longlong
    -> bitnot
    -> cat( $top_i + 1 )
    -> transpose
    -> qsortvec
    -> slice([], [ 0, TOP - 1 ])

)-> slice([ 0 ], [])
 -> inplace
 -> bitnot;

say $result;
say time - $t;

__END__

[
 [   525 837799]
 [   509 626331]
...
 [   445 886953]
 [   445 906175]
 [   445 922524]
 [   445 922525]
]

6.0809600353241