gran problema de asignación de bajo rango denso


9

¿Existe un método razonablemente barato para resolver el problema de asignación grande, denso y de bajo rango , donde ejecuta sobre todas las permutaciones. De ?maxπiAπi,i1 : nπ1:n

Aquí es una matriz de bajo rango . Los tamaños típicos serían (posiblemente mucho más grande), .n × n r n = 10000 r = 15An×nrn=10000  r=15


1
Por ¿te refieres al producto para que estés avanzando a través de la matriz por diferentes ? ππiπ
Bill Barth

π recorre todas las permutaciones.
Arnold Neumaier

¿No debería ser entonces? Aπ(i),i
Jack Poulson

@JackPoulson: y son dos notaciones diferentes para la imagen de bajo la permutación . π i i π\i(i)πiiπ
Arnold Neumaier

¡Interesante pregunta! ¿Está buscando para explotar la estructura de bajo rango sólo para almacenamiento razones --- es decir, para salvar de tener que formar la matriz entera cuando se aplica un algoritmo de asignación tradicional? ¿O estás buscando una manera de explotar el rango bajo para acelerar la búsqueda?
Michael Grant

Respuestas:


3

Como con R 1 , R 2R n × r , tenemos i A π i , i = i ( P π A ) i , i = trace ( P π R 1 R T 2 ) donde P π es la matriz de permutación correspondiente a π .A=R1R2TR1,R2Rn×r

iAπi,i=i(PπA)i,i=trace(PπR1R2T)
Pππ

Para cualquier , la traza se puede calcular como traza ( P π R 1 R T 2 ) = i k ( P π R 1 ) i , k ( R T 2 ) k , i = i , k ( ( P π R 1 ) R 2 ) i , k . (Esta cantidad también se conoce comoπ

trace(PπR1R2T)=ik(PπR1)i,k(R2T)k,i=i,k((PπR1)R2)i,k.
Producto Frobenius , ).PπR1:R2

Esta idea no quita la carga de tener que pasar por todas las permutaciones y la búsqueda de fuerza bruta para obtener el máximo de todos los productos Frobenius, y de hecho tiene la misma complejidad aritmética que computar explícitamente . Sin embargo, tiene requisitos de memoria mucho más bajos, ya que nunca tenga que formar en realidad una .A=R1R2TA

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.