Получить вторые минимальные значения на столбец в 2D массиве


15

Как я могу получить второе минимальное значение из каждого столбца? У меня есть этот массив:

A = [[72 76 44 62 81 31]
     [54 36 82 71 40 45]
     [63 59 84 36 34 51]
     [58 53 59 22 77 64]
     [35 77 60 76 57 44]]

Я хотел бы иметь вывод, как:

A = [54 53 59 36 40 44]

ты что-нибудь пробовал? ?
Меха Парех

второй минимум на столбец ?
Николя Жерве

@NicolasGervais да
г-н Дан

Ответы:


12

Попробуйте это всего одной строкой:

[sorted(i)[1] for i in zip(*A)]

в действии:

In [12]: A = [[72, 76, 44, 62, 81, 31], 
    ...:      [54 ,36 ,82 ,71 ,40, 45], 
    ...:      [63 ,59, 84, 36, 34 ,51], 
    ...:      [58, 53, 59, 22, 77 ,64], 
    ...:      [35 ,77, 60, 76, 57, 44]] 

In [18]: [sorted(i)[1] for i in zip(*A)]                                                                                                                                                                           
Out[18]: [54, 53, 59, 36, 40, 44]

zip(*A) перенесет ваш список списка, чтобы столбцы стали строками.

и если у вас есть повторяющееся значение, например:

In [19]: A = [[72, 76, 44, 62, 81, 31], 
    ...:  [54 ,36 ,82 ,71 ,40, 45], 
    ...:  [63 ,59, 84, 36, 34 ,51], 
    ...:  [35, 53, 59, 22, 77 ,64],   # 35
    ...:  [35 ,77, 50, 76, 57, 44],]  # 35

Если вам нужно пропустить оба параметра 35, вы можете использовать set():

In [29]: [sorted(list(set(i)))[1] for i in zip(*A)]                                                                                                                                                                
Out[29]: [54, 53, 50, 36, 40, 44]

6

Операции над numpyмассивами должны выполняться с помощью numpyфункций, поэтому посмотрите на это:

np.sort(A, axis=0)[1, :]
Out[61]: array([54, 53, 59, 36, 40, 44])

Насколько мне известно, это должно быть лучшим решением, оно сохраняет все numpy, я думаю, что решение lambdaдолжно быть замедлено heapq.nsmallest. Кажется, лучше держать все как можно быстрееnumpy
jamylak

3

Вы можете использовать heapq.nsmallest

from heapq import nsmallest

[nsmallest(2, e)[-1] for e in zip(*A)]

вывод:

[54, 53, 50, 36, 40, 44]

Я добавил простой тест для сравнения производительности различных уже опубликованных решений:

введите описание изображения здесь

from simple_benchmark import BenchmarkBuilder
from heapq import nsmallest


b = BenchmarkBuilder()

@b.add_function()
def MehrdadPedramfar(A):
    return [sorted(i)[1] for i in zip(*A)]

@b.add_function()
def NicolasGervais(A):
    return np.sort(A, axis=0)[1, :]

@b.add_function()
def imcrazeegamerr(A):
    rotated = zip(*A[::-1])

    result = []
    for arr in rotated:
        # sort each 1d array from min to max
        arr = sorted(list(arr))
        # add the second minimum value to result array
        result.append(arr[1])

    return result

@b.add_function()
def Daweo(A):
    return np.apply_along_axis(lambda x:heapq.nsmallest(2,x)[-1], 0, A)

@b.add_function()       
def kederrac(A):
    return [nsmallest(2, e)[-1] for e in zip(*A)]


@b.add_arguments('Number of row/cols (A is  square matrix)')
def argument_provider():
    for exp in range(2, 18):
        size = 2**exp
        yield size, [[randint(0, 1000) for _ in range(size)] for _ in range(size)]

r = b.run()
r.plot()

Использование zipс sortedфункцией - самое быстрое решение для небольших 2-мерных списков, а использование zipс heapq.nsmallestшоу - для лучших в больших 2-мерных списках.


1
Просто дикая мысль: могут ли эти результаты повлиять на тот факт, что вы генерировали числа, которые не являются типичными dtypes? Кроме того, не будет ли встроенный randint возвращать список вместо массива?
Николя Жерве

1

Я надеюсь, что я правильно понял ваш вопрос, но в любом случае вот мое решение, я уверен, что есть более гибкий способ сделать это, но это работает

A = [[72,76,44,62,81,31]
 ,[54,36,82,71,40,45]
 ,[63,59,84,36,34,51]
 ,[58,53,59,22,77,64]
 ,[35,77,50,76,57,44]]

#rotate the array 90deg
rotated = zip(*A[::-1])

result = []
for arr in rotated:
    # sort each 1d array from min to max
    arr = sorted(list(arr))
    # add the second minimum value to result array
    result.append(arr[1])
print(result)

введите описание изображения здесь


0

Предполагая, что Aэто так numpy.array(если это верно, рассмотрите возможность добавления numpyтега к своему вопросу), то вы можете использовать apply_along_axisэтот способ следующим образом:

import heap
import numpy as np
A = np.array([[72, 76, 44, 62, 81, 31],
              [54, 36, 82, 71, 40, 45],
              [63, 59, 84, 36, 34, 51],
              [58, 53, 59, 22, 77, 64],
              [35, 77, 60, 76, 57, 44]])
second_mins = np.apply_along_axis(lambda x:heapq.nsmallest(2,x)[-1], 0, A)
print(second_mins)  # [54 53 59 36 40 44]

Обратите внимание, что я использовал heapq.nsmallest, так как он выполняет столько сортировки, сколько требуется для получения 2 самых маленьких элементов, в отличие от того, sortedкоторый выполняет полную сортировку.


0
>>> A = np.arange(30).reshape(5,6).tolist()
>>> A
[[0, 1, 2, 3, 4, 5], 
 [6, 7, 8, 9, 10, 11], 
 [12, 13, 14, 15, 16, 17], 
 [18, 19, 20, 21, 22, 23],
 [24, 25, 26, 27, 28, 29]]

Обновлено : используйте setдля предотвращения дублирования и транспонирования списка с помощьюzip(*A)

>>> [sorted(set(items))[1] for items in zip(*A)]
[6, 7, 8, 9, 10, 11]

old: второй минимальный элемент в каждой строке

>>> [sorted(set(items))[1] for items in A]
[1, 7, 13, 19, 25]

Разве это не получение второго элемента в каждой строке, а не в столбце?
paxdiablo

@paxdiablo Да, спасибо, что сообщили. обновленный ответ.
Дишин Х Гояни
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.