TSP 1d approximatif avec comparaisons linéaires?

21

O(nlogn)1+O(nc)cO(n)(maxmin)n(c+1)de sa valeur d'origine, puis utilisez le tri radix. Mais les modèles avec arrondi ont une théorie de complexité problématique et cela m'a amené à me demander, qu'en est-il des modèles de calcul plus faibles?

Ainsi, avec quelle précision le TSP unidimensionnel peut-il être approximé, dans un modèle d'arbre de comparaison linéaire de calcul (chaque nœud de comparaison teste le signe d'une fonction linéaire des valeurs d'entrée), par un algorithme dont la complexité temporelle est ? La même méthode d'arrondi permet d'obtenir n'importe quel rapport d'approximation de la forme n ^ {1-o (1)} (en utilisant des recherches binaires pour effectuer l'arrondi, et en arrondissant beaucoup plus grossièrement pour le rendre assez rapide). Mais est-il possible d'atteindre même un rapport d'approximation comme O (n ^ {1- \ epsilon}) pour certains \ epsilon> 0 ?n 1 - o ( 1 ) O ( n 1 - ϵ ) ϵ > 0o(nlogn)n1o(1)O(n1ϵ)ϵ>0

David Eppstein
la source
Je ne connais pas 1D TSP. Pouvez-vous le définir?
Tyson Williams
4
@Tyson Williams: Le problème de chemin du vendeur itinérant 1D est le cas particulier du problème de chemin du vendeur voyageur euclidien où toutes les villes sont sur l'axe des x. Ou formellement, on vous donne n nombres réels a_1,…, a_n, et votre objectif est de sortir une permutation π: {1,…, n} → {1,…, n} telle que ∑_ {i = 1} ^ {n − 1} | a_ {π (i)} - a_ {π (i + 1)} | est minimisé.
Tsuyoshi Ito

Réponses:

10

EDIT (UPDATE): La limite inférieure dans ma réponse ci-dessous a été prouvée (par une preuve différente) dans "Sur la complexité de l'approximation des tournées des vendeurs itinérants euclidiens et des arbres couvrant minimum", par Das et al; Algorithmica 19: 447-460 (1997).


est-il possible d'atteindre même un rapport d'approximation comme pendant un certain temps en utilisant un algorithme basé sur la comparaison?O(n1ϵ)o ( n log n )ϵ>0o(nlogn)

Non, voici une borne inférieure.

Prétendre. Pour tout , chaque algorithme d'approximation basé sur la comparaison nécessite des comparaisons dans le pire des cas.n 1 - ϵ Ω ( ϵ n log n )ϵ>0n1ϵΩ(ϵnlogn)

Par "basé sur la comparaison", j'entends tout algorithme qui interroge uniquement l'entrée avec des requêtes binaires (Vrai / Faux).

Voici une tentative de preuve. Espérons qu'il n'y ait pas d'erreurs. FWIW la borne inférieure semble susceptible de s'étendre aux algorithmes randomisés.


Fixez tout et tout arbitrairement petit mais constant .ϵ > 0nϵ>0

Considérez juste leinstances d'entrée de "permutation" qui sont des permutations de . La solution optimale pour une telle instance a coûté .( x 1 , x 2 , , x n ) [ n ] n - 1n!(x1,x2,,xn)[n]n1

Définissez le coût d'une permutation comme étant. Modélisez l'algorithme en prenant comme entrée une permutation , en sortant une permutation et en payant le coût .c ( π ) = i | π ( i + 1 ) - π ( i ) | π π d ( π , π ) = c ( π π )πc(π)=i|π(i+1)π(i)|ππd(π,π)=c(ππ)

Définissez comme étant le nombre minimum de comparaisons pour tout algorithme basé sur la comparaison pour atteindre le rapport compétitif sur ces instances. Puisque opt est , l'algorithme doit garantir le coût au plus .n 1 - ϵ n - 1 n 2 - ϵCn1ϵn1n2ϵ

Nous allons montrer .CΩ(ϵnlogn)

Définissez comme étant, pour toute sortie , la fraction des entrées possibles pour lesquelles la sortie coûterait au plus . Cette fraction est indépendante de .π π n 2 - ϵ π Pππn2ϵπ

