Incorporação isométrica de L2 em L1

27

Sabe-se que, dado um subconjunto de pontos de (ou seja, dados pontos em com distância euclidiana), é possível incorporá-los isometricamente em .d 2 n R d ( nn2dnRd1(n2)

A isometria é computável em tempo polinomial (possivelmente aleatório)?

Como existem problemas de precisão finita, a pergunta exata é

Dado um conjunto de pontos em e , existe um mapeamento de computável (possivelmente usando aleatoriedade) em polinômio de tempo em logarítmico em modo que para cada tenhamosn R d ϵ > 0 f : X R ( nXnRdϵ>0 n1/ϵx,yXf:XR(n2)n1/ϵx,yX

||f(x)f(y)||1||xy||2(1+ϵ)||f(x)f(y)||1

(Nota: Estou ciente de que um mapeamento com distorção pode ser encontrado com alta probabilidade de polinomial vez em e projetando em linhas aleatórias, mas não tenho certeza se o número de dimensões pode ser reduzido construtivamente para ou mesmo quando é muito maior que e não sei se existe é um método de tempo polinomial para lidar com o caso em que é exponencial em .)n 1 / ϵ O ( ϵ - 2log n ) ( n(1+ϵ)n1/ϵO(ϵ2logn) O(N2)1/εn1/εn(n2)O(n2)1/ϵn1/ϵn

Luca Trevisan
fonte
1
Esta é uma pergunta muito boa. @Luca, você acha que pode ser difícil? (é claro que meu primeiro pensamento foi para olhar sobre 'Hamming atende Euclides', e então eu vi a identidade do autor da pergunta :)
Suresh Venkat
1
Esta referência parece estar relacionada: Pjotr ​​Indyk, "Princípios da incerteza, extratores e incorporação explícita de l2 em l1", Proc. STOC'07.
Martin Schwarz
2
@ David: é o número de pontos, corrigi o local onde usei para a dimensão. Que pontos no espaço euclidiano (de qualquer dimensão) podem ser incorporados isometricamente em é provado aqui: www-math.mit.edu/~goemans/18409-2006/lec1.pdf, mas não . -constructively (teorema de Carathéodory para ir de dimensão finita, mas grande para dimensão com arbitrariamente pequeno erro, e um argumento de compacidade para ir de arbitrariamente pequeno erro para erro zero.)n n ( nnnn(n1(n2)(n2)
Luca Trevisan
1
@ Martin: obrigado pela referência. O artigo de Piotr lida com o problema mais difícil de mapear todos os (não apenas um conjunto fixo de pontos) para . Para esse problema, acredito que seja um problema aberto até alcançar, construtivamente, e distorção . (Piotr recebe e .) n m 1 m = p o l y ( d , 1 / ε ) ( 1 + ε ) m = d O ( log d ) ε = 1 / d2dn1mm=poly(d,1/ϵ)(1+ϵ)m=dO(logd)ϵ=1/d
Luca Trevisan
1
@LucaTrevisan: re: a dureza de incorporar em L1, isso é verdade (é mencionado no Capítulo 1 ou 2 do livro Deza e Laurent - Acho via MAX CUT)
Suresh Venkat

Respostas:

7

Suresh me pediu para reunir meus comentários acima em uma resposta, então aqui está. Não tenho muita certeza de que seja uma resposta para a pergunta original, pois não é óbvio como torná-lo um tempo polinomial quando a dimensão do espaço euclidiano de entrada é constante. Ele tem pelo menos a vantagem de evitar qualquer problema com como a pergunta original pede, porque não envolve nenhuma aproximação e parece polinomial para a constante .d1/ϵd

De qualquer forma: a partir da geometria integral, há uma medida padrão em conjuntos de hiperplanos em espaço euclidiano dimensional que é invariante sob congruências euclidianos. Tem a propriedade de que o comprimento de qualquer curva de comprimento limitado é proporcional à medida dos hiperplanos que cruzam (com multiplicidade, o que significa que se um hiperplano cruza C duas vezes, isso contribui duas vezes para a medida total dos hiperplanos cruzando C ) Em particular, se C é um segmento de linha, a complicação de multiplicidade não surge e podemos normalizar a medida nos hiperplanos que cruzam C para ser exatamente o comprimento de CC CdCCCCCCC. (Os hiperplanos que contêm têm a medida zero, portanto, não se preocupe com a multiplicidade infinita.)C

Agora, dado um conjunto de n pontos no espaço d-dimensional, faça uma coordenada para cada uma das partições dos pontos em dois subconjuntos induzidos por um hiperplano que não passa por nenhum dos pontos. Informe os pontos de um lado do valor da coordenada da partição zero e os pontos do outro lado do valor da coordenada da partição iguais à medida do conjunto de hiperplanos que induzem essa partição.1

Se e q são quaisquer dois dos n pontos, deixa K ser o conjunto de hiperplanos atravessam segmento de linha P Q , e deixar K i ser os subconjuntos de K formados por cada partição hiperplana possível que tem p de um lado e q no de outros. Então K é a união disjunta de K i , e as diferenças de coordenadas entre p e q são apenas as medidas dos subconjuntos K i . Portanto, o 1pqnKpqKiKpqKKipqKi1A distância entre as coordenadas de e q (a soma das medidas de K i ) é a medida de K , que é apenas a distância 2 original entre p e q .pqKiK2pq

Para geômetros computacionais, uma descrição alternativa da mesma construção pode ser útil: use a dualidade projetiva para transformar os pontos de entrada em n hiperplanos e separar os hiperplanos em pontos. A medida da geometria integral em conjuntos de hiperplanos é então transformada em uma medida mais padrão em conjuntos de pontos, a distância entre p e q se dualiza à medida da cunha dupla entre dois hiperplanos, e o arranjo do hiperplano divide essa cunha dupla em células menores . O valor da coordenada para um ponto é a medida de uma das células no arranjo (se o hiperplano duplo estiver abaixo da célula da coordenada) ou zero (se o hiperplano duplo estiver acima da célula). Portanto, o nnpq distância entre p e q é apenas a soma das medidas das células na cunha dupla, que é a mesma que a medida de toda a cunha dupla. Esse ponto de vista duplo também facilita o cálculo da dimensão da incorporação encontrada dessa maneira: é apenas o número de células na disposição do hiperplano, que é O ( n d ) , ou mais precisamente no máximod i = 0 ( n1pqO(nd) .i=0d(ni)

Até agora, isso fornece uma incorporação completamente determinística e exata em . Mas queríamos uma dimensão menor, ( n1O(nd). Aqui é ondeentrao comentário de Luca sobreo teorema de Carathéodory. O conjunto demétricas representativasde1forma um cone poliédrico no(n1(n2)1 espaço dimensional de todas as funções, de pares de pontos não ordenados a números reais, e o argumento geométrico acima diz que a métrica euclidiana pertence a esse cone. Os pontos nos raios extremos do cone sãopseudométricosunidimensionais1(nos quais os pontos são divididos em dois conjuntos, todas as distâncias em um único conjunto são zero e todas as distâncias na divisão são iguais), e Carathéodory diz que qualquer ponto dentro do cone (incluindo aquele com o qual nos preocupamos) pode ser representado como uma combinação convexa de pontos em raios extremos cujo número é no máximo a dimensão do espaço ambiente ( n(n2)1 . Mas uma combinação convexa de no máximo ( n(n2) métricasunidimensionais1é ( n(n2)1métrica.1(n2)

Por fim, como podemos realmente computar o incorporação dimensional? Neste ponto, não temos apenas um ponto no ( n(n2) cone convexo dimensional de1métrica (a métrica de distância com a qual começamos), mas também temos um conjunto deO(nd)pontos extremos do cone (correspondentes às partições da entrada em dois subconjuntos induzidos por hiperplanos) ) de forma que nossa métrica seja uma combinação convexa desses pontos extremos - paradpequeno, essa é uma grande melhoria em relação aos2n-2raios extremos que o cone possui em geral. Agora, tudo o que precisamos fazer é aplicar um algoritmo ganancioso que elimina pontos extremos do nosso conjunto, um por um, até apenas ( n(n2)1O(nd)d2n2 deles são deixados. Em cada etapa, precisamos manter como invariável que nossa métrica ainda está dentro do casco convexo dos pontos extremos restantes, o que é apenas um problema linear de viabilidade de programação. Se fizermos isso, o Carathéodory garantirá que sempre exista um conjunto de ( n(n2) pontos extremos cujo casco convexo contém a métrica de entrada.(n2)

David Eppstein
fonte