π c ( π ) n 2 - ϵ π I P d ( π , I ) n 2 - ϵ d ( π , I ) = c ( π )P également égal à la probabilité que, pour une permutation aléatoire , son coût soit au plus . (Pour voir pourquoi, prenons comme permutation d'identité Alors est la fraction des entrées pour lesquelles au plus , mais .)πc(π)n2ϵπIPd(π,I)n2ϵd(π,I)=c(π)

Lemme 1 .Clog21/P

Preuve. Corrigez tout algorithme qui utilise toujours moins de comparaisons queL'arbre de décision de l'algorithme a une profondeur inférieure à , il y a donc moins de feuilles, et, pour certaines permutations de sortie , l'algorithme donne comme sortie pour plus d'un fraction des entrées. Par définition de , pour au moins une de ces entrées, la sortie donne un coût supérieur à . QEDlog 2 1 / P 1 / P π π P P π n 2 - ϵlog21/Plog21/P1/PππPPπn2ϵ

Lemme 2. .Pexp(Ω(ϵnlogn))

Avant de donner la preuve du lemme 2, notons que les deux lemmes donnent ensemble la revendication:

C  log21P = log2exp(Ω(ϵnlogn)) = Ω(ϵnlogn).

Preuve du lemme 2. Soit une permutation aléatoire. Rappelons que est égal à la probabilité que son coût soit au maximum . Disons que toute paire est une arête avec un coût, donc est la somme des coûts de bord.P c ( π ) n 2 - ϵ ( i , i + 1 ) | π ( i + 1 ) - π ( i ) | c ( π )πPc(π)n2ϵ(i,i+1)|π(i+1)π(i)|c(π)

Supposons que .c(π)n2ϵ

Ensuite, pour tout , au plus des arêtes ont coûté ou plus. Disons que les arêtes de coût inférieur à sont bon marché .n 2 - ϵ / q q qq>0n2ϵ/qqq

Fixez . Substituer et simplifier, au plus des bords ne sont pas bon marché. n 1 - ϵ / 2q=n1ϵ/2n1ϵ/2

Ainsi, au moins des arêtes sont bon marché. Ainsi, il existe un ensemble contenant bords bon marché.S n / 2nn1ϵ/2n/2Sn/2

Prétendre. Pour tout ensemble donné de arêtes, la probabilité que toutes les arêtes de soient bon marché est au plus .n / 2 S exp ( - Ω ( ϵ n log n ) )Sn/2Sexp(Ω(ϵnlogn))

Avant de prouver l'allégation, notez qu'elle implique le lemme comme suit. D'après la revendication et la liaison naïve, la probabilité qu'il existe un tel ensemble est au plus ( nS

(nn/2)exp(Ω(ϵnlogn))  2nexp(Ω(ϵnlogn))
  exp(O(n)Ω(ϵnlogn))  exp(Ω(ϵnlogn)).

Preuve de réclamation. Choisissez par le processus suivant. Choisissez uniformément dans , puis choisissez uniformément dans , puis choisissez uniformément dans , etc.π ( 1 ) [ n ] π ( 2 ) [ n ] - { π ( 1 ) }ππ(1)[n]π(2)[n]{π(1)}[ n ] - { π ( 1 ) ,π(3)[n]{π(1),π(2)}

Considérons une arête à . Considérez le temps juste après que a été choisi, quand est sur le point d'être choisi. Quels que soient les premiers choix (pour pour ), il y a au moins choix pour , et au plus de ceux-ci les choix donneront le bord moins cher que (ce qui le rend bon marché).S π ( i ) π ( i + 1 ) i π ( j ) j i n - i π ( i + 1 ) 2 n 1 - ϵ(i,i+1)Sπ(i)π(i+1)iπ(j)jiniπ(i+1) (i,i+1)n 1 -2n1ϵ/2(i,i+1)n1ϵ/2

Ainsi, conditionnée aux premiers choix, la probabilité que l'arête soit bon marché est au maximum . Ainsi, la probabilité que toutes les arêtes de soient bon marché est au plus Puisque , il y a au moins arêtes dans avec . Ainsi, ce produit est au maximum 2 n 1i n/2S(i,i+1)S2n 1 - ϵ / 22n1ϵ/2nin/2S| S| n/2

(i,i+1)S2n1ϵ/2ni.
|S|n/2S n - i n / 4 ( 2 n 1 - ϵ / 2n/4Snin/4
(2n1ϵ/2n/4)n/4  (8nϵ/2)n/4 = exp(O(n)Ω(ϵnlogn)) = exp(Ω(ϵnlogn)).

QED

Neal Young
la source
6
ps J'ai reçu une demande pour faire ce citable, donc je l'ai mis sur arvix.org ici .
Neal Young